Home
last modified time | relevance | path

Searched refs:rk (Results 1 – 12 of 12) sorted by relevance

/petsc/src/ts/impls/explicit/rk/
H A Dmrk.c21 TS_RK *rk = (TS_RK *)ts->data; in TSReset_RK_MultirateNonsplit() local
22 RKTableau tab = rk->tableau; in TSReset_RK_MultirateNonsplit()
25 PetscCall(VecDestroy(&rk->X0)); in TSReset_RK_MultirateNonsplit()
26 PetscCall(VecDestroyVecs(tab->s, &rk->YdotRHS_slow)); in TSReset_RK_MultirateNonsplit()
32 TS_RK *rk = (TS_RK *)ts->data; in TSInterpolate_RK_MultirateNonsplit() local
33 PetscInt s = rk->tableau->s, p = rk->tableau->p, i, j; in TSInterpolate_RK_MultirateNonsplit()
37 const PetscReal *B = rk->tableau->binterp; in TSInterpolate_RK_MultirateNonsplit()
40 …tscObject)ts), PETSC_ERR_SUP, "TSRK %s does not have an interpolation formula", rk->tableau->name); in TSInterpolate_RK_MultirateNonsplit()
41 t = (itime - rk->ptime) / h; in TSInterpolate_RK_MultirateNonsplit()
47 PetscCall(VecCopy(rk->X0, X)); in TSInterpolate_RK_MultirateNonsplit()
[all …]
H A Drk.c487 TS_RK *rk = (TS_RK *)ts->data; in TSRKGetTableau_RK() local
488 RKTableau tab = rk->tableau; in TSRKGetTableau_RK()
550 TS_RK *rk = (TS_RK *)ts->data; in TSEvaluateStep_RK() local
551 RKTableau tab = rk->tableau; in TSEvaluateStep_RK()
552 PetscScalar *w = rk->work; in TSEvaluateStep_RK()
557 switch (rk->status) { in TSEvaluateStep_RK()
569 if (rk->status == TS_STEP_INCOMPLETE) { in TSEvaluateStep_RK()
571 for (j = 0; j < s; j++) w[j] = h * tab->b[j] / rk->dtratio; in TSEvaluateStep_RK()
572 PetscCall(VecMAXPY(X, s, w, rk->YdotRHS)); in TSEvaluateStep_RK()
577 if (rk->status == TS_STEP_INCOMPLETE) { /*Complete with the embedded method (be)*/ in TSEvaluateStep_RK()
[all …]
/petsc/src/ts/tests/output/
H A Dex2_rk.out3 …TSAdapt basic rk 0:3bs step 0 accepted t=0 + 1.000e-03 dt=1.000e-02 wlte=4.09e-06 wlt…
6 …TSAdapt basic rk 0:3bs step 1 accepted t=0.001 + 1.000e-02 dt=5.597e-02 wlte=0.00416 wlte…
9 …TSAdapt basic rk 0:3bs step 2 accepted t=0.011 + 5.597e-02 dt=5.465e-02 wlte=0.783 wltea=…
12 …TSAdapt basic rk 0:3bs step 3 accepted t=0.0669669 + 5.465e-02 dt=5.325e-02 wlte=0.788 wltea=…
15 …TSAdapt basic rk 0:3bs step 4 accepted t=0.121616 + 5.325e-02 dt=5.205e-02 wlte=0.781 wltea=…
18 …TSAdapt basic rk 0:3bs step 5 accepted t=0.174867 + 5.205e-02 dt=5.102e-02 wlte=0.774 wltea=…
21 …TSAdapt basic rk 0:3bs step 6 accepted t=0.226919 + 5.102e-02 dt=5.014e-02 wlte=0.768 wltea=…
24 …TSAdapt basic rk 0:3bs step 7 accepted t=0.277942 + 5.014e-02 dt=4.938e-02 wlte=0.763 wltea=…
27 …TSAdapt basic rk 0:3bs step 8 accepted t=0.328081 + 4.938e-02 dt=4.872e-02 wlte=0.759 wltea=…
30 …TSAdapt basic rk 0:3bs step 9 accepted t=0.377459 + 4.872e-02 dt=4.816e-02 wlte=0.755 wltea=…
[all …]
H A Dex4_8.out13 type: rk
/petsc/src/snes/tutorials/
H A Dex3k.kokkos.cxx124 Vec rk; in StubFunction() local
128 PetscCall(DMGetGlobalVector(da, &rk)); in StubFunction()
130 PetscCall(KokkosFunction(snes, x, rk, ctx)); in StubFunction()
131 PetscCall(VecAXPY(rk, -1.0, r)); in StubFunction()
132 PetscCall(VecNorm(rk, NORM_2, &norm)); in StubFunction()
133 PetscCall(DMRestoreGlobalVector(da, &rk)); in StubFunction()
/petsc/src/sys/utils/
H A Dsortd.c167 PetscReal rk, rtmp; in PetscSortRealWithArrayInt() local
174 rk = r[k]; in PetscSortRealWithArrayInt()
176 if (rk > r[j]) { in PetscSortRealWithArrayInt()
178 rk = r[k]; in PetscSortRealWithArrayInt()
/petsc/doc/overview/
H A Dintegrator_table.md25 * - rk*
/petsc/src/snes/interface/
H A Dsnes.c5600 PetscReal rk = ared / pred; in KSPPreSolve_SNESEW() local
5601 if (rk < kctx->v4_p1) rtol = 1. - 2. * kctx->v4_p1; in KSPPreSolve_SNESEW()
5602 else if (rk < kctx->v4_p2) rtol = kctx->rtol_last; in KSPPreSolve_SNESEW()
5603 else if (rk < kctx->v4_p3) rtol = kctx->v4_m1 * kctx->rtol_last; in KSPPreSolve_SNESEW()
5609 kctx->rk_last = rk; in KSPPreSolve_SNESEW()
/petsc/src/mat/tests/
H A Dex23.c305 …PetscInt ri, ci, rr[3] = {0, 1, 0}, cr[4] = {1, 2, 0, 1}, rk[3] = {0, 2, 1}, ck[4] =… in main() local
311 r = (PetscInt *)(ri == 0 ? rr : rk); in main()
/petsc/systems/Apple/iOS/PETSc/PETSc.xcodeproj/
H A Dproject.pbxproj608 …D89D1ED9207AA1AC00FACC19 /* rk.c in Sources */ = {isa = PBXBuildFile; fileRef = D89D1BF1207AA1AA00…
1360 …89D1BF1207AA1AA00FACC19 /* rk.c */ = {isa = PBXFileReference; fileEncoding = 4; lastKnownFileType …
2131 D89D1BF1207AA1AA00FACC19 /* rk.c */,
2956 D89D1ED9207AA1AC00FACC19 /* rk.c in Sources */,
/petsc/doc/manual/
H A Dts.md217 They can also be set with the options database option `-ts_type euler, rk, beuler, cn, theta, gl, p…
/petsc/share/petsc/datafiles/meshes/
H A Dtestcase3D.cas521 (spark/rk (0.002))
2016 (dpm/multicomponent-min-rk-h-evap 1e-12)
2017 (dpm/multicomponent-min-rk-h-react 9.999999999999999e-21)
2019 (dpm/multicomponent-remove-particle-at-rk-failure? #f)
2148 (dpm/acceleration-to-faces-for-rk? #t)
2457 (dpm/use-prev-rk-timestep? #t)