| /petsc/src/ts/impls/explicit/rk/ |
| H A D | mrk.c | 58 Vec *Y = rk->Y, *YdotRHS = rk->YdotRHS; in TSStepRefine_RK_MultirateNonsplit() local
77 PetscCall(VecMAXPY(vec_fast, i, w, YdotRHS)); in TSStepRefine_RK_MultirateNonsplit()
83 …Call(TSComputeRHSFunction(ts, t + k * h / rk->dtratio + h / rk->dtratio * c[i], Y[i], YdotRHS[i])); in TSStepRefine_RK_MultirateNonsplit()
101 PetscCall(VecCopy(YdotRHS[i], rk->YdotRHS_slow[i])); in TSStepRefine_RK_MultirateNonsplit()
122 Vec *Y = rk->Y, *YdotRHS = rk->YdotRHS, *YdotRHS_slow = rk->YdotRHS_slow; in TSStep_RK_MultirateNonsplit() local
148 rk->YdotRHS = YdotRHS_slow; in TSStep_RK_MultirateNonsplit()
152 rk->YdotRHS = YdotRHS; in TSStep_RK_MultirateNonsplit()
263 …for (i = 0; i < s; i++) PetscCall(VecGetSubVector(rk->YdotRHS[i], rk->is_slow, &rk->YdotRHS_slow[i… in TSInterpolate_RK_MultirateSplit()
268 …for (i = 0; i < s; i++) PetscCall(VecRestoreSubVector(rk->YdotRHS[i], rk->is_slow, &rk->YdotRHS_sl… in TSInterpolate_RK_MultirateSplit()
312 Vec *YdotRHS = rk->YdotRHS, *YdotRHS_fast = rk->YdotRHS_fast; in TSStepRefine_RK_MultirateSplit() local
[all …]
|
| H A D | rk.h | 25 …Vec *YdotRHS; /* Function evaluations for the non-stiff part and contains all compone… member
|
| H A D | rk.c | 572 PetscCall(VecMAXPY(X, s, w, rk->YdotRHS)); in TSEvaluateStep_RK()
580 PetscCall(VecMAXPY(X, s, w, rk->YdotRHS)); in TSEvaluateStep_RK()
584 PetscCall(VecMAXPY(X, s, w, rk->YdotRHS)); in TSEvaluateStep_RK()
643 Vec *Y = rk->Y, *YdotRHS = rk->YdotRHS; in TSRollBack_RK() local
660 PetscCall(VecMAXPY(ts->vec_sol, s, w, YdotRHS)); in TSRollBack_RK()
782 Vec *Y = rk->Y, *YdotRHS = rk->YdotRHS; in TSStep_RK() local
792 if (FSAL) PetscCall(VecCopy(YdotRHS[s - 1], YdotRHS[0])); in TSStep_RK()
804 PetscCall(VecMAXPY(Y[i], i, w, YdotRHS)); in TSStep_RK()
810 PetscCall(TSComputeRHSFunction(ts, t + h * c[i], Y[i], YdotRHS[i])); in TSStep_RK()
1077 PetscCall(VecMAXPY(X, s, b, rk->YdotRHS)); in TSInterpolate_RK()
[all …]
|
| /petsc/src/ts/impls/eimex/ |
| H A D | eimex.c | 15 Vec YdotRHS; /* g(x) Work vector holding YdotRHS during residual evaluation */ member
184 PetscCall(VecDestroy(&ext->YdotRHS)); in TSReset_EIMEX()
203 static PetscErrorCode TSEIMEXGetVecs(TS ts, DM dm, Vec *Z, Vec *Ydot, Vec *YdotI, Vec *YdotRHS) in TSEIMEXGetVecs() argument
220 if (YdotRHS) { in TSEIMEXGetVecs()
221 if (dm && dm != ts->dm) PetscCall(DMGetNamedGlobalVector(dm, "TSEIMEX_YdotRHS", YdotRHS)); in TSEIMEXGetVecs()
222 else *YdotRHS = ext->YdotRHS; in TSEIMEXGetVecs()
227 static PetscErrorCode TSEIMEXRestoreVecs(TS ts, DM dm, Vec *Z, Vec *Ydot, Vec *YdotI, Vec *YdotRHS) in TSEIMEXRestoreVecs() argument
239 if (YdotRHS) { in TSEIMEXRestoreVecs()
240 if (dm && dm != ts->dm) PetscCall(DMRestoreNamedGlobalVector(dm, "TSEIMEX_YdotRHS", YdotRHS)); in TSEIMEXRestoreVecs()
342 PetscCall(VecDuplicate(ts->vec_sol, &ext->YdotRHS)); in TSSetUp_EIMEX()
|
| /petsc/src/ts/impls/multirate/ |
| H A D | mprk.c | 50 Vec *YdotRHS; member
637 PetscCall(VecMAXPY(X, s, wf, mprk->YdotRHS)); in TSEvaluateStep_MPRK()
644 …Vec *Y = mprk->Y, *YdotRHS = mprk->YdotRHS, *YdotRHS_fast = mprk->YdotRHS_fast, *YdotR… in TSStep_MPRK() local
661 …for (j = 0; j < i; j++) PetscCall(VecGetSubVector(YdotRHS[j], mprk->is_slowbuffer, &YdotRHS_slowbu… in TSStep_MPRK()
665 …for (j = 0; j < i; j++) PetscCall(VecRestoreSubVector(YdotRHS[j], mprk->is_slowbuffer, &YdotRHS_sl… in TSStep_MPRK()
668 … for (j = 0; j < i; j++) PetscCall(VecGetSubVector(YdotRHS[j], mprk->is_slow, &YdotRHS_slow[j])); in TSStep_MPRK()
672 …for (j = 0; j < i; j++) PetscCall(VecRestoreSubVector(YdotRHS[j], mprk->is_slow, &YdotRHS_slow[j])… in TSStep_MPRK()
677 for (j = 0; j < i; j++) PetscCall(VecGetSubVector(YdotRHS[j], mprk->is_fast, &YdotRHS_fast[j])); in TSStep_MPRK()
681 …for (j = 0; j < i; j++) PetscCall(VecRestoreSubVector(YdotRHS[j], mprk->is_fast, &YdotRHS_fast[j])… in TSStep_MPRK()
690 …for (j = 0; j < i; j++) PetscCall(VecGetSubVector(YdotRHS[j], mprk->is_medium, &YdotRHS_medium[j])… in TSStep_MPRK()
[all …]
|
| /petsc/src/ts/impls/arkimex/ |
| H A D | arkimex.h | 27 Vec *YdotRHS; /* Function evaluations for the non-stiff part */ member
|
| H A D | arkimex.c | 1258 PetscCall(VecMAXPY(X, s, w, ark->YdotRHS)); in TSEvaluateStep_ARKIMEX()
1275 PetscCall(VecMAXPY(X, s, w, ark->YdotRHS)); in TSEvaluateStep_ARKIMEX()
1286 PetscCall(VecMAXPY(X, s, w, ark->YdotRHS)); in TSEvaluateStep_ARKIMEX()
1337 …Vec *Y = ark->Y, *YdotI = ark->YdotI, *YdotRHS = ark->YdotRHS, Ydot = ark->Ydot, Ydot0… in TSStep_ARKIMEX() local
1364 if (tab->additive && hasE) PetscCall(VecCopy(YdotRHS[i], ark->YdotRHS_prev[i])); in TSStep_ARKIMEX()
1446 PetscCall(VecMAXPY(Y[i], i, w, YdotRHS)); in TSStep_ARKIMEX()
1458 PetscCall(VecMAXPY(Z, i, w, YdotRHS)); in TSStep_ARKIMEX()
1500 PetscCall(TSComputeRHSFunction(ts, t + h * c[i], Y[i], YdotRHS[i])); in TSStep_ARKIMEX()
1502 PetscCall(VecZeroEntries(YdotRHS[i])); in TSStep_ARKIMEX()
1706 if (hasE) PetscCall(VecMAXPY(X, s, b, ark->YdotRHS)); in TSInterpolate_ARKIMEX()
[all …]
|