125   TS        quadts = ts->quadraturets;  in TSThetaEvaluateCostIntegral()  local
131       PetscCall(TSComputeRHSFunction(quadts, th->ptime0, th->X0, ts->vec_costintegrand));  in TSThetaEvaluateCostIntegral()
132 …    PetscCall(VecAXPY(quadts->vec_sol, th->time_step0 * (1.0 - th->Theta), ts->vec_costintegrand));  in TSThetaEvaluateCostIntegral()
134     PetscCall(TSComputeRHSFunction(quadts, ts->ptime, ts->vec_sol, ts->vec_costintegrand));  in TSThetaEvaluateCostIntegral()
135     PetscCall(VecAXPY(quadts->vec_sol, th->time_step0 * th->Theta, ts->vec_costintegrand));  in TSThetaEvaluateCostIntegral()
137     PetscCall(TSComputeRHSFunction(quadts, th->stage_time, th->X, ts->vec_costintegrand));  in TSThetaEvaluateCostIntegral()
138     PetscCall(VecAXPY(quadts->vec_sol, th->time_step0, ts->vec_costintegrand));  in TSThetaEvaluateCostIntegral()
146   TS        quadts = ts->quadraturets;  in TSForwardCostIntegral_Theta()  local
150   PetscCall(VecCopy(quadts->vec_sol, th->VecCostIntegral0));  in TSForwardCostIntegral_Theta()
274   TS             quadts       = ts->quadraturets;  in TSAdjointStepBEuler_Private()  local
295   if (quadts) {  in TSAdjointStepBEuler_Private()
296     PetscCall(TSGetRHSJacobian(quadts, &quadJ, NULL, NULL, NULL));  in TSAdjointStepBEuler_Private()
297     PetscCall(TSGetRHSJacobianP(quadts, &quadJp, NULL, NULL));  in TSAdjointStepBEuler_Private()
304   if (quadts) PetscCall(TSComputeRHSJacobian(quadts, th->stage_time, ts->vec_sol, quadJ, NULL));  in TSAdjointStepBEuler_Private()
381     if (quadts) PetscCall(TSComputeRHSJacobianP(quadts, th->stage_time, ts->vec_sol, quadJp));  in TSAdjointStepBEuler_Private()
419   TS           quadts       = ts->quadraturets;  in TSAdjointStep_Theta()  local
437   if (quadts) {  in TSAdjointStep_Theta()
438     PetscCall(TSGetRHSJacobian(quadts, &quadJ, NULL, NULL, NULL));  in TSAdjointStep_Theta()
439     PetscCall(TSGetRHSJacobianP(quadts, &quadJp, NULL, NULL));  in TSAdjointStep_Theta()
453   if (quadts) {  in TSAdjointStep_Theta()
455       PetscCall(TSComputeRHSJacobian(quadts, th->stage_time, ts->vec_sol, quadJ, NULL));  in TSAdjointStep_Theta()
457       PetscCall(TSComputeRHSJacobian(quadts, th->stage_time, th->X, quadJ, NULL));  in TSAdjointStep_Theta()
530     if (quadts) PetscCall(TSComputeRHSJacobian(quadts, adjoint_ptime, th->X0, quadJ, NULL));  in TSAdjointStep_Theta()
570       if (quadts) PetscCall(TSComputeRHSJacobianP(quadts, th->stage_time, ts->vec_sol, quadJp));  in TSAdjointStep_Theta()
605       if (quadts) PetscCall(TSComputeRHSJacobianP(quadts, adjoint_ptime, th->X0, quadJp));  in TSAdjointStep_Theta()
640     if (quadts) PetscCall(TSComputeRHSJacobian(quadts, th->stage_time, th->X, quadJ, NULL));  in TSAdjointStep_Theta()
656       if (quadts) PetscCall(TSComputeRHSJacobianP(quadts, th->stage_time, th->X, quadJp));  in TSAdjointStep_Theta()
728   TS        quadts = ts->quadraturets;  in TSRollBack_Theta()  local
731   if (quadts && ts->costintegralfwd) PetscCall(VecCopy(th->VecCostIntegral0, quadts->vec_sol));  in TSRollBack_Theta()
734 …if (quadts && quadts->mat_sensip) PetscCall(MatCopy(th->MatIntegralSensip0, quadts->mat_sensip, SA…  in TSRollBack_Theta()
741   TS           quadts               = ts->quadraturets;  in TSForwardStep_Theta()  local
754 …if (quadts && quadts->mat_sensip) PetscCall(MatCopy(quadts->mat_sensip, th->MatIntegralSensip0, SA…  in TSForwardStep_Theta()
757   if (quadts) {  in TSForwardStep_Theta()
758     PetscCall(TSGetRHSJacobian(quadts, &quadJ, NULL, NULL, NULL));  in TSForwardStep_Theta()
759     PetscCall(TSGetRHSJacobianP(quadts, &quadJp, NULL, NULL));  in TSForwardStep_Theta()
809     if (quadts && quadts->mat_sensip) {  in TSForwardStep_Theta()
810       PetscCall(TSComputeRHSJacobian(quadts, th->ptime0, th->X0, quadJ, NULL));  in TSForwardStep_Theta()
811       PetscCall(TSComputeRHSJacobianP(quadts, th->ptime0, th->X0, quadJp));  in TSForwardStep_Theta()
814 …PetscCall(MatAXPY(quadts->mat_sensip, th->time_step0 * (1. - th->Theta), th->MatIntegralSensipTemp…  in TSForwardStep_Theta()
845   if (quadts && quadts->mat_sensip) {  in TSForwardStep_Theta()
848       PetscCall(TSComputeRHSJacobian(quadts, th->stage_time, th->X, quadJ, NULL));  in TSForwardStep_Theta()
849       PetscCall(TSComputeRHSJacobianP(quadts, th->stage_time, th->X, quadJp));  in TSForwardStep_Theta()
852 …PetscCall(MatAXPY(quadts->mat_sensip, th->time_step0, th->MatIntegralSensipTemp, SAME_NONZERO_PATT…  in TSForwardStep_Theta()
855       PetscCall(TSComputeRHSJacobian(quadts, th->stage_time, ts->vec_sol, quadJ, NULL));  in TSForwardStep_Theta()
856       PetscCall(TSComputeRHSJacobianP(quadts, th->stage_time, ts->vec_sol, quadJp));  in TSForwardStep_Theta()
859 …PetscCall(MatAXPY(quadts->mat_sensip, th->time_step0 * th->Theta, th->MatIntegralSensipTemp, SAME_…  in TSForwardStep_Theta()
991   TS        quadts = ts->quadraturets;  in TSForwardSetUp_Theta()  local
997   if (quadts && quadts->mat_sensip) {  in TSForwardSetUp_Theta()
998     PetscCall(MatDuplicate(quadts->mat_sensip, MAT_DO_NOT_COPY_VALUES, &th->MatIntegralSensipTemp));  in TSForwardSetUp_Theta()
999     PetscCall(MatDuplicate(quadts->mat_sensip, MAT_DO_NOT_COPY_VALUES, &th->MatIntegralSensip0));  in TSForwardSetUp_Theta()
1011   TS        quadts = ts->quadraturets;  in TSForwardReset_Theta()  local
1014   if (quadts && quadts->mat_sensip) {  in TSForwardReset_Theta()
1027   TS        quadts = ts->quadraturets;  in TSSetUp_Theta()  local
1031   if (!th->VecCostIntegral0 && quadts && ts->costintegralfwd) { /* back up cost integral */  in TSSetUp_Theta()
1032     PetscCall(VecDuplicate(quadts->vec_sol, &th->VecCostIntegral0));  in TSSetUp_Theta()