| /petsc/src/snes/tutorials/ |
| H A D | ex9.c | 36 …const PetscReal r = x * x + y * y, r0 = 0.9, psi0 = PetscSqrtReal(1.0 - r0 * r0), dpsi0 = -r0 / ps… in psi() local
37 if (r <= r0) { in psi()
40 return psi0 + dpsi0 * (r - r0); in psi()
|
| /petsc/src/dm/impls/swarm/tests/ |
| H A D | ex5.c | 343 const PetscReal r0 = (p / n) + 1.; in circleMultipleX() local
346 x[0] = r0 * PetscCosReal(th0); in circleMultipleX()
347 x[1] = r0 * PetscSinReal(th0); in circleMultipleX()
355 const PetscReal r0 = (p / n) + 1.; in circleMultipleV() local
357 const PetscReal omega = PetscSqrtReal(1000. / r0) / r0; in circleMultipleV()
359 v[0] = -r0 * omega * PetscSinReal(th0); in circleMultipleV()
360 v[1] = r0 * omega * PetscCosReal(th0); in circleMultipleV()
447 const PetscReal r0 = DMPlex_NormD_Internal(dim, &coords[p * dim]); in ComputeError() local
450 const PetscReal omega = v0 / r0; in ComputeError()
455 const PetscReal xe[3] = {r0 * ct, r0 * st, 0.0}; in ComputeError()
|
| H A D | ex6.c | 879 const PetscReal r0 = (p / n) + 1.; in circleMultipleX() local
882 x[0] = r0 * PetscCosReal(th0); in circleMultipleX()
883 x[1] = r0 * PetscSinReal(th0); in circleMultipleX()
891 const PetscReal r0 = (p / n) + 1.; in circleMultipleV() local
893 const PetscReal omega = PetscSqrtReal(1000. / r0) / r0; in circleMultipleV()
895 v[0] = -r0 * omega * PetscSinReal(th0); in circleMultipleV()
896 v[1] = r0 * omega * PetscCosReal(th0); in circleMultipleV()
985 const PetscReal r0 = DMPlex_NormD_Internal(dim, &coords[p * dim]); in ComputeError() local
988 const PetscReal omega = v0 / r0; in ComputeError()
993 const PetscReal xe[3] = {r0 * ct, r0 * st, 0.0}; in ComputeError()
|
| /petsc/src/ksp/pc/impls/bjacobi/bjkokkos/ |
| H A D | bjkokkos.kokkos.cxx | 93 PetscReal dp, dpold, w, dpest, tau, psi, cm, r0; in BJSolve_TFQMR() local
166 r0 = dp = PetscSqrtReal(PetscRealPart(dpi)); in BJSolve_TFQMR()
256 if (dpest / r0 < rtol) { in BJSolve_TFQMR()
261 if (dpest / r0 > dtol) { in BJSolve_TFQMR()
263 … printf("ERROR block %d diverged: %d it, res=%e, r_0=%e\n", team.league_rank(), it, dpest, r0); }); in BJSolve_TFQMR()
267 if (dpest / r0 > dtol) { in BJSolve_TFQMR()
275 … TFQMR %d:%d it, res=%e, r_0=%e r_res=%e\n", team.league_rank(), it, m, dpest, r0, dpest / r0); }); in BJSolve_TFQMR()
335 PetscReal dp, r0; in BJSolve_BICG() local
398 r0 = dp = PetscSqrtReal(PetscRealPart(dpi)); in BJSolve_BICG()
477 if (dp / r0 < rtol) { in BJSolve_BICG()
[all …]
|
| /petsc/src/ksp/ksp/impls/ibcgs/ |
| H A D | ibcgs.c | 57 …const PetscScalar *PETSC_RESTRICT r0, *PETSC_RESTRICT f0, *PETSC_RESTRICT qn, *PETSC_RESTRICT b, *… in KSPSolve_IBCGS() local
85 PetscCall(VecGetArrayRead(R0, (const PetscScalar **)&r0)); in KSPSolve_IBCGS()
207 phin += r0[i] * sn[i]; in KSPSolve_IBCGS()
208 pin += r0[i] * qn[i]; in KSPSolve_IBCGS()
|
| /petsc/src/dm/impls/swarm/ |
| H A D | data_ex.c | 259 static void _get_tags(PetscInt counter, PetscMPIInt N, PetscMPIInt r0, PetscMPIInt r1, PetscMPIInt … in _get_tags() argument
263 st = (N * r0 + r1 + N * N * counter) % maxtag; in _get_tags()
264 rt = (N * r1 + r0 + N * N * counter) % maxtag; in _get_tags()
333 PetscMPIInt symm_nn = 0, *symm_procs = NULL, r0, n, st, rt, size, *maxtag, iflg; in DMSwarmDataExTopologyFinalize() local
358 r0 = d->rank; in DMSwarmDataExTopologyFinalize()
362 _get_tags(d->instance, size, r0, r1, *maxtag, &st, &rt); in DMSwarmDataExTopologyFinalize()
|
| /petsc/src/ts/tutorials/ |
| H A D | ex11.h | 588 static PetscScalar s, c0, p0, r0, u0, w0, x0, x2, ri, cx, sgn0, wsp0, gasc1, gasc2, gasc3, gasc4; in riemannsolver() local
617 r0 = cvmgm_(rl, rr, &sgn0); in riemannsolver()
625 c0 = PetscSqrtScalar(*gam * p0 / r0); in riemannsolver()
627 w0 = *gam * r0 * p0 * (gasc3 * streng / p0 + 1.); in riemannsolver()
628 rstars = r0 / (1. - r0 * streng / w0); in riemannsolver()
631 rstarr = r0 * PetscPowScalar(d__1, d__2); in riemannsolver()
647 s = p0 / PetscPowScalar(r0, *gam); in riemannsolver()
652 *rx = cvmgp_(rx, &r0, &d__1); in riemannsolver()
|