Home
last modified time | relevance | path

Searched refs:idxn (Results 1 – 16 of 16) sorted by relevance

/petsc/src/mat/tests/
H A Dex219f.F907 PetscInt :: n, m, idxm(1), idxn(1), nl1, nl2, zero, one, i
22 idxn(1) = i
25 PetscCallA(MatSetValues(A, one, idxn, one, idxm, v, INSERT_VALUES, ierr))
33 idxn(1) = 0
35 if ((idxn(1) >= nl1) .and. (idxn(1) <= nl2 - 1)) then
37 PetscCallA(MatSetValues(A, one, idxn, one, idxm, v, INSERT_VALUES, ierr))
42 if ((idxn(1) >= nl1) .and. (idxn(1) <= nl2 - 1)) then
43 PetscCallA(MatGetValues(A, one, idxn, one, idxm, v, ierr))
54 idxn(1) = 0
55 PetscCallA(MatGetValuesLocal(A, one, idxn, one, [zero], values, ierr))
[all …]
H A Dex94.c29 PetscInt i, m, n, j, *idxn, M, N, nzp, rstart, rend; in main() local
82 PetscCall(PetscMalloc2(nzp + 1, &idxn, nzp + 1, &a)); in main()
197 idxn[j] = (PetscInt)(PetscRealPart(rval) * PN); in main()
199 PetscCall(MatSetValues(P, 1, &i, nzp, idxn, a, ADD_VALUES)); in main()
284 idxn[j] = (PetscInt)(PetscRealPart(rval) * PN); in main()
286 PetscCall(MatSetValues(P, 1, &i, nzp, idxn, a, ADD_VALUES)); in main()
428 PetscCall(PetscFree2(idxn, a)); in main()
H A Dex62.c28 PetscInt i, j, *idxn, PM, PN = PETSC_DECIDE, rstart, rend; in main() local
81 PetscCall(PetscMalloc2(nzp, &idxn, nzp, &a)); in main()
156 idxn[j] = (PetscInt)(PetscRealPart(rval) * PN); in main()
158 PetscCall(MatSetValues(P, 1, &i, nzp, idxn, a, ADD_VALUES)); in main()
308 PetscCall(PetscFree2(idxn, a)); in main()
/petsc/src/ksp/ksp/tests/
H A Dex54.c11 PetscInt idxn[3] = {0, 1, 2}; in fill() local
34 PetscCall(MatSetValues(m, 4, idxm1, 3, idxn, values1, INSERT_VALUES)); in fill()
40 PetscCall(MatSetValues(m, 4, idxm2, 3, idxn, values2, INSERT_VALUES)); in fill()
49 PetscCall(MatSetValues(m, 4, idxm, 3, idxn, values, INSERT_VALUES)); in fill()
58 PetscCall(MatSetValues(m, 4, idxm, 3, idxn, values, INSERT_VALUES)); in fill()
/petsc/src/mat/utils/
H A Dmatstash.c258 …luesRow_Private(MatStash *stash, PetscInt row, PetscInt n, const PetscInt idxn[], const PetscScala… in MatStashValuesRow_Private() argument
271 space->idy[k] = idxn[i]; in MatStashValuesRow_Private()
296 …luesCol_Private(MatStash *stash, PetscInt row, PetscInt n, const PetscInt idxn[], const PetscScala… in MatStashValuesCol_Private() argument
309 space->idy[k] = idxn[i]; in MatStashValuesCol_Private()
338 …Blocked_Private(MatStash *stash, PetscInt row, PetscInt n, const PetscInt idxn[], const PetscScala… in MatStashValuesRowBlocked_Private() argument
352 space->idy[l] = idxn[i]; in MatStashValuesRowBlocked_Private()
389 …Blocked_Private(MatStash *stash, PetscInt row, PetscInt n, const PetscInt idxn[], const PetscScala… in MatStashValuesColBlocked_Private() argument
403 space->idy[l] = idxn[i]; in MatStashValuesColBlocked_Private()
/petsc/src/mat/impls/sell/mpi/
H A Dmpisell.c248 …ELL(Mat mat, PetscInt m, const PetscInt idxm[], PetscInt n, const PetscInt idxn[], PetscScalar v[]) in MatGetValues_MPISELL() argument
261 if (idxn[j] < 0) continue; /* negative column */ in MatGetValues_MPISELL()
262 …etscCheck(idxn[j] < mat->cmap->N, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Column too large: co… in MatGetValues_MPISELL()
263 if (idxn[j] >= cstart && idxn[j] < cend) { in MatGetValues_MPISELL()
264 col = idxn[j] - cstart; in MatGetValues_MPISELL()
269 PetscCall(PetscHMapIGetWithDefault(sell->colmap, idxn[j] + 1, 0, &col)); in MatGetValues_MPISELL()
272 col = sell->colmap[idxn[j]] - 1; in MatGetValues_MPISELL()
274 if (col < 0 || sell->garray[col] != idxn[j]) *(v + i * n + j) = 0.0; in MatGetValues_MPISELL()
/petsc/src/ksp/ksp/tutorials/
H A Dex74.c141 PetscInt idxm[1] = {0}, idxn[1] = {0}; in main() local
145 PetscCall(MatSetValuesBlocked(A_baij, 1, idxm, 1, idxn, A, INSERT_VALUES)); in main()
/petsc/src/mat/impls/sbaij/mpi/
H A Dmpisbaij.c650 …AIJ(Mat mat, PetscInt m, const PetscInt idxm[], PetscInt n, const PetscInt idxn[], PetscScalar v[]) in MatGetValues_MPISBAIJ() argument
663 if (idxn[j] < 0) continue; /* negative column */ in MatGetValues_MPISBAIJ()
664 …etscCheck(idxn[j] < mat->cmap->N, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Column too large: co… in MatGetValues_MPISBAIJ()
665 if (idxn[j] >= bscstart && idxn[j] < bscend) { in MatGetValues_MPISBAIJ()
666 col = idxn[j] - bscstart; in MatGetValues_MPISBAIJ()
671 PetscCall(PetscHMapIGetWithDefault(baij->colmap, idxn[j] / bs + 1, 0, &data)); in MatGetValues_MPISBAIJ()
674 data = baij->colmap[idxn[j] / bs] - 1; in MatGetValues_MPISBAIJ()
676 if (data < 0 || baij->garray[data / bs] != idxn[j] / bs) *(v + i * n + j) = 0.0; in MatGetValues_MPISBAIJ()
678 col = data + idxn[j] % bs; in MatGetValues_MPISBAIJ()
/petsc/src/mat/interface/
H A Dmatrix.c1532 …es(Mat mat, PetscInt m, const PetscInt idxm[], PetscInt n, const PetscInt idxn[], const PetscScala… in MatSetValues() argument
1539 PetscAssertPointer(idxn, 5); in MatSetValues()
1554 …, (double)PetscRealPart(v[i * n + j]), (double)PetscImaginaryPart(v[i * n + j]), idxm[i], idxn[j]); in MatSetValues()
1556 …%g at matrix entry (%" PetscInt_FMT ",%" PetscInt_FMT ")", (double)v[i * n + j], idxm[i], idxn[j]); in MatSetValues()
1562 …cCheck(idxn[i] < mat->cmap->N, PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Cannot insert in column %" P… in MatSetValues()
1570 PetscUseTypeMethod(mat, setvalues, m, idxm, n, idxn, v, addv); in MatSetValues()
1788 …at mat, PetscInt m, const MatStencil idxm[], PetscInt n, const MatStencil idxn[], const PetscScala… in MatSetValuesStencil() argument
1792 …PetscInt *starts = mat->stencil.starts, *dxm = (PetscInt *)idxm, *dxn = (PetscInt *)idxn, sdim = d… in MatSetValuesStencil()
1799 PetscAssertPointer(idxn, 5); in MatSetValuesStencil()
1903 …at mat, PetscInt m, const MatStencil idxm[], PetscInt n, const MatStencil idxn[], const PetscScala… in MatSetValuesBlockedStencil() argument
[all …]
/petsc/src/mat/impls/dense/mpi/
H A Dmpidense.c133 …se(Mat mat, PetscInt m, const PetscInt idxm[], PetscInt n, const PetscInt idxn[], const PetscScala… in MatSetValues_MPIDense() argument
146 PetscCall(MatSetValues(A->A, 1, &row, n, idxn, PetscSafePointerPlusOffset(v, i * n), addv)); in MatSetValues_MPIDense()
149 if (idxn[j] < 0) continue; in MatSetValues_MPIDense()
150 … PetscCheck(idxn[j] < mat->cmap->N, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Column too large"); in MatSetValues_MPIDense()
151 …PetscCall(MatSetValues(A->A, 1, &row, 1, &idxn[j], PetscSafePointerPlusOffset(v, i + j * m), addv)… in MatSetValues_MPIDense()
157 …PetscCall(MatStashValuesRow_Private(&mat->stash, idxm[i], n, idxn, PetscSafePointerPlusOffset(v, i… in MatSetValues_MPIDense()
159 …PetscCall(MatStashValuesCol_Private(&mat->stash, idxm[i], n, idxn, PetscSafePointerPlusOffset(v, i… in MatSetValues_MPIDense()
166 …nse(Mat mat, PetscInt m, const PetscInt idxm[], PetscInt n, const PetscInt idxn[], PetscScalar v[]) in MatGetValues_MPIDense() argument
178 if (idxn[j] < 0) continue; /* negative column */ in MatGetValues_MPIDense()
179 … PetscCheck(idxn[j] < mat->cmap->N, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Column too large"); in MatGetValues_MPIDense()
[all …]
/petsc/src/mat/impls/baij/mpi/
H A Dmpibaij.c652 …AIJ(Mat mat, PetscInt m, const PetscInt idxm[], PetscInt n, const PetscInt idxn[], PetscScalar v[]) in MatGetValues_MPIBAIJ() argument
665 if (idxn[j] < 0) continue; /* negative column */ in MatGetValues_MPIBAIJ()
666 …etscCheck(idxn[j] < mat->cmap->N, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Column too large: co… in MatGetValues_MPIBAIJ()
667 if (idxn[j] >= bscstart && idxn[j] < bscend) { in MatGetValues_MPIBAIJ()
668 col = idxn[j] - bscstart; in MatGetValues_MPIBAIJ()
673 PetscCall(PetscHMapIGetWithDefault(baij->colmap, idxn[j] / bs + 1, 0, &data)); in MatGetValues_MPIBAIJ()
676 data = baij->colmap[idxn[j] / bs] - 1; in MatGetValues_MPIBAIJ()
678 if (data < 0 || baij->garray[data / bs] != idxn[j] / bs) *(v + i * n + j) = 0.0; in MatGetValues_MPIBAIJ()
680 col = data + idxn[j] % bs; in MatGetValues_MPIBAIJ()
/petsc/src/ts/impls/implicit/irk/
H A Dirk.c132 PetscInt idxm[1] = {0}, idxn[1] = {0}; in TSIRKCreate_Gauss() local
137 PetscCall(MatSetValuesBlocked(A_baij, 1, idxm, 1, idxn, gauss_A, INSERT_VALUES)); in TSIRKCreate_Gauss()
/petsc/doc/manual/
H A Dmat.md77 MatSetValues(Mat A,PetscInt m,const PetscInt idxm[],PetscInt n,const PetscInt idxn[],const PetscSca…
83 MatSetValues(Mat A,PetscInt m,const PetscInt idxm[],PetscInt n,const PetscInt idxn[],const PetscSca…
87 `m*n` into the matrix. The integer indices `idxm` and `idxn`,
95 `idxn[j]` is located in `values[i*n+j]`. To allow the insertion of
/petsc/src/mat/impls/aij/mpi/
H A Dmpiaij.c725 …AIJ(Mat mat, PetscInt m, const PetscInt idxm[], PetscInt n, const PetscInt idxn[], PetscScalar v[]) in MatGetValues_MPIAIJ() argument
738 if (idxn[j] < 0) continue; /* negative column */ in MatGetValues_MPIAIJ()
739 …etscCheck(idxn[j] < mat->cmap->N, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Column too large: co… in MatGetValues_MPIAIJ()
740 if (idxn[j] >= cstart && idxn[j] < cend) { in MatGetValues_MPIAIJ()
741 col = idxn[j] - cstart; in MatGetValues_MPIAIJ()
746 PetscCall(PetscHMapIGetWithDefault(aij->colmap, idxn[j] + 1, 0, &col)); in MatGetValues_MPIAIJ()
749 col = aij->colmap[idxn[j]] - 1; in MatGetValues_MPIAIJ()
751 if ((col < 0) || (aij->garray[col] != idxn[j])) *(v + i * n + j) = 0.0; in MatGetValues_MPIAIJ()
/petsc/src/mat/impls/hypre/
H A Dmhypre.c2104 …HYPRE(Mat A, PetscInt m, const PetscInt idxm[], PetscInt n, const PetscInt idxn[], PetscScalar v[]) argument
2117 …MatrixGetValues(hA->ij, 1, &hn, (HYPRE_BigInt *)&idxm[i], (HYPRE_BigInt *)idxn, (HYPRE_Complex *)(…
/petsc/src/ksp/pc/impls/patch/
H A Dpcpatch.c2176 …te(Mat mat, PetscInt m, const PetscInt idxm[], PetscInt n, const PetscInt idxn[], const PetscScala… in MatSetValues_PCPatch_Private() argument
2189 …PetscCheck(idxm[i] == idxn[i], PetscObjectComm((PetscObject)mat), PETSC_ERR_ARG_WRONG, "Row and co… in MatSetValues_PCPatch_Private()