1567 …Vec             *Y = ark->Y, Ydot = ark->Ydot, *VecsDeltaLam = ark->VecsDeltaLam, *VecsSensiTemp =…  in TSAdjointStep_ARKIMEX()  local
1601         PetscCall(VecMAXPY(VecsSensiTemp[nadj], s - i - 1, w, &VecsDeltaLam[nadj * s + i + 1]));  in TSAdjointStep_ARKIMEX()
1607         PetscCall(MatMultTranspose(Jim, VecsSensiTemp[nadj], VecsDeltaLam[nadj * s + i]));  in TSAdjointStep_ARKIMEX()
1609 …PetscCall(VecAXPY(VecsDeltaLam[nadj * s + i], ark->scoeff / adjoint_time_step, VecsSensiTemp[nadj]…  in TSAdjointStep_ARKIMEX()
1617         PetscCall(VecSet(VecsDeltaLam[nadj * s + i], 0)); // make sure it is initialized  in TSAdjointStep_ARKIMEX()
1624         PetscCall(VecMAXPY(VecsSensiTemp[nadj], s - i - 1, w, &VecsDeltaLam[nadj * s + i + 1]));  in TSAdjointStep_ARKIMEX()
1630 …PetscCall(MatMultTransposeAdd(Jex, VecsSensiTemp[nadj], VecsDeltaLam[nadj * s + i], VecsDeltaLam[n…  in TSAdjointStep_ARKIMEX()
1640         PetscCall(VecScale(VecsDeltaLam[nadj * s + i], adjoint_time_step));  in TSAdjointStep_ARKIMEX()
1646         PetscCall(VecScale(VecsDeltaLam[nadj * s + i], 1. / At[i * s + i]));  in TSAdjointStep_ARKIMEX()
1647         PetscCall(KSPSolveTranspose(ksp, VecsDeltaLam[nadj * s + i], VecsDeltaLam[nadj * s + i]));  in TSAdjointStep_ARKIMEX()
1655           PetscCall(MatMultTranspose(ts->Jacp, VecsDeltaLam[nadj * s + i], VecsSensiPTemp[nadj]));  in TSAdjointStep_ARKIMEX()
1664     PetscCall(VecMAXPY(ts->vecs_sensi[nadj], s, w, &VecsDeltaLam[nadj * s]));  in TSAdjointStep_ARKIMEX()
1783   PetscCall(VecDestroyVecs(tab->s * ts->numcost, &ark->VecsDeltaLam));  in TSAdjointReset_ARKIMEX()
2054   PetscCall(VecDuplicateVecs(ts->vecs_sensi[0], tab->s * ts->numcost, &ark->VecsDeltaLam));  in TSAdjointSetUp_ARKIMEX()
2406   ark->VecsDeltaLam   = NULL;  in TSCreate_ARKIMEX()