283   PetscReal      adjoint_time_step;  in TSAdjointStepBEuler_Private()  local
301   adjoint_time_step = -ts->time_step; /* always positive since time_step is negative */  in TSAdjointStepBEuler_Private()
308     PetscCall(VecScale(VecsSensiTemp[nadj], 1. / adjoint_time_step)); /* lambda_{n+1}/h */  in TSAdjointStepBEuler_Private()
319   th->shift = 1. / adjoint_time_step;  in TSAdjointStepBEuler_Private()
344       PetscCall(VecScale(VecsSensi2Temp[nadj], 1. / adjoint_time_step));  in TSAdjointStepBEuler_Private()
369       PetscCall(VecScale(VecsSensiTemp[nadj], -adjoint_time_step));  in TSAdjointStepBEuler_Private()
373         PetscCall(VecScale(VecsSensi2Temp[nadj], -adjoint_time_step));  in TSAdjointStepBEuler_Private()
379 …PetscCall(VecAXPBYPCZ(th->Xdot, -1. / adjoint_time_step, 1.0 / adjoint_time_step, 0, th->X0, ts->v…  in TSAdjointStepBEuler_Private()
380 …PetscCall(TSComputeIJacobianP(ts, th->stage_time, ts->vec_sol, th->Xdot, 1. / adjoint_time_step, t…  in TSAdjointStepBEuler_Private()
391       PetscCall(VecAXPY(ts->vecs_sensip[nadj], -adjoint_time_step, VecsDeltaMu[nadj]));  in TSAdjointStepBEuler_Private()
395         PetscCall(VecAXPY(ts->vecs_sensip[nadj], adjoint_time_step, ts->vec_drdp_col));  in TSAdjointStepBEuler_Private()
401         PetscCall(VecAXPY(ts->vecs_sensi2p[nadj], -adjoint_time_step, VecsDeltaMu2[nadj]));  in TSAdjointStepBEuler_Private()
402 …if (ts->vecs_fpu) PetscCall(VecAXPY(ts->vecs_sensi2p[nadj], -adjoint_time_step, ts->vecs_fpu[nadj]…  in TSAdjointStepBEuler_Private()
403 …if (ts->vecs_fpp) PetscCall(VecAXPY(ts->vecs_sensi2p[nadj], -adjoint_time_step, ts->vecs_fpp[nadj]…  in TSAdjointStepBEuler_Private()
426   PetscReal    adjoint_time_step;  in TSAdjointStep_Theta()  local
444   adjoint_time_step = -ts->time_step; /* always positive since time_step is negative */  in TSAdjointStep_Theta()
463     PetscCall(VecScale(VecsSensiTemp[nadj], 1. / (th->Theta * adjoint_time_step)));  in TSAdjointStep_Theta()
474   th->shift = 1. / (th->Theta * adjoint_time_step);  in TSAdjointStep_Theta()
525     th->shift      = 1. / ((th->Theta - 1.) * adjoint_time_step);  in TSAdjointStep_Theta()
567       th->shift = 1.0 / (adjoint_time_step * th->Theta);  in TSAdjointStep_Theta()
569 …bianP(ts, th->stage_time, ts->vec_sol, th->Xdot, -1. / (th->Theta * adjoint_time_step), ts->Jacp, …  in TSAdjointStep_Theta()
573 …      PetscCall(VecAXPY(ts->vecs_sensip[nadj], -adjoint_time_step * th->Theta, VecsDeltaMu[nadj]));  in TSAdjointStep_Theta()
577 …        PetscCall(VecAXPY(ts->vecs_sensip[nadj], adjoint_time_step * th->Theta, ts->vec_drdp_col));  in TSAdjointStep_Theta()
596 …    PetscCall(VecAXPY(ts->vecs_sensi2p[nadj], -adjoint_time_step * th->Theta, VecsDeltaMu2[nadj]));  in TSAdjointStep_Theta()
597 …if (ts->vecs_fpu) PetscCall(VecAXPY(ts->vecs_sensi2p[nadj], -adjoint_time_step * th->Theta, ts->ve…  in TSAdjointStep_Theta()
598 …if (ts->vecs_fpp) PetscCall(VecAXPY(ts->vecs_sensi2p[nadj], -adjoint_time_step * th->Theta, ts->ve…  in TSAdjointStep_Theta()
604 …ianP(ts, adjoint_ptime, th->X0, th->Xdot, 1. / ((th->Theta - 1.0) * adjoint_time_step), ts->Jacp, …  in TSAdjointStep_Theta()
608 …PetscCall(VecAXPY(ts->vecs_sensip[nadj], -adjoint_time_step * (1.0 - th->Theta), VecsDeltaMu[nadj]…  in TSAdjointStep_Theta()
612 …PetscCall(VecAXPY(ts->vecs_sensip[nadj], adjoint_time_step * (1.0 - th->Theta), ts->vec_drdp_col));  in TSAdjointStep_Theta()
629 …PetscCall(VecAXPY(ts->vecs_sensi2p[nadj], -adjoint_time_step * (1.0 - th->Theta), VecsDeltaMu2[nad…  in TSAdjointStep_Theta()
630 …if (ts->vecs_fpu) PetscCall(VecAXPY(ts->vecs_sensi2p[nadj], -adjoint_time_step * (1.0 - th->Theta)…  in TSAdjointStep_Theta()
631 …if (ts->vecs_fpp) PetscCall(VecAXPY(ts->vecs_sensi2p[nadj], -adjoint_time_step * (1.0 - th->Theta)…  in TSAdjointStep_Theta()
643       PetscCall(VecAXPY(ts->vecs_sensi[nadj], -adjoint_time_step, VecsSensiTemp[nadj]));  in TSAdjointStep_Theta()
647         PetscCall(VecAXPY(ts->vecs_sensi[nadj], adjoint_time_step, ts->vec_drdu_col));  in TSAdjointStep_Theta()
653       th->shift = 1.0 / (adjoint_time_step * th->Theta);  in TSAdjointStep_Theta()
659         PetscCall(VecAXPY(ts->vecs_sensip[nadj], -adjoint_time_step, VecsDeltaMu[nadj]));  in TSAdjointStep_Theta()
663           PetscCall(VecAXPY(ts->vecs_sensip[nadj], adjoint_time_step, ts->vec_drdp_col));  in TSAdjointStep_Theta()