Home
last modified time | relevance | path

Searched refs:AinvB (Results 1 – 3 of 3) sorted by relevance

/petsc/src/dm/dt/tests/
H A Dex13.c63 Mat AinvB; in test() local
141 PetscCall(MatDuplicate(mat_moments_T, MAT_DO_NOT_COPY_VALUES, &AinvB)); in test()
143 PetscCall(MatMatSolve(mat_trimmed, mat_moments_T, AinvB)); in test()
144 PetscCall(MatTransposeMatMult(mat_moments_T, AinvB, MAT_INITIAL_MATRIX, PETSC_DETERMINE, &Mm1)); in test()
149 PetscCall(MatDestroy(&AinvB)); in test()
163 Mat AinvB; in test() local
229 PetscCall(MatDuplicate(mat_k_moment_T, MAT_DO_NOT_COPY_VALUES, &AinvB)); in test()
230 PetscCall(MatMatSolve(mat_trimmed, mat_k_moment_T, AinvB)); in test()
231 … PetscCall(MatTransposeMatMult(mat_k_moment_T, AinvB, MAT_INITIAL_MATRIX, PETSC_DETERMINE, &prod)); in test()
237 PetscCall(MatDestroy(&AinvB)); in test()
/petsc/src/ksp/ksp/utils/schurm/
H A Dschurm.c524 Mat AinvB; in MatSchurComplementComputeExplicitOperator() local
532 …cObjectComm((PetscObject)A), AinvBd->rmap->n, A->cmap->n, AinvBd->rmap->N, A->cmap->N, v, &AinvB)); in MatSchurComplementComputeExplicitOperator()
533 PetscCall(MatDenseCUDAReplaceArray(AinvB, v)); in MatSchurComplementComputeExplicitOperator()
541 …cObjectComm((PetscObject)A), AinvBd->rmap->n, A->cmap->n, AinvBd->rmap->N, A->cmap->N, v, &AinvB)); in MatSchurComplementComputeExplicitOperator()
542 PetscCall(MatDenseHIPReplaceArray(AinvB, v)); in MatSchurComplementComputeExplicitOperator()
547 …cObjectComm((PetscObject)A), AinvBd->rmap->n, A->cmap->n, AinvBd->rmap->N, A->cmap->N, v, &AinvB)); in MatSchurComplementComputeExplicitOperator()
548 PetscCall(MatDenseReplaceArray(AinvB, v)); // let MatDestroy() free the data pointer in MatSchurComplementComputeExplicitOperator()
552 …PetscCall(MatHeaderReplace(AinvBd, &AinvB)); // replace the input composed Mat with just A00^-1 A0… in MatSchurComplementComputeExplicitOperator()
/petsc/src/ksp/pc/impls/fieldsplit/
H A Dfieldsplit.c1193 Mat AinvB = NULL; in PCApply_FieldSplit_Schur() local
1274 PetscCall(PetscObjectQuery((PetscObject)jac->schur, "AinvB", (PetscObject *)&AinvB)); in PCApply_FieldSplit_Schur()
1275 if (AinvB) { in PCApply_FieldSplit_Schur()
1276 PetscCall(MatGetSize(AinvB, NULL, &N)); in PCApply_FieldSplit_Schur()
1285 PetscCall(MatDenseGetArrayAndMemType(AinvB, &array, &mtype)); in PCApply_FieldSplit_Schur()
1293 PetscCall(MatDenseRestoreArrayAndMemType(AinvB, &array)); in PCApply_FieldSplit_Schur()
1319 if (!AinvB) { in PCApply_FieldSplit_Schur()
1326 } else PetscCall(MatMultAdd(AinvB, ilinkD->y, ilinkA->y, ilinkA->y)); in PCApply_FieldSplit_Schur()
1397 Mat AinvB = NULL; in PCMatApply_FieldSplit_Schur() local
1461 PetscCall(PetscObjectQuery((PetscObject)jac->schur, "AinvB", (PetscObject *)&AinvB)); in PCMatApply_FieldSplit_Schur()
[all …]