| /petsc/src/mat/tests/ |
| H A D | ex68.c | 11 IS isrow, iscol, identity; in main() local
49 PetscCall(MatGetOrdering(mat, MATORDERINGNATURAL, &isrow, &iscol)); in main()
51 PetscCall(MatPermute(mat, isrow, iscol, &B)); in main()
56 PetscCall(MatReorderForNonzeroDiagonal(mat, 1.e-8, isrow, iscol)); in main()
57 PetscCall(MatPermute(mat, isrow, iscol, &B)); in main()
61 PetscCall(ISView(isrow, viewer)); in main()
66 PetscCall(ISDestroy(&isrow)); in main()
69 PetscCall(MatGetOrdering(mat, MATORDERINGND, &isrow, &iscol)); in main()
70 PetscCall(MatPermute(mat, isrow, iscol, &B)); in main()
75 PetscCall(ISView(isrow, viewer)); in main()
[all …]
|
| H A D | ex35.c | 10 IS isrow; in main() local
44 PetscCall(ISCreateStride(PETSC_COMM_WORLD, m, 0, 1, &isrow)); in main()
45 PetscCall(MatCreateSubMatrices(A, 1, &isrow, &isrow, MAT_INITIAL_MATRIX, &Bsub)); in main()
47 PetscCall(ISDestroy(&isrow)); in main()
52 PetscCall(ISCreateStride(PETSC_COMM_WORLD, m, 0, 2, &isrow)); in main()
53 PetscCall(MatCreateSubMatrices(A, 1, &isrow, &isrow, MAT_INITIAL_MATRIX, &Bsub)); in main()
55 PetscCall(ISDestroy(&isrow)); in main()
60 PetscCall(ISCreateStride(PETSC_COMM_WORLD, m, N - m - 1, 1, &isrow)); in main()
61 PetscCall(MatCreateSubMatrices(A, 1, &isrow, &isrow, MAT_INITIAL_MATRIX, &Bsub)); in main()
63 PetscCall(ISDestroy(&isrow)); in main()
|
| H A D | ex67f.F90 | 14 IS isrow
35 PetscCallA(ISCreateGeneral(PETSC_COMM_SELF, nis, zero, PETSC_COPY_VALUES, isrow, ierr))
37 PetscCallA(MatCreateSubmatrices(A, nis, [isrow], [isrow], MAT_INITIAL_MATRIX, B, ierr))
43 PetscCallA(MatCreateSubmatrices(A, nis, [isrow], [isrow], MAT_REUSE_MATRIX, B, ierr))
49 PetscCallA(ISDestroy(isrow, ierr))
|
| H A D | ex151.c | 33 IS isrow, iscol; in main() local
52 …cCall(ISCreateGeneral(PETSC_COMM_WORLD, rend - rstart, ixrow + rstart, PETSC_USE_POINTER, &isrow)); in main()
54 PetscCall(ISSetPermutation(isrow)); in main()
64 PetscCall(MatPermute(A, isrow, iscol, &B)); in main()
70 PetscCall(ISView(isrow, viewer)); in main()
75 PetscCall(ISDestroy(&isrow)); in main()
|
| H A D | ex181.c | 11 IS isrow; in main() local
50 PetscCall(ISCreateStride(PETSC_COMM_WORLD, 1, rstart, 1, &isrow)); in main()
52 PetscCall(ISCreateStride(PETSC_COMM_WORLD, rend - rstart, rstart, 1, &isrow)); in main()
54 PetscCall(MatCreateSubMatrix(C, isrow, NULL, MAT_INITIAL_MATRIX, &A)); in main()
66 PetscCall(MatCreateSubMatrix(C, isrow, NULL, MAT_REUSE_MATRIX, &A)); in main()
75 PetscCall(ISDestroy(&isrow)); in main()
|
| H A D | ex166.c | 33 IS isrow, iscol; in main() local
52 …cCall(ISCreateGeneral(PETSC_COMM_WORLD, rend - rstart, ixrow + rstart, PETSC_USE_POINTER, &isrow)); in main()
54 PetscCall(ISSetPermutation(isrow)); in main()
64 PetscCall(MatPermute(A, isrow, iscol, &B)); in main()
70 PetscCall(ISView(isrow, viewer)); in main()
77 PetscCall(ISDestroy(&isrow)); in main()
|
| H A D | ex211.c | 6 PetscErrorCode ISGetSeqIS_SameColDist_Private(Mat mat, IS isrow, IS iscol, IS *isrow_d, IS *iscol_d… in ISGetSeqIS_SameColDist_Private() argument
57 PetscCall(ISGetLocalSize(isrow, &m)); in ISGetSeqIS_SameColDist_Private()
60 PetscCall(ISGetIndices(isrow, &is_idx)); in ISGetSeqIS_SameColDist_Private()
62 PetscCall(ISRestoreIndices(isrow, &is_idx)); in ISGetSeqIS_SameColDist_Private()
65 PetscCall(ISGetBlockSize(isrow, &i)); in ISGetSeqIS_SameColDist_Private()
131 IS isrow, iscol; in main() local
165 PetscCall(ISCreateStride(PETSC_COMM_WORLD, (rend - rstart) / 2, rstart, 2, &isrow)); in main()
171 PetscCall(ISGetSeqIS_SameColDist_Private(C, isrow, iscol, &isrow_d, &iscol_d, &iscol_o, &garray)); in main()
178 PetscCall(MatCreateSubMatrix(C, isrow, iscol, MAT_INITIAL_MATRIX, &A)); in main()
179 PetscCall(MatCreateSubMatrix(C, isrow, iscol, MAT_REUSE_MATRIX, &A)); in main()
[all …]
|
| H A D | ex59.c | 11 IS isrow, iscol; in main() local
50 PetscCall(ISCreateStride(PETSC_COMM_WORLD, (rend - rstart) / 2, rstart, 2, &isrow)); in main()
54 PetscCall(MatCreateSubMatrix(C, isrow, iscol, MAT_INITIAL_MATRIX, &A)); in main()
55 PetscCall(MatCreateSubMatrix(C, isrow, iscol, MAT_REUSE_MATRIX, &A)); in main()
64 PetscCall(ISDestroy(&isrow)); in main()
|
| H A D | ex87.c | 13 IS isrow, iscol, irow[n], icol[n]; in main() local
33 PetscCall(ISCreateStride(PETSC_COMM_SELF, issize, 0, 1, &isrow)); in main()
34 irow[0] = irow[1] = isrow; in main()
42 icol[0] = icol[1] = isrow; in main()
53 PetscCall(ISDestroy(&isrow)); in main()
|
| H A D | ex194.c | 11 IS isrow, iscol; in main() local
47 …ll(ISCreateStride(PETSC_COMM_WORLD, rend - rstart > 0 ? rend - rstart - 1 : 0, rstart, 1, &isrow)); in main()
49 PetscCall(MatCreateSubMatrix(C, isrow, NULL, MAT_INITIAL_MATRIX, &A)); in main()
61 PetscCall(MatCreateSubMatrix(C, isrow, NULL, MAT_REUSE_MATRIX, &A)); in main()
66 PetscCall(ISDestroy(&isrow)); in main()
|
| H A D | ex11.c | 10 IS isrow; in main() local
47 PetscCall(ISCreateStride(PETSC_COMM_SELF, (m * n) / 2, 0, 2, &isrow)); in main()
52 PetscCall(MatZeroRowsIS(C, isrow, five, 0, 0)); in main()
56 PetscCall(ISDestroy(&isrow)); in main()
|
| H A D | ex57.c | 20 IS isrow, iscol; in main() local
43 PetscCall(ISCreateStride(PETSC_COMM_SELF, m, start, 1, &isrow)); in main()
45 PetscCall(MatCreateSubMatrices(A, 1, &isrow, &iscol, MAT_INITIAL_MATRIX, &B)); in main()
59 PetscCall(ISDestroy(&isrow)); in main()
|
| H A D | ex17.c | 10 IS isrow, row, col; in main() local
51 PetscCall(ISCreateStride(PETSC_COMM_SELF, (m * n) / 2, 0, 2, &isrow)); in main()
52 PetscCall(MatZeroRowsIS(C, isrow, five, 0, 0)); in main()
78 PetscCall(ISDestroy(&isrow)); in main()
|
| /petsc/src/mat/tutorials/ |
| H A D | ex1.c | 19 IS isrow, iscol; /* row and column permutations */ in main() local
63 PetscCall(MatGetOrdering(A, rtype, &isrow, &iscol)); in main()
64 PetscCall(ISView(isrow, PETSC_VIEWER_STDOUT_WORLD)); in main()
71 PetscCall(ISDestroy(&isrow)); in main()
|
| /petsc/src/mat/interface/ftn-custom/ |
| H A D | zmatrixf.c | 35 PETSC_EXTERN void matcreatesubmatrices_(Mat *mat, PetscInt *n, IS *isrow, IS *iscol, MatReuse *scal… in matcreatesubmatrices_() argument
40 *ierr = MatCreateSubMatrices(*mat, *n, isrow, iscol, *scall, &lsmat); in matcreatesubmatrices_()
44 *ierr = MatCreateSubMatrices(*mat, *n, isrow, iscol, *scall, &lsmat); in matcreatesubmatrices_()
48 PETSC_EXTERN void matcreatesubmatricesmpi_(Mat *mat, PetscInt *n, IS *isrow, IS *iscol, MatReuse *s… in matcreatesubmatricesmpi_() argument
53 *ierr = MatCreateSubMatricesMPI(*mat, *n, isrow, iscol, *scall, &lsmat); in matcreatesubmatricesmpi_()
59 *ierr = MatCreateSubMatricesMPI(*mat, *n, isrow, iscol, *scall, &lsmat); in matcreatesubmatricesmpi_()
|
| /petsc/src/mat/impls/submat/ |
| H A D | submat.c | 4 IS isrow, iscol; /* rows and columns in submatrix, only used to check consistency */ member
152 PetscCall(ISDestroy(&Na->isrow)); in MatDestroy_SubMatrix()
185 PetscErrorCode MatCreateSubMatrixVirtual(Mat A, IS isrow, IS iscol, Mat *newmat) in MatCreateSubMatrixVirtual() argument
194 PetscValidHeaderSpecific(isrow, IS_CLASSID, 2); in MatCreateSubMatrixVirtual()
200 PetscCall(ISGetLocalSize(isrow, &m)); in MatCreateSubMatrixVirtual()
208 PetscCall(PetscObjectReference((PetscObject)isrow)); in MatCreateSubMatrixVirtual()
210 Na->isrow = isrow; in MatCreateSubMatrixVirtual()
236 PetscCall(VecScatterCreate(Na->lwork, isrow, left, NULL, &Na->lrestrict)); in MatCreateSubMatrixVirtual()
277 PetscErrorCode MatSubMatrixVirtualUpdate(Mat N, Mat A, IS isrow, IS iscol) in MatSubMatrixVirtualUpdate() argument
285 PetscValidHeaderSpecific(isrow, IS_CLASSID, 3); in MatSubMatrixVirtualUpdate()
[all …]
|
| /petsc/src/mat/impls/baij/seq/ |
| H A D | baijsolvtran1.c | 7 IS iscol = a->col, isrow = a->row; in MatSolveTranspose_SeqBAIJ_1() local
20 PetscCall(ISGetIndices(isrow, &rout)); in MatSolveTranspose_SeqBAIJ_1()
51 PetscCall(ISRestoreIndices(isrow, &rout)); in MatSolveTranspose_SeqBAIJ_1()
63 IS iscol = a->col, isrow = a->row; in MatSolveTranspose_SeqBAIJ_1_inplace() local
76 PetscCall(ISGetIndices(isrow, &rout)); in MatSolveTranspose_SeqBAIJ_1_inplace()
106 PetscCall(ISRestoreIndices(isrow, &rout)); in MatSolveTranspose_SeqBAIJ_1_inplace()
|
| H A D | baijsolvtrann.c | 7 IS iscol = a->col, isrow = a->row; in MatSolveTranspose_SeqBAIJ_N_inplace() local
21 PetscCall(ISGetIndices(isrow, &rout)); in MatSolveTranspose_SeqBAIJ_N_inplace()
61 PetscCall(ISRestoreIndices(isrow, &rout)); in MatSolveTranspose_SeqBAIJ_N_inplace()
72 IS iscol = a->col, isrow = a->row; in MatSolveTranspose_SeqBAIJ_N() local
87 PetscCall(ISGetIndices(isrow, &rout)); in MatSolveTranspose_SeqBAIJ_N()
127 PetscCall(ISRestoreIndices(isrow, &rout)); in MatSolveTranspose_SeqBAIJ_N()
|
| H A D | baijsolvtran2.c | 7 IS iscol = a->col, isrow = a->row; in MatSolveTranspose_SeqBAIJ_2_inplace() local
20 PetscCall(ISGetIndices(isrow, &rout)); in MatSolveTranspose_SeqBAIJ_2_inplace()
82 PetscCall(ISRestoreIndices(isrow, &rout)); in MatSolveTranspose_SeqBAIJ_2_inplace()
93 IS iscol = a->col, isrow = a->row; in MatSolveTranspose_SeqBAIJ_2() local
107 PetscCall(ISGetIndices(isrow, &rout)); in MatSolveTranspose_SeqBAIJ_2()
167 PetscCall(ISRestoreIndices(isrow, &rout)); in MatSolveTranspose_SeqBAIJ_2()
|
| H A D | baijsolv.c | 7 IS iscol = a->col, isrow = a->row; in MatSolve_SeqBAIJ_N_inplace() local
22 PetscCall(ISGetIndices(isrow, &rout)); in MatSolve_SeqBAIJ_N_inplace()
55 PetscCall(ISRestoreIndices(isrow, &rout)); in MatSolve_SeqBAIJ_N_inplace()
66 IS iscol = a->col, isrow = a->row; in MatSolve_SeqBAIJ_7_inplace() local
79 PetscCall(ISGetIndices(isrow, &rout)); in MatSolve_SeqBAIJ_7_inplace()
175 PetscCall(ISRestoreIndices(isrow, &rout)); in MatSolve_SeqBAIJ_7_inplace()
186 IS iscol = a->col, isrow = a->row; in MatSolve_SeqBAIJ_7() local
199 PetscCall(ISGetIndices(isrow, &rout)); in MatSolve_SeqBAIJ_7()
294 PetscCall(ISRestoreIndices(isrow, &rout)); in MatSolve_SeqBAIJ_7()
305 IS iscol = a->col, isrow = a->row; in MatSolve_SeqBAIJ_6_inplace() local
[all …]
|
| H A D | baijfact2.c | 118 IS isrow = b->row, isicol = b->icol; in MatLUFactorNumeric_SeqBAIJ_N() local
129 PetscCall(ISGetIndices(isrow, &r)); in MatLUFactorNumeric_SeqBAIJ_N()
204 PetscCall(ISRestoreIndices(isrow, &r)); in MatLUFactorNumeric_SeqBAIJ_N()
206 PetscCall(ISIdentity(isrow, &row_identity)); in MatLUFactorNumeric_SeqBAIJ_N()
252 static PetscErrorCode MatILUFactorSymbolic_SeqBAIJ_ilu0(Mat fact, Mat A, IS isrow, IS iscol, const … in MatILUFactorSymbolic_SeqBAIJ_ilu0() argument
311 PetscErrorCode MatILUFactorSymbolic_SeqBAIJ(Mat fact, Mat A, IS isrow, IS iscol, const MatFactorInf… in MatILUFactorSymbolic_SeqBAIJ() argument
344 PetscCall(ISIdentity(isrow, &row_identity)); in MatILUFactorSymbolic_SeqBAIJ()
351 PetscCall(MatILUFactorSymbolic_SeqBAIJ_ilu0(fact, A, isrow, iscol, info)); in MatILUFactorSymbolic_SeqBAIJ()
360 b->row = isrow; in MatILUFactorSymbolic_SeqBAIJ()
363 PetscCall(PetscObjectReference((PetscObject)isrow)); in MatILUFactorSymbolic_SeqBAIJ()
[all …]
|
| H A D | baijsolvtran3.c | 7 IS iscol = a->col, isrow = a->row; in MatSolveTranspose_SeqBAIJ_3_inplace() local
20 PetscCall(ISGetIndices(isrow, &rout)); in MatSolveTranspose_SeqBAIJ_3_inplace()
90 PetscCall(ISRestoreIndices(isrow, &rout)); in MatSolveTranspose_SeqBAIJ_3_inplace()
101 IS iscol = a->col, isrow = a->row; in MatSolveTranspose_SeqBAIJ_3() local
115 PetscCall(ISGetIndices(isrow, &rout)); in MatSolveTranspose_SeqBAIJ_3()
183 PetscCall(ISRestoreIndices(isrow, &rout)); in MatSolveTranspose_SeqBAIJ_3()
|
| /petsc/src/ksp/ksp/tests/ |
| H A D | ex16f.F90 | 30 IS isrow
73 PetscCallA(ISCreateStride(PETSC_COMM_SELF, ifive, izero, ione, isrow, ierr))
74 PetscCallA(MatZeroRowsIS(A, isrow, five, PETSC_NULL_VEC, PETSC_NULL_VEC, ierr))
98 PetscCallA(ISDestroy(isrow, ierr))
|
| /petsc/src/mat/impls/localref/ |
| H A D | mlocalref.c | 215 PetscErrorCode MatCreateLocalRef(Mat A, IS isrow, IS iscol, Mat *newmat) in MatCreateLocalRef() argument
224 PetscValidHeaderSpecific(isrow, IS_CLASSID, 2); in MatCreateLocalRef()
231 PetscCall(ISGetLocalSize(isrow, &m)); in MatCreateLocalRef()
262 PetscCall(ISL2GCompose(isrow, A->rmap->mapping, &rltog)); in MatCreateLocalRef()
263 if (isrow == iscol && A->rmap->mapping == A->cmap->mapping) { in MatCreateLocalRef()
271 PetscCall(PetscObjectTypeCompare((PetscObject)isrow, ISBLOCK, &lr->rowisblock)); in MatCreateLocalRef()
278 PetscCall(ISGetBlockSize(isrow, &rbs)); in MatCreateLocalRef()
290 PetscCall(ISL2GComposeBlock(isrow, A->rmap->mapping, &rltog)); in MatCreateLocalRef()
291 if (isrow == iscol && A->rmap->mapping == A->cmap->mapping) { in MatCreateLocalRef()
|
| /petsc/src/mat/impls/aij/seq/ |
| H A D | aijfact.c | 54 PetscErrorCode MatLUFactorSymbolic_SeqAIJ(Mat B, Mat A, IS isrow, IS iscol, const MatFactorInfo *in… in MatLUFactorSymbolic_SeqAIJ() argument
74 PetscCall(ISGetIndices(isrow, &r)); in MatLUFactorSymbolic_SeqAIJ()
141 PetscCall(ISRestoreIndices(isrow, &r)); in MatLUFactorSymbolic_SeqAIJ()
161 b->row = isrow; in MatLUFactorSymbolic_SeqAIJ()
163 PetscCall(PetscObjectReference((PetscObject)isrow)); in MatLUFactorSymbolic_SeqAIJ()
222 IS isrow = b->row, isicol = b->icol; in MatLUFactorNumeric_SeqAIJ() local
260 PetscCall(ISGetIndices(isrow, &r)); in MatLUFactorNumeric_SeqAIJ()
359 PetscCall(ISRestoreIndices(isrow, &r)); in MatLUFactorNumeric_SeqAIJ()
361 PetscCall(ISIdentity(isrow, &row_identity)); in MatLUFactorNumeric_SeqAIJ()
401 IS isrow = b->row, isicol = b->icol; in MatLUFactorNumeric_SeqAIJ_inplace() local
[all …]
|