Lines Matching refs:YdotI
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()
218 else *YdotI = ext->YdotI; in TSEIMEXGetVecs()
227 static PetscErrorCode TSEIMEXRestoreVecs(TS ts, DM dm, Vec *Z, Vec *Ydot, Vec *YdotI, Vec *YdotRHS) in TSEIMEXRestoreVecs() argument
236 if (YdotI) { in TSEIMEXRestoreVecs()
237 if (dm && dm != ts->dm) PetscCall(DMRestoreNamedGlobalVector(dm, "TSEIMEX_YdotI", YdotI)); in TSEIMEXRestoreVecs()
341 PetscCall(VecDuplicate(ts->vec_sol, &ext->YdotI)); in TSSetUp_EIMEX()