132 …    PetscCall(VecAXPY(quadts->vec_sol, th->time_step0 * (1.0 - th->Theta), 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()
215     if (th->extrapolate && !ts->steprestart) PetscCall(VecAXPY(th->X, 1 / th->shift, th->Xdot));  in TSStep_Theta()
234         PetscCall(VecAXPY(ts->vec_sol, ts->time_step, th->Xdot));  in TSStep_Theta()
312       PetscCall(VecAXPY(VecsSensiTemp[nadj], 1., ts->vec_drdu_col));  in TSAdjointStepBEuler_Private()
345       PetscCall(VecAXPY(VecsSensi2Temp[nadj], -1., ts->vecs_fuu[nadj]));  in TSAdjointStepBEuler_Private()
346       if (ts->vecs_fup) PetscCall(VecAXPY(VecsSensi2Temp[nadj], -1., ts->vecs_fup[nadj]));  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()
467       PetscCall(VecAXPY(VecsSensiTemp[nadj], 1., ts->vec_drdu_col));  in TSAdjointStep_Theta()
508       PetscCall(VecAXPY(VecsSensi2Temp[nadj], -1., ts->vecs_fuu[nadj]));  in TSAdjointStep_Theta()
509       if (ts->vecs_fup) PetscCall(VecAXPY(VecsSensi2Temp[nadj], -1., ts->vecs_fup[nadj]));  in TSAdjointStep_Theta()
536         PetscCall(VecAXPY(ts->vecs_sensi[nadj], -1., ts->vec_drdu_col));  in TSAdjointStep_Theta()
557         PetscCall(VecAXPY(ts->vecs_sensi2[nadj], 1., ts->vecs_fuu[nadj]));  in TSAdjointStep_Theta()
558         if (ts->vecs_fup) PetscCall(VecAXPY(ts->vecs_sensi2[nadj], 1., ts->vecs_fup[nadj]));  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()
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()
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()