| /petsc/include/petsc/private/kernels/ |
| H A D | blockinvert.h | 27 #define PetscKernel_A_gets_inverse_A_4_nopivot(mat) \ argument
28 …mat[0]; mat[0] = d = 1.0 / di; mat[4] *= -d; mat[8] *= -d; mat[12] *= -d; mat[1] *= d; mat[2] *= d…
|
| /petsc/src/mat/interface/ |
| H A D | matrix.c | 165 PetscErrorCode MatFactorGetErrorZeroPivot(Mat mat, PetscReal *pivot, PetscInt *row) in MatFactorGetErrorZeroPivot() argument
168 PetscValidHeaderSpecific(mat, MAT_CLASSID, 1); in MatFactorGetErrorZeroPivot()
171 *pivot = mat->factorerror_zeropivot_value; in MatFactorGetErrorZeroPivot()
172 *row = mat->factorerror_zeropivot_row; in MatFactorGetErrorZeroPivot()
195 PetscErrorCode MatFactorGetError(Mat mat, MatFactorError *err) in MatFactorGetError() argument
198 PetscValidHeaderSpecific(mat, MAT_CLASSID, 1); in MatFactorGetError()
200 *err = mat->factorerrortype; in MatFactorGetError()
220 PetscErrorCode MatFactorClearError(Mat mat) in MatFactorClearError() argument
223 PetscValidHeaderSpecific(mat, MAT_CLASSID, 1); in MatFactorClearError()
224 mat->factorerrortype = MAT_FACTOR_NOERROR; in MatFactorClearError()
[all …]
|
| H A D | matproduct.c | 141 static PetscErrorCode MatProductNumeric_ABC_Unsafe(Mat mat) in MatProductNumeric_ABC_Unsafe() argument
143 Mat_Product *product = mat->product; in MatProductNumeric_ABC_Unsafe()
151 PetscCall((*mat->ops->matmultnumeric)(A, BC, mat)); in MatProductNumeric_ABC_Unsafe()
156 static PetscErrorCode MatProductSymbolic_ABC_Unsafe(Mat mat) in MatProductSymbolic_ABC_Unsafe() argument
158 Mat_Product *product = mat->product; in MatProductSymbolic_ABC_Unsafe()
163 …PetscCall(PetscInfo(mat, "for A %s, B %s, C %s is used\n", ((PetscObject)product->A)->type_name, (… in MatProductSymbolic_ABC_Unsafe()
173 PetscCall(MatProductSetType(mat, MATPRODUCT_AB)); in MatProductSymbolic_ABC_Unsafe()
174 PetscCall(MatProductSetAlgorithm(mat, MATPRODUCTALGORITHMDEFAULT)); in MatProductSymbolic_ABC_Unsafe()
177 PetscCall(MatProductSetFromOptions(mat)); in MatProductSymbolic_ABC_Unsafe()
178 PetscCall(MatProductSymbolic(mat)); in MatProductSymbolic_ABC_Unsafe()
[all …]
|
| H A D | matreg.c | 27 PetscErrorCode MatGetRootType_Private(Mat mat, MatType *rootType) in MatGetRootType_Private() argument
34 PetscValidHeaderSpecific(mat, MAT_CLASSID, 1); in MatGetRootType_Private()
35 PetscCall(MatGetType(mat, &inType)); in MatGetRootType_Private()
64 PetscErrorCode MatGetMPIMatType_Private(Mat mat, MatType *MPIType) in MatGetMPIMatType_Private() argument
71 PetscValidHeaderSpecific(mat, MAT_CLASSID, 1); in MatGetMPIMatType_Private()
72 PetscCall(MatGetType(mat, &inType)); in MatGetMPIMatType_Private()
107 PetscErrorCode MatSetType(Mat mat, MatType matype) in MatSetType() argument
114 PetscValidHeaderSpecific(mat, MAT_CLASSID, 1); in MatSetType()
117 PetscCall(PetscObjectTypeCompare((PetscObject)mat, matype, &sametype)); in MatSetType()
124 …if (((PetscObject)mat)->type_name) PetscCall(PetscStrbeginswith(((PetscObject)mat)->type_name, "mp… in MatSetType()
[all …]
|
| /petsc/src/mat/impls/shell/ftn-custom/ |
| H A D | zshellf.c | 46 static PetscErrorCode ourmult(Mat mat, Vec x, Vec y) in ourmult() argument
48 …, Vec *, Vec *, PetscErrorCode *))(((PetscObject)mat)->fortran_func_pointers[FORTRAN_MATOP_MULT]))… in ourmult()
52 static PetscErrorCode ourmultadd(Mat mat, Vec x, Vec y, Vec z) in ourmultadd() argument
54 …Vec *, Vec *, PetscErrorCode *))(((PetscObject)mat)->fortran_func_pointers[FORTRAN_MATOP_MULT_ADD]… in ourmultadd()
58 static PetscErrorCode ourmulttranspose(Mat mat, Vec x, Vec y) in ourmulttranspose() argument
60 … *, Vec *, PetscErrorCode *))(((PetscObject)mat)->fortran_func_pointers[FORTRAN_MATOP_MULT_TRANSPO… in ourmulttranspose()
64 static PetscErrorCode ourmulthermitiantranspose(Mat mat, Vec x, Vec y) in ourmulthermitiantranspose() argument
66 … Vec *, Vec *, PetscErrorCode *))(((PetscObject)mat)->fortran_func_pointers[FORTRAN_MATOP_MULT_HT]… in ourmulthermitiantranspose()
70 static PetscErrorCode ourmulttransposeadd(Mat mat, Vec x, Vec y, Vec z) in ourmulttransposeadd() argument
72 …, Vec *, PetscErrorCode *))(((PetscObject)mat)->fortran_func_pointers[FORTRAN_MATOP_MULT_TRANSPOSE… in ourmulttransposeadd()
[all …]
|
| /petsc/src/mat/impls/dense/seq/hdf5/ |
| H A D | densehdf5.c | 10 PetscErrorCode MatLoad_Dense_HDF5(Mat mat, PetscViewer viewer) in MatLoad_Dense_HDF5() argument
30 …SETERRQ(PetscObjectComm((PetscObject)mat), PETSC_ERR_SUP, "PetscViewerFormat %s not supported for … in MatLoad_Dense_HDF5()
36 …PetscCheck(((PetscObject)mat)->name, PetscObjectComm((PetscObject)mat), PETSC_ERR_SUP, "Mat name m… in MatLoad_Dense_HDF5()
47 PetscCall(PetscObjectGetComm((PetscObject)mat, &comm)); in MatLoad_Dense_HDF5()
50 PetscCall(PetscObjectGetName((PetscObject)mat, &mat_name)); in MatLoad_Dense_HDF5()
53 PetscCall(PetscLayoutCreate(PetscObjectComm((PetscObject)mat), &vmap)); in MatLoad_Dense_HDF5()
54 if (mat->rmap->n >= 0 && mat->cmap->N < 0) { in MatLoad_Dense_HDF5()
56 PetscCall(PetscViewerHDF5ReadSizes(viewer, mat_name, &mat->cmap->N, NULL)); in MatLoad_Dense_HDF5()
58 vmap->bs = mat->cmap->N; in MatLoad_Dense_HDF5()
59 vmap->n = (mat->rmap->n < 0 || mat->cmap->N < 0) ? -1 : mat->rmap->n * mat->cmap->N; in MatLoad_Dense_HDF5()
[all …]
|
| /petsc/src/mat/impls/composite/ |
| H A D | mcomposite.c | 7 Mat mat; member
31 static PetscErrorCode MatDestroy_Composite(Mat mat) in MatDestroy_Composite() argument
38 PetscCall(MatShellGetContext(mat, &shell)); in MatDestroy_Composite()
41 PetscCall(MatDestroy(&next->mat)); in MatDestroy_Composite()
58 PetscCall(PetscObjectComposeFunction((PetscObject)mat, "MatCompositeAddMat_C", NULL)); in MatDestroy_Composite()
59 PetscCall(PetscObjectComposeFunction((PetscObject)mat, "MatCompositeSetType_C", NULL)); in MatDestroy_Composite()
60 PetscCall(PetscObjectComposeFunction((PetscObject)mat, "MatCompositeGetType_C", NULL)); in MatDestroy_Composite()
61 PetscCall(PetscObjectComposeFunction((PetscObject)mat, "MatCompositeSetMergeType_C", NULL)); in MatDestroy_Composite()
62 PetscCall(PetscObjectComposeFunction((PetscObject)mat, "MatCompositeSetMatStructure_C", NULL)); in MatDestroy_Composite()
63 PetscCall(PetscObjectComposeFunction((PetscObject)mat, "MatCompositeGetMatStructure_C", NULL)); in MatDestroy_Composite()
[all …]
|
| /petsc/src/mat/tutorials/ |
| H A D | ex8.c | 29 PetscErrorCode RegisterMatScaleUserImpl(Mat mat) in RegisterMatScaleUserImpl() argument
34 PetscCallMPI(MPI_Comm_size(PetscObjectComm((PetscObject)mat), &size)); in RegisterMatScaleUserImpl()
36 …PetscCall(PetscObjectComposeFunction((PetscObject)mat, "MatScaleUserImpl_C", MatScaleUserImpl_SeqA… in RegisterMatScaleUserImpl()
39 PetscCall(MatMPIAIJGetSeqAIJ(mat, &AA, &AB, NULL)); in RegisterMatScaleUserImpl()
40 …PetscCall(PetscObjectComposeFunction((PetscObject)mat, "MatScaleUserImpl_C", MatScaleUserImpl_MPIA… in RegisterMatScaleUserImpl()
50 PetscErrorCode DeRegisterMatScaleUserImpl(Mat mat) in DeRegisterMatScaleUserImpl() argument
55 PetscCallMPI(MPI_Comm_size(PetscObjectComm((PetscObject)mat), &size)); in DeRegisterMatScaleUserImpl()
57 PetscCall(PetscObjectComposeFunction((PetscObject)mat, "MatScaleUserImpl_C", NULL)); in DeRegisterMatScaleUserImpl()
60 PetscCall(MatMPIAIJGetSeqAIJ(mat, &AA, &AB, NULL)); in DeRegisterMatScaleUserImpl()
61 PetscCall(PetscObjectComposeFunction((PetscObject)mat, "MatScaleUserImpl_C", NULL)); in DeRegisterMatScaleUserImpl()
[all …]
|
| /petsc/src/mat/tests/ |
| H A D | ex156.c | 5 PetscErrorCode Assemble(Mat mat) in Assemble() argument
20 PetscCall(MatSetValues(mat, 4, idx, 4, idx, vals, ADD_VALUES)); in Assemble()
26 PetscCall(MatSetValues(mat, 4, idx, 4, idx, vals, ADD_VALUES)); in Assemble()
33 PetscCall(MatSetValues(mat, 4, idx, 4, idx, vals, ADD_VALUES)); in Assemble()
39 PetscCall(MatSetValues(mat, 4, idx, 4, idx, vals, ADD_VALUES)); in Assemble()
41 PetscCall(MatAssemblyBegin(mat, MAT_FINAL_ASSEMBLY)); in Assemble()
42 PetscCall(MatAssemblyEnd(mat, MAT_FINAL_ASSEMBLY)); in Assemble()
48 Mat mat; in main() local
54 PetscCall(MatCreate(PETSC_COMM_WORLD, &mat)); in main()
55 if (rank == 0) PetscCall(MatSetSizes(mat, 6, 6, PETSC_DETERMINE, PETSC_DETERMINE)); in main()
[all …]
|
| H A D | ex2.c | 5 static PetscErrorCode TransposeAXPY(Mat C, PetscScalar alpha, Mat mat, PetscErrorCode (*f)(Mat, Mat… in TransposeAXPY() argument
11 PetscCall(MatGetType(mat, &mtype)); in TransposeAXPY()
17 PetscCall(MatDuplicate(mat, MAT_COPY_VALUES, &C)); in TransposeAXPY()
20 PetscCall(MatAXPY(E, alpha, mat, SAME_NONZERO_PATTERN)); in TransposeAXPY()
32 PetscCall(MatTranspose(mat, MAT_INITIAL_MATRIX, &D)); in TransposeAXPY()
34 PetscCall(MatHermitianTranspose(mat, MAT_INITIAL_MATRIX, &D)); in TransposeAXPY()
37 PetscCall(MatDuplicate(mat, MAT_COPY_VALUES, &C)); in TransposeAXPY()
48 PetscCall(MatDuplicate(mat, MAT_COPY_VALUES, &C)); in TransposeAXPY()
51 PetscCall(f(mat, &F)); in TransposeAXPY()
64 PetscCall(MatConvert(mat, MATSHELL, MAT_INITIAL_MATRIX, &C)); in TransposeAXPY()
[all …]
|
| H A D | ex49.c | 7 Mat mat, tmat = 0; in main() local
34 PetscCall(MatCreate(PETSC_COMM_WORLD, &mat)); in main()
35 PetscCall(MatSetSizes(mat, PETSC_DECIDE, PETSC_DECIDE, m, n)); in main()
36 PetscCall(MatSetFromOptions(mat)); in main()
37 PetscCall(MatSetUp(mat)); in main()
38 PetscCall(MatGetOwnershipRange(mat, &rstart, &rend)); in main()
42 PetscCall(MatSetValues(mat, 1, &i, 1, &j, &v, INSERT_VALUES)); in main()
45 PetscCall(MatAssemblyBegin(mat, MAT_FINAL_ASSEMBLY)); in main()
46 PetscCall(MatAssemblyEnd(mat, MAT_FINAL_ASSEMBLY)); in main()
49 PetscCall(MatGetInfo(mat, MAT_GLOBAL_SUM, &info)); in main()
[all …]
|
| H A D | ex117.c | 10 Mat mat, fact, B; in main() local
26 PetscCall(MatCreateSeqSBAIJ(PETSC_COMM_SELF, 2, 6, 6, 0, nnz, &mat)); in main()
34 PetscCall(MatSetValues(mat, 2, ind1, 2, ind1, temp, INSERT_VALUES)); in main()
41 PetscCall(MatSetValues(mat, 2, ind1, 2, ind2, temp, INSERT_VALUES)); in main()
48 PetscCall(MatSetValues(mat, 2, ind1, 2, ind1, temp, INSERT_VALUES)); in main()
55 PetscCall(MatSetValues(mat, 2, ind1, 2, ind1, temp, INSERT_VALUES)); in main()
57 PetscCall(MatAssemblyBegin(mat, MAT_FINAL_ASSEMBLY)); in main()
58 PetscCall(MatAssemblyEnd(mat, MAT_FINAL_ASSEMBLY)); in main()
60 PetscCall(MatDuplicate(mat, MAT_SHARE_NONZERO_PATTERN, &B)); in main()
67 PetscCall(MatSetValues(mat, 2, ind1, 2, ind1, temp, INSERT_VALUES)); in main()
[all …]
|
| H A D | ex68.c | 7 Mat mat, B, C; in main() local
19 PetscCall(MatCreate(PETSC_COMM_WORLD, &mat)); in main()
20 PetscCall(MatSetSizes(mat, PETSC_DECIDE, PETSC_DECIDE, 4, 4)); in main()
21 PetscCall(MatSetFromOptions(mat)); in main()
22 PetscCall(MatSetUp(mat)); in main()
28 PetscCall(MatSetValues(mat, 1, &i, 1, &j, &v, INSERT_VALUES)); in main()
32 PetscCall(MatSetValues(mat, 1, &i, 1, &j, &v, INSERT_VALUES)); in main()
36 PetscCall(MatSetValues(mat, 1, &i, 1, &j, &v, INSERT_VALUES)); in main()
40 PetscCall(MatSetValues(mat, 1, &i, 1, &j, &v, INSERT_VALUES)); in main()
41 PetscCall(MatAssemblyBegin(mat, MAT_FINAL_ASSEMBLY)); in main()
[all …]
|
| H A D | ex1.c | 10 Mat mat, RHS, SOLU; in createMatsAndVecs() local
40 PetscCall(MatCreate(PETSC_COMM_WORLD, &mat)); in createMatsAndVecs()
41 PetscCall(MatSetSizes(mat, PETSC_DECIDE, PETSC_DECIDE, m, n)); in createMatsAndVecs()
42 PetscCall(MatSetType(mat, MATDENSE)); in createMatsAndVecs()
43 PetscCall(MatSetFromOptions(mat)); in createMatsAndVecs()
44 PetscCall(MatSetUp(mat)); in createMatsAndVecs()
45 PetscCall(MatGetOwnershipRange(mat, &rstart, &rend)); in createMatsAndVecs()
46 PetscCall(MatGetOwnershipRangeColumn(mat, &cstart, &cend)); in createMatsAndVecs()
51 PetscCall(MatSetValues(mat, 1, &i, 1, &i, &value, INSERT_VALUES)); in createMatsAndVecs()
55 PetscCall(MatSetValues(mat, 1, &i, 1, &j, &value, INSERT_VALUES)); in createMatsAndVecs()
[all …]
|
| H A D | ex107.c | 7 Mat mat; in main() local
24 PetscCall(MatCreate(PETSC_COMM_WORLD, &mat)); in main()
25 PetscCall(MatSetSizes(mat, PETSC_DECIDE, PETSC_DECIDE, m, n)); in main()
26 PetscCall(MatSetFromOptions(mat)); in main()
27 if (struct_only) PetscCall(MatSetOption(mat, MAT_STRUCTURE_ONLY, PETSC_TRUE)); in main()
28 PetscCall(MatSetUp(mat)); in main()
29 PetscCall(MatGetOwnershipRange(mat, &rstart, &rend)); in main()
33 PetscCall(MatSetValues(mat, 1, &i, 1, &j, &v, INSERT_VALUES)); in main()
36 PetscCall(MatAssemblyBegin(mat, MAT_FINAL_ASSEMBLY)); in main()
37 PetscCall(MatAssemblyEnd(mat, MAT_FINAL_ASSEMBLY)); in main()
[all …]
|
| H A D | ex4.c | 7 Mat mat, submat, submat1, *submatrices; in main() local
19 PetscCall(MatCreate(PETSC_COMM_WORLD, &mat)); in main()
20 PetscCall(MatSetSizes(mat, PETSC_DECIDE, PETSC_DECIDE, m, n)); in main()
21 PetscCall(MatSetFromOptions(mat)); in main()
22 PetscCall(MatSetUp(mat)); in main()
23 PetscCall(MatGetOwnershipRange(mat, &rstart, &rend)); in main()
27 PetscCall(MatSetValues(mat, 1, &tmp, 1, &i, &value, INSERT_VALUES)); in main()
29 PetscCall(MatAssemblyBegin(mat, MAT_FINAL_ASSEMBLY)); in main()
30 PetscCall(MatAssemblyEnd(mat, MAT_FINAL_ASSEMBLY)); in main()
32 PetscCall(MatView(mat, PETSC_VIEWER_STDOUT_WORLD)); in main()
[all …]
|
| /petsc/src/mat/interface/ftn-custom/ |
| H A D | zmatrixf90.c | 104 PETSC_EXTERN void matgetghosts_(Mat *mat, F90Array1d *ptr, int *ierr PETSC_F90_2PTR_PROTO(ptrd)) in matgetghosts_() argument
109 *ierr = MatGetGhosts(*mat, &N, &ghosts); in matgetghosts_()
113 PETSC_EXTERN void matdensegetarray2d_(Mat *mat, F90Array2d *ptr, int *ierr PETSC_F90_2PTR_PROTO(ptr… in matdensegetarray2d_() argument
117 *ierr = MatDenseGetArray(*mat, &fa); in matdensegetarray2d_()
119 *ierr = MatGetLocalSize(*mat, &m, NULL); in matdensegetarray2d_()
121 *ierr = MatGetSize(*mat, NULL, &N); in matdensegetarray2d_()
123 *ierr = MatDenseGetLDA(*mat, &lda); in matdensegetarray2d_()
126 …*ierr = PetscError(((PetscObject)*mat)->comm, __LINE__, PETSC_FUNCTION_NAME, __FILE__, PETSC_ERR_A… in matdensegetarray2d_()
131 PETSC_EXTERN void matdenserestorearray2d_(Mat *mat, F90Array2d *ptr, int *ierr PETSC_F90_2PTR_PROTO… in matdenserestorearray2d_() argument
138 *ierr = MatDenseRestoreArray(*mat, &fa); in matdenserestorearray2d_()
[all …]
|
| /petsc/src/mat/impls/aij/seq/hdf5/ |
| H A D | aijhdf5.c | 6 PetscErrorCode MatLoad_AIJ_HDF5(Mat mat, PetscViewer viewer) in MatLoad_AIJ_HDF5() argument
16 PetscInt bs = mat->rmap->bs; in MatLoad_AIJ_HDF5()
34 …SETERRQ(PetscObjectComm((PetscObject)mat), PETSC_ERR_SUP, "PetscViewerFormat %s not supported for … in MatLoad_AIJ_HDF5()
37 PetscCall(PetscObjectGetComm((PetscObject)mat, &comm)); in MatLoad_AIJ_HDF5()
40 PetscCall(PetscObjectGetName((PetscObject)mat, &mat_name)); in MatLoad_AIJ_HDF5()
58 if (flg) PetscCall(MatSetBlockSize(mat, bs)); in MatLoad_AIJ_HDF5()
65 if (format == PETSC_VIEWER_HDF5_MAT && mat->symmetric != PETSC_BOOL3_TRUE) { in MatLoad_AIJ_HDF5()
68 …PetscCheck(!mat->preallocated, comm, PETSC_ERR_SUP, "Not for preallocated matrix - we would need t… in MatLoad_AIJ_HDF5()
69 tmp = mat->rmap; in MatLoad_AIJ_HDF5()
70 mat->rmap = mat->cmap; in MatLoad_AIJ_HDF5()
[all …]
|
| /petsc/src/ksp/pc/tests/ |
| H A D | ex3.c | 10 Mat mat; /* matrix */ in main() local
31 PetscCall(MatCreate(PETSC_COMM_SELF, &mat)); in main()
32 PetscCall(MatSetType(mat, MATSEQAIJ)); in main()
33 PetscCall(MatSetSizes(mat, n, n, n, n)); in main()
34 PetscCall(MatSetFromOptions(mat)); in main()
35 PetscCall(MatSeqAIJSetPreallocation(mat, 3, NULL)); in main()
36 PetscCall(MatSeqBAIJSetPreallocation(mat, 1, 3, NULL)); in main()
37 PetscCall(MatSeqSBAIJSetPreallocation(mat, 1, 3, NULL)); in main()
38 PetscCall(MatSeqSELLSetPreallocation(mat, 3, NULL)); in main()
46 PetscCall(MatSetValues(mat, 1, &i, 3, col, value, INSERT_VALUES)); in main()
[all …]
|
| /petsc/src/mat/impls/mffd/ |
| H A D | mffd.c | 69 static PetscErrorCode MatMFFDSetType_MFFD(Mat mat, MatMFFDType ftype) in MatMFFDSetType_MFFD() argument
76 PetscValidHeaderSpecific(mat, MAT_CLASSID, 1); in MatMFFDSetType_MFFD()
78 PetscCall(MatShellGetContext(mat, &ctx)); in MatMFFDSetType_MFFD()
88 …PetscCheck(r, PetscObjectComm((PetscObject)mat), PETSC_ERR_ARG_UNKNOWN_TYPE, "Unknown MatMFFD type… in MatMFFDSetType_MFFD()
117 PetscErrorCode MatMFFDSetType(Mat mat, MatMFFDType ftype) in MatMFFDSetType() argument
120 PetscValidHeaderSpecific(mat, MAT_CLASSID, 1); in MatMFFDSetType()
122 PetscTryMethod(mat, "MatMFFDSetType_C", (Mat, MatMFFDType), (mat, ftype)); in MatMFFDSetType()
129 static PetscErrorCode MatMFFDSetFunctioniBase_MFFD(Mat mat, FCN1 func) in MatMFFDSetFunctioniBase_MFFD() argument
134 PetscCall(MatShellGetContext(mat, &ctx)); in MatMFFDSetFunctioniBase_MFFD()
140 static PetscErrorCode MatMFFDSetFunctioni_MFFD(Mat mat, FCN2 funci) in MatMFFDSetFunctioni_MFFD() argument
[all …]
|
| /petsc/src/mat/tests/output/ |
| H A D | ex182_15.out | 2 [1] <mat:seqaij> MatAssemblyEnd_SeqAIJ(): Matrix size: 5 X 10; storage space: 0 unneeded, 0 used
3 [1] <mat:seqaij> MatAssemblyEnd_SeqAIJ(): Matrix size: 5 X 5; storage space: 0 unneeded, 0 used
4 [1] <mat:seqaij> MatAssemblyEnd_SeqAIJ(): Matrix size: 5 X 5; storage space: 0 unneeded, 5 used
5 [1] <mat:seqaij> MatAssemblyEnd_SeqAIJ(): Maximum nonzeros in any row is 0
6 [1] <mat:seqaij> MatAssemblyEnd_SeqAIJ(): Maximum nonzeros in any row is 0
7 [1] <mat:seqaij> MatAssemblyEnd_SeqAIJ(): Maximum nonzeros in any row is 1
8 [1] <mat:seqaij> MatAssemblyEnd_SeqAIJ(): Number of mallocs during MatSetValues() is 0
9 [1] <mat:seqaij> MatAssemblyEnd_SeqAIJ(): Number of mallocs during MatSetValues() is 0
10 [1] <mat:seqaij> MatAssemblyEnd_SeqAIJ(): Number of mallocs during MatSetValues() is 0
11 [1] <mat:seqaij> MatCheckCompressedRow(): Found the ratio (num_zerorows 0)/(num_localrows 5) < 0.6.…
[all …]
|
| H A D | ex182_12.out | 2 [1] <mat:seqaij> MatAssemblyEnd_SeqAIJ(): Matrix size: 5 X 10; storage space: 0 unneeded, 0 used
3 [1] <mat:seqaij> MatAssemblyEnd_SeqAIJ(): Matrix size: 5 X 5; storage space: 0 unneeded, 0 used
4 [1] <mat:seqaij> MatAssemblyEnd_SeqAIJ(): Matrix size: 5 X 5; storage space: 0 unneeded, 5 used
5 [1] <mat:seqaij> MatAssemblyEnd_SeqAIJ(): Maximum nonzeros in any row is 0
6 [1] <mat:seqaij> MatAssemblyEnd_SeqAIJ(): Maximum nonzeros in any row is 0
7 [1] <mat:seqaij> MatAssemblyEnd_SeqAIJ(): Maximum nonzeros in any row is 1
8 [1] <mat:seqaij> MatAssemblyEnd_SeqAIJ(): Number of mallocs during MatSetValues() is 0
9 [1] <mat:seqaij> MatAssemblyEnd_SeqAIJ(): Number of mallocs during MatSetValues() is 0
10 [1] <mat:seqaij> MatAssemblyEnd_SeqAIJ(): Number of mallocs during MatSetValues() is 0
11 [1] <mat:seqaij> MatCheckCompressedRow(): Found the ratio (num_zerorows 0)/(num_localrows 5) < 0.6.…
[all …]
|
| /petsc/src/mat/graphops/order/ |
| H A D | sorder.c | 11 PETSC_INTERN PetscErrorCode MatGetOrdering_Natural(Mat mat, MatOrderingType type, IS *irow, IS *ico… in MatGetOrdering_Natural() argument
18 PetscCall(PetscObjectGetComm((PetscObject)mat, &comm)); in MatGetOrdering_Natural()
19 PetscCall(MatGetRowIJ(mat, 0, PETSC_FALSE, PETSC_TRUE, &n, NULL, NULL, &done)); in MatGetOrdering_Natural()
20 PetscCall(MatRestoreRowIJ(mat, 0, PETSC_FALSE, PETSC_TRUE, NULL, NULL, NULL, &done)); in MatGetOrdering_Natural()
35 PetscCall(MatGetOwnershipRange(mat, &start, &end)); in MatGetOrdering_Natural()
49 PETSC_INTERN PetscErrorCode MatGetOrdering_RowLength(Mat mat, MatOrderingType type, IS *irow, IS *i… in MatGetOrdering_RowLength() argument
56 PetscCall(MatGetRowIJ(mat, 0, PETSC_FALSE, PETSC_TRUE, &n, &ia, &ja, &done)); in MatGetOrdering_RowLength()
57 PetscCheck(done, PetscObjectComm((PetscObject)mat), PETSC_ERR_SUP, "Cannot get rows for matrix"); in MatGetOrdering_RowLength()
64 PetscCall(MatRestoreRowIJ(mat, 0, PETSC_FALSE, PETSC_TRUE, NULL, &ia, &ja, &done)); in MatGetOrdering_RowLength()
149 PetscErrorCode MatGetOrdering(Mat mat, MatOrderingType type, IS *rperm, IS *cperm) in MatGetOrdering() argument
[all …]
|
| /petsc/src/mat/utils/ |
| H A D | convert.c | 7 PetscErrorCode MatConvert_Basic(Mat mat, MatType newtype, MatReuse reuse, Mat *newmat) in MatConvert_Basic() argument
16 if (!mat->ops->getrow) { /* missing get row, use matvecs */ in MatConvert_Basic()
17 PetscCall(MatConvert_Shell(mat, newtype, reuse, newmat)); in MatConvert_Basic()
20 PetscCall(PetscObjectTypeCompare((PetscObject)mat, MATSEQSBAIJ, &isSBAIJ)); in MatConvert_Basic()
21 if (!isSBAIJ) PetscCall(PetscObjectTypeCompare((PetscObject)mat, MATMPISBAIJ, &isSBAIJ)); in MatConvert_Basic()
22 …PetscCheck(!isSBAIJ, PetscObjectComm((PetscObject)mat), PETSC_ERR_SUP, "Cannot convert from SBAIJ … in MatConvert_Basic()
28 PetscCall(MatGetSize(mat, &m, &n)); in MatConvert_Basic()
29 PetscCall(MatGetLocalSize(mat, &lm, &ln)); in MatConvert_Basic()
30 PetscCall(MatCreate(PetscObjectComm((PetscObject)mat), &M)); in MatConvert_Basic()
32 PetscCall(MatSetBlockSizesFromMats(M, mat, mat)); in MatConvert_Basic()
[all …]
|
| /petsc/src/tao/matrix/ |
| H A D | submatfree.c | 28 PetscErrorCode MatCreateSubMatrixFree(Mat mat, IS Rows, IS Cols, Mat *J) in MatCreateSubMatrixFree() argument
30 MPI_Comm comm = PetscObjectComm((PetscObject)mat); in MatCreateSubMatrixFree()
36 ctx->A = mat; in MatCreateSubMatrixFree()
37 PetscCall(MatGetSize(mat, &m, &n)); in MatCreateSubMatrixFree()
38 PetscCall(MatGetLocalSize(mat, &mloc, &nloc)); in MatCreateSubMatrixFree()
39 PetscCall(MatCreateVecs(mat, NULL, &ctx->VC)); in MatCreateSubMatrixFree()
41 PetscCall(PetscObjectReference((PetscObject)mat)); in MatCreateSubMatrixFree()
67 PetscErrorCode MatSMFResetRowColumn(Mat mat, IS Rows, IS Cols) in MatSMFResetRowColumn() argument
72 PetscCall(MatShellGetContext(mat, &ctx)); in MatSMFResetRowColumn()
82 PetscErrorCode MatMult_SMF(Mat mat, Vec a, Vec y) in MatMult_SMF() argument
[all …]
|