impls/hpddm/pchpddm.cxx

1770   Mat                                       *sub, A, P, N, C = nullptr, uaux = nullptr, weighted, subA[2], S;
1855         PetscCall((*create)(dm, &uis, &uaux, &usetup, &uctx));
1864       if (!uaux) {
1865         PetscCall(PetscObjectQuery((PetscObject)pc, "_PCHPDDM_Neumann_Mat", (PetscObject *)&uaux));
1866         PetscCall(PetscObjectReference((PetscObject)uaux));
1873       if (!uaux) {
1874         PetscCall(PetscObjectQuery((PetscObject)P, "_PCHPDDM_Neumann_Mat", (PetscObject *)&uaux));
1875         PetscCall(PetscObjectReference((PetscObject)uaux));
1878     PetscCall(PCHPDDMSetAuxiliaryMat(pc, uis, uaux, usetup, uctx));
1879     PetscCall(MatDestroy(&uaux));
2062               uaux = data->aux;
2064               PetscCall(PetscObjectReference((PetscObject)uaux));
2086               PetscCall(PCHPDDMSetAuxiliaryMat(std::get<0>(*ctx)[1], uis, uaux, nullptr, nullptr)); /* transfer ownership of the auxiliary inputs from the inner (PCKSP) to the inner-most (PCHPDDM) PC */
2091                 PetscCall(PCHPDDMCheckMatStructure_Private(pc, sub[0], uaux));
2096               PetscCall(PetscObjectDereference((PetscObject)uaux));
2337             PetscCall(MatConvert(P, MATMPIBAIJ, MAT_INITIAL_MATRIX, &uaux));
2348       } else if (!uaux || overlap != -1) {
2388               PetscCall(MatCreateSubMatrices(uaux ? uaux : P, 1, ov + !flg, ov + 1, MAT_INITIAL_MATRIX, &a)); /* submatrix from above, either square (!flg) or rectangular (flg) */
2410                 if (uaux) PetscCall(MatConvert(a[0], MATSEQSBAIJ, MAT_INPLACE_MATRIX, a)); /* initial Pmat was MATSBAIJ, convert back to the same format since the rectangular A_12 submatrix has been created */
2431                 if (uaux) { /* initial Pmat was MATSBAIJ, convert back to the same format since this submatrix is square */
2524                 if (uaux) PetscCall(MatDestroy(&uaux));
2537               PetscCall(MatCreateSubMatrices(uaux ? uaux : P, 1, is, is, MAT_INITIAL_MATRIX, &sub));
2545             if (uaux) {
2546               PetscCall(MatDestroy(&uaux));
2552         PetscCall(MatCreateSubMatrices(uaux, 1, is, is, MAT_INITIAL_MATRIX, &sub));
2553         PetscCall(MatDestroy(&uaux));