| /petsc/src/vec/is/is/utils/ |
| H A D | isblock.c | 25 PetscInt isz, len, i, j, ival, bbs; in ISCompressIndicesGeneral() local
60 isz = 0; in ISCompressIndicesGeneral()
72 PetscCall(PetscHMapISet(gid1_lid1, ival + 1, isz + 1)); in ISCompressIndicesGeneral()
73 isz++; in ISCompressIndicesGeneral()
77 if (!PetscBTLookupSet(table, ival)) nidx[isz++] = ival; in ISCompressIndicesGeneral()
83 PetscCall(PetscMalloc1(isz, &nidx)); in ISCompressIndicesGeneral()
89 …PetscCheck(tt-- <= isz, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "index greater than array-dim"); in ISCompressIndicesGeneral()
94 PetscCheck(j == isz, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "table error: jj != isz"); in ISCompressIndicesGeneral()
95 …PetscCall(ISCreateGeneral(PetscObjectComm((PetscObject)is_in[i]), isz, nidx, PETSC_OWN_POINTER, is… in ISCompressIndicesGeneral()
97 …PetscCall(ISCreateGeneral(PetscObjectComm((PetscObject)is_in[i]), isz, nidx, PETSC_COPY_VALUES, is… in ISCompressIndicesGeneral()
|
| /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
98 isz = nis; in MatIncreaseOverlap_MPISBAIJ()
104 if (!PetscBTLookupSet(table, brow)) nidx[isz++] = brow; in MatIncreaseOverlap_MPISBAIJ()
111 PetscCall(ISCreateGeneral(PETSC_COMM_SELF, isz, nidx, PETSC_COPY_VALUES, is_new + i)); in MatIncreaseOverlap_MPISBAIJ()
164 PetscInt idx, isz, col, *n, *data1, **data1_start, *data2, *data2_i, *data, *data_i, len; in MatIncreaseOverlap_MPISBAIJ_Once() local
412 isz = data[1 + i]; in MatIncreaseOverlap_MPISBAIJ_Once()
415 if (!PetscBTLookupSet(table_i, col)) data_i[isz++] = col; in MatIncreaseOverlap_MPISBAIJ_Once()
417 data[1 + i] = isz; in MatIncreaseOverlap_MPISBAIJ_Once()
470 …PetscInt row, mbs, Mbs, *nidx_i, col, col_max, isz, isz0, *ai, *aj, *bi, *bj, *garray, rstart… in MatIncreaseOverlap_MPISBAIJ_Local() local
492 isz = 0; in MatIncreaseOverlap_MPISBAIJ_Local()
[all …]
|
| /petsc/src/mat/impls/baij/mpi/ |
| H A D | baijov.c | 64 PetscInt *ctr, *tmp, *isz, *isz1, **xdata, **rbuf2, *d_p; in MatIncreaseOverlap_MPIBAIJ_Once() local
166 …PetscCall(PetscCalloc5(imax, &table, imax, &data, imax, &isz, Mbs * imax, &d_p, (Mbs / PETSC_BITS_… in MatIncreaseOverlap_MPIBAIJ_Once()
185 isz_i = isz[i]; in MatIncreaseOverlap_MPIBAIJ_Once()
206 isz[i] = isz_i; in MatIncreaseOverlap_MPIBAIJ_Once()
228 PetscCall(MatIncreaseOverlap_MPIBAIJ_Local(C, imax, table, isz, data)); in MatIncreaseOverlap_MPIBAIJ_Once()
289 isz_i = isz[is_no]; in MatIncreaseOverlap_MPIBAIJ_Once()
296 isz[is_no] = isz_i; in MatIncreaseOverlap_MPIBAIJ_Once()
303 PetscCall(ISCreateGeneral(iscomms[i], isz[i], data[i], PETSC_COPY_VALUES, is + i)); in MatIncreaseOverlap_MPIBAIJ_Once()
320 PetscCall(PetscFree5(table, data, isz, d_p, t_p)); in MatIncreaseOverlap_MPIBAIJ_Once()
343 …IncreaseOverlap_MPIBAIJ_Local(Mat C, PetscInt imax, PetscBT *table, PetscInt *isz, PetscInt **data) in MatIncreaseOverlap_MPIBAIJ_Local() argument
[all …]
|
| /petsc/src/mat/impls/aij/mpi/ |
| H A D | mpiov.c | 234 PetscInt *isz, isz_i, i, j, is_id, data_size; in MatIncreaseOverlap_MPIAIJ_Receive_Scalable() local
241 PetscCall(PetscMalloc1(nidx, &isz)); in MatIncreaseOverlap_MPIAIJ_Receive_Scalable()
245 isz[i] = lsize; in MatIncreaseOverlap_MPIAIJ_Receive_Scalable()
251 …aycpy(PetscSafePointerPlusOffset(indices_temp, i * (max_lsize + nrecvs)), indices_i_temp, isz[i])); in MatIncreaseOverlap_MPIAIJ_Receive_Scalable()
260 isz_i = isz[is_id]; in MatIncreaseOverlap_MPIAIJ_Receive_Scalable()
265 isz[is_id] = isz_i; in MatIncreaseOverlap_MPIAIJ_Receive_Scalable()
270 isz_i = isz[i]; in MatIncreaseOverlap_MPIAIJ_Receive_Scalable()
275 PetscCall(PetscFree(isz)); in MatIncreaseOverlap_MPIAIJ_Receive_Scalable()
506 PetscInt *ctr, *tmp, *isz, **xdata, **rbuf2; in MatIncreaseOverlap_MPIAIJ_Once() local
615 PetscCall(PetscCalloc4(imax, &table, imax, &data, imax, &isz, M_BPB_imax, &t_p)); in MatIncreaseOverlap_MPIAIJ_Once()
[all …]
|
| /petsc/src/mat/impls/sbaij/seq/ |
| H A D | sbaij2.c | 11 PetscInt brow, i, j, k, l, mbs, n, *nidx, isz, bcol, bcol_max, start, end, *ai, *aj, bs; in MatIncreaseOverlap_SeqSBAIJ() local
26 isz = 0; in MatIncreaseOverlap_SeqSBAIJ()
39 nidx[isz++] = brow; in MatIncreaseOverlap_SeqSBAIJ()
50 for (l = k; l < isz; l++) PetscCall(PetscBTSet(table_in, nidx[l])); in MatIncreaseOverlap_SeqSBAIJ()
52 n = isz; /* length of the updated is[i] */ in MatIncreaseOverlap_SeqSBAIJ()
60 nidx[isz++] = bcol; in MatIncreaseOverlap_SeqSBAIJ()
71 if (!PetscBTLookupSet(table_out, brow)) nidx[isz++] = brow; in MatIncreaseOverlap_SeqSBAIJ()
78 PetscCall(ISCreateBlock(PETSC_COMM_SELF, bs, isz, nidx, PETSC_COPY_VALUES, is + i)); in MatIncreaseOverlap_SeqSBAIJ()
|
| /petsc/src/mat/impls/aij/seq/ |
| H A D | aij.c | 2737 PetscInt row, i, j, k, l, ll, m, n, *nidx, isz, val; in MatIncreaseOverlap_SeqAIJ() local
2754 isz = 0; in MatIncreaseOverlap_SeqAIJ()
2764 if (!PetscBTLookupSet(table, idx[j] / bs)) nidx[isz++] = idx[j] / bs; in MatIncreaseOverlap_SeqAIJ()
2771 n = isz; in MatIncreaseOverlap_SeqAIJ()
2779 if (!PetscBTLookupSet(table, val)) nidx[isz++] = val; in MatIncreaseOverlap_SeqAIJ()
2784 PetscCall(ISCreateBlock(PETSC_COMM_SELF, bs, isz, nidx, PETSC_COPY_VALUES, is + i)); in MatIncreaseOverlap_SeqAIJ()
2788 if (!PetscBTLookupSet(table, idx[j])) nidx[isz++] = idx[j]; in MatIncreaseOverlap_SeqAIJ()
2795 n = isz; in MatIncreaseOverlap_SeqAIJ()
2802 if (!PetscBTLookupSet(table, val)) nidx[isz++] = val; in MatIncreaseOverlap_SeqAIJ()
2806 PetscCall(ISCreateGeneral(PETSC_COMM_SELF, isz, nidx, PETSC_COPY_VALUES, is + i)); in MatIncreaseOverlap_SeqAIJ()
|
| /petsc/src/mat/impls/baij/seq/ |
| H A D | baij2.c | 16 PetscInt row, i, j, k, l, m, n, *nidx, isz, val, ival; in MatIncreaseOverlap_SeqBAIJ() local
34 isz = 0; in MatIncreaseOverlap_SeqBAIJ()
45 if (!PetscBTLookupSet(table, ival)) nidx[isz++] = ival; in MatIncreaseOverlap_SeqBAIJ()
52 n = isz; in MatIncreaseOverlap_SeqBAIJ()
59 if (!PetscBTLookupSet(table, val)) nidx[isz++] = val; in MatIncreaseOverlap_SeqBAIJ()
63 PetscCall(ISCreateBlock(PETSC_COMM_SELF, bs, isz, nidx, PETSC_COPY_VALUES, is + i)); in MatIncreaseOverlap_SeqBAIJ()
|
| /petsc/src/ksp/pc/impls/fieldsplit/ |
| H A D | fieldsplit.c | 2365 IS iszl, isz; in PCFieldSplitRestrictIS_FieldSplit() local
2379 PetscCall(ISCreateGeneral(comm, localsize, indcz, PETSC_OWN_POINTER, &isz)); in PCFieldSplitRestrictIS_FieldSplit()
2380 PetscCall(PCFieldSplitRestrictIS(subpc, isz)); in PCFieldSplitRestrictIS_FieldSplit()
2381 PetscCall(ISDestroy(&isz)); in PCFieldSplitRestrictIS_FieldSplit()
|