| /petsc/src/mat/impls/baij/mpi/ |
| H A D | mmbaij.c | 9 Mat_MPIBAIJ *baij = (Mat_MPIBAIJ *)mat->data; in MatSetUpMultiply_MPIBAIJ() local
10 Mat_SeqBAIJ *B = (Mat_SeqBAIJ *)baij->B->data; in MatSetUpMultiply_MPIBAIJ()
20 PetscInt Nbs = baij->Nbs, *indices; in MatSetUpMultiply_MPIBAIJ()
61 PetscCall(PetscLayoutDestroy(&baij->B->cmap)); in MatSetUpMultiply_MPIBAIJ()
62 …teFromSizes(PetscObjectComm((PetscObject)baij->B), ec * mat->rmap->bs, ec * mat->rmap->bs, mat->rm… in MatSetUpMultiply_MPIBAIJ()
90 PetscCall(PetscLayoutDestroy(&baij->B->cmap)); in MatSetUpMultiply_MPIBAIJ()
91 …teFromSizes(PetscObjectComm((PetscObject)baij->B), ec * mat->rmap->bs, ec * mat->rmap->bs, mat->rm… in MatSetUpMultiply_MPIBAIJ()
96 PetscCall(VecCreateSeq(PETSC_COMM_SELF, ec * bs, &baij->lvec)); in MatSetUpMultiply_MPIBAIJ()
108 PetscCall(VecScatterCreate(gvec, from, baij->lvec, to, &baij->Mvctx)); in MatSetUpMultiply_MPIBAIJ()
109 PetscCall(VecScatterViewFromOptions(baij->Mvctx, (PetscObject)mat, "-matmult_vecscatter_view")); in MatSetUpMultiply_MPIBAIJ()
[all …]
|
| H A D | mpibaij.c | 9 Mat_MPIBAIJ *baij = (Mat_MPIBAIJ *)mat->data; in MatDestroy_MPIBAIJ() local
15 PetscCall(MatDestroy(&baij->A)); in MatDestroy_MPIBAIJ()
16 PetscCall(MatDestroy(&baij->B)); in MatDestroy_MPIBAIJ()
18 PetscCall(PetscHMapIDestroy(&baij->colmap)); in MatDestroy_MPIBAIJ()
20 PetscCall(PetscFree(baij->colmap)); in MatDestroy_MPIBAIJ()
22 PetscCall(PetscFree(baij->garray)); in MatDestroy_MPIBAIJ()
23 PetscCall(VecDestroy(&baij->lvec)); in MatDestroy_MPIBAIJ()
24 PetscCall(VecScatterDestroy(&baij->Mvctx)); in MatDestroy_MPIBAIJ()
25 PetscCall(PetscFree2(baij->rowvalues, baij->rowindices)); in MatDestroy_MPIBAIJ()
26 PetscCall(PetscFree(baij->barray)); in MatDestroy_MPIBAIJ()
[all …]
|
| /petsc/src/mat/impls/sbaij/mpi/ |
| H A D | mmsbaij.c | 183 Mat_MPISBAIJ *baij = (Mat_MPISBAIJ *)A->data; in MatDisAssemble_MPISBAIJ() local
184 Mat B = baij->B, Bnew; in MatDisAssemble_MPISBAIJ()
186 PetscInt i, j, mbs = Bbaij->mbs, n = A->cmap->N, col, *garray = baij->garray; in MatDisAssemble_MPISBAIJ()
187 PetscInt k, bs = A->rmap->bs, bs2 = baij->bs2, *rvals, *nz, m = A->rmap->n; in MatDisAssemble_MPISBAIJ()
199 PetscCall(VecDestroy(&baij->lvec)); in MatDisAssemble_MPISBAIJ()
200 PetscCall(VecScatterDestroy(&baij->Mvctx)); in MatDisAssemble_MPISBAIJ()
202 PetscCall(VecDestroy(&baij->slvec0)); in MatDisAssemble_MPISBAIJ()
203 PetscCall(VecDestroy(&baij->slvec0b)); in MatDisAssemble_MPISBAIJ()
204 PetscCall(VecDestroy(&baij->slvec1)); in MatDisAssemble_MPISBAIJ()
205 PetscCall(VecDestroy(&baij->slvec1a)); in MatDisAssemble_MPISBAIJ()
[all …]
|
| H A D | mpisbaij.c | 8 Mat_MPISBAIJ *baij = (Mat_MPISBAIJ *)mat->data; in MatDestroy_MPISBAIJ() local
14 PetscCall(MatDestroy(&baij->A)); in MatDestroy_MPISBAIJ()
15 PetscCall(MatDestroy(&baij->B)); in MatDestroy_MPISBAIJ()
17 PetscCall(PetscHMapIDestroy(&baij->colmap)); in MatDestroy_MPISBAIJ()
19 PetscCall(PetscFree(baij->colmap)); in MatDestroy_MPISBAIJ()
21 PetscCall(PetscFree(baij->garray)); in MatDestroy_MPISBAIJ()
22 PetscCall(VecDestroy(&baij->lvec)); in MatDestroy_MPISBAIJ()
23 PetscCall(VecScatterDestroy(&baij->Mvctx)); in MatDestroy_MPISBAIJ()
24 PetscCall(VecDestroy(&baij->slvec0)); in MatDestroy_MPISBAIJ()
25 PetscCall(VecDestroy(&baij->slvec0b)); in MatDestroy_MPISBAIJ()
[all …]
|
| /petsc/src/mat/impls/baij/seq/ |
| H A D | baij.c | 2315 Mat_SeqBAIJ *baij = (Mat_SeqBAIJ *)A->data; in MatZeroRows_SeqBAIJ() local
2317 PetscInt bs = A->rmap->bs, bs2 = baij->bs2, *sizes, row, bs_max; in MatZeroRows_SeqBAIJ()
2341 if (baij->keepnonzeropattern) { in MatZeroRows_SeqBAIJ()
2352 count = (baij->i[row / bs + 1] - baij->i[row / bs]) * bs; in MatZeroRows_SeqBAIJ()
2353 aa = baij->a + baij->i[row / bs] * bs2 + (row % bs); in MatZeroRows_SeqBAIJ()
2354 if (sizes[i] == bs && !baij->keepnonzeropattern) { in MatZeroRows_SeqBAIJ()
2356 if (baij->ilen[row / bs] > 0) { in MatZeroRows_SeqBAIJ()
2357 baij->ilen[row / bs] = 1; in MatZeroRows_SeqBAIJ()
2358 baij->j[baij->i[row / bs]] = row / bs; in MatZeroRows_SeqBAIJ()
2365 baij->ilen[row / bs] = 0; in MatZeroRows_SeqBAIJ()
[all …]
|
| /petsc/src/mat/impls/sbaij/seq/ |
| H A D | sbaij.c | 862 Mat_SeqSBAIJ *baij = (Mat_SeqSBAIJ *)mat->data; in MatSeqSBAIJSetColumnIndices_SeqSBAIJ() local
866 nz = baij->maxnz; in MatSeqSBAIJSetColumnIndices_SeqSBAIJ()
868 for (i = 0; i < nz; i++) baij->j[i] = indices[i]; in MatSeqSBAIJSetColumnIndices_SeqSBAIJ()
870 baij->nz = nz; in MatSeqSBAIJSetColumnIndices_SeqSBAIJ()
871 for (i = 0; i < n; i++) baij->ilen[i] = baij->imax[i]; in MatSeqSBAIJSetColumnIndices_SeqSBAIJ()
1056 Mat_SeqSBAIJ *baij = (Mat_SeqSBAIJ *)A->data; in MatZeroRowsColumns_SeqSBAIJ() local
1058 PetscInt bs = A->rmap->bs, bs2 = baij->bs2, row, col; in MatZeroRowsColumns_SeqSBAIJ()
1082 for (j = baij->i[row]; j < baij->i[row + 1]; j++) { in MatZeroRowsColumns_SeqSBAIJ()
1084 col = bs * baij->j[j] + k; in MatZeroRowsColumns_SeqSBAIJ()
1086 aa = baij->a + j * bs2 + (i % bs) + bs * k; in MatZeroRowsColumns_SeqSBAIJ()
[all …]
|
| /petsc/src/mat/impls/aij/mpi/ |
| H A D | fdmpiaij.c | 433 Mat_MPIBAIJ *baij = (Mat_MPIBAIJ *)mat->data; in MatFDColoringSetUp_MPIXAIJ() local
435 A = baij->A; in MatFDColoringSetUp_MPIXAIJ()
438 B = baij->B; in MatFDColoringSetUp_MPIXAIJ()
442 if (!baij->colmap) PetscCall(MatCreateColmap_MPIBAIJ_Private(mat)); in MatFDColoringSetUp_MPIXAIJ()
443 colmap = baij->colmap; in MatFDColoringSetUp_MPIXAIJ()
450 for (PetscInt i = 0; i < baij->B->cmap->n / bs; i++) { in MatFDColoringSetUp_MPIXAIJ()
451 for (j = 0; j < bs; j++) garray[i * bs + j] = bs * baij->garray[i] + j; in MatFDColoringSetUp_MPIXAIJ()
|
| H A D | mpiaij.c | 2629 Mat_SeqAIJ *aaij = (Mat_SeqAIJ *)maij->A->data, *baij = (Mat_SeqAIJ *)maij->B->data; in MatMPIAIJGetNumberNonzeros() local
2635 *nz = aaij->i[A->rmap->n] + baij->i[A->rmap->n]; in MatMPIAIJGetNumberNonzeros()
|
| /petsc/src/binding/petsc4py/src/lib-petsc/ |
| H A D | custom.h | 204 PetscErrorCode MatHasPreallocationAIJ(Mat A,PetscBool *aij,PetscBool *baij,PetscBool *sbaij,PetscBo… in MatHasPreallocationAIJ() argument
212 PetscAssertPointer(baij,3); in MatHasPreallocationAIJ()
215 *aij = *baij = *sbaij = *is = PETSC_FALSE; in MatHasPreallocationAIJ()
221 if (f) {*baij = PETSC_TRUE; goto done;} in MatHasPreallocationAIJ()
|
| /petsc/src/mat/impls/maij/ |
| H A D | maij.c | 198 const Mat baij = b->AIJ; in MatMult_MatMultAdd_SeqMAIJ_Template() local
199 const Mat_SeqAIJ *a = (Mat_SeqAIJ *)baij->data; in MatMult_MatMultAdd_SeqMAIJ_Template()
200 const PetscInt m = baij->rmap->n; in MatMult_MatMultAdd_SeqMAIJ_Template()
259 const Mat baij = b->AIJ; in MatMultTranspose_MatMultTransposeAdd_SeqMAIJ_Template() local
260 const Mat_SeqAIJ *a = (Mat_SeqAIJ *)baij->data; in MatMultTranspose_MatMultTransposeAdd_SeqMAIJ_Template()
261 const PetscInt m = baij->rmap->n; in MatMultTranspose_MatMultTransposeAdd_SeqMAIJ_Template()
|
| /petsc/src/binding/petsc4py/src/petsc4py/PETSc/ |
| H A D | petscmat.pxi | 783 cdef PetscBool aij=PETSC_FALSE, baij=PETSC_FALSE, sbaij=PETSC_FALSE, aijis=PETSC_FALSE
784 CHKERR(MatHasPreallocationAIJ(A, &aij, &baij, &sbaij, &aijis))
788 if baij == PETSC_TRUE or sbaij == PETSC_TRUE:
822 if baij == PETSC_TRUE:
834 cdef PetscBool aij=PETSC_FALSE, baij=PETSC_FALSE, sbaij=PETSC_FALSE, aijis=PETSC_FALSE
835 CHKERR(MatHasPreallocationAIJ(A, &aij, &baij, &sbaij, &aijis))
839 if baij == PETSC_TRUE or sbaij == PETSC_TRUE:
875 if baij == PETSC_TRUE:
|
| /petsc/doc/developers/ |
| H A D | testing.md | 433 args: -mat_type baij -matload_block_size {{2 3}shared output}
447 args: -f ${DATAFILESPATH}/matrices/small -mat_type baij
1029 args: -f ${DATAFILESPATH}/matrices/small -mat_type {{aij baij sbaij}} -matload_block_size {{2 3}}
1040 ok 5 In mat...tests: "./ex138 -f ${DATAFILESPATH}/matrices/small -mat_type baij -matload_block_size…
1041 ok 6 In mat...tests: "Diff of ./ex138 -f ${DATAFILESPATH}/matrices/small -mat_type baij -matload_bl…
1080 ok 11 ./ex3 -f /matrices/small -mat_type baij -matload_block_size 2
1082 ok 13 ./ex3 -f /matrices/small -mat_type baij -matload_block_size 3
|
| /petsc/src/dm/impls/da/ |
| H A D | fdda.c | 617 …PetscBool aij = PETSC_FALSE, baij = PETSC_FALSE, sbaij = PETSC_FALSE, sell = PETSC_FALSE, is = P… in DMCreateMatrix_DA() local
677 PetscCall(PetscObjectHasFunction((PetscObject)A, "MatMPIBAIJSetPreallocation_C", &baij)); in DMCreateMatrix_DA()
678 …if (!baij) PetscCall(PetscObjectHasFunction((PetscObject)A, "MatSeqBAIJSetPreallocation_C", &baij)… in DMCreateMatrix_DA()
679 if (!baij) { in DMCreateMatrix_DA()
718 } else if (baij) { in DMCreateMatrix_DA()
|
| /petsc/src/mat/impls/aij/seq/ |
| H A D | matmatmult.c | 1697 PetscBool baij; in MatProductSetFromOptions_SeqXBAIJ_SeqDense_AB() local
1700 PetscCall(PetscObjectTypeCompare((PetscObject)A, MATSEQBAIJ, &baij)); in MatProductSetFromOptions_SeqXBAIJ_SeqDense_AB()
1701 if (!baij) { /* A is seqsbaij */ in MatProductSetFromOptions_SeqXBAIJ_SeqDense_AB()
|
| /petsc/doc/manual/ |
| H A D | profiling.md | 163 mpiexec -n 4 ./ex10 -f0 medium -f1 arco6 -ksp_gmres_classicalgramschmidt -log_view -mat_type baij \
235 mpiexec -n 4 ./ex10 -f0 medium -f1 arco6 -ksp_gmres_classicalgramschmidt -log_view -mat_type baij \
|
| H A D | ksp.md | 2548 * - ``baij``
2568 * - ``baij``
2621 The same holds for `baij` and `sbaij`. For codes that are intended
|
| H A D | mat.md | 44 and the options database option `-dm_mat_type <aij or baij or aijcusparse etc>` Matrices can be cre…
|
| /petsc/src/benchmarks/results/ |
| H A D | performance_arco1 | 3 …rco1 -pc_type ilu -ksp_gmres_unmodifiedgramschmidt -optionsleft -mat_type baij -matload_ignore_inf…
|
| /petsc/systems/Apple/iOS/PETSc/PETSc.xcodeproj/ |
| H A D | project.pbxproj | 296 …D89D1DA1207AA1AB00FACC19 /* baij.c in Sources */ = {isa = PBXBuildFile; fileRef = D89D1AB9207AA1A9…
1048 …D1AB9207AA1A900FACC19 /* baij.c */ = {isa = PBXFileReference; fileEncoding = 4; lastKnownFileType …
1819 D89D1AB9207AA1A900FACC19 /* baij.c */,
2709 D89D1DA1207AA1AB00FACC19 /* baij.c in Sources */,
|
| /petsc/doc/faq/ |
| H A D | index.md | 394 compiler chokes on `src/mat/impls/baij/seq/baijsolvtrannat.c` you can run the following
398 $ make -f gmakefile PCC_FLAGS="-O1" $PETSC_ARCH/obj/src/mat/impls/baij/seq/baijsolvtrannat.o
|