Home
last modified time | relevance | path

Searched refs:stage_time (Results 1 – 17 of 17) sorted by relevance

/petsc/src/ts/impls/explicit/ssp/
H A Dssp.c65 PetscReal stage_time = t0 + dt * (i / (s - 1.)); in TSSSPStep_RK_2() local
66 PetscCall(TSPreStage(ts, stage_time)); in TSSSPStep_RK_2()
67 PetscCall(TSComputeRHSFunction(ts, stage_time, work[0], F)); in TSSSPStep_RK_2()
90 PetscReal c, stage_time; in TSSSPStep_RK_3() local
102 stage_time = t0 + c * dt; in TSSSPStep_RK_3()
103 PetscCall(TSPreStage(ts, stage_time)); in TSSSPStep_RK_3()
104 PetscCall(TSComputeRHSFunction(ts, stage_time, work[0], F)); in TSSSPStep_RK_3()
110 stage_time = t0 + c * dt; in TSSSPStep_RK_3()
111 PetscCall(TSPreStage(ts, stage_time)); in TSSSPStep_RK_3()
112 PetscCall(TSComputeRHSFunction(ts, stage_time, work[0], F)); in TSSSPStep_RK_3()
[all …]
/petsc/src/ts/impls/implicit/alpha/
H A Dalpha1.c20 PetscReal stage_time; member
64 th->stage_time = t + Alpha_f * dt; in TSAlpha_StageTime()
130 th->stage_time = ts->ptime; in TSAlpha_Restart()
134 th->stage_time += ts->time_step; in TSAlpha_Restart()
136 PetscCall(TSPreStage(ts, th->stage_time)); in TSAlpha_Restart()
139 PetscCall(TSPostStage(ts, th->stage_time, 0, &X1)); in TSAlpha_Restart()
140 PetscCall(TSAdaptCheckStage(ts->adapt, ts, th->stage_time, X1, &stageok)); in TSAlpha_Restart()
144 th->stage_time += ts->time_step; in TSAlpha_Restart()
146 PetscCall(TSPreStage(ts, th->stage_time)); in TSAlpha_Restart()
149 PetscCall(TSPostStage(ts, th->stage_time, 0, &X2)); in TSAlpha_Restart()
[all …]
H A Dalpha2.c20 PetscReal stage_time; member
97 th->stage_time = t + Alpha_f * dt; in TSAlpha_StageTime()
174 th->stage_time = ts->ptime; in TSAlpha_Restart()
178 th->stage_time += ts->time_step; in TSAlpha_Restart()
181 PetscCall(TSPreStage(ts, th->stage_time)); in TSAlpha_Restart()
185 PetscCall(TSPostStage(ts, th->stage_time, 0, &X1)); in TSAlpha_Restart()
186 PetscCall(TSAdaptCheckStage(ts->adapt, ts, th->stage_time, X1, &stageok)); in TSAlpha_Restart()
190 th->stage_time += ts->time_step; in TSAlpha_Restart()
193 PetscCall(TSPreStage(ts, th->stage_time)); in TSAlpha_Restart()
197 PetscCall(TSPostStage(ts, th->stage_time, 0, &X2)); in TSAlpha_Restart()
[all …]
/petsc/src/ts/impls/implicit/theta/
H A Dtheta.c11 PetscReal stage_time; member
137 PetscCall(TSComputeRHSFunction(quadts, th->stage_time, th->X, ts->vec_costintegrand)); in TSThetaEvaluateCostIntegral()
213 th->stage_time = ts->ptime + (th->endpoint ? (PetscReal)1 : th->Theta) * ts->time_step; in TSStep_Theta()
222 PetscCall(TSPreStage(ts, th->stage_time)); in TSStep_Theta()
224 PetscCall(TSPostStage(ts, th->stage_time, 0, &th->X)); in TSStep_Theta()
225 PetscCall(TSAdaptCheckStage(ts->adapt, ts, th->stage_time, th->X, &stageok)); in TSStep_Theta()
300 th->stage_time = ts->ptime; in TSAdjointStepBEuler_Private()
304 if (quadts) PetscCall(TSComputeRHSJacobian(quadts, th->stage_time, ts->vec_sol, quadJ, NULL)); in TSAdjointStepBEuler_Private()
339 …PetscCall(TSComputeIHessianProductFunctionUU(ts, th->stage_time, ts->vec_sol, VecsDeltaLam, ts->ve… in TSAdjointStepBEuler_Private()
341 …PetscCall(TSComputeIHessianProductFunctionUP(ts, th->stage_time, ts->vec_sol, VecsDeltaLam, ts->ve… in TSAdjointStepBEuler_Private()
[all …]
/petsc/src/ts/impls/implicit/discgrad/
H A Dtsdiscgrad.c20 PetscReal stage_time; member
235 dg->stage_time = ts->ptime + 0.5 * ts->time_step; in TSStep_DiscGrad()
238 PetscCall(TSPreStage(ts, dg->stage_time)); in TSStep_DiscGrad()
240 PetscCall(TSPostStage(ts, dg->stage_time, 0, &dg->X)); in TSStep_DiscGrad()
241 PetscCall(TSAdaptCheckStage(adapt, ts, dg->stage_time, dg->X, &stageok)); in TSStep_DiscGrad()
318 PetscCall((*dg->Sfunc)(ts, dg->stage_time, x, S, dg->funcCtx)); in SNESTSFormFunction_DiscGrad()
353 PetscCall((*dg->Gfunc)(ts, dg->stage_time, Xquad, den, dg->funcCtx)); in SNESTSFormFunction_DiscGrad()
362 PetscCall((*dg->Ffunc)(ts, dg->stage_time, Xp, &F, dg->funcCtx)); in SNESTSFormFunction_DiscGrad()
363 PetscCall((*dg->Ffunc)(ts, dg->stage_time, X0, &F0, dg->funcCtx)); in SNESTSFormFunction_DiscGrad()
364 PetscCall((*dg->Gfunc)(ts, dg->stage_time, x, G, dg->funcCtx)); in SNESTSFormFunction_DiscGrad()
[all …]
/petsc/src/ts/impls/explicit/rk/
H A Drk.c680 PetscReal stage_time, h = ts->time_step; in TSForwardStep_RK() local
688 stage_time = ts->ptime + h * c[i]; in TSForwardStep_RK()
702 PetscCall(TSComputeRHSJacobian(ts, stage_time, Y[i], J, J)); in TSForwardStep_RK()
705 PetscCall(TSComputeRHSJacobianP(ts, stage_time, Y[i], ts->Jacprhs)); /* get f_p */ in TSForwardStep_RK()
800 rk->stage_time = t + h * c[i]; in TSStep_RK()
801 PetscCall(TSPreStage(ts, rk->stage_time)); in TSStep_RK()
805 PetscCall(TSPostStage(ts, rk->stage_time, i, Y)); in TSStep_RK()
807 PetscCall(TSAdaptCheckStage(adapt, ts, rk->stage_time, Y[i], &stageok)); in TSStep_RK()
897 rk->stage_time = t + h * (1.0 - c[i]); in TSAdjointStep_RK()
899 …if (quadts) PetscCall(TSComputeRHSJacobian(quadts, rk->stage_time, Y[i], Jquad, Jquad)); /* get r_… in TSAdjointStep_RK()
[all …]
H A Dmrk.c138 rk->stage_time = t + h * c[i]; in TSStep_RK_MultirateNonsplit()
139 PetscCall(TSPreStage(ts, rk->stage_time)); in TSStep_RK_MultirateNonsplit()
143 PetscCall(TSPostStage(ts, rk->stage_time, i, Y)); in TSStep_RK_MultirateNonsplit()
334 subrk_fast->stage_time = t + h / rk->dtratio * c[i]; in TSStepRefine_RK_MultirateSplit()
335 PetscCall(TSPreStage(subts_fast, subrk_fast->stage_time)); in TSStepRefine_RK_MultirateSplit()
341 PetscCall(TSPostStage(subts_fast, subrk_fast->stage_time, i, Y)); in TSStepRefine_RK_MultirateSplit()
386 rk->stage_time = t + h * c[i]; in TSStep_RK_MultirateSplit()
387 PetscCall(TSPreStage(ts, rk->stage_time)); in TSStep_RK_MultirateSplit()
396 PetscCall(TSPostStage(ts, rk->stage_time, i, Y)); in TSStep_RK_MultirateSplit()
H A Drk.h36 PetscReal stage_time; member
/petsc/src/ts/impls/mimex/
H A Dmimex.c11 PetscReal stage_time; member
76 PetscCall(TSComputeIFunction(ts, mimex->stage_time, x, Xdot, y, PETSC_TRUE)); in SNESTSFormFunction_Mimex()
137 PetscCall(TSComputeIJacobian(ts, mimex->stage_time, x, Xdot, shift, A, B, PETSC_TRUE)); in SNESTSFormJacobian_Mimex()
162 mimex->stage_time = ts->ptime + ts->time_step; in TSStep_Mimex_Split()
222 mimex->stage_time = ts->ptime + ts->time_step; in TSStep_Mimex_Implicit()
/petsc/src/ts/impls/arkimex/
H A Dfsarkimex.c48 PetscCall(TSComputeIFunction(ark->subts_fast, ark->stage_time, Y, Ydot, F, ark->imex)); in SNESTSFormFunction_ARKIMEX_FastSlowSplit()
69 PetscCall(TSComputeIJacobian(ark->subts_fast, ark->stage_time, Y, Ydot, shift, A, B, ark->imex)); in SNESTSFormJacobian_ARKIMEX_FastSlowSplit()
293 ark->stage_time = t + h * ct[i]; in TSStep_ARKIMEX_FastSlowSplit()
294 PetscCall(TSPreStage(ts, ark->stage_time)); in TSStep_ARKIMEX_FastSlowSplit()
336 PetscCall(TSAdaptCheckStage(adapt, ts, ark->stage_time, Y[i], &stageok)); in TSStep_ARKIMEX_FastSlowSplit()
358 …PetscCall(TSComputeIFunction(ark->subts_fast, ark->stage_time, Y[i], Ydot_fast, YdotI_fast[i], ark… in TSStep_ARKIMEX_FastSlowSplit()
365 …if (fasthasE) PetscCall(TSComputeRHSFunction(ark->subts_fast, ark->stage_time, Y[i], YdotRHS_fast[… in TSStep_ARKIMEX_FastSlowSplit()
374 PetscCall(TSComputeRHSFunction(ark->subts_slow, ark->stage_time, Y[i], YdotRHS_slow[i])); in TSStep_ARKIMEX_FastSlowSplit()
376 PetscCall(TSPostStage(ts, ark->stage_time, i, Y)); in TSStep_ARKIMEX_FastSlowSplit()
H A Darkimex.h37 PetscReal stage_time; member
H A Darkimex.c1437 ark->stage_time = t + h * ct[i]; in TSStep_ARKIMEX()
1438 PetscCall(TSPreStage(ts, ark->stage_time)); in TSStep_ARKIMEX()
1474 PetscCall(TSAdaptCheckStage(adapt, ts, ark->stage_time, Y[i], &stageok)); in TSStep_ARKIMEX()
1506 PetscCall(TSPostStage(ts, ark->stage_time, i, Y)); in TSStep_ARKIMEX()
1581 ark->stage_time = t - adjoint_time_step * (1.0 - ct[i]); in TSAdjointStep_ARKIMEX()
1591 …PetscCall(TSComputeIJacobianP(ts, ark->stage_time, Y[i], Ydot, ark->scoeff / adjoint_time_step, ts… in TSAdjointStep_ARKIMEX()
1891 PetscCall(TSARKIMEXComputeAlgebraicIS(ts, ark->stage_time, Z, &alg_is)); in SNESTSFormFunction_ARKIMEX()
1896 PetscCall(TSComputeIFunction(ts, ark->stage_time, Z, X, F, ark->imex)); in SNESTSFormFunction_ARKIMEX()
1901 PetscCall(TSComputeIFunction(ts, ark->stage_time, X, Ydot, F, ark->imex)); in SNESTSFormFunction_ARKIMEX()
1933 PetscCall(TSComputeIJacobian(ts, ark->stage_time, Z, X, shift, A, B, ark->imex)); in SNESTSFormJacobian_ARKIMEX()
[all …]
/petsc/src/ts/impls/implicit/glle/
H A Dglle.h65 PetscReal stage_time; /* time at current stage */ member
H A Dglle.c848 gl->stage_time = ts->ptime + ts->time_step; in TSSolve_GLLE()
903 gl->stage_time = ts->ptime + PetscRealPart(c[i]) * h; in TSSolve_GLLE()
943 gl->stage_time = ts->ptime + ts->time_step; in TSSolve_GLLE()
951 PetscCall((*gl->Accept)(ts, ts->max_time - gl->stage_time, h, enorm, &accept)); in TSSolve_GLLE()
954 …MT " (t=%g) not accepted, rejections=%" PetscInt_FMT "\n", k, (double)gl->stage_time, rejections)); in TSSolve_GLLE()
967 …MT " (t=%g) not accepted after %" PetscInt_FMT " failures", k, (double)gl->stage_time, rejections); in TSSolve_GLLE()
1061 PetscCall(TSComputeIFunction(ts, gl->stage_time, x, Ydot, f, PETSC_FALSE)); in SNESTSFormFunction_GLLE()
1079 …PetscCall(TSComputeIJacobian(ts, gl->stage_time, x, gl->Ydot[gl->stage], gl->scoeff / ts->time_ste… in SNESTSFormJacobian_GLLE()
/petsc/src/ts/impls/glee/
H A Dglee.c60 PetscReal stage_time; member
526 glee->stage_time = t + h * c[i]; in TSStep_GLEE()
527 PetscCall(TSPreStage(ts, glee->stage_time)); in TSStep_GLEE()
554 PetscCall(TSAdaptCheckStage(adapt, ts, glee->stage_time, YStage[i], &accept)); in TSStep_GLEE()
556 PetscCall(TSPostStage(ts, glee->stage_time, i, YStage)); in TSStep_GLEE()
684 PetscCall(TSComputeIFunction(ts, glee->stage_time, X, Ydot, F, PETSC_FALSE)); in SNESTSFormFunction_GLEE()
705 PetscCall(TSComputeIJacobian(ts, glee->stage_time, X, Ydot, shift, A, B, PETSC_FALSE)); in SNESTSFormJacobian_GLEE()
/petsc/src/ts/impls/rosw/
H A Drosw.c60 PetscReal stage_time; member
1150 ros->stage_time = ts->ptime + h * ASum[i]; in TSStep_RosW()
1151 PetscCall(TSPreStage(ts, ros->stage_time)); in TSStep_RosW()
1188 PetscCall(TSComputeIFunction(ts, ros->stage_time, Zstage, Ydot, Y[i], PETSC_FALSE)); in TSStep_RosW()
1198 … PetscCall(TSComputeIJacobian(ts, ros->stage_time, ts->vec_sol, Ydot, 0, J, Jp, PETSC_FALSE)); in TSStep_RosW()
1205 PetscCall(TSPostStage(ts, ros->stage_time, i, Y)); in TSStep_RosW()
1207 PetscCall(TSAdaptCheckStage(adapt, ts, ros->stage_time, Y[i], &stageok)); in TSStep_RosW()
1433 PetscCall(TSComputeIFunction(ts, ros->stage_time, Ystage, Ydot, F, PETSC_FALSE)); in SNESTSFormFunction_RosW()
1452 PetscCall(TSComputeIJacobian(ts, ros->stage_time, Ystage, Ydot, shift, A, B, PETSC_TRUE)); in SNESTSFormJacobian_RosW()
/petsc/src/ts/impls/multirate/
H A Dmprk.c61 PetscReal stage_time; member
655 mprk->stage_time = t + h * cf[i]; in TSStep_MPRK()
656 PetscCall(TSPreStage(ts, mprk->stage_time)); in TSStep_MPRK()
703 PetscCall(TSPostStage(ts, mprk->stage_time, i, Y)); in TSStep_MPRK()
800 mprk->stage_time = t + h * cf[i]; in TSStep_MPRKSPLIT()
801 PetscCall(TSPreStage(ts, mprk->stage_time)); in TSStep_MPRKSPLIT()