| /petsc/src/mat/impls/sbaij/mpi/ |
| H A D | sbaijov.c | 13 …PetscInt i, N = C->cmap->N, bs = C->rmap->bs, M = C->rmap->N, Mbs = M / bs, *nidx, isz, iov; in MatIncreaseOverlap_MPISBAIJ() local
37 PetscCall(PetscMalloc1(Mbs + 1, &nidx)); in MatIncreaseOverlap_MPISBAIJ()
38 PetscCall(PetscBTCreate(Mbs, &table)); /* for column search */ in MatIncreaseOverlap_MPISBAIJ()
42 PetscCall(ISCreateStride(PETSC_COMM_SELF, Mbs, 0, 1, &is_row[0])); in MatIncreaseOverlap_MPISBAIJ()
90 PetscCall(PetscBTMemzero(Mbs, table)); in MatIncreaseOverlap_MPISBAIJ()
101 for (brow = 0; brow < Mbs; brow++) { in MatIncreaseOverlap_MPISBAIJ()
165 PetscInt Mbs, i, j, k, *odata1, *odata2; in MatIncreaseOverlap_MPISBAIJ_Once() local
183 Mbs = c->Mbs; in MatIncreaseOverlap_MPISBAIJ_Once()
193 PetscCall(PetscMalloc2(len, &table, (Mbs / PETSC_BITS_PER_BYTE + 1) * len, &t_p)); in MatIncreaseOverlap_MPISBAIJ_Once()
194 for (i = 0; i < len; i++) table[i] = t_p + (Mbs / PETSC_BITS_PER_BYTE + 1) * i; in MatIncreaseOverlap_MPISBAIJ_Once()
[all …]
|
| H A D | mpisbaij.c | 587 …] < baij->Mbs, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Block indexed row too large %" PetscInt… in MatSetValuesBlocked_MPISBAIJ()
1901 PetscInt i, mbs, Mbs; in MatMPISBAIJSetPreallocation_MPISBAIJ() local
1918 Mbs = B->rmap->N / bs; in MatMPISBAIJSetPreallocation_MPISBAIJ()
1924 b->Mbs = Mbs; in MatMPISBAIJSetPreallocation_MPISBAIJ()
2425 a->Mbs = oldmat->Mbs; in MatDuplicate_MPISBAIJ()
2513 PetscInt i, bs, mbs, *bi, *bj, brow, j, ncols, krow, kcol, col, row, Mbs, bcol; in MatGetRowMaxAbs_MPISBAIJ() local
2530 Mbs = a->Mbs; in MatGetRowMaxAbs_MPISBAIJ()
2539 PetscCall(PetscCalloc1(bs * Mbs, &work)); in MatGetRowMaxAbs_MPISBAIJ()
|
| /petsc/src/mat/tests/ |
| H A D | ex92.c | 11 …nt bs = 1, mbs = 10, ov = 1, i, j, k, *rows, *cols, nd = 2, *idx, rstart, rend, sz, M, N, Mbs; in main() local
45 Mbs = M / bs; in main()
54 for (i = 0; i < Mbs; i++) { in main()
66 cols[0] = bs * (PetscInt)(PetscRealPart(rval) * Mbs); in main()
118 idx[j * bs] = bs * (PetscInt)(PetscRealPart(rval) * Mbs); in main()
|
| H A D | ex54.c | 8 … bs = 1, m = 11, ov = 1, i, j, k, *rows, *cols, nd = 5, *idx, rstart, rend, sz, mm, nn, M, N, Mbs; in main() local
52 Mbs = M / bs; in main()
63 cols[0] = bs * (int)(PetscRealPart(rval) * Mbs); in main()
109 idx[j * bs] = bs * (int)(PetscRealPart(rval) * Mbs); in main()
|
| /petsc/src/mat/impls/baij/mpi/ |
| H A D | baijov.c | 63 PetscInt Mbs, **rbuf, row, msz, **outdat, **ptr; in MatIncreaseOverlap_MPIBAIJ_Once() local
75 Mbs = c->Mbs; in MatIncreaseOverlap_MPIBAIJ_Once()
166 …PetscCall(PetscCalloc5(imax, &table, imax, &data, imax, &isz, Mbs * imax, &d_p, (Mbs / PETSC_BITS_… in MatIncreaseOverlap_MPIBAIJ_Once()
169 table[i] = t_p + (Mbs / PETSC_BITS_PER_BYTE + 1) * i; in MatIncreaseOverlap_MPIBAIJ_Once()
170 data[i] = d_p + (Mbs)*i; in MatIncreaseOverlap_MPIBAIJ_Once()
411 PetscInt val, max1, max2, Mbs, no_malloc = 0, *tmp, new_estimate, ctr; in MatIncreaseOverlap_MPIBAIJ_Receive() local
416 Mbs = c->Mbs; in MatIncreaseOverlap_MPIBAIJ_Receive()
432 if (c->Mbs) max1 = ct * (a->nz + b->nz) / c->Mbs; in MatIncreaseOverlap_MPIBAIJ_Receive()
439 PetscCall(PetscBTCreate(Mbs, &xtable)); in MatIncreaseOverlap_MPIBAIJ_Receive()
450 PetscCall(PetscBTMemzero(Mbs, xtable)); in MatIncreaseOverlap_MPIBAIJ_Receive()
[all …]
|
| H A D | mpibaij.h | 11 PetscInt Mbs, Nbs; /* number block rows/cols in matrix; M/bs, N/bs */ \
|
| H A D | mpibaij.c | 428 …] < baij->Mbs, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Block indexed row too large %" PetscInt… in MatSetValuesBlocked_MPIBAIJ()
577 …(im[i] < baij->Mbs, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Row too large: row %" PetscInt_FMT… in MatSetValuesBlocked_MPIBAIJ_HT()
2680 b->Mbs = B->rmap->N / bs; in MatMPIBAIJSetPreallocation_MPIBAIJ()
3174 a->Mbs = oldmat->Mbs; in MatDuplicate_MPIBAIJ()
3451 …(im[i] < baij->Mbs, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Row too large, row %" PetscInt_FMT… in matmpibaijsetvaluesblocked_()
|
| /petsc/src/mat/impls/aij/mpi/ |
| H A D | mpihashmat.h | 183 a->Mbs = A->rmap->N / bs; in MatSetUp_MPI_Hash()
|
| /petsc/src/mat/impls/aij/mpi/mumps/ |
| H A D | mumps.c | 1524 PetscInt *pjcns_w, Mbs = 0; in MatConvertToTriples_nest_xaij() local
1619 Mbs += size / bs; in MatConvertToTriples_nest_xaij()
1690 if (Mbs != A->rmap->N) { in MatConvertToTriples_nest_xaij()
1698 PetscCall(PetscMUMPSIntCast(Mbs, &mumps->id.nblk)); in MatConvertToTriples_nest_xaij()
1700 PetscCall(PetscMalloc1(Mbs + 1, &mumps->id.blkptr)); in MatConvertToTriples_nest_xaij()
|