Home
last modified time | relevance | path

Searched refs:cbs (Results 1 – 25 of 30) sorted by relevance

12

/petsc/src/mat/utils/ !
H A Dmatio.c8 PetscInt rbs, cbs; in MatView_Binary_BlockSizes() local
14 PetscCall(MatGetBlockSizes(mat, &rbs, &cbs)); in MatView_Binary_BlockSizes()
18 …if (rbs != cbs) PetscCall(PetscFPrintf(PETSC_COMM_SELF, info, "-matload_block_size %" PetscInt_FMT… in MatView_Binary_BlockSizes()
26 PetscInt rbs, cbs, bs[2], n = 2; in MatLoad_Binary_BlockSizes() local
31 PetscCall(MatGetBlockSizes(mat, &rbs, &cbs)); in MatLoad_Binary_BlockSizes()
33 bs[1] = cbs; in MatLoad_Binary_BlockSizes()
42 if (bs[1] > 0) cbs = bs[1]; in MatLoad_Binary_BlockSizes()
43 PetscCall(MatSetBlockSizes(mat, rbs, cbs)); in MatLoad_Binary_BlockSizes()
H A Dgcreate.c6 PetscErrorCode MatSetBlockSizes_Default(Mat mat, PetscInt rbs, PetscInt cbs) in MatSetBlockSizes_Default() argument
11 … == cbs, PetscObjectComm((PetscObject)mat), PETSC_ERR_SUP, "Cannot change column block size %" Pet… in MatSetBlockSizes_Default()
355 PetscInt cbs; in MatXAIJSetPreallocation() local
364 PetscCall(MatGetBlockSizes(A, &bs, &cbs)); in MatXAIJSetPreallocation()
379 if (bs == cbs && bs == 1) { in MatXAIJSetPreallocation()
391 if (dnnz) sdnnz[i] = dnnz[i / bs] * cbs; in MatXAIJSetPreallocation()
392 if (onnz) sonnz[i] = onnz[i / bs] * cbs; in MatXAIJSetPreallocation()
/petsc/src/mat/impls/localref/ !
H A Dmlocalref.c53 PetscInt rbs, cbs, buf[4096], *irowm, *icolm; in MatSetValuesBlockedLocal_LocalRef_Scalar() local
56 PetscCall(MatGetBlockSizes(A, &rbs, &cbs)); in MatSetValuesBlockedLocal_LocalRef_Scalar()
57 IndexSpaceGet(buf, nrow * rbs, ncol * cbs, irowm, icolm); in MatSetValuesBlockedLocal_LocalRef_Scalar()
59 BlockIndicesExpand(ncol, icol, cbs, icolm); in MatSetValuesBlockedLocal_LocalRef_Scalar()
61 PetscCall(ISLocalToGlobalMappingApplyBlock(A->cmap->mapping, ncol * cbs, icolm, icolm)); in MatSetValuesBlockedLocal_LocalRef_Scalar()
62 PetscCall((*lr->SetValues)(lr->Top, nrow * rbs, irowm, ncol * cbs, icolm, y, addv)); in MatSetValuesBlockedLocal_LocalRef_Scalar()
63 IndexSpaceRestore(buf, nrow * rbs, ncol * cbs, irowm, icolm); in MatSetValuesBlockedLocal_LocalRef_Scalar()
252 PetscInt arbs, acbs, rbs, cbs; in MatCreateLocalRef() local
279 PetscCall(ISGetBlockSize(iscol, &cbs)); in MatCreateLocalRef()
282 PetscCall(PetscLayoutSetBlockSize(B->cmap, cbs)); in MatCreateLocalRef()
[all …]
/petsc/src/mat/tests/ !
H A Dex46.c11 PetscInt bs[4], l2gbs[4], rbs, cbs, l2grbs, l2gcbs, i; in main() local
41 PetscCall(MatGetBlockSizes(A, &rbs, &cbs)); in main()
46 …t_FMT " %" PetscInt_FMT " (l2g %" PetscInt_FMT " %" PetscInt_FMT ")\n", rbs, cbs, l2grbs, l2gcbs)); in main()
H A Dex23.c25 PetscInt rst, ren, cst, cen, nr, nc, rbs = 1, cbs = 1; in main() local
47 PetscCall(PetscOptionsGetInt(NULL, NULL, "-cbs", &cbs, NULL)); in main()
95 PetscCall(MatSetBlockSizes(A, rbs, cbs)); in main()
161 PetscCall(MatSetBlockSizes(B, rbs, cbs)); in main()
919 PetscCall(MatSetValuesBlocked(A, lm / rbs, idxs1, ln / cbs, idxs2, vals, ADD_VALUES)); in main()
920 PetscCall(MatSetValuesBlocked(B, lm / rbs, idxs1, ln / cbs, idxs2, vals, ADD_VALUES)); in main()
/petsc/src/dm/dt/fe/tests/ !
H A Dex1.c11 PetscInt cbs; /* Number of cells in an integration block */ member
20 options->cbs = 8; in ProcessOptions()
26 …-cbs", "The number of cells in an integration block", "ex1.c", options->cbs, &options->cbs, NULL)); in ProcessOptions()
180 static PetscErrorCode TestIntegration(DM dm, PetscInt cbs, PetscInt its) in TestIntegration() argument
188 PetscInt depth, cStart, cEnd, cell, chunkSize = cbs, Nch = 0, Nf, f, totDim, i, k; in TestIntegration()
252 static PetscErrorCode TestIntegration2(DM dm, PetscInt cbs, PetscInt its) in TestIntegration2() argument
304 PetscCall(TestIntegration(dm, ctx.cbs, ctx.its)); in main()
305 PetscCall(TestIntegration2(dm, ctx.cbs, ctx.its)); in main()
H A Dex2.c233 static PetscErrorCode TestIntegration(DM dm, PetscInt cbs, PetscInt its) in TestIntegration() argument
241 PetscInt depth, cStart, cEnd, cell, chunkSize = cbs, Nch = 0, Nf, f, totDim, i, k; in TestIntegration()
/petsc/src/binding/petsc4py/src/petsc4py/PETSc/ !
H A Dpetscmat.pxi766 cdef PetscInt rbs = 0, cbs = 0, m = 0, n = 0, M = 0, N = 0
767 Mat_Sizes(size, bsize, &rbs, &cbs, &m, &n, &M, &N)
769 if cbs == PETSC_DECIDE: cbs = rbs
771 Sys_Layout(ccomm, cbs, &n, &N)
776 CHKERR(MatSetBlockSizes(mat, rbs, cbs))
924 cdef PetscInt rbs=1, cbs=1
925 if blocked: CHKERR(MatGetBlockSizes(A, &rbs, &cbs))
927 if cbs < 1: cbs = 1
936 if ni*nj*rbs*cbs != nv: raise ValueError(
951 cdef PetscInt rbs=1, cbs=1
[all …]
H A DMat.pyx580 cdef PetscInt rbs = 0, cbs = 0, m = 0, n = 0, M = 0, N = 0
581 Mat_Sizes(size, bsize, &rbs, &cbs, &m, &n, &M, &N)
584 if cbs != PETSC_DECIDE:
585 CHKERR(MatSetBlockSizes(self.mat, rbs, cbs))
625 cdef PetscInt cbs = asInt(col_bsize)
626 CHKERR(MatSetBlockSizes(self.mat, rbs, cbs))
1059 cdef PetscInt rbs = 0, cbs = 0, m = 0, n = 0, M = 0, N = 0
1060 Mat_Sizes(size, bsize, &rbs, &cbs, &m, &n, &M, &N)
1062 if cbs == PETSC_DECIDE: cbs = rbs
1064 Sys_Layout(ccomm, cbs, &n, &N)
[all …]
H A DDMSwarm.pyx1061 cdef PetscInt cbs = 0
1066 CHKERR(DMSwarmGetFieldInfo(self.dm, ccoord, &cbs, NULL))
1067 cdef object moments = oarray_r(empty_r(asInt(cbs) + 2), NULL, &mom)
H A Dlibpetsc4py.pyx855 cdef PetscInt rbs = -1, cbs = -1
857 CHKERR(PetscLayoutGetBlockSize(mat.cmap, &cbs))
859 if cbs == -1: cbs = rbs
861 CHKERR(PetscLayoutSetBlockSize(mat.cmap, cbs))
/petsc/src/binding/petsc4py/test/ !
H A Dtest_mat_aij.py54 rbs, cbs = BS
56 cbs = cbs or 1
58 rbs = cbs = BS or 1
61 self.vals = N.array(range(1, 1 + len(self.adjy) * rbs * cbs), dtype=sdt)
62 self.vals.shape = (-1, rbs, cbs)
66 colsz = (m * n * cbs, None)
256 rbs, cbs = self.A.getBlockSizes()
257 if rbs != cbs:
420 rbs, cbs = bsizes
422 rbs = cbs = bsizes
[all …]
/petsc/src/ksp/pc/impls/mg/ !
H A Dgdsw.c223 PetscInt NG, ns, n, i, c, rbs, cbs[2]; in PCMGGDSWCreateCoarseSpace_Private() local
253 PetscCall(ISGetMinMax(GidxMult, &cbs[0], &cbs[1])); in PCMGGDSWCreateCoarseSpace_Private()
254 PetscCall(PetscGlobalMinMaxInt(PetscObjectComm((PetscObject)A), cbs, cbs)); in PCMGGDSWCreateCoarseSpace_Private()
297 PetscCall(MatSetBlockSizes(*cspace, rbs, cbs[0] == cbs[1] ? cbs[0] : 1)); in PCMGGDSWCreateCoarseSpace_Private()
/petsc/src/mat/impls/shell/ !
H A Dshellcnv.c10 PetscInt cst, cen, Nbs, mbs, nbs, rbs, cbs; in MatConvert_Shell() local
25 PetscCall(MatGetBlockSizes(mat, &rbs, &cbs)); in MatConvert_Shell()
27 nbs = n / cbs; in MatConvert_Shell()
28 Nbs = N / cbs; in MatConvert_Shell()
29 cst = cst / cbs; in MatConvert_Shell()
/petsc/src/mat/impls/is/ !
H A Dmatis.c640 PetscInt rbs, cbs, lc, dr, dc, oc, str, stc, nnz, i, jd, jo, cum; in MatConvert_XAIJ_IS() local
651 PetscCall(MatGetBlockSizes(A, &rbs, &cbs)); in MatConvert_XAIJ_IS()
653 if (reuse != MAT_REUSE_MATRIX && cong && rbs == cbs) { in MatConvert_XAIJ_IS()
677 PetscCall(ISLocalToGlobalMappingGetBlockSize(cl2g, &cbs)); in MatConvert_XAIJ_IS()
680 PetscCall(ISCreateBlock(comm, cbs, nc / cbs, cidx, PETSC_USE_POINTER, &cols)); in MatConvert_XAIJ_IS()
710 cbs = 1; /* We cannot guarantee the off-process matrix will respect the column block size */ in MatConvert_XAIJ_IS()
753 PetscCall(PetscMalloc1((dc + oc) / cbs, &aux)); in MatConvert_XAIJ_IS()
754 for (i = 0; i < dc / cbs; i++) aux[i] = i + stc / cbs; in MatConvert_XAIJ_IS()
755 for (i = 0; i < oc / cbs; i++) aux[i + dc / cbs] = garray[i]; in MatConvert_XAIJ_IS()
756 PetscCall(ISCreateBlock(comm, cbs, (dc + oc) / cbs, aux, PETSC_OWN_POINTER, &is)); in MatConvert_XAIJ_IS()
[all …]
/petsc/src/mat/tests/output/ !
H A Dex23_3.out248 rows=11, cols=10, rbs=1, cbs=2
254 rows=0, cols=0, rbs=1, cbs=2
260 rows=0, cols=0, rbs=1, cbs=2
266 rows=0, cols=0, rbs=1, cbs=2
272 rows=0, cols=0, rbs=1, cbs=2
H A Dex23_6.out49 rows=12, cols=12, rbs=6, cbs=3
89 rows=12, cols=12, rbs=6, cbs=3
129 rows=12, cols=12, rbs=6, cbs=3
169 rows=12, cols=12, rbs=6, cbs=3
H A Dex23_5.out314 rows=12, cols=12, rbs=2, cbs=1
320 rows=0, cols=0, rbs=2, cbs=1
326 rows=0, cols=0, rbs=2, cbs=1
332 rows=0, cols=0, rbs=2, cbs=1
338 rows=0, cols=0, rbs=2, cbs=1
344 rows=0, cols=0, rbs=2, cbs=1
/petsc/src/ksp/ksp/tutorials/output/ !
H A Dex43_6.out58 rows=27, cols=27, rbs=3, cbs=1
65 rows=27, cols=27, rbs=3, cbs=1
72 rows=108, cols=108, rbs=3, cbs=1
/petsc/src/dm/impls/swarm/ !
H A Dswarmpic.c1314 PetscInt *cbs, *bs, *uOff, *uOff_x; in DMProjectFieldLocal_Swarm() local
1345 PetscCall(PetscMalloc2(Nfc, &coordVals, Nfc, &cbs)); in DMProjectFieldLocal_Swarm()
1346 …for (PetscInt i = 0; i < Nfc; ++i) PetscCall(DMSwarmGetField(dm, coordFields[i], &cbs[i], NULL, (v… in DMProjectFieldLocal_Swarm()
1381 … for (PetscInt b = 0; b < cbs[i]; ++b, ++off) xr[off] = coordVals[i][pindices[p] * cbs[i] + b]; in DMProjectFieldLocal_Swarm()
1408 …for (PetscInt i = 0; i < Nfc; ++i) PetscCall(DMSwarmRestoreField(dm, coordFields[i], &cbs[i], NULL… in DMProjectFieldLocal_Swarm()
1414 PetscCall(PetscFree2(coordVals, cbs)); in DMProjectFieldLocal_Swarm()
/petsc/src/mat/interface/ !
H A Dmatrix.c1118 PetscInt rows, cols, rbs, cbs; in MatView() local
1165 PetscCall(MatGetBlockSizes(mat, &rbs, &cbs)); in MatView()
1166 if (rbs != 1 || cbs != 1) { in MatView()
1167 …cbs) PetscCall(PetscViewerASCIIPrintf(viewer, "rows=%" PetscInt_FMT ", cols=%" PetscInt_FMT ", rbs… in MatView()
2086 PetscInt rbs, cbs, M, N, i; in MatSetValuesBlocked() local
2087 PetscCall(MatGetBlockSizes(mat, &rbs, &cbs)); in MatSetValuesBlocked()
2091 …cbs < N, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Column block %" PetscInt_FMT " contains an in… in MatSetValuesBlocked()
2102 PetscInt i, j, bs, cbs; in MatSetValuesBlocked() local
2104 PetscCall(MatGetBlockSizes(mat, &bs, &cbs)); in MatSetValuesBlocked()
2105 if ((m * bs + n * cbs) <= (PetscInt)PETSC_STATIC_ARRAY_LENGTH(buf)) { in MatSetValuesBlocked()
[all …]
/petsc/src/mat/impls/sbaij/mpi/ !
H A Dmpisbaij.c2803 PetscInt m, N, i, rstart, nnz, Ii, bs, cbs; in MatCreateMPIMatConcatenateSeqMat_MPISBAIJ() local
2815 PetscCall(MatGetBlockSizes(inmat, &bs, &cbs)); in MatCreateMPIMatConcatenateSeqMat_MPISBAIJ()
2817 Nbs = N / cbs; in MatCreateMPIMatConcatenateSeqMat_MPISBAIJ()
2818 if (n == PETSC_DECIDE) PetscCall(PetscSplitOwnershipBlock(comm, cbs, &n, &N)); in MatCreateMPIMatConcatenateSeqMat_MPISBAIJ()
2819 nbs = n / cbs; in MatCreateMPIMatConcatenateSeqMat_MPISBAIJ()
2845 PetscCall(MatSetBlockSizes(*outmat, bs, cbs)); in MatCreateMPIMatConcatenateSeqMat_MPISBAIJ()
2853 PetscCall(MatGetBlockSizes(inmat, &bs, &cbs)); in MatCreateMPIMatConcatenateSeqMat_MPISBAIJ()
/petsc/src/mat/impls/aij/mpi/ !
H A Dmpiaij.c191 static PetscErrorCode MatSetBlockSizes_MPIAIJ(Mat M, PetscInt rbs, PetscInt cbs) in MatSetBlockSizes_MPIAIJ() argument
197 PetscCall(MatSetBlockSizes(mat->A, rbs, cbs)); in MatSetBlockSizes_MPIAIJ()
3118 PetscInt cbs; in ISGetSeqIS_Private() local
3119 PetscCall(ISGetBlockSize(iscol, &cbs)); in ISGetSeqIS_Private()
3121 PetscCall(ISSetBlockSize(iscol_local, cbs)); in ISGetSeqIS_Private()
3524 PetscInt i, m, n, rstart, row, rend, nz, j, bs, cbs; in MatCreateSubMatrix_MPIAIJ_SameRowDist() local
3599 PetscCall(ISGetBlockSize(iscol, &cbs)); in MatCreateSubMatrix_MPIAIJ_SameRowDist()
3600 PetscCall(ISSetBlockSize(iscol_sub, cbs)); in MatCreateSubMatrix_MPIAIJ_SameRowDist()
3668 PetscCall(ISGetBlockSize(iscol, &cbs)); in MatCreateSubMatrix_MPIAIJ_SameRowDist()
3672 PetscCall(MatSetBlockSizes(M, bs, cbs)); in MatCreateSubMatrix_MPIAIJ_SameRowDist()
[all …]
/petsc/src/mat/impls/baij/mpi/ !
H A Dmpibaij.c3579 PetscInt m, N, i, rstart, nnz, Ii, bs, cbs; in MatCreateMPIMatConcatenateSeqMat_MPIBAIJ() local
3591 PetscCall(MatGetBlockSizes(inmat, &bs, &cbs)); in MatCreateMPIMatConcatenateSeqMat_MPIBAIJ()
3593 Nbs = N / cbs; in MatCreateMPIMatConcatenateSeqMat_MPIBAIJ()
3594 if (n == PETSC_DECIDE) PetscCall(PetscSplitOwnershipBlock(comm, cbs, &n, &N)); in MatCreateMPIMatConcatenateSeqMat_MPIBAIJ()
3595 nbs = n / cbs; in MatCreateMPIMatConcatenateSeqMat_MPIBAIJ()
3619 PetscCall(MatSetBlockSizes(*outmat, bs, cbs)); in MatCreateMPIMatConcatenateSeqMat_MPIBAIJ()
3628 PetscCall(MatGetBlockSizes(inmat, &bs, &cbs)); in MatCreateMPIMatConcatenateSeqMat_MPIBAIJ()
/petsc/src/mat/impls/aij/seq/ !
H A Daij.c2442 PetscInt rbs, cbs; in MatCreateSubMatrix_SeqAIJ() local
2446 PetscCall(ISGetBlockSize(iscol, &cbs)); in MatCreateSubMatrix_SeqAIJ()
2447 PetscCall(MatSetBlockSizes(C, rbs, cbs)); in MatCreateSubMatrix_SeqAIJ()
2501 PetscInt rbs, cbs; in MatCreateSubMatrix_SeqAIJ() local
2505 PetscCall(ISGetBlockSize(iscol, &cbs)); in MatCreateSubMatrix_SeqAIJ()
2506 if (rbs > 1 || cbs > 1) PetscCall(MatSetBlockSizes(C, rbs, cbs)); in MatCreateSubMatrix_SeqAIJ()

12