860     PetscCall(VecDuplicateVecs(ts->vecs_sensi[0], s * ts->numcost, &rk->VecsDeltaLam2));  in TSAdjointSetUp_RK()
881 …Vec             *VecsDeltaLam2 = rk->VecsDeltaLam2, VecDeltaMu2 = rk->VecDeltaMu2, *VecsSensi2Temp…  in TSAdjointStep_RK()  local
991           PetscCall(VecMAXPY(VecsSensi2Temp[nadj], s - i - 1, w, &VecsDeltaLam2[nadj * s + i + 1]));  in TSAdjointStep_RK()
992           PetscCall(MatMultTranspose(J, VecsSensi2Temp[nadj], VecsDeltaLam2[nadj * s + i]));  in TSAdjointStep_RK()
993           PetscCall(VecScale(VecsDeltaLam2[nadj * s + i], -h * b[i]));  in TSAdjointStep_RK()
994           PetscCall(VecAXPY(VecsDeltaLam2[nadj * s + i], -h * b[i], ts->vecs_guu[nadj]));  in TSAdjointStep_RK()
995 …if (ts->vecs_sensip) PetscCall(VecAXPY(VecsDeltaLam2[nadj * s + i], -h * b[i], ts->vecs_gup[nadj])…  in TSAdjointStep_RK()
998           PetscCall(VecSet(VecsDeltaLam2[nadj * s + i], 0));  in TSAdjointStep_RK()
999           PetscCall(VecMAXPY(VecsSensi2Temp[nadj], s - i - 1, w, &VecsDeltaLam2[nadj * s + i + 1]));  in TSAdjointStep_RK()
1000           PetscCall(MatMultTranspose(J, VecsSensi2Temp[nadj], VecsDeltaLam2[nadj * s + i]));  in TSAdjointStep_RK()
1001           PetscCall(VecScale(VecsDeltaLam2[nadj * s + i], -h));  in TSAdjointStep_RK()
1002           PetscCall(VecAXPY(VecsDeltaLam2[nadj * s + i], -h, ts->vecs_guu[nadj]));  in TSAdjointStep_RK()
1003 …      if (ts->vecs_sensip) PetscCall(VecAXPY(VecsDeltaLam2[nadj * s + i], -h, ts->vecs_gup[nadj]));  in TSAdjointStep_RK()
1025     if (ts->vecs_sensi2) PetscCall(VecMAXPY(ts->vecs_sensi2[nadj], s, w, &VecsDeltaLam2[nadj * s]));  in TSAdjointStep_RK()
1040   PetscCall(VecDestroyVecs(tab->s * ts->numcost, &rk->VecsDeltaLam2));  in TSAdjointReset_RK()