| /petsc/src/ts/utils/dmplexlandau/ |
| H A D | land_tensors.h | 71 PETSC_DEVICE_FUNC_DECL static void LandauTensor2D(const PetscReal x[], const PetscReal rp, const Pe… in LandauTensor2D() argument
78 rp2 = PetscSqr(rp); in LandauTensor2D()
86 s = mask * 2 * r * rp / l; /* mask for vectorization */ in LandauTensor2D()
99 Ud[0][1] = Ud[1][0] = Uk[0][1] = pi4pow * (zmzp) * (r * i2func - rp * i3func); in LandauTensor2D()
100 Uk[1][1] = Ud[1][1] = -pi4pow * ((r2prp2)*i2func - 2 * r * rp * i3func) * mask; in LandauTensor2D()
101 Uk[0][0] = -pi4pow * (zmzp2 * i3func + r * rp * i1func); in LandauTensor2D()
102 …Uk[1][0] = pi4pow * (zmzp) * (r * i3func - rp * i2func); /* 48 mults + 21 + 21 = 90 mul… in LandauTensor2D()
|
| /petsc/src/ksp/ksp/impls/bcgs/fbcgsr/ |
| H A D | fbcgsr.c | 24 PetscScalar *PETSC_RESTRICT rp, *PETSC_RESTRICT r, *PETSC_RESTRICT p; in KSPSolve_FBCGSR() local
41 PetscCall(VecGetArray(RP, (PetscScalar **)&rp)); in KSPSolve_FBCGSR()
112 tau += r[j] * rp[j]; /* tau <- (r,rp) */ in KSPSolve_FBCGSR()
113 sigma += v[j] * rp[j]; /* sigma <- (v,rp) */ in KSPSolve_FBCGSR()
143 xi4 += t[j] * rp[j]; /* xi4 <- (t,rp) */ in KSPSolve_FBCGSR()
|
| /petsc/share/petsc/datafiles/meshes/ |
| H A D | testcase3D.cas | 12 (rp-seg? #t)
14 (rp-inviscid? #f)
15 (rp-lam? #f)
16 (rp-atm? #f)
17 (rp-hvac? #f)
18 (rp-les? #f)
19 (rp-les-sgske? #f)
20 (rp-sa? #f)
21 (rp-des? #f)
22 (rp-ke? #f)
[all …]
|
| /petsc/src/mat/impls/sbaij/mpi/ |
| H A D | mpisbaij.c | 204 rp = aj + ai[brow]; \
215 if (rp[t] > bcol) high = t; \
219 if (rp[_i] > bcol) break; \
220 if (rp[_i] == bcol) { \
229 …MatSeqXAIJReallocateAIJ(A, a->mbs, bs2, nrow, brow, bcol, rmax, aa, ai, aj, rp, ap, aimax, a->none…
232 PetscCall(PetscArraymove(rp + _i + 1, rp + _i, N - _i + 1)); \
235 rp[_i] = bcol; \
244 rp = bj + bi[brow]; \
255 if (rp[t] > bcol) high = t; \
259 if (rp[_i] > bcol) break; \
[all …]
|
| /petsc/src/mat/impls/sbaij/seq/ |
| H A D | sbaij.c | 513 PetscInt *rp, k, low, high, t, row, nrow, i, col, l, *aj = a->j; in MatGetValues_SeqSBAIJ() local
527 rp = aj + ai[brow]; in MatGetValues_SeqSBAIJ()
544 if (rp[t] > bcol) high = t; in MatGetValues_SeqSBAIJ()
548 if (rp[i] > bcol) break; in MatGetValues_SeqSBAIJ()
549 if (rp[i] == bcol) { 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
592 rp = aj + ai[row]; in MatSetValuesBlocked_SeqSBAIJ()
615 if (rp[t] > col) high = t; in MatSetValuesBlocked_SeqSBAIJ()
619 if (rp[i] > col) break; in MatSetValuesBlocked_SeqSBAIJ()
620 if (rp[i] == col) { in MatSetValuesBlocked_SeqSBAIJ()
[all …]
|
| H A D | sbaijfact2.c | 1678 const PetscInt mbs = a->mbs, *ai = a->i, *aj = a->j, *rp, *vj, *adiag = a->diag; in MatSolve_SeqSBAIJ_1() local
1688 PetscCall(ISGetIndices(isrow, &rp)); in MatSolve_SeqSBAIJ_1()
1691 for (k = 0; k < mbs; k++) t[k] = b[rp[k]]; in MatSolve_SeqSBAIJ_1()
1707 x[rp[k]] = t[k]; in MatSolve_SeqSBAIJ_1()
1710 PetscCall(ISRestoreIndices(isrow, &rp)); in MatSolve_SeqSBAIJ_1()
1721 const PetscInt mbs = a->mbs, *ai = a->i, *aj = a->j, *rp, *vj; in MatSolve_SeqSBAIJ_1_inplace() local
1731 PetscCall(ISGetIndices(isrow, &rp)); in MatSolve_SeqSBAIJ_1_inplace()
1734 for (k = 0; k < mbs; k++) t[k] = b[rp[k]]; in MatSolve_SeqSBAIJ_1_inplace()
1750 x[rp[k]] = t[k]; in MatSolve_SeqSBAIJ_1_inplace()
1753 PetscCall(ISRestoreIndices(isrow, &rp)); in MatSolve_SeqSBAIJ_1_inplace()
[all …]
|
| /petsc/src/mat/impls/baij/seq/ |
| H A D | baij.c | 1292 PetscInt *rp, k, low, high, t, ii, jj, row, nrow, i, col, l, N, m = *mm, n = *nn; in matsetvaluesblocked4_() local
1303 rp = aj + ai[row]; in matsetvaluesblocked4_()
1316 if (rp[t] > col) high = t; in matsetvaluesblocked4_()
1320 if (rp[i] > col) break; in matsetvaluesblocked4_()
1321 if (rp[i] == col) { in matsetvaluesblocked4_()
1333 rp[ii + 1] = rp[ii]; in matsetvaluesblocked4_()
1337 rp[i] = col; in matsetvaluesblocked4_()
1360 PetscInt *rp, k, low, high, t, row, nrow, i, col, l, N, n = *nn, m = *mm; in matsetvalues4_() local
1370 rp = aj + ai[brow]; in matsetvalues4_()
1386 if (rp[t] > bcol) high = t; in matsetvalues4_()
[all …]
|
| /petsc/src/mat/impls/blockmat/seq/ |
| H A D | blockmat.c | 217 PetscInt *rp, k, low, high, t, row, nrow, i, col, l, lastcol = -1; in MatSetValues_BlockMat() local
231 rp = aj + ai[brow]; in MatSetValues_BlockMat()
252 if (rp[t] > bcol) high = t; in MatSetValues_BlockMat()
256 if (rp[i] > bcol) break; in MatSetValues_BlockMat()
257 if (rp[i] == bcol) goto noinsert1; in MatSetValues_BlockMat()
262 …MatSeqXAIJReallocateAIJ(A, a->mbs, 1, nrow, brow, bcol, rmax, aa, ai, aj, rp, ap, imax, nonew, Mat… in MatSetValues_BlockMat()
267 rp[ii + 1] = rp[ii]; in MatSetValues_BlockMat()
271 rp[i] = bcol; in MatSetValues_BlockMat()
|
| /petsc/src/mat/impls/baij/mpi/ |
| H A D | mpibaij.c | 163 rp = PetscSafePointerPlusOffset(aj, ai[brow]); \
174 if (rp[t] > bcol) high = t; \
178 if (rp[_i] > bcol) break; \
179 if (rp[_i] == bcol) { \
188 …MatSeqXAIJReallocateAIJ(A, a->mbs, bs2, nrow, brow, bcol, rmax, aa, ai, aj, rp, ap, aimax, a->none…
191 PetscCall(PetscArraymove(rp + _i + 1, rp + _i, N - _i + 1)); \
194 rp[_i] = bcol; \
203 rp = PetscSafePointerPlusOffset(bj, bi[brow]); \
214 if (rp[t] > bcol) high = t; \
218 if (rp[_i] > bcol) break; \
[all …]
|
| /petsc/src/mat/impls/aij/seq/ |
| H A D | aij.c | 365 const PetscInt *rp, *ai = a->i, *ailen = a->ilen, *aj = a->j; in MatSeqAIJSetValuesLocalFast() local
373 rp = aj + ai[row]; in MatSeqAIJSetValuesLocalFast()
387 if (rp[t] > col) high = t; in MatSeqAIJSetValuesLocalFast()
391 if (rp[i] == col) { in MatSeqAIJSetValuesLocalFast()
405 PetscInt *rp, k, low, high, t, ii, row, nrow, i, col, l, rmax, N; in MatSetValues_SeqAIJ() local
418 rp = PetscSafePointerPlusOffset(aj, ai[row]); in MatSetValues_SeqAIJ()
436 if (rp[t] > col) high = t; in MatSetValues_SeqAIJ()
440 if (rp[i] > col) break; in MatSetValues_SeqAIJ()
441 if (rp[i] == col) { in MatSetValues_SeqAIJ()
456 …MatSeqXAIJReallocateAIJ_structure_only(A, A->rmap->n, 1, nrow, row, col, rmax, ai, aj, rp, imax, n… in MatSetValues_SeqAIJ()
[all …]
|
| /petsc/src/dm/impls/plex/transform/interface/ |
| H A D | plextransform.c | 1004 PetscInt rt, Nct, n, off, rp; in DMPlexTransformGetTargetPoint() local
1015 PetscCall(DMLabelGetStratumPointIndex(trType, rt, p, &rp)); in DMPlexTransformGetTargetPoint()
1016 …PetscCheck(rp >= 0, PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Cell type %s point %" PetscInt_FMT " do… in DMPlexTransformGetTargetPoint()
1019 rp = p - ctS; in DMPlexTransformGetTargetPoint()
1028 newp += rp * rsize[n] + r; in DMPlexTransformGetTargetPoint()
1062 PetscInt rt = -1, rtTmp, Nct, n, rp = 0, rO = 0, pO; in DMPlexTransformGetSourcePoint() local
1140 rp = c - ctS; in DMPlexTransformGetSourcePoint()
1144 rp = tmp / rsize[n]; in DMPlexTransformGetSourcePoint()
1151 pO = rp + ctS; in DMPlexTransformGetSourcePoint()
2101 PetscSFNode *rp, *rtmp; in DMPlexTransformCreateSF() local
[all …]
|
| /petsc/src/ksp/ksp/utils/dm/ |
| H A D | dmproject.c | 1117 const PetscInt rp = rpidx[rv]; in DMSwarmRemap_Colella_Internal() local
1118 …PetscCheck(rp >= 0 && rp < Np, PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Particle index %" PetscInt_F… in DMSwarmRemap_Colella_Internal()
1119 rw[rp] += wp * Wx[0] * Wx[1] * Wv[0] * Wv[1]; in DMSwarmRemap_Colella_Internal()
1120 …tscInt_FMT "\n", (double)(wp * Wx[0] * Wx[1] * Wv[0] * Wv[1]), (double)PetscRealPart(rw[rp]), rp)); in DMSwarmRemap_Colella_Internal()
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexinterpolate.c | 1115 PetscSFNode cpmin = {-1, -1}, rp = {-1, -1}; in DMPlexAddSharedFace_Private() local
1128 PetscCall(DMPlexMapToGlobalPoint(dm, p, &rp, NULL)); in DMPlexAddSharedFace_Private()
1131 …FMT ") Min Cone Point (%" PetscInt_FMT ", %" PetscInt_FMT ")\n", rank, rp.rank, rp.index, cpmin.ra… in DMPlexAddSharedFace_Private()
1133 if (isShared && (rp.rank == cpmin.rank && rp.index == cpmin.index)) { in DMPlexAddSharedFace_Private()
|
| H A D | plexgeometry.c | 4066 const PetscReal rp = (x[2] + 1) * 0.5 * (ro - ri) + ri; in coordMap_shell() local
4070 xp[0] = rp * PetscCosReal(thetap) * PetscCosReal(phip); in coordMap_shell()
4071 xp[1] = rp * PetscCosReal(thetap) * PetscSinReal(phip); in coordMap_shell()
4072 xp[2] = rp * PetscSinReal(thetap); in coordMap_shell()
|
| /petsc/src/binding/petsc4py/src/petsc4py/PETSc/ |
| H A D | Mat.pyx | 4700 rp : IS
4712 cdef IS rp = IS(), cp = IS()
4713 CHKERR(MatGetOrdering(self.mat, cval, &rp.iset, &cp.iset))
4714 return (rp, cp)
4741 cdef PetscIS rp = isrow.iset, cp = iscol.iset
4742 CHKERR(MatReorderForNonzeroDiagonal(self.mat, rval, rp, cp))
|
| /petsc/doc/ |
| H A D | petsc.bib | 5859 @InProceedings{ rp, inproceedings
|