| /petsc/src/mat/tests/ |
| H A D | ex77.c | 9 Mat sA; /* symmetric part of the matrices */ in main() local
30 PetscCall(MatCreateSeqSBAIJ(PETSC_COMM_WORLD, bs, n, n, nz, NULL, &sA)); in main()
31 PetscCall(MatSetOption(sA, MAT_NEW_NONZERO_ALLOCATION_ERR, PETSC_FALSE)); in main()
35 PetscCall(MatGetOwnershipRange(sA, &i, &j)); in main()
50 PetscCall(MatSetValues(sA, 1, &i, 3, col, value, INSERT_VALUES)); in main()
61 PetscCall(MatSetValues(sA, 1, &i, 3, col, value, INSERT_VALUES)); in main()
72 PetscCall(MatSetValues(sA, 1, &i, 3, col, value, INSERT_VALUES)); in main()
82 PetscCall(MatSetValues(sA, 1, &Ii, 1, &J, &neg_one, INSERT_VALUES)); in main()
87 PetscCall(MatSetValues(sA, 1, &Ii, 1, &J, &neg_one, INSERT_VALUES)); in main()
92 PetscCall(MatSetValues(sA, 1, &Ii, 1, &J, &neg_one, INSERT_VALUES)); in main()
[all …]
|
| H A D | ex74.c | 10 Mat sA, sB, sFactor, B, C; /* symmetric matrices */ in main() local
35 PetscCall(MatCreate(PETSC_COMM_SELF, &sA)); in main()
36 PetscCall(MatSetSizes(sA, n, n, PETSC_DETERMINE, PETSC_DETERMINE)); in main()
37 PetscCall(MatSetType(sA, MATSEQSBAIJ)); in main()
38 PetscCall(MatSetFromOptions(sA)); in main()
39 PetscCall(MatGetType(sA, &type)); in main()
40 PetscCall(PetscObjectTypeCompare((PetscObject)sA, MATSEQSBAIJ, &doIcc)); in main()
41 PetscCall(MatSeqSBAIJSetPreallocation(sA, bs, nz, NULL)); in main()
42 PetscCall(MatSetOption(sA, MAT_IGNORE_LOWER_TRIANGULAR, PETSC_TRUE)); in main()
46 PetscCall(MatGetOwnershipRange(sA, &i, &j)); in main()
[all …]
|
| H A D | ex75.c | 8 Mat A, sA, sB; in main() local
130 PetscCall(MatConvert(A, MATSBAIJ, MAT_INITIAL_MATRIX, &sA)); in main()
133 PetscCall(MatGetSize(sA, &i, &j)); in main()
142 PetscCall(MatGetLocalSize(sA, &i, &j)); in main()
169 PetscCall(MatNorm(sA, NORM_FROBENIUS, &r2)); in main()
173 PetscCall(MatNorm(sA, NORM_INFINITY, &r2)); in main()
177 PetscCall(MatNorm(sA, NORM_1, &r2)); in main()
182 PetscCall(MatGetOwnershipRange(sA, &rstart, &rend)); in main()
193 PetscCall(MatDiagonalScale(sA, x, x)); in main()
194 PetscCall(MatMultEqual(A, sA, 10, &flg)); in main()
[all …]
|
| H A D | ex128.c | 14 Mat C, sC, sA; in main() local
83 PetscCall(MatGetFactor(sC, MATSOLVERPETSC, MAT_FACTOR_CHOLESKY, &sA)); in main()
84 PetscCall(MatCholeskyFactorSymbolic(sA, sC, row, &info)); in main()
89 PetscCall(MatGetFactor(sC, MATSOLVERPETSC, MAT_FACTOR_ICC, &sA)); in main()
90 PetscCall(MatICCFactorSymbolic(sA, sC, row, &info)); in main()
92 PetscCall(MatCholeskyFactorNumeric(sA, sC, &info)); in main()
99 PetscCall(MatForwardSolve(sA, b, ytmp)); in main()
101 PetscCall(MatBackwardSolve(sA, ytmp, y)); in main()
108 PetscCall(MatSolve(sA, b, y)); in main()
110 PetscCall(MatDestroy(&sA)); in main()
|
| H A D | ex92.c | 9 Mat A, Atrans, sA, *submatA, *submatsA; in main() local
95 PetscCall(MatConvert(A, MATSBAIJ, MAT_INITIAL_MATRIX, &sA)); in main()
100 PetscCall(MatView(sA, PETSC_VIEWER_STDOUT_WORLD)); in main()
104 PetscCall(MatMultEqual(A, sA, 10, &flg)); in main()
147 PetscCall(MatIncreaseOverlap(sA, nd, is2, ov)); in main()
179 PetscCall(MatCreateSubMatrices(sA, nd, is1, is1, MAT_INITIAL_MATRIX, &submatsA)); in main()
181 PetscCall(MatMultEqual(A, sA, 10, &flg)); in main()
186 PetscCall(MatCreateSubMatrices(sA, nd, is1, is1, MAT_REUSE_MATRIX, &submatsA)); in main()
187 PetscCall(MatMultEqual(A, sA, 10, &flg)); in main()
206 PetscCall(MatDestroy(&sA)); in main()
|
| H A D | ex76.c | 9 Mat sA, sC; /* symmetric part of the matrices */ in main() local
154 PetscCall(MatConvert(A, MATSEQSBAIJ, MAT_INITIAL_MATRIX, &sA)); in main()
155 PetscCall(MatMultEqual(A, sA, 20, &equal)); in main()
160 PetscCall(MatGetOwnershipRange(sA, &i, &j)); in main()
258 PetscCall(MatGetFactor(sA, MATSOLVERPETSC, MAT_FACTOR_CHOLESKY, &sC)); in main()
259 PetscCall(MatCholeskyFactorSymbolic(sC, sA, perm, &factinfo)); in main()
264 PetscCall(MatGetFactor(sA, MATSOLVERPETSC, MAT_FACTOR_ICC, &sC)); in main()
265 PetscCall(MatICCFactorSymbolic(sC, sA, perm, &factinfo)); in main()
267 PetscCall(MatCholeskyFactorNumeric(sC, sA, &factinfo)); in main()
280 PetscCall(MatMult(sA, x, b)); in main()
[all …]
|
| H A D | ex91.c | 7 Mat A, Atrans, sA, *submatA, *submatsA; in main() local
80 PetscCall(MatConvert(A, MATSEQSBAIJ, MAT_INITIAL_MATRIX, &sA)); in main()
91 PetscCall(MatMult(sA, xx, s2)); in main()
124 PetscCall(MatIncreaseOverlap(sA, nd, is2, ov)); in main()
137 PetscCall(MatCreateSubMatrices(sA, nd, is2, is2, MAT_INITIAL_MATRIX, &submatsA)); in main()
161 PetscCall(MatCreateSubMatrices(sA, nd, is2, is2, MAT_REUSE_MATRIX, &submatsA)); in main()
198 PetscCall(MatDestroy(&sA)); in main()
|
| H A D | ex221.c | 341 Mat sA, sS, dA, dS, At, St; in main() local
349 PetscCall(MatCreateSubMatrix(A, r, c, MAT_INITIAL_MATRIX, &sA)); in main()
351 PetscCall(MatMultAddEqual(sA, sS, 10, &flg)); in main()
353 PetscCall(MatMultTransposeAddEqual(sA, sS, 10, &flg)); in main()
355 PetscCall(MatConvert(sA, MATDENSE, MAT_INITIAL_MATRIX, &dA)); in main()
360 PetscCall(MatDestroy(&sA)); in main()
366 PetscCall(MatCreateSubMatrix(At, c, r, MAT_INITIAL_MATRIX, &sA)); in main()
368 PetscCall(MatMultAddEqual(sA, sS, 10, &flg)); in main()
370 PetscCall(MatMultTransposeAddEqual(sA, sS, 10, &flg)); in main()
372 PetscCall(MatConvert(sA, MATDENSE, MAT_INITIAL_MATRIX, &dA)); in main()
[all …]
|
| /petsc/src/ksp/ksp/tests/ |
| H A D | ex9.c | 16 Mat sA, nA; in replace_submats() local
19 PetscCall(MatCreateSubMatrix(A, r[i], c[j], MAT_INITIAL_MATRIX, &sA)); in replace_submats()
20 PetscCall(MatDuplicate(sA, MAT_COPY_VALUES, &nA)); in replace_submats()
21 PetscCall(MatGetOptionsPrefix(sA, &prefix)); in replace_submats()
25 PetscCall(MatDestroy(&sA)); in replace_submats()
37 Mat M, A, P, sA[2][2], sP[2][2]; in main() local
53 PetscCall(MatCreateSubMatrix(M, f[i], f[j], MAT_INITIAL_MATRIX, &sA[i][j])); in main()
57 PetscCall(MatCreateNest(PetscObjectComm((PetscObject)M), 2, f, 2, f, &sA[0][0], &A)); in main()
83 PetscCall(MatDestroy(&sA[i][j])); in main()
|
| H A D | ex32.c | 67 Mat sA; in main() local
72 PetscCall(MatConvert(A, MATSBAIJ, MAT_INITIAL_MATRIX, &sA)); in main()
74 A = sA; in main()
|
| /petsc/src/ksp/pc/impls/mpi/ |
| H A D | pcmpi.c | 98 Mat sA; in PCMPICreate() local
108 PetscCall(PCGetOperators(pc, &sA, &sA)); in PCMPICreate()
109 PetscCall(MatGetSize(sA, &N[0], &N[1])); in PCMPICreate()
161 Mat sA; in PCMPISetMat() local
186 PetscCall(PCGetOperators(pc, &sA, &sA)); in PCMPISetMat()
187 PetscCall(MatGetSize(sA, &matproperties[0], &matproperties[1])); in PCMPISetMat()
188 PetscCall(MatGetBlockSize(sA, &bs)); in PCMPISetMat()
190 PetscCall(MatIsSymmetricKnown(sA, &isset, &issymmetric)); in PCMPISetMat()
192 PetscCall(MatIsHermitianKnown(sA, &isset, &ishermitian)); in PCMPISetMat()
194 PetscCall(MatIsSPDKnown(sA, &isset, &isspd)); in PCMPISetMat()
[all …]
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexpreallocate.c | 239 …PetscBool useCone, PetscBool useClosure, PetscBool useAnchors, PetscSection *sA, PetscInt **colIdx) in DMPlexCreateAdjacencySection_Static() argument
696 *sA = sectionAdj; in DMPlexCreateAdjacencySection_Static()
|
| H A D | plextree.c | 84 PetscInt size, i, sA = -1, sB, sOrientB, sConeSize; in DMPlexReferenceTreeGetChildSymmetry_Default() local
94 sA = supp[i]; in DMPlexReferenceTreeGetChildSymmetry_Default()
95 if (sA == parent) continue; in DMPlexReferenceTreeGetChildSymmetry_Default()
96 PetscCall(DMPlexGetTreeParent(dm, sA, &sParent, NULL)); in DMPlexReferenceTreeGetChildSymmetry_Default()
102 …PetscCall(DMPlexReferenceTreeGetChildSymmetry_Default(dm, parent, parentOrientA, 0, sA, parentOrie… in DMPlexReferenceTreeGetChildSymmetry_Default()
103 PetscCall(DMPlexGetConeSize(dm, sA, &sConeSize)); in DMPlexReferenceTreeGetChildSymmetry_Default()
104 PetscCall(DMPlexGetCone(dm, sA, &coneA)); in DMPlexReferenceTreeGetChildSymmetry_Default()
106 PetscCall(DMPlexGetConeOrientation(dm, sA, &oA)); in DMPlexReferenceTreeGetChildSymmetry_Default()
|
| /petsc/src/dm/impls/forest/p4est/ |
| H A D | pforest.h | 1943 PetscInt size, i, sA = -1, sB, sOrientB, sConeSize; in DMReferenceTreeGetChildSymmetry_pforest() local
1953 sA = supp[i]; in DMReferenceTreeGetChildSymmetry_pforest()
1954 if (sA == parent) continue; in DMReferenceTreeGetChildSymmetry_pforest()
1955 PetscCall(DMPlexGetTreeParent(dm, sA, &sParent, NULL)); in DMReferenceTreeGetChildSymmetry_pforest()
1961 …PetscCall(DMReferenceTreeGetChildSymmetry_pforest(dm, parent, parentOrientA, 0, sA, parentOrientB,… in DMReferenceTreeGetChildSymmetry_pforest()
1962 PetscCall(DMPlexGetConeSize(dm, sA, &sConeSize)); in DMReferenceTreeGetChildSymmetry_pforest()
1963 PetscCall(DMPlexGetCone(dm, sA, &coneA)); in DMReferenceTreeGetChildSymmetry_pforest()
1965 PetscCall(DMPlexGetConeOrientation(dm, sA, &oA)); in DMReferenceTreeGetChildSymmetry_pforest()
|
| /petsc/src/mat/impls/is/ |
| H A D | matis.c | 3541 Mat sA = (*submat)[i]; in MatCreateSubMatrices_IS()
3543 PetscCall(MatConvert(sA,ltype,MAT_INPLACE_MATRIX,&sA)); in MatCreateSubMatrices_IS()
3544 (*submat)[i] = sA; in MatCreateSubMatrices_IS()
|
| /petsc/src/ksp/pc/impls/bddc/ |
| H A D | bddcprivate.c | 9796 static PetscErrorCode MatAIJExtractRows(Mat A, IS rows, Mat *sA) in MatAIJExtractRows() argument
9807 PetscCall(MatCreateSubMatrix(A, rows, NULL, MAT_INITIAL_MATRIX, sA)); in MatAIJExtractRows()
9861 …scCall(MatCreateMPIAIJWithArrays(comm, ni, A->cmap->n, PETSC_DECIDE, A->cmap->N, di, dj, daa, sA)); in MatAIJExtractRows()
|