134     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()
563     th->stage_time = ts->ptime; /* recover the old value */  in TSAdjointStep_Theta()
678   PetscReal dt = t - ts->ptime;  in TSInterpolate_Theta()
706     PetscReal   h = ts->time_step, h_prev = ts->ptime - ts->ptime_prev;  in TSEvaluateWLTE_Theta()