| /petsc/src/ksp/ksp/tests/ |
| H A D | ex22.c | 7 Mat A11, A12, A21, A22, A, tmp[2][2]; in test_solve() local 28 PetscCall(MatCreate(PETSC_COMM_WORLD, &A11)); in test_solve() 29 PetscCall(MatSetSizes(A11, PETSC_DECIDE, PETSC_DECIDE, n, n)); in test_solve() 30 PetscCall(MatSetType(A11, MATAIJ)); in test_solve() 31 PetscCall(MatSeqAIJSetPreallocation(A11, n, NULL)); in test_solve() 32 PetscCall(MatMPIAIJSetPreallocation(A11, np, NULL, np, NULL)); in test_solve() 33 PetscCall(MatDiagonalSet(A11, diag, INSERT_VALUES)); in test_solve() 57 tmp[0][0] = A11; in test_solve() 106 PetscCall(MatDestroy(&A11)); in test_solve() 118 Mat A11, A12, A21, A; in test_solve_matgetvecs() local [all …]
|
| H A D | ex11.c | 11 Mat A11, A22, A12, A21; in replace_submats() local 16 PetscCall(MatCreateSubMatrix(A, isu, isu, MAT_INITIAL_MATRIX, &A11)); in replace_submats() 20 PetscCall(MatDuplicate(A11, MAT_COPY_VALUES, &nA11)); in replace_submats() 24 PetscCall(MatGetOptionsPrefix(A11, &prefix)); in replace_submats() 32 PetscCall(MatDestroy(&A11)); in replace_submats() 43 PetscErrorCode LSCLoadTestOperators(Mat *A11, Mat *A12, Mat *A21, Mat *A22, Vec *b1, Vec *b2) in LSCLoadTestOperators() argument 50 PetscCall(MatCreate(PETSC_COMM_WORLD, A11)); in LSCLoadTestOperators() 54 PetscCall(MatSetOptionsPrefix(*A11, "a11_")); in LSCLoadTestOperators() 56 PetscCall(MatSetFromOptions(*A11)); in LSCLoadTestOperators() 64 PetscCall(MatLoad(*A11, viewer)); in LSCLoadTestOperators() [all …]
|
| /petsc/src/mat/tutorials/ |
| H A D | ex2.c | 13 Mat A, A11, A12, A21, A22; in main() local 61 PetscCall(MatCreate(PETSC_COMM_SELF, &A11)); in main() 62 PetscCall(MatSetSizes(A11, size1, size1, size1, size1)); in main() 63 PetscCall(MatSetType(A11, MATSEQDENSE)); in main() 64 PetscCall(MatSeqDenseSetPreallocation(A11, b)); in main() 65 PetscCall(MatDenseSetLDA(A11, size)); in main() 66 PetscCall(MatAssemblyBegin(A11, MAT_FINAL_ASSEMBLY)); in main() 67 PetscCall(MatAssemblyEnd(A11, MAT_FINAL_ASSEMBLY)); in main() 103 PetscCall(MatMult(A11, X1, Z1)); in main() 133 PetscCall(MatMult(A11, X1, Z1)); in main() [all …]
|
| /petsc/src/binding/petsc4py/test/ |
| H A D | test_mat_schur.py | 18 A11 = PETSc.Mat().create(comm=COMM) 19 A11.setType(TYPE) 20 A11.setSizes([[3, 3*comm_size], [3, 3*comm_size]]) 21 S = PETSc.Mat().createSchurComplement(A00, A00, A01, A10, A11) 33 self.assertEqual(A11_dup.id, A11.id)
|
| /petsc/src/ksp/ksp/tutorials/output/ |
| H A D | ex87_1_system-stokes.out | 43 Schur complement A11 - A10 inv(A00) A01 44 A11 79 Schur complement A11 - A10 inv(A00) A01 80 A11 92 Schur complement A11 - A10 inv(A00) A01 93 A11 131 Schur complement A11 - A10 inv(A00) A01 132 A11 144 Schur complement A11 - A10 inv(A00) A01 145 A11
|
| H A D | ex87_1_system-elasticity.out | 46 Schur complement A11 - A10 inv(A00) A01 47 A11 82 Schur complement A11 - A10 inv(A00) A01 83 A11 95 Schur complement A11 - A10 inv(A00) A01 96 A11 134 Schur complement A11 - A10 inv(A00) A01 135 A11 147 Schur complement A11 - A10 inv(A00) A01 148 A11
|
| /petsc/src/ksp/ksp/utils/schurm/ |
| H A D | schurm.c | 175 PetscErrorCode MatCreateSchurComplement(Mat A00, Mat Ap00, Mat A01, Mat A10, Mat A11, Mat *S) in MatCreateSchurComplement() argument 181 PetscCall(MatSchurComplementSetSubMatrices(*S, A00, Ap00, A01, A10, A11)); in MatCreateSchurComplement() 217 PetscErrorCode MatSchurComplementSetSubMatrices(Mat S, Mat A00, Mat Ap00, Mat A01, Mat A10, Mat A11) in MatSchurComplementSetSubMatrices() argument 238 if (A11) { in MatSchurComplementSetSubMatrices() 239 PetscValidHeaderSpecific(A11, MAT_CLASSID, 6); in MatSchurComplementSetSubMatrices() 240 PetscCheckSameComm(A00, 2, A11, 6); in MatSchurComplementSetSubMatrices() 241 …= A11->rmap->n, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "Local rows of A10 %" PetscInt_FMT " do not eq… in MatSchurComplementSetSubMatrices() 249 PetscCall(PetscObjectReference((PetscObject)A11)); in MatSchurComplementSetSubMatrices() 254 Na->D = A11; in MatSchurComplementSetSubMatrices() 359 …cErrorCode MatSchurComplementUpdateSubMatrices(Mat S, Mat A00, Mat Ap00, Mat A01, Mat A10, Mat A11) in MatSchurComplementUpdateSubMatrices() argument [all …]
|
| /petsc/src/mat/tests/ |
| H A D | ex159.c | 119 Mat A00, A11, A0a0a, A0a0b; in main() local 121 PetscCall(MatGetLocalSubMatrix(A, isl1, isl1, &A11)); in main() 131 PetscCall(MatSetValueLocal(A11, 0, 0, 1000 * (rank + 1) + 1, ADD_VALUES)); in main() 132 PetscCall(MatSetValueLocal(A11, 1, 2, 1000 * (rank + 1) + 6, ADD_VALUES)); in main() 137 PetscCall(MatRestoreLocalSubMatrix(A, isl1, isl1, &A11)); in main()
|
| /petsc/src/ksp/ksp/tests/output/ |
| H A D | ex11_test_fs-0_fc_pc_fieldsplit_off_diag_use_amat-0_fc_pc_fieldsplit_diag_use_amat-1.out | 83 Preconditioner for the Schur complement formed from A11 113 KSP solver for S = A11 - A10 inv(A00) A01 157 Schur complement A11 - A10 inv(A00) A01 158 A11 258 Preconditioner for the Schur complement formed from A11 288 KSP solver for S = A11 - A10 inv(A00) A01 332 Schur complement A11 - A10 inv(A00) A01 333 A11
|
| H A D | ex11_test_fs-0_fc_pc_fieldsplit_off_diag_use_amat-1_fc_pc_fieldsplit_diag_use_amat-1.out | 83 Preconditioner for the Schur complement formed from A11 113 KSP solver for S = A11 - A10 inv(A00) A01 157 Schur complement A11 - A10 inv(A00) A01 158 A11 258 Preconditioner for the Schur complement formed from A11 288 KSP solver for S = A11 - A10 inv(A00) A01 332 Schur complement A11 - A10 inv(A00) A01 333 A11
|
| H A D | ex11_test_fs-0_fc_pc_fieldsplit_off_diag_use_amat-0_fc_pc_fieldsplit_diag_use_amat-0.out | 83 Preconditioner for the Schur complement formed from A11 113 KSP solver for S = A11 - A10 inv(A00) A01 157 Schur complement A11 - A10 inv(A00) A01 158 A11 258 Preconditioner for the Schur complement formed from A11 288 KSP solver for S = A11 - A10 inv(A00) A01 332 Schur complement A11 - A10 inv(A00) A01 333 A11
|
| H A D | ex11_test_fs-0_fc_pc_fieldsplit_off_diag_use_amat-1_fc_pc_fieldsplit_diag_use_amat-0.out | 83 Preconditioner for the Schur complement formed from A11 113 KSP solver for S = A11 - A10 inv(A00) A01 157 Schur complement A11 - A10 inv(A00) A01 158 A11 258 Preconditioner for the Schur complement formed from A11 288 KSP solver for S = A11 - A10 inv(A00) A01 332 Schur complement A11 - A10 inv(A00) A01 333 A11
|
| H A D | ex11_test_fs-1_fc_pc_fieldsplit_off_diag_use_amat-1_fc_pc_fieldsplit_diag_use_amat-0.out | 83 Preconditioner for the Schur complement formed from A11 113 KSP solver for S = A11 - A10 inv(A00) A01 157 Schur complement A11 - A10 inv(A00) A01 158 A11 266 Preconditioner for the Schur complement formed from A11 296 KSP solver for S = A11 - A10 inv(A00) A01 340 Schur complement A11 - A10 inv(A00) A01 341 A11
|
| H A D | ex11_test_fs-1_fc_pc_fieldsplit_off_diag_use_amat-1_fc_pc_fieldsplit_diag_use_amat-1.out | 83 Preconditioner for the Schur complement formed from A11 119 KSP solver for S = A11 - A10 inv(A00) A01 163 Schur complement A11 - A10 inv(A00) A01 164 A11 272 Preconditioner for the Schur complement formed from A11 308 KSP solver for S = A11 - A10 inv(A00) A01 352 Schur complement A11 - A10 inv(A00) A01 353 A11
|
| H A D | ex11_test_fs-1_fc_pc_fieldsplit_off_diag_use_amat-0_fc_pc_fieldsplit_diag_use_amat-1.out | 83 Preconditioner for the Schur complement formed from A11 119 KSP solver for S = A11 - A10 inv(A00) A01 163 Schur complement A11 - A10 inv(A00) A01 164 A11 272 Preconditioner for the Schur complement formed from A11 308 KSP solver for S = A11 - A10 inv(A00) A01 352 Schur complement A11 - A10 inv(A00) A01 353 A11
|
| H A D | ex11_test_fs-1_fc_pc_fieldsplit_off_diag_use_amat-0_fc_pc_fieldsplit_diag_use_amat-0.out | 83 Preconditioner for the Schur complement formed from A11 113 KSP solver for S = A11 - A10 inv(A00) A01 157 Schur complement A11 - A10 inv(A00) A01 158 A11 266 Preconditioner for the Schur complement formed from A11 296 KSP solver for S = A11 - A10 inv(A00) A01 340 Schur complement A11 - A10 inv(A00) A01 341 A11
|
| H A D | ex11_2_test_fs-0_fc_pc_fieldsplit_off_diag_use_amat-0_fc_pc_fieldsplit_diag_use_amat-0.out | 65 Preconditioner for the Schur complement formed from A11 113 KSP solver for S = A11 - A10 inv(A00) A01 175 Schur complement A11 - A10 inv(A00) A01 176 A11 272 Preconditioner for the Schur complement formed from A11 320 KSP solver for S = A11 - A10 inv(A00) A01 382 Schur complement A11 - A10 inv(A00) A01 383 A11
|
| H A D | ex11_2_test_fs-0_fc_pc_fieldsplit_off_diag_use_amat-1_fc_pc_fieldsplit_diag_use_amat-0.out | 65 Preconditioner for the Schur complement formed from A11 113 KSP solver for S = A11 - A10 inv(A00) A01 175 Schur complement A11 - A10 inv(A00) A01 176 A11 272 Preconditioner for the Schur complement formed from A11 320 KSP solver for S = A11 - A10 inv(A00) A01 382 Schur complement A11 - A10 inv(A00) A01 383 A11
|
| H A D | ex11_2_test_fs-0_fc_pc_fieldsplit_off_diag_use_amat-1_fc_pc_fieldsplit_diag_use_amat-1.out | 65 Preconditioner for the Schur complement formed from A11 113 KSP solver for S = A11 - A10 inv(A00) A01 175 Schur complement A11 - A10 inv(A00) A01 176 A11 272 Preconditioner for the Schur complement formed from A11 320 KSP solver for S = A11 - A10 inv(A00) A01 382 Schur complement A11 - A10 inv(A00) A01 383 A11
|
| H A D | ex11_2_test_fs-0_fc_pc_fieldsplit_off_diag_use_amat-0_fc_pc_fieldsplit_diag_use_amat-1.out | 65 Preconditioner for the Schur complement formed from A11 113 KSP solver for S = A11 - A10 inv(A00) A01 175 Schur complement A11 - A10 inv(A00) A01 176 A11 272 Preconditioner for the Schur complement formed from A11 320 KSP solver for S = A11 - A10 inv(A00) A01 382 Schur complement A11 - A10 inv(A00) A01 383 A11
|
| H A D | ex11_2_test_fs-1_fc_pc_fieldsplit_off_diag_use_amat-1_fc_pc_fieldsplit_diag_use_amat-1.out | 65 Preconditioner for the Schur complement formed from A11 119 KSP solver for S = A11 - A10 inv(A00) A01 181 Schur complement A11 - A10 inv(A00) A01 182 A11 286 Preconditioner for the Schur complement formed from A11 340 KSP solver for S = A11 - A10 inv(A00) A01 402 Schur complement A11 - A10 inv(A00) A01 403 A11
|
| H A D | ex11_2_test_fs-1_fc_pc_fieldsplit_off_diag_use_amat-0_fc_pc_fieldsplit_diag_use_amat-1.out | 65 Preconditioner for the Schur complement formed from A11 119 KSP solver for S = A11 - A10 inv(A00) A01 181 Schur complement A11 - A10 inv(A00) A01 182 A11 286 Preconditioner for the Schur complement formed from A11 340 KSP solver for S = A11 - A10 inv(A00) A01 402 Schur complement A11 - A10 inv(A00) A01 403 A11
|
| H A D | ex11_2_test_fs-1_fc_pc_fieldsplit_off_diag_use_amat-1_fc_pc_fieldsplit_diag_use_amat-0.out | 65 Preconditioner for the Schur complement formed from A11 113 KSP solver for S = A11 - A10 inv(A00) A01 175 Schur complement A11 - A10 inv(A00) A01 176 A11 280 Preconditioner for the Schur complement formed from A11 328 KSP solver for S = A11 - A10 inv(A00) A01 390 Schur complement A11 - A10 inv(A00) A01 391 A11
|
| H A D | ex11_2_test_fs-1_fc_pc_fieldsplit_off_diag_use_amat-0_fc_pc_fieldsplit_diag_use_amat-0.out | 65 Preconditioner for the Schur complement formed from A11 113 KSP solver for S = A11 - A10 inv(A00) A01 175 Schur complement A11 - A10 inv(A00) A01 176 A11 280 Preconditioner for the Schur complement formed from A11 328 KSP solver for S = A11 - A10 inv(A00) A01 390 Schur complement A11 - A10 inv(A00) A01 391 A11
|
| /petsc/src/ksp/pc/impls/hpddm/ |
| H A D | pchpddm.cxx | 1907 Mat A11 = nullptr; in PCSetUp_HPDDM() local 1927 … PetscCall(MatSchurComplementGetSubMatrices(P, nullptr, nullptr, nullptr, nullptr, &A11)); in PCSetUp_HPDDM() 1928 PetscCall(MatGetOwnershipRanges(A11, &ranges)); in PCSetUp_HPDDM() 1929 PetscCallMPI(MPI_Comm_size(PetscObjectComm((PetscObject)A11), &size)); in PCSetUp_HPDDM() 1974 …PetscCall(ISCreateStride(PetscObjectComm((PetscObject)data->is), A11->rmap->n, A11->rmap->rstart, … in PCSetUp_HPDDM() 2004 PetscCall(MatCreateVecs(A11, &d, nullptr)); in PCSetUp_HPDDM() 2005 PetscCall(MatGetDiagonal(A11, d)); in PCSetUp_HPDDM() 2006 …PetscCall(PetscObjectTypeCompareAny((PetscObject)A11, &flg, MATDIAGONAL, MATCONSTANTDIAGONAL, "")); in PCSetUp_HPDDM() 2009 PetscCall(MatMultEqual(A11, D, 20, &flg)); in PCSetUp_HPDDM() 2014 if ((PetscDefined(USE_DEBUG) || (data->Neumann != PETSC_BOOL3_TRUE && !flg)) && A11) { in PCSetUp_HPDDM() [all …]
|