Lines Matching refs:aux
37 PetscCall(MatDestroy(&data->aux)); in PCReset_HPDDM()
95 PetscCall(MatDestroy(&data->aux)); in PCHPDDMSetAuxiliaryMat_Private()
96 data->aux = A; in PCHPDDMSetAuxiliaryMat_Private()
122 Mat *splitting[2] = {}, aux; in PCHPDDMSetAuxiliaryMatNormal_Private() local
155 …etscCall(MatTransposeMatMult(conjugate, *splitting[0], MAT_INITIAL_MATRIX, PETSC_DETERMINE, &aux)); in PCHPDDMSetAuxiliaryMatNormal_Private()
158 PetscCall(MatNorm(aux, NORM_FROBENIUS, &norm)); in PCHPDDMSetAuxiliaryMatNormal_Private()
159 PetscCall(MatSetOption(aux, MAT_NEW_NONZERO_ALLOCATION_ERR, PETSC_FALSE)); in PCHPDDMSetAuxiliaryMatNormal_Private()
175 PetscCall(MatDiagonalSet(aux, d, ADD_VALUES)); in PCHPDDMSetAuxiliaryMatNormal_Private()
179 if (!diagonal) PetscCall(MatShift(aux, PETSC_SMALL * norm)); in PCHPDDMSetAuxiliaryMatNormal_Private()
180 PetscCall(MatEliminateZeros(aux, PETSC_TRUE)); in PCHPDDMSetAuxiliaryMatNormal_Private()
191 if (flg) PetscCall(MatCreateNormal(*splitting[0], &aux)); in PCHPDDMSetAuxiliaryMatNormal_Private()
192 else PetscCall(MatCreateNormalHermitian(*splitting[0], &aux)); in PCHPDDMSetAuxiliaryMatNormal_Private()
195 PetscCall(PCHPDDMSetAuxiliaryMat(pc, *is, aux, nullptr, nullptr)); in PCHPDDMSetAuxiliaryMatNormal_Private()
198 PetscCall(MatDestroy(&aux)); in PCHPDDMSetAuxiliaryMatNormal_Private()
618 Mat aux[2]; in PCView_HPDDM() local
632 if (data->aux) { in PCView_HPDDM()
633 PetscCall(MatGetSize(data->aux, &m, &n)); in PCView_HPDDM()
634 PetscCall(MatCreate(PetscObjectComm((PetscObject)pc), aux)); in PCView_HPDDM()
635 PetscCall(MatSetSizes(aux[0], m, n, PETSC_DETERMINE, PETSC_DETERMINE)); in PCView_HPDDM()
636 PetscCall(PetscObjectBaseTypeCompare((PetscObject)data->aux, MATSEQAIJ, &flg)); in PCView_HPDDM()
637 if (flg) PetscCall(MatSetType(aux[0], MATMPIAIJ)); in PCView_HPDDM()
639 PetscCall(PetscObjectBaseTypeCompare((PetscObject)data->aux, MATSEQBAIJ, &flg)); in PCView_HPDDM()
640 if (flg) PetscCall(MatSetType(aux[0], MATMPIBAIJ)); in PCView_HPDDM()
642 PetscCall(PetscObjectBaseTypeCompare((PetscObject)data->aux, MATSEQSBAIJ, &flg)); in PCView_HPDDM()
643 …any of the following: MATSEQAIJ, MATSEQBAIJ, or MATSEQSBAIJ", ((PetscObject)data->aux)->type_name); in PCView_HPDDM()
644 PetscCall(MatSetType(aux[0], MATMPISBAIJ)); in PCView_HPDDM()
647 PetscCall(MatSetBlockSizesFromMats(aux[0], data->aux, data->aux)); in PCView_HPDDM()
648 PetscCall(MatAssemblyBegin(aux[0], MAT_FINAL_ASSEMBLY)); in PCView_HPDDM()
649 PetscCall(MatAssemblyEnd(aux[0], MAT_FINAL_ASSEMBLY)); in PCView_HPDDM()
650 PetscCall(MatGetDiagonalBlock(aux[0], aux + 1)); in PCView_HPDDM()
651 PetscCall(MatCopy(data->aux, aux[1], DIFFERENT_NONZERO_PATTERN)); in PCView_HPDDM()
653 PetscCall(MatView(aux[0], subviewer)); in PCView_HPDDM()
655 PetscCall(MatDestroy(aux)); in PCView_HPDDM()
1263 …PetscCallBack("PCHPDDM Neumann callback", (*data->setup)(data->aux, t, x, xt, s, data->is, data->s… in PCHPDDMSetUpNeumannOverlap_Private()
1531 Mat aux; in PCHPDDMAlgebraicAuxiliaryMat_Private() local
1536 PetscCall(MatDuplicate((*sub)[0], MAT_COPY_VALUES, &aux)); in PCHPDDMAlgebraicAuxiliaryMat_Private()
1537 PetscCall(MatSetOption(aux, MAT_NEW_NONZERO_ALLOCATION_ERR, PETSC_FALSE)); in PCHPDDMAlgebraicAuxiliaryMat_Private()
1541 PetscCall(MatCreateVecs(aux, sum, sum + 1)); in PCHPDDMAlgebraicAuxiliaryMat_Private()
1543 PetscCall(MatGetRowSum(aux, sum[1])); in PCHPDDMAlgebraicAuxiliaryMat_Private()
1547 PetscCall(MatDiagonalSet(aux, sum[0], ADD_VALUES)); in PCHPDDMAlgebraicAuxiliaryMat_Private()
1564 PetscCall(MatCreateDense(PETSC_COMM_SELF, aux->cmap->n, bs, aux->cmap->n, bs, ptr, &ones)); in PCHPDDMAlgebraicAuxiliaryMat_Private()
1566 PetscCall(MatMatMult(aux, ones, MAT_INITIAL_MATRIX, PETSC_CURRENT, sum + 1)); in PCHPDDMAlgebraicAuxiliaryMat_Private()
1578 …for (PetscInt n = 0; n < aux->cmap->n / bs; ++n) PetscCall(MatSetValuesBlocked(aux, 1, &n, 1, &n, … in PCHPDDMAlgebraicAuxiliaryMat_Private()
1579 PetscCall(MatAssemblyBegin(aux, MAT_FINAL_ASSEMBLY)); in PCHPDDMAlgebraicAuxiliaryMat_Private()
1580 PetscCall(MatAssemblyEnd(aux, MAT_FINAL_ASSEMBLY)); in PCHPDDMAlgebraicAuxiliaryMat_Private()
1584 PetscCall(MatSetOption(aux, MAT_NEW_NONZERO_ALLOCATION_ERR, PETSC_TRUE)); in PCHPDDMAlgebraicAuxiliaryMat_Private()
1586 PetscCall(PetscObjectCompose((PetscObject)(*sub)[0], "_PCHPDDM_Neumann_Mat", (PetscObject)aux)); in PCHPDDMAlgebraicAuxiliaryMat_Private()
1891 PetscCall(MatDestroy(&data->aux)); in PCSetUp_HPDDM()
1898 PetscCall(MatDestroy(&data->aux)); in PCSetUp_HPDDM()
2001 if (!data->aux) { in PCSetUp_HPDDM()
2018 … build an auxiliary Mat (which was%s initially provided by the user)\n", data->aux ? "" : " not")); in PCSetUp_HPDDM()
2019 PetscCall(MatDestroy(&data->aux)); in PCSetUp_HPDDM()
2023 …if (!data->aux) { /* if A11 is near zero, e.g., Stokes equation, or diagonal, build an auxiliary (… in PCSetUp_HPDDM()
2028 PetscCall(MatDestroy(&data->aux)); in PCSetUp_HPDDM()
2058 PetscCall(MatCreateDiagonal(v, &data->aux)); in PCSetUp_HPDDM()
2062 uaux = data->aux; in PCSetUp_HPDDM()
2119 PetscCall(MatGetDiagonalBlock(A11, &data->aux)); in PCSetUp_HPDDM()
2120 PetscCall(PetscObjectReference((PetscObject)data->aux)); in PCSetUp_HPDDM()
2212 PetscCall(PetscObjectTypeCompare((PetscObject)data->aux, MATSEQAIJ, &flg)); in PCSetUp_HPDDM()
2253 if (data->aux) PetscCall(MatDuplicate(data->aux, MAT_COPY_VALUES, &daux)); in PCSetUp_HPDDM()
2321 if ((ctx || data->N > 1) && (data->aux || ismatis || algebraic)) { in PCSetUp_HPDDM()
2345 …etscCall(PetscObjectQuery((PetscObject)sub[0], "_PCHPDDM_Neumann_Mat", (PetscObject *)&data->aux)); in PCSetUp_HPDDM()
2350 if (PetscBool3ToBool(data->Neumann)) sub = &data->aux; in PCSetUp_HPDDM()
2526 …(MatCreateShell(PETSC_COMM_SELF, P->rmap->n, n[1] - n[0], P->rmap->n, n[1] - n[0], h, &data->aux)); in PCSetUp_HPDDM()
2529 … PetscCall(MatShellSetOperation(data->aux, MATOP_MULT, (PetscErrorCodeFn *)MatMult_Harmonic)); in PCSetUp_HPDDM()
2530 …PetscCall(MatShellSetOperation(data->aux, MATOP_MULT_TRANSPOSE, (PetscErrorCodeFn *)MatMultTranspo… in PCSetUp_HPDDM()
2531 …PetscCall(MatShellSetMatProductOperation(data->aux, MATPRODUCT_AB, nullptr, MatProduct_AB_Harmonic… in PCSetUp_HPDDM()
2532 …PetscCall(MatShellSetMatProductOperation(data->aux, MATPRODUCT_AtB, nullptr, MatProduct_AtB_Harmon… in PCSetUp_HPDDM()
2533 …PetscCall(MatShellSetOperation(data->aux, MATOP_DESTROY, (PetscErrorCodeFn *)MatDestroy_Harmonic)); in PCSetUp_HPDDM()
2598 …rivate(unsorted, data->is, &uis, PetscBool3ToBool(data->Neumann) ? sub[0] : data->aux, &C, &perm)); in PCSetUp_HPDDM()
2628 PetscCall(MatAXPY(D, 1.0, data->aux, structure)); in PCSetUp_HPDDM()
2646 std::swap(C, data->aux); in PCSetUp_HPDDM()
2671 …PetscCall(PCHPDDMPermute_Private(unsorted, data->is, &uis, data->aux, &C, nullptr)); /* permute si… in PCSetUp_HPDDM()
2672 std::swap(C, data->aux); in PCSetUp_HPDDM()
2719 b[3] = data->aux; in PCSetUp_HPDDM()
2814 …s[0]->P->structure(loc, data->is, !ctx ? sub[0] : nullptr, ismatis ? C : data->aux, data->levels)); in PCSetUp_HPDDM()
2818 if (data->deflation || overlap != -1) weighted = data->aux; in PCSetUp_HPDDM()
2856 …matis ? C : (algebraic && !block && overlap == -1 ? sub[0] : (!ctx ? data->aux : S)), weighted, da… in PCSetUp_HPDDM()
2867 else N = data->aux; in PCSetUp_HPDDM()
3002 std::swap(C, data->aux); in PCSetUp_HPDDM()
3007 if (algebraic) PetscCall(MatDestroy(&data->aux)); in PCSetUp_HPDDM()
3019 …aux && daux) || (!data->aux && !daux), PETSC_COMM_SELF, PETSC_ERR_PLIB, "A Mat pointer is NULL but… in PCSetUp_HPDDM()
3020 if (data->aux) { in PCSetUp_HPDDM()
3021 PetscCall(MatMultEqual(data->aux, daux, 20, &flg)); in PCSetUp_HPDDM()