599   TS               quadts = ts->quadraturets;  in TSForwardCostIntegral_RK()  local
610 …PetscCall(TSComputeRHSFunction(quadts, rk->ptime + rk->time_step * c[i], Y[i], ts->vec_costintegra…  in TSForwardCostIntegral_RK()
611     PetscCall(VecAXPY(quadts->vec_sol, rk->time_step * b[i], ts->vec_costintegrand));  in TSForwardCostIntegral_RK()
620   TS               quadts = ts->quadraturets;  in TSAdjointCostIntegral_RK()  local
629 …PetscCall(TSComputeRHSFunction(quadts, ts->ptime + ts->time_step * (1.0 - c[i]), Y[i], ts->vec_cos…  in TSAdjointCostIntegral_RK()
630     PetscCall(VecAXPY(quadts->vec_sol, -ts->time_step * b[i], ts->vec_costintegrand));  in TSAdjointCostIntegral_RK()
638   TS               quadts = ts->quadraturets;  in TSRollBack_RK()  local
661   if (quadts && ts->costintegralfwd) {  in TSRollBack_RK()
664       PetscCall(TSComputeRHSFunction(quadts, rk->ptime + h * c[j], Y[j], ts->vec_costintegrand));  in TSRollBack_RK()
665       PetscCall(VecAXPY(quadts->vec_sol, -h * b[j], ts->vec_costintegrand));  in TSRollBack_RK()
874   TS               quadts = ts->quadraturets;  in TSAdjointStep_RK()  local
891   if (quadts) PetscCall(TSGetRHSJacobian(quadts, &Jquad, NULL, NULL, NULL));  in TSAdjointStep_RK()
899 …if (quadts) PetscCall(TSComputeRHSJacobian(quadts, rk->stage_time, Y[i], Jquad, Jquad)); /* get r_…  in TSAdjointStep_RK()
902 …if (quadts) PetscCall(TSComputeRHSJacobianP(quadts, rk->stage_time, Y[i], quadts->Jacprhs)); /* ge…  in TSAdjointStep_RK()
919         if (quadts) {  in TSAdjointStep_RK()
939           if (quadts) {  in TSAdjointStep_RK()
940             PetscCall(MatDenseGetColumn(quadts->Jacprhs, nadj, &xarr));  in TSAdjointStep_RK()
944             PetscCall(MatDenseRestoreColumn(quadts->Jacprhs, &xarr));  in TSAdjointStep_RK()
959       if (quadts) {  in TSAdjointStep_RK()
961 …PetscCall(TSComputeRHSHessianProductFunctionUU(quadts, rk->stage_time, Y[i], NULL, ts->vec_sensip_…  in TSAdjointStep_RK()
966         if (quadts) {  in TSAdjointStep_RK()
968 …PetscCall(TSComputeRHSHessianProductFunctionUP(quadts, rk->stage_time, Y[i], NULL, ts->vec_sensip_…  in TSAdjointStep_RK()
976         if (b[i] && quadts) {  in TSAdjointStep_RK()
978 …PetscCall(TSComputeRHSHessianProductFunctionPU(quadts, rk->stage_time, Y[i], NULL, ts->vec_sensip_…  in TSAdjointStep_RK()
980 …PetscCall(TSComputeRHSHessianProductFunctionPP(quadts, rk->stage_time, Y[i], NULL, ts->vec_dir, ts…  in TSAdjointStep_RK()
1146   TS quadts = ts->quadraturets;  in TSSetUp_RK()  local
1152   if (quadts && ts->costintegralfwd) {  in TSSetUp_RK()
1154     PetscCall(TSGetRHSJacobian(quadts, &Jquad, NULL, NULL, NULL));  in TSSetUp_RK()