Home
last modified time | relevance | path

Searched refs:ptime (Results 1 – 25 of 50) sorted by relevance

12

/petsc/src/ts/impls/explicit/euler/
H A Deuler.c18 PetscCall(TSPreStage(ts, ts->ptime)); in TSStep_Euler()
19 PetscCall(TSComputeRHSFunction(ts, ts->ptime, solution, update)); in TSStep_Euler()
21 PetscCall(TSPostStage(ts, ts->ptime, 0, &solution)); in TSStep_Euler()
22 PetscCall(TSAdaptCheckStage(ts->adapt, ts, ts->ptime, solution, &stageok)); in TSStep_Euler()
27 PetscCall(TSFunctionDomainError(ts, ts->ptime + ts->time_step, update, &stageok)); in TSStep_Euler()
40 ts->ptime += ts->time_step; in TSStep_Euler()
90 PetscReal alpha = (ts->ptime - t) / ts->time_step; in TSInterpolate_Euler()
/petsc/src/ts/impls/symplectic/basicsymplectic/
H A Dbasicsymplectic.c211 PetscReal ptime, next_time_step = ts->time_step; in TSStep_BasicSymplectic() local
218 PetscCall(TSGetTime(ts, &ptime)); in TSStep_BasicSymplectic()
219 PetscCall(TSSetTime(subts_p, ptime)); in TSStep_BasicSymplectic()
220 PetscCall(TSSetTime(subts_q, ptime)); in TSStep_BasicSymplectic()
224 PetscCall(TSPreStage(ts, ptime)); in TSStep_BasicSymplectic()
229 PetscCall(TSComputeRHSFunction(subts_q, ptime, p, q_update)); in TSStep_BasicSymplectic()
234 ptime = ptime + scheme->d[iter] * ts->time_step; in TSStep_BasicSymplectic()
235 PetscCall(TSComputeRHSFunction(subts_p, ptime, q, p_update)); in TSStep_BasicSymplectic()
240 PetscCall(TSPostStage(ts, ptime, 0, &solution)); in TSStep_BasicSymplectic()
241 PetscCall(TSAdaptCheckStage(ts->adapt, ts, ptime, solution, &stageok)); in TSStep_BasicSymplectic()
[all …]
/petsc/src/ts/interface/
H A Dtsmon.c28 PetscErrorCode TSMonitor(TS ts, PetscInt step, PetscReal ptime, Vec u) in TSMonitor() argument
38 PetscCall(DMSetOutputSequenceNumber(dm, step, ptime)); in TSMonitor()
41 for (i = 0; i < n; i++) PetscCall((*ts->monitor[i])(ts, step, ptime, u, ts->monitorcontext[i])); in TSMonitor()
192 PetscErrorCode TSMonitorDefault(TS ts, PetscInt step, PetscReal ptime, Vec v, PetscViewerAndFormat … in TSMonitorDefault() argument
205 …time %g between steps %" PetscInt_FMT " and %" PetscInt_FMT "\n", (double)ptime, ts->steps - 1, ts… in TSMonitorDefault()
207 … PetscInt_FMT " TS dt %g time %g%s", step, (double)ts->time_step, (double)ptime, ts->steprollback … in TSMonitorDefault()
219 PetscCall(PetscRealView(1, &ptime, viewer)); in TSMonitorDefault()
221 PetscCall(PetscRealView(0, &ptime, viewer)); in TSMonitorDefault()
285 PetscErrorCode TSMonitorWallClockTime(TS ts, PetscInt step, PetscReal ptime, Vec v, PetscViewerAndF… in TSMonitorWallClockTime() argument
306 …cInt_FMT " ksp %" PetscInt_FMT "\n", step, (double)ts->time_step, (double)ptime, ts->steprollback … in TSMonitorWallClockTime()
[all …]
H A Dts.c125 …PetscCall(PetscOptionsReal("-ts_init_time", "Initial time", "TSSetTime", ts->ptime, &ts->ptime, NU… in TSSetFromOptions()
2003 PetscCallSAWs(SAWs_Register, (dir, &ts->ptime, 1, SAWs_READ, SAWs_DOUBLE)); in TSView()
3459 …ptime, PetscObjectComm((PetscObject)ts), PETSC_ERR_ARG_OUTOFRANGE, "Requested time %g not in last … in TSInterpolate()
3488 PetscReal ptime; in TSStep() local
3514 if (!ts->steps) ts->ptime_prev = ts->ptime; in TSStep()
3515 ptime = ts->ptime; in TSStep()
3525 ts->ptime_prev = ptime; in TSStep()
3985 …PetscCallBack("TS callback resize setup", (*ts->resizesetup)(ts, ts->steps, ts->ptime, ts->vec_sol… in TSResize()
4080 …PetscBool is_close = PetscIsCloseAtTol(ts->ptime, ts->eval_times->time_points[i], ts->eval_times->… in TSSolve()
4081 if (ts->ptime <= ts->eval_times->time_points[i] || is_close) { in TSSolve()
[all …]
H A Dtscreate.c43 t->ptime = 0.0; in TSCreate()
/petsc/src/sys/mpiuni/
H A Dmpitime.c24 double ptime; in MPI_Wtime() local
58 ptime = (double)SecInTick * dTime; in MPI_Wtime()
59 return ptime; in MPI_Wtime()
/petsc/src/ts/utils/
H A Ddmdats.c48 static PetscErrorCode TSComputeIFunction_DMDA(TS ts, PetscReal ptime, Vec X, Vec Xdot, Vec F, Petsc… in TSComputeIFunction_DMDA() argument
75 PetscCall((*dmdats->ifunctionlocal)(&info, ptime, x, xdot, f, dmdats->ifunctionlocalctx)); in TSComputeIFunction_DMDA()
85 PetscCall((*dmdats->ifunctionlocal)(&info, ptime, x, xdot, f, dmdats->ifunctionlocalctx)); in TSComputeIFunction_DMDA()
103 static PetscErrorCode TSComputeIJacobian_DMDA(TS ts, PetscReal ptime, Vec X, Vec Xdot, PetscReal sh… in TSComputeIJacobian_DMDA() argument
126 …PetscCall((*dmdats->ijacobianlocal)(&info, ptime, x, xdot, shift, A, B, dmdats->ijacobianlocalctx)… in TSComputeIJacobian_DMDA()
140 static PetscErrorCode TSComputeRHSFunction_DMDA(TS ts, PetscReal ptime, Vec X, Vec F, PetscCtx ctx) in TSComputeRHSFunction_DMDA() argument
163 PetscCall((*dmdats->rhsfunctionlocal)(&info, ptime, x, f, dmdats->rhsfunctionlocalctx)); in TSComputeRHSFunction_DMDA()
173 PetscCall((*dmdats->rhsfunctionlocal)(&info, ptime, x, f, dmdats->rhsfunctionlocalctx)); in TSComputeRHSFunction_DMDA()
189 static PetscErrorCode TSComputeRHSJacobian_DMDA(TS ts, PetscReal ptime, Vec X, Mat A, Mat B, PetscC… in TSComputeRHSJacobian_DMDA() argument
208 PetscCall((*dmdats->rhsjacobianlocal)(&info, ptime, x, A, B, dmdats->rhsjacobianlocalctx)); in TSComputeRHSJacobian_DMDA()
[all …]
H A Ddmnetworkts.c109 PetscErrorCode TSMonitorLGCtxNetworkSolution(TS ts, PetscInt step, PetscReal ptime, Vec u, void *dc… in TSMonitorLGCtxNetworkSolution() argument
159 PetscCall(PetscDrawLGAddCommonPoint(ctx->lg[i], ptime, xv + offset)); in TSMonitorLGCtxNetworkSolution()
170 PetscCall(PetscDrawLGAddCommonPoint(ctx->lg[i], ptime, xv + offset)); in TSMonitorLGCtxNetworkSolution()
/petsc/src/ts/impls/mimex/
H A Dmimex.c113 PetscCall(TSComputeRHSFunction(ts, ts->ptime, Xstar, G)); in SNESTSFormFunction_Mimex()
160 PetscCall(TSPreStage(ts, ts->ptime)); in TSStep_Mimex_Split()
162 mimex->stage_time = ts->ptime + ts->time_step; in TSStep_Mimex_Split()
186 PetscCall(TSComputeRHSFunction(ts, ts->ptime, sol, update)); in TSStep_Mimex_Split()
207 PetscCall(TSPostStage(ts, ts->ptime, 0, &sol)); in TSStep_Mimex_Split()
208 ts->ptime += ts->time_step; in TSStep_Mimex_Split()
220 PetscCall(TSPreStage(ts, ts->ptime)); in TSStep_Mimex_Implicit()
222 mimex->stage_time = ts->ptime + ts->time_step; in TSStep_Mimex_Implicit()
223 ts->ptime += ts->time_step; in TSStep_Mimex_Implicit()
227 PetscCall(TSPostStage(ts, ts->ptime, 0, &sol)); in TSStep_Mimex_Implicit()
[all …]
/petsc/src/ts/impls/explicit/rk/
H A Dmrk.c41 t = (itime - rk->ptime) / h; in TSInterpolate_RK_MultirateNonsplit()
64 PetscReal t = ts->ptime, h = ts->time_step; in TSStepRefine_RK_MultirateNonsplit()
96 ts->ptime = t + k * h / rk->dtratio; in TSStepRefine_RK_MultirateNonsplit()
107 ts->ptime = t; in TSStepRefine_RK_MultirateNonsplit()
130 PetscReal next_time_step = ts->time_step, t = ts->ptime, h = ts->time_step; in TSStep_RK_MultirateNonsplit()
159 rk->ptime = t; in TSStep_RK_MultirateNonsplit()
163 ts->ptime = t + ts->time_step; in TSStep_RK_MultirateNonsplit()
249 t = (itime - ts->ptime) / h; in TSInterpolate_RK_MultirateSplit()
252 h = ts->ptime - ts->ptime_prev; in TSInterpolate_RK_MultirateSplit()
253 t = (itime - ts->ptime) / h + 1; /* In the interval [0,1] */ in TSInterpolate_RK_MultirateSplit()
[all …]
H A Drk.h38 PetscReal ptime; member
H A Drk.c563 h = ts->ptime - ts->ptime_prev; in TSEvaluateStep_RK()
610 …PetscCall(TSComputeRHSFunction(quadts, rk->ptime + rk->time_step * c[i], Y[i], ts->vec_costintegra… in TSForwardCostIntegral_RK()
629 …PetscCall(TSComputeRHSFunction(quadts, ts->ptime + ts->time_step * (1.0 - c[i]), Y[i], ts->vec_cos… in TSAdjointCostIntegral_RK()
654 h = ts->ptime - ts->ptime_prev; in TSRollBack_RK()
664 PetscCall(TSComputeRHSFunction(quadts, rk->ptime + h * c[j], Y[j], ts->vec_costintegrand)); in TSRollBack_RK()
688 stage_time = ts->ptime + h * c[i]; in TSForwardStep_RK()
797 PetscReal t = ts->ptime; in TSStep_RK()
828 rk->ptime = ts->ptime; in TSStep_RK()
832 ts->ptime += ts->time_step; in TSStep_RK()
884 PetscReal t = ts->ptime; in TSAdjointStep_RK()
[all …]
/petsc/src/ts/impls/pseudo/
H A Dposindep.c73 PetscCall(TSComputeIFunction(ts, ts->ptime, solution, pseudo->xdot, pseudo->func, PETSC_FALSE)); in TSPseudoComputeFunction()
95 PetscCall(TSPreStage(ts, ts->ptime + ts->time_step)); in TSStep_Pseudo()
102 PetscCall(TSPostStage(ts, ts->ptime + ts->time_step, 0, &ts->vec_sol)); in TSStep_Pseudo()
104 PetscCall(TSAdaptCheckStage(adapt, ts, ts->ptime + ts->time_step, ts->vec_sol, &accept)); in TSStep_Pseudo()
114 ts->ptime += ts->time_step; in TSStep_Pseudo()
229 PetscCall(TSComputeIFunction(ts, ts->ptime + ts->time_step, X, pseudo->xdot, Y, PETSC_FALSE)); in SNESTSFormFunction_Pseudo()
250 …PetscCall(TSComputeIJacobian(ts, ts->ptime + ts->time_step, X, Xdot, 1. / ts->time_step, AA, BB, P… in SNESTSFormJacobian_Pseudo()
264 static PetscErrorCode TSPseudoMonitorDefault(TS ts, PetscInt step, PetscReal ptime, Vec v, void *du… in TSPseudoMonitorDefault() argument
272 …Int_FMT " dt %g time %g fnorm %g\n", step, (double)ts->time_step, (double)ptime, (double)pseudo->f… in TSPseudoMonitorDefault()
/petsc/src/ts/impls/eimex/
H A Deimex.c60 ext->ctime = ts->ptime + h * i; in TSStage_EIMEX()
141 ts->ptime += ts->time_step; in TSStep_EIMEX()
154 const PetscReal h = ts->ptime - ts->ptime_prev; in TSInterpolate_EIMEX()
157 t = (itime - ts->ptime + h) / h; in TSInterpolate_EIMEX()
161 PetscCall(TSComputeIFunction(ts, ts->ptime - h, Y0, Ydot, YdotI, PETSC_FALSE)); in TSInterpolate_EIMEX()
167 PetscCall(TSComputeIFunction(ts, ts->ptime, Y1, Ydot, YdotI, PETSC_FALSE)); in TSInterpolate_EIMEX()
289 PetscCall(TSComputeIJacobian(ts, ts->ptime, X, Ydot, ext->shift, A, B, PETSC_TRUE)); in SNESTSFormJacobian_EIMEX()
/petsc/src/ts/interface/sensitivity/
H A Dtssen.c1237 static PetscErrorCode TSAdjointMonitorSensi(TS ts, PetscInt step, PetscReal ptime, Vec v, PetscInt … in TSAdjointMonitorSensi() argument
1426 PetscErrorCode TSAdjointMonitorDrawSensi(TS ts, PetscInt step, PetscReal ptime, Vec u, PetscInt num… in TSAdjointMonitorDrawSensi() argument
1438 PetscCall(PetscSNPrintf(time, 32, "Timestep %" PetscInt_FMT " Time %g", step, (double)ptime)); in TSAdjointMonitorDrawSensi()
1518 ts->ptime_prev = ts->ptime; in TSAdjointStep()
1590 PetscCall(TSTrajectoryGet(ts->trajectory, ts, ts->steps, &ts->ptime)); in TSAdjointSolve()
1591 …PetscCall(TSAdjointMonitor(ts, ts->steps, ts->ptime, ts->vec_sol, ts->numcost, ts->vecs_sensi, ts-… in TSAdjointSolve()
1597 PetscCall(TSTrajectoryGet(ts->trajectory, ts, ts->steps, &ts->ptime)); in TSAdjointSolve()
1598 …PetscCall(TSAdjointMonitor(ts, ts->steps, ts->ptime, ts->vec_sol, ts->numcost, ts->vecs_sensi, ts-… in TSAdjointSolve()
1600 ts->solvetime = ts->ptime; in TSAdjointSolve()
1632 PetscErrorCode TSAdjointMonitor(TS ts, PetscInt step, PetscReal ptime, Vec u, PetscInt numcost, Vec… in TSAdjointMonitor() argument
[all …]
/petsc/src/ts/impls/arkimex/
H A Dfsarkimex.c88 h_prev = ts->ptime - ts->ptime_prev; in TSExtrapolate_ARKIMEX_FastSlowSplit()
137 h = ts->ptime - ts->ptime_prev; in TSEvaluateStep_ARKIMEX_FastSlowSplit()
260 PetscCall(TSSetTime(ts_start, ts->ptime)); in TSStep_ARKIMEX_FastSlowSplit()
262 PetscCall(TSSetMaxTime(ts_start, ts->ptime + ts->time_step)); in TSStep_ARKIMEX_FastSlowSplit()
271 PetscCall(TSGetTime(ts_start, &ts->ptime)); in TSStep_ARKIMEX_FastSlowSplit()
290 PetscReal t = ts->ptime; in TSStep_ARKIMEX_FastSlowSplit()
392 ts->ptime += ts->time_step; in TSStep_ARKIMEX_FastSlowSplit()
/petsc/src/ts/impls/explicit/ssp/
H A Dssp.c195 PetscCall((*ssp->onestep)(ts, ts->ptime, ts->time_step, sol)); in TSStep_SSP()
196 PetscCall(TSPostStage(ts, ts->ptime, 0, &sol)); in TSStep_SSP()
197 PetscCall(TSAdaptCheckStage(ts->adapt, ts, ts->ptime + ts->time_step, sol, &stageok)); in TSStep_SSP()
209 ts->ptime += ts->time_step; in TSStep_SSP()
/petsc/src/ts/impls/implicit/alpha/
H A Dalpha1.c57 PetscReal t = ts->ptime; in TSAlpha_StageTime()
130 th->stage_time = ts->ptime; in TSAlpha_Restart()
220 ts->ptime += ts->time_step; in TSStep_Alpha()
256 PetscReal h = ts->time_step, h_prev = ts->ptime - ts->ptime_prev; in TSEvaluateWLTE_Alpha()
277 PetscReal dt = t - ts->ptime; in TSInterpolate_Alpha()
/petsc/src/ts/impls/implicit/theta/
H A Dtheta.c134 PetscCall(TSComputeRHSFunction(quadts, ts->ptime, ts->vec_sol, ts->vec_costintegrand)); in TSThetaEvaluateCostIntegral()
161 th->ptime0 = ts->ptime + ts->time_step; in TSAdjointCostIntegral_Theta()
213 th->stage_time = ts->ptime + (th->endpoint ? (PetscReal)1 : th->Theta) * ts->time_step; in TSStep_Theta()
219 PetscCall(TSComputeIFunction(ts, ts->ptime, th->X0, th->Xdot, th->affine, PETSC_FALSE)); in TSStep_Theta()
235 … PetscCall(TSAdaptCheckStage(ts->adapt, ts, ts->ptime + ts->time_step, ts->vec_sol, &stageok)); in TSStep_Theta()
253 th->ptime0 = ts->ptime; in TSStep_Theta()
256 ts->ptime += ts->time_step; in TSStep_Theta()
300 th->stage_time = ts->ptime; in TSAdjointStepBEuler_Private()
442 th->stage_time = th->endpoint ? ts->ptime : (ts->ptime + (1. - th->Theta) * ts->time_step); in TSAdjointStep_Theta()
443 adjoint_ptime = ts->ptime + ts->time_step; in TSAdjointStep_Theta()
[all …]
/petsc/src/ts/impls/implicit/irk/
H A Dirk.c285 h = ts->ptime - ts->ptime_prev; in TSEvaluateStep_IRK()
351 ts->ptime += ts->time_step; in TSStep_IRK()
380 t = (itime - ts->ptime) / h; in TSInterpolate_IRK()
383 h = ts->ptime - ts->ptime_prev; in TSInterpolate_IRK()
384 t = (itime - ts->ptime) / h + 1; /* In the interval [0,1] */ in TSInterpolate_IRK()
478 …PetscCall(TSComputeIFunction(ts, ts->ptime + ts->time_step * c[i], Y[i], Ydot, YdotI[i], PETSC_FAL… in SNESTSFormFunction_IRK()
514 …PetscCall(TSComputeIJacobian(ts, ts->ptime + ts->time_step * c[nstages - 1], Y[nstages - 1], Ydot,… in SNESTSFormJacobian_IRK()
/petsc/src/ts/tutorials/
H A Dex7.c232 PetscErrorCode MyTSMonitor(TS ts, PetscInt step, PetscReal ptime, Vec v, PetscCtx ctx) in MyTSMonitor() argument
241 …PetscCall(PetscPrintf(comm, "timestep %" PetscInt_FMT " time %g norm %g\n", step, (double)ptime, (… in MyTSMonitor()
H A Dex12.c247 PetscErrorCode MyTSMonitor(TS ts, PetscInt step, PetscReal ptime, Vec v, PetscCtx ctx) in MyTSMonitor() argument
256 …PetscCall(PetscPrintf(comm, "timestep %" PetscInt_FMT " time %g norm %g\n", step, (double)ptime, (… in MyTSMonitor()
/petsc/src/ts/impls/bdf/
H A Dbdf.c225 PetscCall(TSBDF_Advance(ts, ts->ptime, ts->vec_sol)); in TSBDF_Restart()
230 bdf->time[0] = ts->ptime + ts->time_step / 2; in TSBDF_Restart()
260 PetscCall(TSBDF_Advance(ts, ts->ptime, ts->vec_sol)); in TSStep_BDF()
270 bdf->time[0] = ts->ptime + ts->time_step; in TSStep_BDF()
290 ts->ptime += ts->time_step; in TSStep_BDF()
/petsc/src/ts/impls/glee/
H A Dglee.c449 h = ts->ptime - ts->ptime_prev; in TSEvaluateStep_GLEE()
517 t = ts->ptime; in TSStep_GLEE()
569 ts->ptime += ts->time_step; in TSStep_GLEE()
575 …cCall(PetscObjectComposedDataSetReal((PetscObject)ts->vec_sol, explicit_stage_time_id, ts->ptime)); in TSStep_GLEE()
604 t = (itime - ts->ptime) / h; in TSInterpolate_GLEE()
607 h = ts->ptime - ts->ptime_prev; in TSInterpolate_GLEE()
608 t = (itime - ts->ptime) / h + 1; /* In the interval [0,1] */ in TSInterpolate_GLEE()
/petsc/src/ts/adapt/interface/
H A Dtsadapt.c960 PetscReal t = ts->ptime + ts->time_step, tend, tmax, h1, hmax; in TSAdaptChoose()
1005 (double)ts->ptime, (double)h, (double)*next_h)); in TSAdaptChoose()
1008 … (double)ts->ptime, (double)h, (double)*next_h, (double)wlte, (double)wltea, (double)wlter)); in TSAdaptChoose()
1087 … (double)ts->ptime, (double)ts->time_step, ts->num_snes_failures)); in TSAdaptCheckStage()
1116 (double)ts->ptime, (double)dt, (double)new_dt)); in TSAdaptCheckStage()

12