52     else *X0 = ts->vec_sol;  in TSThetaGetX0AndXdot()
132 …    PetscCall(VecAXPY(quadts->vec_sol, th->time_step0 * (1.0 - th->Theta), ts->vec_costintegrand));  in TSThetaEvaluateCostIntegral()
134     PetscCall(TSComputeRHSFunction(quadts, ts->ptime, ts->vec_sol, ts->vec_costintegrand));  in TSThetaEvaluateCostIntegral()
135     PetscCall(VecAXPY(quadts->vec_sol, th->time_step0 * th->Theta, ts->vec_costintegrand));  in TSThetaEvaluateCostIntegral()
138     PetscCall(VecAXPY(quadts->vec_sol, th->time_step0, ts->vec_costintegrand));  in TSThetaEvaluateCostIntegral()
150   PetscCall(VecCopy(quadts->vec_sol, th->VecCostIntegral0));  in TSForwardCostIntegral_Theta()
207     PetscCall(VecCopy(ts->vec_sol, th->X0));  in TSStep_Theta()
217       if (!th->affine) PetscCall(VecDuplicate(ts->vec_sol, &th->affine));  in TSStep_Theta()
229       PetscCall(VecCopy(th->X, ts->vec_sol));  in TSStep_Theta()
232 …if (th->Theta == 1.0) PetscCall(VecCopy(th->X, ts->vec_sol));              /* BEULER, stage alread…  in TSStep_Theta()
234         PetscCall(VecAXPY(ts->vec_sol, ts->time_step, th->Xdot));  in TSStep_Theta()
235 …     PetscCall(TSAdaptCheckStage(ts->adapt, ts, ts->ptime + ts->time_step, ts->vec_sol, &stageok));  in TSStep_Theta()
237           PetscCall(VecCopy(th->X0, ts->vec_sol));  in TSStep_Theta()
247       PetscCall(VecCopy(th->X0, ts->vec_sol));  in TSStep_Theta()
304   if (quadts) PetscCall(TSComputeRHSJacobian(quadts, th->stage_time, ts->vec_sol, quadJ, NULL));  in TSAdjointStepBEuler_Private()
320   PetscCall(TSComputeSNESJacobian(ts, ts->vec_sol, J, Jpre));  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()
363     PetscCall(TSComputeSNESJacobian(ts, ts->vec_sol, J, Jpre));  in TSAdjointStepBEuler_Private()
379 …l(VecAXPBYPCZ(th->Xdot, -1. / adjoint_time_step, 1.0 / adjoint_time_step, 0, th->X0, ts->vec_sol));  in TSAdjointStepBEuler_Private()
380 …PetscCall(TSComputeIJacobianP(ts, th->stage_time, ts->vec_sol, th->Xdot, 1. / adjoint_time_step, t…  in TSAdjointStepBEuler_Private()
381     if (quadts) PetscCall(TSComputeRHSJacobianP(quadts, th->stage_time, ts->vec_sol, quadJp));  in TSAdjointStepBEuler_Private()
384 …PetscCall(TSComputeIHessianProductFunctionPU(ts, th->stage_time, ts->vec_sol, VecsDeltaLam, ts->ve…  in TSAdjointStepBEuler_Private()
386 …PetscCall(TSComputeIHessianProductFunctionPP(ts, th->stage_time, ts->vec_sol, VecsDeltaLam, ts->ve…  in TSAdjointStepBEuler_Private()
448 …ecAXPBYPCZ(th->X0, 1.0 / (1.0 - th->Theta), th->Theta / (th->Theta - 1.0), 0, th->X, ts->vec_sol));  in TSAdjointStep_Theta()
455       PetscCall(TSComputeRHSJacobian(quadts, th->stage_time, ts->vec_sol, quadJ, NULL));  in TSAdjointStep_Theta()
476     PetscCall(TSComputeSNESJacobian(ts, ts->vec_sol, J, Jpre));  in TSAdjointStep_Theta()
500 …PetscCall(TSComputeIHessianProductFunctionUU(ts, th->stage_time, ts->vec_sol, VecsDeltaLam, ts->ve…  in TSAdjointStep_Theta()
504 …PetscCall(TSComputeIHessianProductFunctionUP(ts, th->stage_time, ts->vec_sol, VecsDeltaLam, ts->ve…  in TSAdjointStep_Theta()
568       PetscCall(VecAXPBYPCZ(th->Xdot, -th->shift, th->shift, 0, th->X0, ts->vec_sol));  in TSAdjointStep_Theta()
569 …PetscCall(TSComputeIJacobianP(ts, th->stage_time, ts->vec_sol, th->Xdot, -1. / (th->Theta * adjoin…  in TSAdjointStep_Theta()
570       if (quadts) PetscCall(TSComputeRHSJacobianP(quadts, th->stage_time, ts->vec_sol, quadJp));  in TSAdjointStep_Theta()
587 …PetscCall(TSComputeIHessianProductFunctionPU(ts, th->stage_time, ts->vec_sol, VecsDeltaLam, ts->ve…  in TSAdjointStep_Theta()
592 …PetscCall(TSComputeIHessianProductFunctionPP(ts, th->stage_time, ts->vec_sol, VecsDeltaLam, ts->ve…  in TSAdjointStep_Theta()
681   PetscCall(VecCopy(ts->vec_sol, th->X));  in TSInterpolate_Theta()
690   Vec       X  = ts->vec_sol;      /* X = solution */  in TSEvaluateWLTE_Theta()
731   if (quadts && ts->costintegralfwd) PetscCall(VecCopy(th->VecCostIntegral0, quadts->vec_sol));  in TSRollBack_Theta()
775       PetscCall(VecAXPBYPCZ(th->Xdot, -th->shift, th->shift, 0, th->X0, ts->vec_sol));  in TSForwardStep_Theta()
776 …PetscCall(TSComputeIJacobianP(ts, th->stage_time, ts->vec_sol, th->Xdot, th->shift, ts->Jacp, PETS…  in TSForwardStep_Theta()
797 …PetscCall(TSComputeIJacobian(ts, th->stage_time, ts->vec_sol, th->Xdot, th->shift, J, Jpre, PETSC_…  in TSForwardStep_Theta()
855       PetscCall(TSComputeRHSJacobian(quadts, th->stage_time, ts->vec_sol, quadJ, NULL));  in TSForwardStep_Theta()
856       PetscCall(TSComputeRHSJacobianP(quadts, th->stage_time, ts->vec_sol, quadJp));  in TSForwardStep_Theta()
1004   PetscCall(VecDuplicate(ts->vec_sol, &th->VecDeltaFwdSensipCol));  in TSForwardSetUp_Theta()
1032     PetscCall(VecDuplicate(quadts->vec_sol, &th->VecCostIntegral0));  in TSSetUp_Theta()
1034   if (!th->X) PetscCall(VecDuplicate(ts->vec_sol, &th->X));  in TSSetUp_Theta()
1035   if (!th->Xdot) PetscCall(VecDuplicate(ts->vec_sol, &th->Xdot));  in TSSetUp_Theta()
1036   if (!th->X0) PetscCall(VecDuplicate(ts->vec_sol, &th->X0));  in TSSetUp_Theta()
1037   if (th->endpoint) PetscCall(VecDuplicate(ts->vec_sol, &th->affine));  in TSSetUp_Theta()
1050     if (!th->vec_sol_prev) PetscCall(VecDuplicate(ts->vec_sol, &th->vec_sol_prev));  in TSSetUp_Theta()
1051     if (!th->vec_lte_work) PetscCall(VecDuplicate(ts->vec_sol, &th->vec_lte_work));  in TSSetUp_Theta()
1178       th->Stages[1] = ts->vec_sol; /* stiffly accurate */  in TSGetStages_Theta()