| /petsc/src/mat/tests/ |
| H A D | ex79f.F90 | 9 Mat A, Ad, Ao
29 PetscCallA(MatMPIAIJGetSeqAIJ(A, Ad, Ao, icol, ierr))
38 PetscCallA(MatGetRowIJ(Ad, one, bb, bb, n, ia, ja, done, ierr))
39 PetscCallA(MatSeqAIJGetArray(Ad, aa, ierr))
46 PetscCallA(MatRestoreRowIJ(Ad, one, bb, bb, n, ia, ja, done, ierr))
47 PetscCallA(MatSeqAIJRestoreArray(Ad, aa, ierr))
59 PetscCallA(MatMPIAIJRestoreSeqAIJ(A, Ad, Ao, icol, ierr))
65 PetscCallA(MatGetDiagonalBlock(A, Ad, ierr))
66 PetscCallA(MatView(Ad, PETSC_VIEWER_STDOUT_WORLD, ierr))
|
| H A D | ex248.c | 7 Mat A, B, C, K, Ad, Bd; in main() local
14 PetscCall(MatCreateDense(PETSC_COMM_SELF, m, n, m, n, NULL, &Ad)); in main()
16 PetscCall(MatSetRandom(Ad, NULL)); in main()
18 PetscCall(MatFilter(Ad, 0.2, PETSC_FALSE, PETSC_FALSE)); in main()
20 PetscCall(MatConvert(Ad, MATAIJ, MAT_INITIAL_MATRIX, &A)); in main()
45 PetscCall(MatDestroy(&Ad)); in main()
|
| H A D | ex268.c | 23 Mat A, Ad, F, Fd, X, Xd, B; in main() local
58 PetscCall(MatConvert(A, MATDENSE, MAT_INITIAL_MATRIX, &Ad)); in main()
59 PetscCall(MatMultEqual(A, Ad, 10, &flg)); in main()
64 PetscCall(MatViewFromOptions(Ad, NULL, "-Ad")); in main()
68 PetscCall(MatGetFactor(Ad, MATSOLVERPETSC, i == 0 ? MAT_FACTOR_LU : MAT_FACTOR_CHOLESKY, &Fd)); in main()
75 PetscCall(MatLUFactorSymbolic(Fd, Ad, NULL, NULL, NULL)); in main()
76 PetscCall(MatLUFactorNumeric(Fd, Ad, NULL)); in main()
80 PetscCall(MatCholeskyFactorSymbolic(Fd, Ad, NULL, NULL)); in main()
81 PetscCall(MatCholeskyFactorNumeric(Fd, Ad, NULL)); in main()
125 PetscCall(MatDestroy(&Ad)); in main()
|
| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex21.c | 27 Mat A, Ae = NULL, Ad = NULL; in main() local
71 PetscCall(MatCreateDense(PETSC_COMM_WORLD, n, n, PETSC_DECIDE, PETSC_DECIDE, NULL, &Ad)); in main()
73 …for (j = 0; j < n; j++) PetscCall(MatSetValue(Ad, i, j, RBF(sdim, coords + i * sdim, coords + j * … in main()
75 PetscCall(MatAssemblyBegin(Ad, MAT_FINAL_ASSEMBLY)); in main()
76 PetscCall(MatAssemblyEnd(Ad, MAT_FINAL_ASSEMBLY)); in main()
95 PetscCall(MatAXPY(Ae, -1.0, Ad, SAME_NONZERO_PATTERN)); in main()
107 PetscCall(MatMult(Ad, u, b)); in main()
108 PetscCall(MatViewFromOptions(Ad, NULL, "-Ad_view")); in main()
110 PetscCall(KSPSetOperators(ksp, Ad, A)); in main()
127 PetscCall(MatDiagonalSet(Ad, x, ADD_VALUES)); in main()
[all …]
|
| /petsc/src/ksp/ksp/tests/ |
| H A D | ex34.c | 7 Mat A, Ad, B; in main() local
21 PetscCall(MatCreateDense(PETSC_COMM_WORLD, PETSC_DECIDE, PETSC_DECIDE, M, N, NULL, &Ad)); in main()
22 PetscCall(MatSetRandom(Ad, NULL)); in main()
23 PetscCall(MatConvert(Ad, flg ? mtype : MATAIJ, MAT_INITIAL_MATRIX, &A)); in main()
49 PetscCall(MatDestroy(&Ad)); in main()
|
| /petsc/src/mat/impls/baij/mpi/ftn-custom/ |
| H A D | zmpibaijf.c | 10 PETSC_EXTERN void matmpibaijgetseqbaij_(Mat *A, Mat *Ad, Mat *Ao, PetscInt *ic, size_t *iic, PetscE… in matmpibaijgetseqbaij_() argument
13 *ierr = MatMPIBAIJGetSeqBAIJ(*A, Ad, Ao, &i); in matmpibaijgetseqbaij_()
|
| /petsc/src/mat/graphops/order/ |
| H A D | sorder.c | 173 Mat Ad, Ao; in MatGetOrdering() local
179 PetscCall(MatMPIAIJGetSeqAIJ(mat, &Ad, &Ao, &colmap)); in MatGetOrdering()
180 PetscCall(MatGetOrdering(Ad, type, &lrowperm, &lcolperm)); in MatGetOrdering()
|
| /petsc/src/binding/petsc4py/demo/legacy/kspsolve/ |
| H A D | petsc-cg.tex | 11 & \quad\quad\quad q \Leftarrow Ad \\
|
| /petsc/src/mat/impls/aij/mpi/kokkos/ |
| H A D | mpiaijkok.kokkos.cxx | 97 Mat Ad, Ao; in MatMPIAIJGetLocalMatMerge_MPIAIJKokkos() local
101 PetscCall(MatMPIAIJGetSeqAIJ(mat, &Ad, &Ao, &cmap)); in MatMPIAIJGetLocalMatMerge_MPIAIJKokkos()
102 PetscCall(MatSeqAIJKokkosMergeMats(Ad, Ao, reuse, C)); in MatMPIAIJGetLocalMatMerge_MPIAIJKokkos()
105 PetscCall(MatGetLocalSize(Ad, NULL, &dn)); in MatMPIAIJGetLocalMatMerge_MPIAIJKokkos()
111 …PetscCall(ISCreateGeneral(PetscObjectComm((PetscObject)Ad), dn + on, gidx, PETSC_OWN_POINTER, glob… in MatMPIAIJGetLocalMatMerge_MPIAIJKokkos()
1054 KokkosCsrMatrix Adt, Aot, Ad, Ao, Bd, Bo; in MatProductSymbolic_MPIAIJKokkos_AtB() local
1062 PetscCall(MatSeqAIJKokkosGetKokkosCsrMatrix(ampi->A, &Ad)); in MatProductSymbolic_MPIAIJKokkos_AtB()
1173 KokkosCsrMatrix Ad, Ao, Bd, Bo; in MatProductSymbolic_MPIAIJKokkos_AB() local
1176 PetscCall(MatSeqAIJKokkosGetKokkosCsrMatrix(ampi->A, &Ad)); in MatProductSymbolic_MPIAIJKokkos_AB()
1201 PetscCallCXX(KokkosSparse::spgemm_symbolic(mm->kh1, Ad, false, Bd, false, mm->C1)); in MatProductSymbolic_MPIAIJKokkos_AB()
[all …]
|
| /petsc/src/mat/impls/aij/mpi/mpihipsparse/ |
| H A D | mpiaijhipsparse.hip.cxx | 242 Mat Ad, Ao; in MatMPIAIJGetLocalMatMerge_MPIAIJHIPSPARSE() local
246 PetscCall(MatMPIAIJGetSeqAIJ(A, &Ad, &Ao, &cmap)); in MatMPIAIJGetLocalMatMerge_MPIAIJHIPSPARSE()
247 PetscCall(MatSeqAIJHIPSPARSEMergeMats(Ad, Ao, scall, A_loc)); in MatMPIAIJGetLocalMatMerge_MPIAIJHIPSPARSE()
251 PetscCall(MatGetLocalSize(Ad, NULL, &dn)); in MatMPIAIJGetLocalMatMerge_MPIAIJHIPSPARSE()
257 …PetscCall(ISCreateGeneral(PetscObjectComm((PetscObject)Ad), dn + on, gidx, PETSC_OWN_POINTER, glob… in MatMPIAIJGetLocalMatMerge_MPIAIJHIPSPARSE()
|
| /petsc/src/mat/impls/aij/mpi/mpicusparse/ |
| H A D | mpiaijcusparse.cu | 242 Mat Ad, Ao; in MatMPIAIJGetLocalMatMerge_MPIAIJCUSPARSE() local
246 PetscCall(MatMPIAIJGetSeqAIJ(A, &Ad, &Ao, &cmap)); in MatMPIAIJGetLocalMatMerge_MPIAIJCUSPARSE()
247 PetscCall(MatSeqAIJCUSPARSEMergeMats(Ad, Ao, scall, A_loc)); in MatMPIAIJGetLocalMatMerge_MPIAIJCUSPARSE()
251 PetscCall(MatGetLocalSize(Ad, NULL, &dn)); in MatMPIAIJGetLocalMatMerge_MPIAIJCUSPARSE()
257 …PetscCall(ISCreateGeneral(PetscObjectComm((PetscObject)Ad), dn + on, gidx, PETSC_OWN_POINTER, glob… in MatMPIAIJGetLocalMatMerge_MPIAIJCUSPARSE()
|
| /petsc/src/ksp/pc/impls/asm/ |
| H A D | asm.c | 1365 Mat Ad = NULL, adj; in PCASMCreateSubdomains() local
1381 if (hasop) PetscCall(MatGetDiagonalBlock(A, &Ad)); in PCASMCreateSubdomains()
1382 if (Ad) { in PCASMCreateSubdomains()
1383 PetscCall(PetscObjectBaseTypeCompare((PetscObject)Ad, MATSEQBAIJ, &isbaij)); in PCASMCreateSubdomains()
1384 if (!isbaij) PetscCall(PetscObjectBaseTypeCompare((PetscObject)Ad, MATSEQSBAIJ, &isbaij)); in PCASMCreateSubdomains()
1386 if (Ad && n > 1) { in PCASMCreateSubdomains()
1397 PetscCall(MatGetRowIJ(Ad, 0, PETSC_TRUE, isbaij, &na, &ia, &ja, &done)); in PCASMCreateSubdomains()
1442 PetscCall(MatRestoreRowIJ(Ad, 0, PETSC_TRUE, isbaij, &na, &ia, &ja, &done)); in PCASMCreateSubdomains()
|
| /petsc/src/ksp/pc/impls/gasm/ |
| H A D | gasm.c | 1325 Mat Ad = NULL, adj; in PCGASMCreateLocalSubdomains() local
1339 if (hasop) PetscCall(MatGetDiagonalBlock(A, &Ad)); in PCGASMCreateLocalSubdomains()
1340 if (Ad) { in PCGASMCreateLocalSubdomains()
1341 PetscCall(PetscObjectBaseTypeCompare((PetscObject)Ad, MATSEQBAIJ, &isbaij)); in PCGASMCreateLocalSubdomains()
1342 if (!isbaij) PetscCall(PetscObjectBaseTypeCompare((PetscObject)Ad, MATSEQSBAIJ, &isbaij)); in PCGASMCreateLocalSubdomains()
1344 if (Ad && nloc > 1) { in PCGASMCreateLocalSubdomains()
1355 PetscCall(MatGetRowIJ(Ad, 0, PETSC_TRUE, isbaij, &na, &ia, &ja, &done)); in PCGASMCreateLocalSubdomains()
1398 PetscCall(MatRestoreRowIJ(Ad, 0, PETSC_TRUE, isbaij, &na, &ia, &ja, &done)); in PCGASMCreateLocalSubdomains()
|
| /petsc/src/mat/impls/is/ |
| H A D | matis.c | 418 Mat Ad, Ao; in MatMPIXAIJComputeLocalToGlobalMapping_Private() local
600 PetscCall(MatMPIAIJGetSeqAIJ(A, &Ad, &Ao, &garray)); in MatMPIXAIJComputeLocalToGlobalMapping_Private()
602 PetscCall(MatMPIBAIJGetSeqBAIJ(A, &Ad, &Ao, &garray)); in MatMPIXAIJComputeLocalToGlobalMapping_Private()
607 PetscCall(MatGetLocalSize(Ad, NULL, &dc)); in MatMPIXAIJComputeLocalToGlobalMapping_Private()
629 Mat lA, Ad, Ao, B = NULL; in MatConvert_XAIJ_IS() local
711 PetscCall(MatMPIAIJGetSeqAIJ(A, &Ad, &Ao, &garray)); in MatConvert_XAIJ_IS()
713 PetscCall(MatMPIBAIJGetSeqBAIJ(A, &Ad, &Ao, &garray)); in MatConvert_XAIJ_IS()
714 PetscCall(MatConvert(Ad, MATSEQAIJ, MAT_INITIAL_MATRIX, &Ad)); in MatConvert_XAIJ_IS()
717 PetscCall(MatSeqAIJGetArray(Ad, &dd)); in MatConvert_XAIJ_IS()
723 PetscCall(MatGetLocalSize(Ad, &dr, &dc)); in MatConvert_XAIJ_IS()
[all …]
|
| /petsc/src/mat/impls/aij/mpi/ |
| H A D | mpiaij.c | 4230 Mat_SeqAIJ *Ad = (Mat_SeqAIJ *)Aij->A->data; in MatUpdateMPIAIJWithArrays() local
4233 const PetscInt *Adi = Ad->i; in MatUpdateMPIAIJWithArrays()
4299 Mat_SeqAIJ *Ad = (Mat_SeqAIJ *)Aij->A->data; in MatUpdateMPIAIJWithArray() local
4302 const PetscInt *Adi = Ad->i, *Adj = Ao->i; in MatUpdateMPIAIJWithArray()
4556 PetscErrorCode MatMPIAIJGetSeqAIJ(Mat A, Mat *Ad, Mat *Ao, const PetscInt *colmap[]) in MatMPIAIJGetSeqAIJ() argument
4564 if (Ad) *Ad = a->A; in MatMPIAIJGetSeqAIJ()
5249 Mat Ao, Ad; in MatMPIAIJGetLocalMatMerge() local
5255 PetscCall(MatMPIAIJGetSeqAIJ(A, &Ad, &Ao, &cmap)); in MatMPIAIJGetLocalMatMerge()
5259 PetscCall(PetscObjectReference((PetscObject)Ad)); in MatMPIAIJGetLocalMatMerge()
5260 *A_loc = Ad; in MatMPIAIJGetLocalMatMerge()
[all …]
|
| /petsc/src/mat/impls/aij/mpi/mumps/ |
| H A D | mumps.c | 1179 Mat Ad, Ao; in MatConvertToTriples_mpiaij_mpiaij() local
1184 PetscCall(MatMPIAIJGetSeqAIJ(A, &Ad, &Ao, &garray)); in MatConvertToTriples_mpiaij_mpiaij()
1185 PetscCall(MatSeqAIJGetArrayRead(Ad, &av)); in MatConvertToTriples_mpiaij_mpiaij()
1188 aa = (Mat_SeqAIJ *)Ad->data; in MatConvertToTriples_mpiaij_mpiaij()
1239 PetscCall(MatSeqAIJRestoreArrayRead(Ad, &av)); in MatConvertToTriples_mpiaij_mpiaij()
1321 Mat Ad, Ao; in MatConvertToTriples_mpiaij_mpisbaij() local
1333 PetscCall(MatMPIAIJGetSeqAIJ(A, &Ad, &Ao, &garray)); in MatConvertToTriples_mpiaij_mpisbaij()
1334 PetscCall(MatSeqAIJGetArrayRead(Ad, &av)); in MatConvertToTriples_mpiaij_mpisbaij()
1337 aa = (Mat_SeqAIJ *)Ad->data; in MatConvertToTriples_mpiaij_mpisbaij()
1343 PetscCall(MatGetDiagonalMarkers_SeqAIJ(Ad, &adiag, NULL)); in MatConvertToTriples_mpiaij_mpisbaij()
[all …]
|
| /petsc/src/mat/impls/sell/mpi/ |
| H A D | mpisell.c | 1539 PetscErrorCode MatMPISELLGetSeqSELL(Mat A, Mat *Ad, Mat *Ao, const PetscInt *colmap[]) in MatMPISELLGetSeqSELL() argument
1547 if (Ad) *Ad = a->A; in MatMPISELLGetSeqSELL()
|
| /petsc/src/mat/impls/baij/mpi/ |
| H A D | mpibaij.c | 3370 PetscErrorCode MatMPIBAIJGetSeqBAIJ(Mat A, Mat *Ad, Mat *Ao, const PetscInt *colmap[]) in MatMPIBAIJGetSeqBAIJ() argument
3378 if (Ad) *Ad = a->A; in MatMPIBAIJGetSeqBAIJ()
|