857   PetscCall(VecDuplicateVecs(ts->vecs_sensi[0], ts->numcost, &rk->VecsSensiTemp));  in TSAdjointSetUp_RK()
880 … *VecsDeltaLam = rk->VecsDeltaLam, VecDeltaMu = rk->VecDeltaMu, *VecsSensiTemp = rk->VecsSensiTemp;  in TSAdjointStep_RK()  local
915 …PetscCall(VecCopy(ts->vecs_sensi[nadj], VecsSensiTemp[nadj])); /* VecDeltaLam is an vec array of s…  in TSAdjointStep_RK()
916         PetscCall(VecMAXPY(VecsSensiTemp[nadj], s - i - 1, w, &VecsDeltaLam[nadj * s + i + 1]));  in TSAdjointStep_RK()
917 …PetscCall(MatMultTranspose(J, VecsSensiTemp[nadj], VecsDeltaLam[nadj * s + i])); /* VecsSensiTemp …  in TSAdjointStep_RK()
928         PetscCall(VecSet(VecsSensiTemp[nadj], 0));  in TSAdjointStep_RK()
929         PetscCall(VecMAXPY(VecsSensiTemp[nadj], s - i - 1, w, &VecsDeltaLam[nadj * s + i + 1]));  in TSAdjointStep_RK()
930         PetscCall(MatMultTranspose(J, VecsSensiTemp[nadj], VecsDeltaLam[nadj * s + i]));  in TSAdjointStep_RK()
936         PetscCall(MatMultTranspose(ts->Jacprhs, VecsSensiTemp[nadj], VecDeltaMu));  in TSAdjointStep_RK()
958 …PetscCall(TSComputeRHSHessianProductFunctionUU(ts, rk->stage_time, Y[i], VecsSensiTemp, ts->vec_se…  in TSAdjointStep_RK()
965 …PetscCall(TSComputeRHSHessianProductFunctionUP(ts, rk->stage_time, Y[i], VecsSensiTemp, ts->vec_di…  in TSAdjointStep_RK()
973 …PetscCall(TSComputeRHSHessianProductFunctionPU(ts, rk->stage_time, Y[i], VecsSensiTemp, ts->vec_se…  in TSAdjointStep_RK()
975 …PetscCall(TSComputeRHSHessianProductFunctionPP(ts, rk->stage_time, Y[i], VecsSensiTemp, ts->vec_di…  in TSAdjointStep_RK()
1038   PetscCall(VecDestroyVecs(ts->numcost, &rk->VecsSensiTemp));  in TSAdjointReset_RK()