| /petsc/src/mat/impls/aij/seq/ |
| H A D | fdaij.c | 22 c->bcols = 1; in MatFDColoringCreate_SeqXAIJ()
26 PetscInt nz, brows, bcols; in MatFDColoringCreate_SeqXAIJ() local
37 bcols = (PetscInt)(0.5 * mem / (m * sizeof(PetscScalar))); in MatFDColoringCreate_SeqXAIJ()
38 if (!bcols) bcols = 1; in MatFDColoringCreate_SeqXAIJ()
39 brows = 1000 / bcols; in MatFDColoringCreate_SeqXAIJ()
40 if (bcols > nis) bcols = nis; in MatFDColoringCreate_SeqXAIJ()
43 c->bcols = bcols; in MatFDColoringCreate_SeqXAIJ()
64 PetscInt i, j, nrows, nbcols, brows = c->brows, bcols = c->bcols, mbs = c->m, nis = c->ncolors; in MatFDColoringSetUpBlocked_AIJ_Private() local
69 PetscCall(PetscMalloc2(bcols + 1, &color_start, bcols, &row_start)); in MatFDColoringSetUpBlocked_AIJ_Private()
71 PetscCall(PetscMalloc1(bcols * mat->rmap->n, &c->dy)); in MatFDColoringSetUpBlocked_AIJ_Private()
[all …]
|
| /petsc/src/mat/impls/aij/mpi/ |
| H A D | fdmpiaij.c | 247 if (coloring->bcols > 1) { /* use blocked insertion of Jentry */ in MatFDColoringApply_AIJ()
248 PetscInt i, m = J->rmap->n, nbcols, bcols = coloring->bcols; in MatFDColoringApply_AIJ() local
252 for (k = 0; k < ncolors; k += bcols) { in MatFDColoringApply_AIJ()
259 if (k + bcols > ncolors) bcols = ncolors - k; in MatFDColoringApply_AIJ()
260 for (i = 0; i < bcols; i++) { in MatFDColoringApply_AIJ()
416 PetscInt bcols = c->bcols; in MatFDColoringSetUp_MPIXAIJ() local
634 if (bcols > 1) { /* reorder Jentry for faster MatFDColoringApply() */ in MatFDColoringSetUp_MPIXAIJ()
654 …brows %" PetscInt_FMT " and bcols %" PetscInt_FMT " are used.\n", c->ncolors, c->brows, c->bcols)); in MatFDColoringSetUp_MPIXAIJ()
671 c->bcols = 1; in MatFDColoringCreate_MPIXAIJ()
677 PetscInt nz, brows, bcols; in MatFDColoringCreate_MPIXAIJ() local
[all …]
|
| H A D | mpiaij.c | 1575 PetscInt *acols = dnnz, *bcols = onnz; /* Repurpose now-unneeded arrays */ in MatPermute_MPIAIJ() local
1584 for (; j < PetscMin(rowlen, j0 + m); j++) bcols[j - j0] = gcdest[bj[bi[i] + j]]; in MatPermute_MPIAIJ()
1585 PetscCall(MatSetValues(Aperm, 1, &rdest[i], j - j0, bcols, ba + bi[i] + j0, INSERT_VALUES)); in MatPermute_MPIAIJ()
|
| /petsc/src/mat/matfd/ |
| H A D | fdmatrix.c | 195 PetscErrorCode MatFDColoringSetBlockSize(MatFDColoring matfd, PetscInt brows, PetscInt bcols) in MatFDColoringSetBlockSize() argument
200 PetscValidLogicalCollectiveInt(matfd, bcols, 3); in MatFDColoringSetBlockSize()
202 if (bcols != PETSC_DEFAULT) matfd->bcols = bcols; in MatFDColoringSetBlockSize()
346 …cols", "Number of block columns", "MatFDColoringSetBlockSize", matfd->bcols, &matfd->bcols, &flg)); in MatFDColoringSetFromOptions()
347 if (flg && matfd->bcols > matfd->ncolors) { in MatFDColoringSetFromOptions()
349 matfd->bcols = matfd->ncolors; in MatFDColoringSetFromOptions()
|
| /petsc/src/tao/interface/ |
| H A D | taosolver_hj.c | 169 const PetscInt *bcols; in TaoTestHessian() local
181 PetscCall(MatGetRow(B, row, &bncols, &bcols, &bvals)); in TaoTestHessian()
185 ccols[cncols] = bcols[j]; in TaoTestHessian()
191 PetscCall(MatRestoreRow(B, row, &bncols, &bcols, &bvals)); in TaoTestHessian()
|
| /petsc/src/snes/interface/ |
| H A D | snes.c | 2859 const PetscInt *bcols; in SNESTestJacobian() local
2873 PetscCall(MatGetRow(B, row, &bncols, &bcols, &bvals)); in SNESTestJacobian()
2877 ccols[cncols] = bcols[j]; in SNESTestJacobian()
2883 PetscCall(MatRestoreRow(B, row, &bncols, &bcols, &bvals)); in SNESTestJacobian()
|
| /petsc/include/petsc/private/ |
| H A D | matimpl.h | 669 …PetscInt brows, bcols; /* number of block rows or columns for speedup inserting… member
|