Lines Matching refs:cbs
640 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()
759 PetscCall(ISCreateBlock(comm, cbs, 0, NULL, PETSC_OWN_POINTER, &is)); in MatConvert_XAIJ_IS()
769 PetscCall(MatSetBlockSizes(B, rbs, cbs)); in MatConvert_XAIJ_IS()
832 PetscInt i, j, nr, nc, rbs, cbs; in MatConvert_Nest_IS() local
1047 PetscCall(MatGetBlockSizes(A, &rbs, &cbs)); in MatConvert_Nest_IS()
1048 PetscCall(MatSetBlockSizes(B, rbs, cbs)); in MatConvert_Nest_IS()
1325 PetscInt buf[2 * MATIS_MAX_ENTRIES_INSERTION], *rows_l = NULL, *cols_l = NULL, rbs, cbs; in MatSetValuesBlockedLocal_SubMat_IS() local
1330 PetscCall(ISLocalToGlobalMappingGetBlockSize(A->cmap->mapping, &cbs)); in MatSetValuesBlockedLocal_SubMat_IS()
1331 IndexSpaceGet(buf, m * rbs, n * cbs, rows_l, cols_l); in MatSetValuesBlockedLocal_SubMat_IS()
1333 BlockIndicesExpand(n, cols, cbs, cols_l); in MatSetValuesBlockedLocal_SubMat_IS()
1335 PetscCall(ISLocalToGlobalMappingApply(A->cmap->mapping, n * cbs, cols_l, cols_l)); in MatSetValuesBlockedLocal_SubMat_IS()
1336 PetscCall(MatSetValuesLocal_IS(A, m * rbs, rows_l, n * cbs, cols_l, values, addv)); in MatSetValuesBlockedLocal_SubMat_IS()
1337 IndexSpaceRestore(buf, m * rbs, n * cbs, rows_l, cols_l); in MatSetValuesBlockedLocal_SubMat_IS()
1420 PetscInt ll, newloc, irbs, icbs, arbs, acbs, rbs, cbs; in MatCreateSubMatrix_IS() local
1429 cbs = acbs == icbs ? icbs : 1; in MatCreateSubMatrix_IS()
1436 PetscCall(MatSetBlockSizes(*newmat, rbs, cbs)); in MatCreateSubMatrix_IS()
1492 PetscCall(ISLocalToGlobalMappingSetBlockSize(cl2g, cbs)); in MatCreateSubMatrix_IS()
1840 PetscInt rbs, cbs, rows, cols, lrows, lcols; in MatConvert_IS_XAIJ() local
1851 PetscInt rbs, cbs; in MatConvert_IS_XAIJ() local
1854 PetscCall(ISLocalToGlobalMappingGetBlockSize(matis->cmapping, &cbs)); in MatConvert_IS_XAIJ()
1879 nw = nw / cbs; in MatConvert_IS_XAIJ()
1885 PetscCall(ISCreateBlock(PETSC_COMM_SELF, cbs, nw, cidxs, PETSC_USE_POINTER, &cols)); in MatConvert_IS_XAIJ()
1926 PetscCall(ISLocalToGlobalMappingGetBlockSize(matis->cmapping, &cbs)); in MatConvert_IS_XAIJ()
1952 PetscCall(MatSetBlockSizes(MT, rbs, cbs)); in MatConvert_IS_XAIJ()
2003 … == cbs, PetscObjectComm((PetscObject)mat), PETSC_ERR_SUP, "Cannot reuse matrix. Wrong col block s… in MatConvert_IS_XAIJ()
2037 PetscInt rbs, cbs, m, n, M, N; in MatDuplicate_IS() local
2042 PetscCall(ISLocalToGlobalMappingGetBlockSize(mat->cmap->mapping, &cbs)); in MatDuplicate_IS()
2047 PetscCall(MatSetBlockSize(B, rbs == cbs ? rbs : 1)); in MatDuplicate_IS()
2548 PetscInt nr, rbs, nc, cbs; in MatISSetUpScatters_Private() local
2555 PetscCall(ISLocalToGlobalMappingGetBlockSize(is->cmapping, &cbs)); in MatISSetUpScatters_Private()
2575 …PetscCall(ISCreateBlock(PetscObjectComm((PetscObject)A), cbs, nc / cbs, garray, PETSC_USE_POINTER,… in MatISSetUpScatters_Private()
2665 PetscInt nr, rbs, nc, cbs; in MatSetLocalToGlobalMapping_IS() local
2735 PetscCall(ISLocalToGlobalMappingGetBlockSize(is->cmapping, &cbs)); in MatSetLocalToGlobalMapping_IS()
2738 if (nr == nc && cbs == rbs) { in MatSetLocalToGlobalMapping_IS()
2755 PetscCall(PetscLayoutSetBlockSize(A->cmap, cbs)); in MatSetLocalToGlobalMapping_IS()
2767 PetscCall(MatSetBlockSizes(is->A, rbs, cbs)); in MatSetLocalToGlobalMapping_IS()
3249 PetscInt nrg, rbs, cbs; in MatGetLocalSubMatrix_IS() local
3254 PetscCall(ISGetBlockSize(col, &cbs)); in MatGetLocalSubMatrix_IS()
3286 if (ncg % cbs) ncg = cbs * (ncg / cbs + 1); in MatGetLocalSubMatrix_IS()
3295 PetscCall(ISLocalToGlobalMappingSetBlockSize(cl2g, cbs)); in MatGetLocalSubMatrix_IS()
3463 PetscInt rbs, cbs, bs; in MatISGetAssembled_Private() local
3468 PetscCall(ISLocalToGlobalMappingGetBlockSize(a->cmapping, &cbs)); in MatISGetAssembled_Private()
3470 bs = rbs == cbs ? rbs : 1; in MatISGetAssembled_Private()
3604 static PetscErrorCode MatSetBlockSizes_IS(Mat A, PetscInt rbs, PetscInt cbs) in MatSetBlockSizes_IS() argument
3609 if (a->A) PetscCall(MatSetBlockSizes(a->A, rbs, cbs)); in MatSetBlockSizes_IS()