| /petsc/src/dm/impls/plex/tests/ |
| H A D | dmplexcomputecellgeometryfem.F90 | 106 function quad_v(rs)
107 PetscReal, intent(in) :: rs(2)
112 r = rs(1)
113 s = rs(2)
119 function quad_J(rs)
120 PetscReal, intent(in) :: rs(2)
126 r = rs(1)
127 s = rs(2)
|
| /petsc/src/mat/graphops/color/utils/ |
| H A D | bipartite.c | 7 PetscInt i, j, rs, re, cs, ce, cn; in MatColoringCreateBipartiteGraph() local
16 PetscCall(MatGetOwnershipRange(m, &rs, &re)); in MatColoringCreateBipartiteGraph()
20 for (i = rs; i < re; i++) { in MatColoringCreateBipartiteGraph()
28 for (i = rs; i < re; i++) { in MatColoringCreateBipartiteGraph()
|
| /petsc/src/mat/utils/ |
| H A D | getcolv.c | 29 PetscInt i, j, nz, N, Rs, Re, rs, re; in MatGetColumnVector() local
40 PetscCall(VecGetOwnershipRange(yy, &rs, &re)); in MatGetColumnVector()
41 …rs && Re == re, PETSC_COMM_SELF, PETSC_ERR_ARG_INCOMP, "Matrix %" PetscInt_FMT " %" PetscInt_FMT "… in MatGetColumnVector()
56 if (idx[j] == col) y[i - rs] = v[j]; in MatGetColumnVector()
|
| /petsc/src/mat/impls/aij/seq/lusol/ |
| H A D | lusol.c | 220 int i, rs, re; in MatLUFactorNumeric_LUSOL() local
256 rs = a->i[i]; in MatLUFactorNumeric_LUSOL()
259 while (rs < re) { in MatLUFactorNumeric_LUSOL()
260 if (a->a[rs] != 0.0) { in MatLUFactorNumeric_LUSOL()
262 lusol->indr[nz] = a->j[rs] + 1; in MatLUFactorNumeric_LUSOL()
263 lusol->data[nz] = a->a[rs]; in MatLUFactorNumeric_LUSOL()
266 rs++; in MatLUFactorNumeric_LUSOL()
|
| /petsc/src/ksp/ksp/impls/gmres/pgmres/ |
| H A D | pgmres.c | 238 PetscScalar *hh, *cc, *ss, *rs; in KSPPGMRESUpdateHessenberg() local
247 rs = RS(0); /* right-hand side of least squares system */ in KSPPGMRESUpdateHessenberg()
253 hapbnd = PetscMin(PetscAbsScalar(hh[it + 1] / rs[it]), pgmres->haptol); in KSPPGMRESUpdateHessenberg()
292 rs[it + 1] = -ss[it] * rs[it]; in KSPPGMRESUpdateHessenberg()
293 rs[it] = PetscConj(cc[it]) * rs[it]; in KSPPGMRESUpdateHessenberg()
294 *res = PetscAbsScalar(rs[it + 1]); in KSPPGMRESUpdateHessenberg()
|
| /petsc/src/dm/dt/interface/ftn-custom/ |
| H A D | zdsf.c | 22 PETSC_EXTERN void petscdssetriemannsolver_(PetscDS *prob, PetscInt *f, void (*rs)(PetscInt *, Petsc… in petscdssetriemannsolver_()
24 …tscObject)*prob, PETSC_FORTRAN_CALLBACK_CLASS, &riemannsolver, (PetscFortranCallbackFn *)rs, NULL); in petscdssetriemannsolver_()
|
| /petsc/src/mat/impls/aij/seq/ |
| H A D | aijfact.c | 233 PetscReal rs; in MatLUFactorNumeric_SeqAIJ() local
248 rs = -PetscAbsScalar(d) - PetscRealPart(d); in MatLUFactorNumeric_SeqAIJ()
251 for (j = 0; j < nz; j++) rs += PetscAbsScalar(v[j]); in MatLUFactorNumeric_SeqAIJ()
252 if (rs > sctx.shift_top) sctx.shift_top = rs; in MatLUFactorNumeric_SeqAIJ()
309 rs = 0.0; in MatLUFactorNumeric_SeqAIJ()
316 rs += PetscAbsScalar(pv[j]); in MatLUFactorNumeric_SeqAIJ()
325 rs += PetscAbsScalar(pv[j]); in MatLUFactorNumeric_SeqAIJ()
328 sctx.rs = rs; in MatLUFactorNumeric_SeqAIJ()
409 PetscReal rs = 0.0; in MatLUFactorNumeric_SeqAIJ_inplace() local
427 rs = -PetscAbsScalar(d) - PetscRealPart(d); in MatLUFactorNumeric_SeqAIJ_inplace()
[all …]
|
| H A D | inode.c | 1211 PetscReal rs; in MatLUFactorNumeric_SeqAIJ_Inode() local
1227 rs = -PetscAbsScalar(d) - PetscRealPart(d); in MatLUFactorNumeric_SeqAIJ_Inode()
1230 for (j = 0; j < nz; j++) rs += PetscAbsScalar(v[j]); in MatLUFactorNumeric_SeqAIJ_Inode()
1231 if (rs > sctx.shift_top) sctx.shift_top = rs; in MatLUFactorNumeric_SeqAIJ_Inode()
1334 rs = 0.0; in MatLUFactorNumeric_SeqAIJ_Inode()
1341 rs += PetscAbsScalar(pv[j]); in MatLUFactorNumeric_SeqAIJ_Inode()
1350 rs += PetscAbsScalar(pv[j]); in MatLUFactorNumeric_SeqAIJ_Inode()
1354 sctx.rs = rs; in MatLUFactorNumeric_SeqAIJ_Inode()
1426 rs = 0.0; in MatLUFactorNumeric_SeqAIJ_Inode()
1434 rs += PetscAbsScalar(pc1[j]); in MatLUFactorNumeric_SeqAIJ_Inode()
[all …]
|
| /petsc/src/ksp/ksp/impls/gmres/pipefgmres/ |
| H A D | pipefgmres.c | 380 PetscScalar *hh, *cc, *ss, *rs; in KSPPIPEFGMRESUpdateHessenberg() local
389 rs = RS(0); /* right-hand side of least squares system */ in KSPPIPEFGMRESUpdateHessenberg()
395 hapbnd = PetscMin(PetscAbsScalar(hh[it + 1] / rs[it]), pipefgmres->haptol); in KSPPIPEFGMRESUpdateHessenberg()
434 rs[it + 1] = -ss[it] * rs[it]; in KSPPIPEFGMRESUpdateHessenberg()
435 rs[it] = PetscConj(cc[it]) * rs[it]; in KSPPIPEFGMRESUpdateHessenberg()
436 *res = PetscAbsScalar(rs[it + 1]); in KSPPIPEFGMRESUpdateHessenberg()
|
| /petsc/src/mat/impls/baij/seq/ |
| H A D | baijfact.c | 463 PetscReal rs; in MatLUFactorNumeric_SeqBAIJ_1() local
476 rs = -PetscAbsScalar(d) - PetscRealPart(d); in MatLUFactorNumeric_SeqBAIJ_1()
479 for (j = 0; j < nz; j++) rs += PetscAbsScalar(v[j]); in MatLUFactorNumeric_SeqBAIJ_1()
480 if (rs > sctx.shift_top) sctx.shift_top = rs; in MatLUFactorNumeric_SeqBAIJ_1()
537 rs = 0.0; in MatLUFactorNumeric_SeqBAIJ_1()
544 rs += PetscAbsScalar(pv[j]); in MatLUFactorNumeric_SeqBAIJ_1()
553 rs += PetscAbsScalar(pv[j]); in MatLUFactorNumeric_SeqBAIJ_1()
556 sctx.rs = rs; in MatLUFactorNumeric_SeqBAIJ_1()
757 PetscReal rs; in MatCholeskyFactorNumeric_SeqBAIJ_N() local
833 rs = 0.0; in MatCholeskyFactorNumeric_SeqBAIJ_N()
[all …]
|
| /petsc/src/mat/impls/sbaij/seq/ |
| H A D | sbaijfact.c | 1197 PetscReal rs; in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace() local
1283 rs = 0.0; in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1289 rs += PetscAbsScalar(rtmp[*bcol]); in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1294 sctx.rs = rs; in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1355 PetscReal rs; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering() local
1446 rs = 0.0; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1453 rs += PetscAbsScalar(ba[j]); in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1462 sctx.rs = rs; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1506 PetscReal rs; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace() local
1579 rs = 0.0; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace()
[all …]
|
| /petsc/src/mat/tests/ |
| H A D | ex70.c | 205 Vec r, l, rs, ls; in main() local
516 PetscCall(MatCreateVecs(T2, &rs, &ls)); in main()
517 PetscCall(VecCopy(r, rs)); in main()
521 PetscCall(MatDiagonalScale(T2, rs, rs)); in main()
524 PetscCall(MatDiagonalScale(T2, ls, rs)); in main()
535 PetscCall(VecDestroy(&rs)); in main()
|
| /petsc/src/binding/petsc4py/test/ |
| H A D | test_vec.py | 247 rs, re = A.getOwnershipRange()
250 for i in range(rs, re):
261 for i in range(rs, re):
|
| H A D | test_mat_aij.py | 282 rs, re = self.A.getOwnershipRange()
284 rows = N.array(range(rs, re), dtype=PETSc.IntType)
308 rs, re = self.A.getOwnershipRange()
310 rows = N.array(range(rs, re), dtype=PETSc.IntType)
|
| /petsc/src/ksp/ksp/impls/gmres/ |
| H A D | gmres.c | 37 PetscInt hh, hes, rs, cc; in KSPSetUp_GMRES() local
45 rs = (max_k + 2); in KSPSetUp_GMRES()
48 …PetscCall(PetscCalloc5(hh, &gmres->hh_origin, hes, &gmres->hes_origin, rs, &gmres->rs_origin, cc, … in KSPSetUp_GMRES()
|
| /petsc/src/binding/petsc4py/src/petsc4py/PETSc/ |
| H A D | petscmat.pxi | 1023 cdef PetscInt rs=0, re=ni-1
1028 CHKERR(MatGetOwnershipRange(A, &rs, &re))
1029 rs //= rbs; re //= rbs
1030 nm = re - rs
1053 irow = m[k] if m!=NULL else rs+k
|
| /petsc/src/mat/impls/hypre/ |
| H A D | mhypre.c | 1618 PetscInt i, rs, re, cs, ce, bs; local
1624 rs = A->rmap->rstart;
1629 PetscCallHYPRE(HYPRE_IJMatrixCreate(hA->comm, rs, re - 1, cs, ce - 1, &hA->ij));
1634 …rs, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Inconsistent local rows: IJMatrix [%" PetscHYPRE_BigInt_FMT …
|
| /petsc/include/petsc/private/ |
| H A D | matimpl.h | 769 PetscReal rs; /* active row sum of abs(off-diagonals) */ member
791 PetscReal _rs = sctx->rs; in MatPivotCheck_nz()
809 PetscReal _rs = sctx->rs; in MatPivotCheck_pd()
|
| /petsc/src/mat/impls/baij/mpi/ |
| H A D | mpibaij.c | 1131 PetscInt header[4], M, N, m, rs, cs, bs, cnt, i, j, ja, jb, k, l; in MatView_MPIBAIJ_Binary() local
1143 rs = mat->rmap->rstart; in MatView_MPIBAIJ_Binary()
1161 PetscCall(PetscViewerBinaryWriteAll(viewer, rowlens, m, rs, M, PETSC_INT)); in MatView_MPIBAIJ_Binary()
3230 PetscInt *rowidxs, *colidxs, rs, cs, ce; in MatLoad_MPIBAIJ_Binary() local
3262 PetscCall(PetscLayoutGetRange(mat->rmap, &rs, NULL)); in MatLoad_MPIBAIJ_Binary()
3317 PetscInt row = rs + i, s = rowidxs[i], e = rowidxs[i + 1]; in MatLoad_MPIBAIJ_Binary()
|
| /petsc/src/mat/impls/aij/mpi/ |
| H A D | mpiaij.c | 1191 PetscInt header[4], M, N, m, rs, cs, cnt, i, ja, jb; in MatView_MPIAIJ_Binary() local
1204 rs = mat->rmap->rstart; in MatView_MPIAIJ_Binary()
1220 PetscCall(PetscViewerBinaryWriteAll(viewer, rowlens, m, rs, M, PETSC_INT)); in MatView_MPIAIJ_Binary()
7408 const PetscInt rs = C->rmap->rstart; in MatProductSymbolic_MPIAIJBACKEND() local
7414 if (gr < rs || gr >= re) ncoo_o += nz; /* this row is offproc */ in MatProductSymbolic_MPIAIJBACKEND()
7457 const PetscInt rs = C->rmap->rstart; in MatProductSymbolic_MPIAIJBACKEND() local
7464 if (gr < rs || gr >= re) { /* this is an offproc row */ in MatProductSymbolic_MPIAIJBACKEND()
7523 const PetscInt rs = C->rmap->rstart; in MatProductSymbolic_MPIAIJBACKEND() local
7531 const PetscInt gr = i + rs; in MatProductSymbolic_MPIAIJBACKEND()
7548 if (gr >= rs && gr < re) { /* local rows */ in MatProductSymbolic_MPIAIJBACKEND()
|