Home
last modified time | relevance | path

Searched refs:nrow (Results 1 – 25 of 31) sorted by relevance

12

/petsc/src/mat/graphops/order/
H A Dspqmd.c9 PetscInt i, *deg, *marker, *rchset, *nbrhd, *qsize, *qlink, nofsub, *iperm, nrow, *perm; in MatGetOrdering_QMD() local
14 PetscCall(MatGetRowIJ(mat, 1, PETSC_TRUE, PETSC_TRUE, &nrow, &ia, &ja, &done)); in MatGetOrdering_QMD()
17 PetscCall(PetscMalloc1(nrow, &perm)); in MatGetOrdering_QMD()
18 PetscCall(PetscMalloc5(nrow, &iperm, nrow, &deg, nrow, &marker, nrow, &rchset, nrow, &nbrhd)); in MatGetOrdering_QMD()
19 PetscCall(PetscMalloc2(nrow, &qsize, nrow, &qlink)); in MatGetOrdering_QMD()
21 …PetscCall(SPARSEPACKgenqmd(&nrow, ia, ja, perm, iperm, deg, marker, rchset, nbrhd, qsize, qlink, &… in MatGetOrdering_QMD()
26 for (i = 0; i < nrow; i++) perm[i]--; in MatGetOrdering_QMD()
27 PetscCall(ISCreateGeneral(PETSC_COMM_SELF, nrow, perm, PETSC_COPY_VALUES, row)); in MatGetOrdering_QMD()
28 PetscCall(ISCreateGeneral(PETSC_COMM_SELF, nrow, perm, PETSC_OWN_POINTER, col)); in MatGetOrdering_QMD()
H A Dwbm.c37 PetscInt *perm, nrow, ncol, nnz, liw, *iw, ldw; in MatGetOrdering_WBM() local
44 PetscCall(MatGetRowIJ(mat, 1, PETSC_TRUE, PETSC_TRUE, &nrow, &ia, &ja, &done)); in MatGetOrdering_WBM()
45 ncol = nrow; in MatGetOrdering_WBM()
46 nnz = ia[nrow]; in MatGetOrdering_WBM()
51 liw = 4 * nrow + ncol; in MatGetOrdering_WBM()
55 liw = 2 * nrow + 2 * ncol; in MatGetOrdering_WBM()
59 liw = 8 * nrow + 2 * ncol + nnz; in MatGetOrdering_WBM()
63 liw = 3 * nrow + 2 * ncol; in MatGetOrdering_WBM()
67 liw = 3 * nrow + 2 * ncol; in MatGetOrdering_WBM()
68 ldw = nrow + 2 * ncol + nnz; in MatGetOrdering_WBM()
[all …]
H A Dsp1wd.c9 PetscInt i, *mask, *xls, nblks, *xblk, *ls, nrow, *perm; in MatGetOrdering_1WD() local
14 PetscCall(MatGetRowIJ(mat, 1, PETSC_TRUE, PETSC_TRUE, &nrow, &ia, &ja, &done)); in MatGetOrdering_1WD()
17 PetscCall(PetscMalloc5(nrow, &mask, nrow + 1, &xls, nrow, &ls, nrow + 1, &xblk, nrow, &perm)); in MatGetOrdering_1WD()
18 PetscCall(SPARSEPACKgen1wd(&nrow, ia, ja, mask, &nblks, xblk, perm, xls, ls)); in MatGetOrdering_1WD()
21 for (i = 0; i < nrow; i++) perm[i]--; in MatGetOrdering_1WD()
23 PetscCall(ISCreateGeneral(PETSC_COMM_SELF, nrow, perm, PETSC_COPY_VALUES, row)); in MatGetOrdering_1WD()
24 PetscCall(ISCreateGeneral(PETSC_COMM_SELF, nrow, perm, PETSC_COPY_VALUES, col)); in MatGetOrdering_1WD()
H A Dspnd.c9 PetscInt i, *mask, *xls, *ls, nrow, *perm; in MatGetOrdering_ND() local
15 PetscCall(MatGetRowIJ(mat, 1, PETSC_TRUE, PETSC_TRUE, &nrow, &ia, &ja, &done)); in MatGetOrdering_ND()
18 PetscCall(MatGetRowIJ(B, 1, PETSC_TRUE, PETSC_TRUE, &nrow, &ia, &ja, &done)); in MatGetOrdering_ND()
21 PetscCall(PetscMalloc4(nrow, &mask, nrow, &perm, nrow + 1, &xls, nrow, &ls)); in MatGetOrdering_ND()
22 PetscCall(SPARSEPACKgennd(&nrow, ia, ja, mask, perm, xls, ls)); in MatGetOrdering_ND()
31 for (i = 0; i < nrow; i++) perm[i]--; in MatGetOrdering_ND()
33 PetscCall(ISCreateGeneral(PETSC_COMM_SELF, nrow, perm, PETSC_COPY_VALUES, row)); in MatGetOrdering_ND()
34 PetscCall(ISCreateGeneral(PETSC_COMM_SELF, nrow, perm, PETSC_COPY_VALUES, col)); in MatGetOrdering_ND()
H A Dsprcm.c9 PetscInt i, *mask, *xls, nrow, *perm; in MatGetOrdering_RCM() local
14 PetscCall(MatGetRowIJ(mat, 1, PETSC_TRUE, PETSC_TRUE, &nrow, &ia, &ja, &done)); in MatGetOrdering_RCM()
17 PetscCall(PetscMalloc3(nrow, &mask, nrow, &perm, 2 * nrow, &xls)); in MatGetOrdering_RCM()
18 PetscCall(SPARSEPACKgenrcm(&nrow, ia, ja, perm, mask, xls)); in MatGetOrdering_RCM()
22 for (i = 0; i < nrow; i++) perm[i]--; in MatGetOrdering_RCM()
23 PetscCall(ISCreateGeneral(PETSC_COMM_SELF, nrow, perm, PETSC_COPY_VALUES, row)); in MatGetOrdering_RCM()
24 PetscCall(ISCreateGeneral(PETSC_COMM_SELF, nrow, perm, PETSC_COPY_VALUES, col)); in MatGetOrdering_RCM()
/petsc/src/mat/graphops/order/metisnd/
H A Dmetisnd.c10 PetscInt i, j, iptr, ival, nrow, *xadj, *adjncy, *perm, *iperm; in MatGetOrdering_METISND() local
18 PetscCall(MatGetRowIJ(mat, 0, PETSC_TRUE, PETSC_TRUE, &nrow, &ia, &ja, &done)); in MatGetOrdering_METISND()
21 PetscCall(MatGetRowIJ(B, 0, PETSC_TRUE, PETSC_TRUE, &nrow, &ia, &ja, &done)); in MatGetOrdering_METISND()
45 PetscCall(PetscMalloc4(nrow + 1, &xadj, ia[nrow], &adjncy, nrow, &perm, nrow, &iperm)); in MatGetOrdering_METISND()
50 for (j = 0; j < nrow; j++) { in MatGetOrdering_METISND()
57 …status = METIS_NodeND(&nrow, (idx_t *)xadj, (idx_t *)adjncy, NULL, options, (idx_t *)perm, (idx_t … in MatGetOrdering_METISND()
78 PetscCall(ISCreateGeneral(PETSC_COMM_SELF, nrow, perm, PETSC_COPY_VALUES, row)); in MatGetOrdering_METISND()
79 PetscCall(ISCreateGeneral(PETSC_COMM_SELF, nrow, perm, PETSC_COPY_VALUES, col)); in MatGetOrdering_METISND()
/petsc/src/mat/impls/localref/
H A Dmlocalref.c12 #define IndexSpaceGet(buf, nrow, ncol, irowm, icolm) \ argument
14 if (nrow + ncol > (PetscInt)PETSC_STATIC_ARRAY_LENGTH(buf)) { \
15 PetscCall(PetscMalloc2(nrow, &irowm, ncol, &icolm)); \
18 icolm = &buf[nrow]; \
22 #define IndexSpaceRestore(buf, nrow, ncol, irowm, icolm) \ argument
24 … if (nrow + ncol > (PetscInt)PETSC_STATIC_ARRAY_LENGTH(buf)) PetscCall(PetscFree2(irowm, icolm)); \
35 static PetscErrorCode MatSetValuesBlockedLocal_LocalRef_Block(Mat A, PetscInt nrow, const PetscInt … in MatSetValuesBlockedLocal_LocalRef_Block() argument
41 if (!nrow || !ncol) PetscFunctionReturn(PETSC_SUCCESS); in MatSetValuesBlockedLocal_LocalRef_Block()
42 IndexSpaceGet(buf, nrow, ncol, irowm, icolm); in MatSetValuesBlockedLocal_LocalRef_Block()
43 PetscCall(ISLocalToGlobalMappingApplyBlock(A->rmap->mapping, nrow, irow, irowm)); in MatSetValuesBlockedLocal_LocalRef_Block()
[all …]
/petsc/src/mat/graphops/order/amd/
H A Damd.c21 PetscInt nrow, *perm; in MatGetOrdering_AMD() local
33 PetscCall(MatGetRowIJ(mat, 0, PETSC_FALSE, PETSC_TRUE, &nrow, &ia, &ja, &done)); in MatGetOrdering_AMD()
51 PetscCall(PetscMalloc1(nrow, &perm)); in MatGetOrdering_AMD()
52 status = amd_AMD_order(nrow, ia, ja, perm, Control, Info); in MatGetOrdering_AMD()
71 PetscCall(ISCreateGeneral(PETSC_COMM_SELF, nrow, perm, PETSC_COPY_VALUES, row)); in MatGetOrdering_AMD()
72 PetscCall(ISCreateGeneral(PETSC_COMM_SELF, nrow, perm, PETSC_OWN_POINTER, col)); in MatGetOrdering_AMD()
/petsc/src/mat/tests/
H A Dex240.c9 PetscInt N, mx = 5, my = 4, i, j, nc, nrow, n, ncols, rstart, *colors, *map; in main() local
49 PetscCall(MatGetLocalSize(A, &nrow, NULL)); in main()
51 PetscCall(PetscCalloc1(nrow * nc, &cm)); in main()
52 for (i = 0; i < nrow; i++) { in main()
56 cm[i + nrow * colors[cols[j]]] = vals[j]; in main()
63 for (i = 0; i < nrow; i++) { in main()
64 …all(PetscSynchronizedPrintf(MPI_COMM_WORLD, "%12.4e ", (double)PetscAbsScalar(cm[i + nrow * j]))); in main()
H A Dex73.c135 PetscInt rstart, i, *nzd, *nzo, nzl, nzmax = 0, *ncols, nrow, j; in main() local
168 nrow = 2 * (i + rstart); in main()
169 PetscCall(MatSetValues(J, 1, &nrow, 2 * nzl, ncols, nvals, INSERT_VALUES)); in main()
170 nrow = 2 * (i + rstart) + 1; in main()
171 PetscCall(MatSetValues(J, 1, &nrow, 2 * nzl, ncols, nvals, INSERT_VALUES)); in main()
/petsc/src/mat/impls/dense/mpi/
H A Dmmdense.c58 PetscInt *nrow, *ncol, *w1, *w3, *w4, start, inrqr; in MatCreateSubMatrices_MPIDense_Local() local
87 …PetscCall(PetscMalloc5(ismax, (PetscInt ***)&irow, ismax, (PetscInt ***)&icol, ismax, &nrow, ismax… in MatCreateSubMatrices_MPIDense_Local()
91 PetscCall(ISGetLocalSize(isrow[i], &nrow[i])); in MatCreateSubMatrices_MPIDense_Local()
108 jmax = nrow[i]; in MatCreateSubMatrices_MPIDense_Local()
187 jmax = nrow[i]; in MatCreateSubMatrices_MPIDense_Local()
277 …PetscCheck(!(submats[i]->rmap->n != nrow[i]) && !(submats[i]->cmap->n != ncol[i]), PETSC_COMM_SELF… in MatCreateSubMatrices_MPIDense_Local()
285 PetscCall(MatSetSizes(submats[i], nrow[i], ncol[i], nrow[i], ncol[i])); in MatCreateSubMatrices_MPIDense_Local()
301 m = nrow[i]; in MatCreateSubMatrices_MPIDense_Local()
326 jmax = nrow[i]; in MatCreateSubMatrices_MPIDense_Local()
349 m = nrow[is_no]; in MatCreateSubMatrices_MPIDense_Local()
[all …]
/petsc/src/ksp/ksp/tutorials/
H A Dex61f.F9085 PetscInt nrow
97 nrow = n*n
98 ncol = nrow
119 …PetscCallA(MatCreateSeqAIJ(PETSC_COMM_SELF, nrow, ncol, nz_per_row, PETSC_NULL_INTEGER_ARRAY, col_…
120 …PetscCallA(VecCreateSeqWithArray(PETSC_COMM_SELF, 1, nrow, PETSC_NULL_SCALAR_ARRAY, col_f_vecb, ie…
/petsc/src/mat/impls/htool/
H A Dhtool.cxx165 PetscInt nrow, m, i; in MatCreateSubMatrices_Htool() local
173 PetscCall(ISGetLocalSize(irow[i], &nrow)); in MatCreateSubMatrices_Htool()
177 …if (scall != MAT_REUSE_MATRIX) PetscCall(MatCreateDense(PETSC_COMM_SELF, nrow, m, nrow, m, nullptr… in MatCreateSubMatrices_Htool()
184 it = std::lower_bound(idxr, idxr + nrow, A->rmap->rstart); in MatCreateSubMatrices_Htool()
185 if (it != idxr + nrow && *it == A->rmap->rstart) { /* rmap->rstart in IS? */ in MatCreateSubMatrices_Htool()
186 …if (std::distance(idxr, it) + A->rmap->n <= nrow) { /* long enough IS to store the local diagonal … in MatCreateSubMatrices_Htool()
198 …scInt k = 0; k < A->rmap->n; ++k) PetscCall(PetscArraycpy(ptr + (m + k) * nrow + m, copy + k * A->… in MatCreateSubMatrices_Htool()
201 … a->wrapper->copy_submatrix(nrow, m, idxr, idxc, ptr); /* vertical block B from above */ in MatCreateSubMatrices_Htool()
205 PetscCall(MatDenseSetLDA(B, nrow)); in MatCreateSubMatrices_Htool()
206 … PetscCall(MatCreateDense(PETSC_COMM_SELF, m, A->rmap->n, m, A->rmap->n, ptr + m * nrow, &BT)); in MatCreateSubMatrices_Htool()
[all …]
/petsc/src/dm/impls/plex/
H A Dplexmetric.c1205 PetscInt nrow; in DMPlexMetricEnforceSPD() local
1207 if (isotropic) nrow = 1; in DMPlexMetricEnforceSPD()
1208 else nrow = dim; in DMPlexMetricEnforceSPD()
1216 PetscCall(DMPlexMetricModify_Private(nrow, h_min, h_max, a_max, vmet, vdet)); in DMPlexMetricEnforceSPD()
1336 PetscInt nrow; in DMPlexMetricNormalize() local
1338 if (isotropic) nrow = 1; in DMPlexMetricNormalize()
1339 else nrow = dim; in DMPlexMetricNormalize()
1349 if (restrictSizes) PetscCall(DMPlexMetricModify_Private(nrow, h_min, h_max, a_max, Mp, detM)); in DMPlexMetricNormalize()
1636 PetscInt dim, vStart, vEnd, nrow; in DMPlexMetricIntersection() local
1640 if (isotropic) nrow = 1; in DMPlexMetricIntersection()
[all …]
/petsc/src/dm/impls/da/hypre/
H A Dmhyp.c25 static PetscErrorCode MatSetValuesLocal_HYPREStruct_3d(Mat mat, PetscInt nrow, const PetscInt irow[… in MatSetValuesLocal_HYPREStruct_3d() argument
34 for (i = 0; i < nrow; i++) { in MatSetValuesLocal_HYPREStruct_3d()
64 static PetscErrorCode MatZeroRowsLocal_HYPREStruct_3d(Mat mat, PetscInt nrow, const PetscInt irow[]… in MatZeroRowsLocal_HYPREStruct_3d() argument
75 for (i = 0; i < nrow; i++) { in MatZeroRowsLocal_HYPREStruct_3d()
339 static PetscErrorCode MatSetValuesLocal_HYPRESStruct_3d(Mat mat, PetscInt nrow, const PetscInt irow… in MatSetValuesLocal_HYPRESStruct_3d() argument
360 for (i = 0; i < nrow; i++) { in MatSetValuesLocal_HYPRESStruct_3d()
399 for (i = 0; i < nrow; i++) { in MatSetValuesLocal_HYPRESStruct_3d()
442 static PetscErrorCode MatZeroRowsLocal_HYPRESStruct_3d(Mat mat, PetscInt nrow, const PetscInt irow[… in MatZeroRowsLocal_HYPRESStruct_3d() argument
472 for (i = 0; i < nrow; i++) { in MatZeroRowsLocal_HYPRESStruct_3d()
483 for (i = 0; i < nrow; i++) { in MatZeroRowsLocal_HYPRESStruct_3d()
/petsc/src/mat/impls/baij/mpi/
H A Dbaijov.c619 PetscInt *nrow, *ncol, start; local
655 …PetscCall(PetscMalloc5(ismax, (PetscInt ***)&irow, ismax, (PetscInt ***)&icol, ismax, &nrow, ismax…
676 PetscCall(ISGetLocalSize(isrow[i], &nrow[i]));
677 if (colflag && nrow[i] == c->Mbs) {
691 …PetscCheck(subc->mbs == nrow[i] && subc->nbs == ncol[i], PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "Cann…
758 jmax = nrow[i];
844 jmax = nrow[i];
1002 for (PetscInt i = 0; i < ismax; i++) jcnt += nrow[i];
1005 …for (PetscInt i = 1; i < ismax; i++) lens[i] = PetscSafePointerPlusOffset(lens[i - 1], nrow[i - 1]…
1010 jmax = nrow[i];
[all …]
/petsc/src/mat/impls/aij/mpi/
H A Dmpiov.c1221 PetscInt nrow, ncol, start; in MatCreateSubMatrices_MPIAIJ_SingleIS_Local() local
1259 PetscCall(ISGetLocalSize(isrow[0], &nrow)); in MatCreateSubMatrices_MPIAIJ_SingleIS_Local()
1279 PetscCall(PetscMalloc1(nrow, &row2proc)); in MatCreateSubMatrices_MPIAIJ_SingleIS_Local()
1284 for (PetscInt j = 0; j < nrow; j++) { in MatCreateSubMatrices_MPIAIJ_SingleIS_Local()
1346 for (PetscInt j = 0; j < nrow; j++) { /* parse the indices of each IS */ in MatCreateSubMatrices_MPIAIJ_SingleIS_Local()
1496 PetscCall(PetscCalloc1(nrow, &lens)); in MatCreateSubMatrices_MPIAIJ_SingleIS_Local()
1499 for (PetscInt j = 0; j < nrow; j++) { in MatCreateSubMatrices_MPIAIJ_SingleIS_Local()
1538 PetscCall(PetscHMapICreateWithSize(nrow, &rmap)); in MatCreateSubMatrices_MPIAIJ_SingleIS_Local()
1539 for (PetscInt j = 0; j < nrow; j++) { in MatCreateSubMatrices_MPIAIJ_SingleIS_Local()
1549 for (PetscInt j = 0; j < nrow; j++) rmap[irow[j]] = j; in MatCreateSubMatrices_MPIAIJ_SingleIS_Local()
[all …]
/petsc/src/mat/impls/sbaij/mpi/
H A Dmpisbaij.c207 nrow = ailen[brow]; \
212 high = nrow; \
229 …MatSeqXAIJReallocateAIJ(A, a->mbs, bs2, nrow, brow, bcol, rmax, aa, ai, aj, rp, ap, aimax, a->none…
230 N = nrow++ - 1; \
238 ailen[brow] = nrow; \
247 nrow = bilen[brow]; \
252 high = nrow; \
269 …MatSeqXAIJReallocateAIJ(B, b->mbs, bs2, nrow, brow, bcol, rmax, ba, bi, bj, rp, ap, bimax, b->none…
270 N = nrow++ - 1; \
278 bilen[brow] = nrow; \
[all …]
/petsc/src/mat/impls/sbaij/seq/
H A Dsbaij.c513 PetscInt *rp, k, low, high, t, row, nrow, i, col, l, *aj = a->j; in MatGetValues_SeqSBAIJ() local
529 nrow = ailen[brow]; in MatGetValues_SeqSBAIJ()
540 high = nrow; in MatGetValues_SeqSBAIJ()
578 PetscInt *rp, k, low, high, t, ii, jj, row, nrow, i, col, l, rmax, N, lastcol = -1; in MatSetValuesBlocked_SeqSBAIJ() local
595 nrow = ailen[row]; in MatSetValuesBlocked_SeqSBAIJ()
597 high = nrow; in MatSetValuesBlocked_SeqSBAIJ()
610 else high = nrow; in MatSetValuesBlocked_SeqSBAIJ()
648 …MatSeqXAIJReallocateAIJ(A, a->mbs, bs2, nrow, row, col, rmax, aa, ai, aj, rp, ap, imax, nonew, Mat… in MatSetValuesBlocked_SeqSBAIJ()
649 N = nrow++ - 1; in MatSetValuesBlocked_SeqSBAIJ()
669 ailen[row] = nrow; in MatSetValuesBlocked_SeqSBAIJ()
[all …]
/petsc/src/mat/impls/baij/seq/
H A Dbaij.c1292 PetscInt *rp, k, low, high, t, ii, jj, row, nrow, i, col, l, N, m = *mm, n = *nn; in matsetvaluesblocked4_() local
1305 nrow = ailen[row]; in matsetvaluesblocked4_()
1307 high = nrow; in matsetvaluesblocked4_()
1311 else high = nrow; in matsetvaluesblocked4_()
1329 N = nrow++ - 1; in matsetvaluesblocked4_()
1345 ailen[row] = nrow; in matsetvaluesblocked4_()
1360 PetscInt *rp, k, low, high, t, row, nrow, i, col, l, N, n = *nn, m = *mm; in matsetvalues4_() local
1372 nrow = ailen[brow]; in matsetvalues4_()
1374 high = nrow; in matsetvalues4_()
1382 else high = nrow; in matsetvalues4_()
[all …]
/petsc/src/mat/impls/blockmat/seq/
H A Dblockmat.c217 PetscInt *rp, k, low, high, t, row, nrow, i, col, l, lastcol = -1; in MatSetValues_BlockMat() local
233 nrow = ailen[brow]; in MatSetValues_BlockMat()
235 high = nrow; in MatSetValues_BlockMat()
248 else high = nrow; in MatSetValues_BlockMat()
262 …MatSeqXAIJReallocateAIJ(A, a->mbs, 1, nrow, brow, bcol, rmax, aa, ai, aj, rp, ap, imax, nonew, Mat… in MatSetValues_BlockMat()
263 N = nrow++ - 1; in MatSetValues_BlockMat()
279 ailen[brow] = nrow; in MatSetValues_BlockMat()
/petsc/src/mat/impls/aij/seq/spqr/
H A Daijspqr.c84 C->nrow = (size_t)AT->cmap->n; in MatWrapCholmod_SPQR_seqaij()
169 …Call(PetscArraycpy(&v[j * lda], &(((PetscScalar *)Y_handle->x)[j * Y_handle->d]), Y_handle->nrow)); in MatMatSolve_SPQR()
227 …Call(PetscArraycpy(&v[j * lda], &(((PetscScalar *)Y_handle->x)[j * Y_handle->d]), Y_handle->nrow)); in MatMatSolveTranspose_SPQR()
/petsc/src/mat/impls/aij/seq/cholmod/
H A Daijcholmod.c36 C->nrow = (size_t)A->cmap->n; in MatWrapCholmod_seqaij()
/petsc/src/mat/impls/sell/seq/
H A Dsell.c1439 PetscInt i, shift, row_in_slice, row, nrow, *cp, lastcol, j, k; in MatAssemblyEnd_SeqSELL() local
1459 nrow = a->rlen[row]; /* number of nonzeros in row */ in MatAssemblyEnd_SeqSELL()
1465 if (nrow > 0) { /* nonempty row */ in MatAssemblyEnd_SeqSELL()
1467 …lastcol = cp[a->sliceheight * (nrow - 1) + row_in_slice]; /* use the index from the last nonzero a… in MatAssemblyEnd_SeqSELL()
1479 for (k = nrow; k < (a->sliidx[i + 1] - shift) / a->sliceheight; ++k) { in MatAssemblyEnd_SeqSELL()
1536 PetscInt shift, i, k, l, low, high, t, ii, row, col, nrow; in MatSetValues_SeqSELL() local
1552 nrow = a->rlen[row]; in MatSetValues_SeqSELL()
1554 high = nrow; in MatSetValues_SeqSELL()
1569 else high = nrow; in MatSetValues_SeqSELL()
1592 …MatSeqXSELLReallocateSELL(A, A->rmap->n, 1, nrow, a->sliidx, a->sliceheight, row / a->sliceheight,… in MatSetValues_SeqSELL()
[all …]
/petsc/src/mat/interface/
H A Dmatrix.c2214 PetscErrorCode MatGetValuesLocal(Mat mat, PetscInt nrow, const PetscInt irow[], PetscInt ncol, cons… in MatGetValuesLocal() argument
2220 if (!nrow || !ncol) PetscFunctionReturn(PETSC_SUCCESS); /* no values to retrieve */ in MatGetValuesLocal()
2229 if (mat->ops->getvalueslocal) PetscUseTypeMethod(mat, getvalueslocal, nrow, irow, ncol, icol, y); in MatGetValuesLocal()
2232 if ((nrow + ncol) <= (PetscInt)PETSC_STATIC_ARRAY_LENGTH(buf)) { in MatGetValuesLocal()
2234 icolm = buf + nrow; in MatGetValuesLocal()
2236 PetscCall(PetscMalloc2(nrow, &bufr, ncol, &bufc)); in MatGetValuesLocal()
2242 PetscCall(ISLocalToGlobalMappingApply(mat->rmap->mapping, nrow, irow, irowm)); in MatGetValuesLocal()
2244 PetscCall(MatGetValues(mat, nrow, irowm, ncol, icolm, y)); in MatGetValuesLocal()
2455 PetscErrorCode MatSetValuesLocal(Mat mat, PetscInt nrow, const PetscInt irow[], PetscInt ncol, cons… in MatSetValuesLocal() argument
2461 if (!nrow || !ncol) PetscFunctionReturn(PETSC_SUCCESS); /* no values to insert */ in MatSetValuesLocal()
[all …]

12