| /petsc/src/ts/impls/eimex/ ! |
| H A D | eimex.c | 16 Vec YdotI; /* xdot-f(x) Work vector holding YdotI = F(t,x,xdot) when xdot =0 */ member
153 Vec Y0 = ext->VecSolPrev, Y1 = ext->Y, Ydot = ext->Ydot, YdotI = ext->YdotI; in TSInterpolate_EIMEX() local
161 PetscCall(TSComputeIFunction(ts, ts->ptime - h, Y0, Ydot, YdotI, PETSC_FALSE)); in TSInterpolate_EIMEX()
165 PetscCall(VecAXPBYPCZ(X, a, b, 0.0, Y0, YdotI)); in TSInterpolate_EIMEX()
167 PetscCall(TSComputeIFunction(ts, ts->ptime, Y1, Ydot, YdotI, PETSC_FALSE)); in TSInterpolate_EIMEX()
170 PetscCall(VecAXPBYPCZ(X, a, b, 1.0, Y1, YdotI)); in TSInterpolate_EIMEX()
185 PetscCall(VecDestroy(&ext->YdotI)); in TSReset_EIMEX()
203 static PetscErrorCode TSEIMEXGetVecs(TS ts, DM dm, Vec *Z, Vec *Ydot, Vec *YdotI, Vec *YdotRHS) in TSEIMEXGetVecs() argument
216 if (YdotI) { in TSEIMEXGetVecs()
217 if (dm && dm != ts->dm) PetscCall(DMGetNamedGlobalVector(dm, "TSEIMEX_YdotI", YdotI)); in TSEIMEXGetVecs()
[all …]
|
| /petsc/src/ts/impls/implicit/irk/ ! |
| H A D | irk.c | 39 Vec *YdotI; /* Work vectors to hold the residual evaluation */ member
273 Vec *YdotI = irk->YdotI; in TSEvaluateStep_IRK() local
293 PetscCall(VecMAXPY(U, irk->nstages, w, YdotI)); in TSEvaluateStep_IRK()
335 PetscCall(VecZeroEntries(irk->YdotI[i])); in TSStep_IRK()
336 …for (j = 0; j < nstages; j++) PetscCall(VecAXPY(irk->YdotI[i], A_inv[i + j * nstages] / ts->time_s… in TSStep_IRK()
337 PetscCall(VecAXPY(irk->YdotI[i], -A_inv_rowsum[i] / ts->time_step, irk->U)); in TSStep_IRK()
394 PetscCall(VecMAXPY(U, nstages, bt, irk->YdotI)); in TSInterpolate_IRK()
420 PetscCall(VecDestroyVecs(irk->nstages, &irk->YdotI)); in TSReset_IRK()
464 Vec U, *YdotI = irk->YdotI, Ydot = irk->Ydot, *Y = irk->Y; in SNESTSFormFunction_IRK() local
478 …PetscCall(TSComputeIFunction(ts, ts->ptime + ts->time_step * c[i], Y[i], Ydot, YdotI[i], PETSC_FAL… in SNESTSFormFunction_IRK()
[all …]
|
| /petsc/src/ts/impls/arkimex/ ! |
| H A D | arkimex.h | 26 Vec *YdotI; /* Time derivatives for the stiff part */ member
|
| H A D | arkimex.c | 1253 PetscCall(VecMAXPY(X, s, w, ark->YdotI)); in TSEvaluateStep_ARKIMEX()
1270 PetscCall(VecMAXPY(X, s, w, ark->YdotI)); in TSEvaluateStep_ARKIMEX()
1281 PetscCall(VecMAXPY(X, tab->s, w, ark->YdotI)); in TSEvaluateStep_ARKIMEX()
1337 …Vec *Y = ark->Y, *YdotI = ark->YdotI, *YdotRHS = ark->YdotRHS, Ydot = ark->Ydot, Ydot0… in TSStep_ARKIMEX() local
1358 PetscCall(VecCopy(YdotI[s - 1], Ydot0)); in TSStep_ARKIMEX()
1363 PetscCall(VecCopy(YdotI[i], ark->YdotI_prev[i])); in TSStep_ARKIMEX()
1419 PetscCall(VecCopy(ark_start->YdotI[ark_start->tableau->s - 1], Ydot0)); in TSStep_ARKIMEX()
1443 PetscCall(VecMAXPY(Y[i], i, w, YdotI)); in TSStep_ARKIMEX()
1455 PetscCall(VecMAXPY(Z, i, w, YdotI)); in TSStep_ARKIMEX()
1486 PetscCall(VecCopy(Ydot0, YdotI[0])); /* YdotI = YdotI(tn-1) */ in TSStep_ARKIMEX()
[all …]
|
| H A D | fsarkimex.c | 276 PetscCall(VecCopy(ark_start->YdotI[ark_start->tableau->s - 1], Ydot0_fast)); in TSStep_ARKIMEX_FastSlowSplit()
|