15   PetscReal    Theta;  member
130     if (th->Theta != 1.0) {  in TSThetaEvaluateCostIntegral()
132 …    PetscCall(VecAXPY(quadts->vec_sol, th->time_step0 * (1.0 - th->Theta), ts->vec_costintegrand));  in TSThetaEvaluateCostIntegral()
135     PetscCall(VecAXPY(quadts->vec_sol, th->time_step0 * th->Theta, ts->vec_costintegrand));  in TSThetaEvaluateCostIntegral()
212     th->shift      = 1 / (th->Theta * ts->time_step);  in TSStep_Theta()
213     th->stage_time = ts->ptime + (th->endpoint ? (PetscReal)1 : th->Theta) * ts->time_step;  in TSStep_Theta()
220       PetscCall(VecScale(th->affine, (th->Theta - 1) / th->Theta));  in TSStep_Theta()
232 …if (th->Theta == 1.0) PetscCall(VecCopy(th->X, ts->vec_sol));              /* BEULER, stage alread…  in TSStep_Theta()
430   if (th->Theta == 1.) {  in TSAdjointStep_Theta()
442   th->stage_time    = th->endpoint ? ts->ptime : (ts->ptime + (1. - th->Theta) * ts->time_step);  in TSAdjointStep_Theta()
448 …PetscCall(VecAXPBYPCZ(th->X0, 1.0 / (1.0 - th->Theta), th->Theta / (th->Theta - 1.0), 0, th->X, ts…  in TSAdjointStep_Theta()
463     PetscCall(VecScale(VecsSensiTemp[nadj], 1. / (th->Theta * adjoint_time_step)));  in TSAdjointStep_Theta()
474   th->shift = 1. / (th->Theta * adjoint_time_step);  in TSAdjointStep_Theta()
525     th->shift      = 1. / ((th->Theta - 1.) * adjoint_time_step);  in TSAdjointStep_Theta()
567       th->shift = 1.0 / (adjoint_time_step * th->Theta);  in TSAdjointStep_Theta()
569 …PetscCall(TSComputeIJacobianP(ts, th->stage_time, ts->vec_sol, th->Xdot, -1. / (th->Theta * adjoin…  in TSAdjointStep_Theta()
573 …      PetscCall(VecAXPY(ts->vecs_sensip[nadj], -adjoint_time_step * th->Theta, VecsDeltaMu[nadj]));  in TSAdjointStep_Theta()
577 …        PetscCall(VecAXPY(ts->vecs_sensip[nadj], adjoint_time_step * th->Theta, ts->vec_drdp_col));  in TSAdjointStep_Theta()
596 …    PetscCall(VecAXPY(ts->vecs_sensi2p[nadj], -adjoint_time_step * th->Theta, VecsDeltaMu2[nadj]));  in TSAdjointStep_Theta()
597 …s_fpu) PetscCall(VecAXPY(ts->vecs_sensi2p[nadj], -adjoint_time_step * th->Theta, ts->vecs_fpu[nadj…  in TSAdjointStep_Theta()
598 …s_fpp) PetscCall(VecAXPY(ts->vecs_sensi2p[nadj], -adjoint_time_step * th->Theta, ts->vecs_fpp[nadj…  in TSAdjointStep_Theta()
604 …PetscCall(TSComputeIJacobianP(ts, adjoint_ptime, th->X0, th->Xdot, 1. / ((th->Theta - 1.0) * adjoi…  in TSAdjointStep_Theta()
608 …PetscCall(VecAXPY(ts->vecs_sensip[nadj], -adjoint_time_step * (1.0 - th->Theta), VecsDeltaMu[nadj]…  in TSAdjointStep_Theta()
612 …PetscCall(VecAXPY(ts->vecs_sensip[nadj], adjoint_time_step * (1.0 - th->Theta), ts->vec_drdp_col));  in TSAdjointStep_Theta()
629 …PetscCall(VecAXPY(ts->vecs_sensi2p[nadj], -adjoint_time_step * (1.0 - th->Theta), VecsDeltaMu2[nad…  in TSAdjointStep_Theta()
630 …PetscCall(VecAXPY(ts->vecs_sensi2p[nadj], -adjoint_time_step * (1.0 - th->Theta), ts->vecs_fpu[nad…  in TSAdjointStep_Theta()
631 …PetscCall(VecAXPY(ts->vecs_sensi2p[nadj], -adjoint_time_step * (1.0 - th->Theta), ts->vecs_fpp[nad…  in TSAdjointStep_Theta()
653       th->shift = 1.0 / (adjoint_time_step * th->Theta);  in TSAdjointStep_Theta()
682   if (th->endpoint) dt *= th->Theta;  in TSInterpolate_Theta()
764     th->shift = 1. / ((th->Theta - 1.) * th->time_step0);  in TSForwardStep_Theta()
767     PetscCall(MatScale(MatDeltaFwdSensip, (th->Theta - 1.) / th->Theta));  in TSForwardStep_Theta()
773 …PetscCall(MatAXPY(MatDeltaFwdSensip, (th->Theta - 1.) / th->Theta, ts->Jacp, SUBSET_NONZERO_PATTER…  in TSForwardStep_Theta()
814 …PetscCall(MatAXPY(quadts->mat_sensip, th->time_step0 * (1. - th->Theta), th->MatIntegralSensipTemp…  in TSForwardStep_Theta()
853 …PetscCall(MatAXPY(ts->mat_sensip, (1. - th->Theta) / th->Theta, MatDeltaFwdSensip, SAME_NONZERO_PA…  in TSForwardStep_Theta()
859 …PetscCall(MatAXPY(quadts->mat_sensip, th->time_step0 * th->Theta, th->MatIntegralSensipTemp, SAME_…  in TSForwardStep_Theta()
862 …if (!th->endpoint) PetscCall(MatAXPY(ts->mat_sensip, 1. / th->Theta, MatDeltaFwdSensip, SAME_NONZE…  in TSForwardStep_Theta()
873     if (!th->endpoint && th->Theta != 1.0) *ns = 1; /* midpoint form */  in TSForwardGetStages_Theta()
877     if (!th->endpoint && th->Theta != 1.0) {  in TSForwardGetStages_Theta()
1039   th->order = (th->Theta == 0.5) ? 2 : 1;  in TSSetUp_Theta()
1040   th->shift = 1 / (th->Theta * ts->time_step);  in TSSetUp_Theta()
1055   ts->stifflyaccurate = (!th->endpoint && th->Theta != 1.0) ? PETSC_FALSE : PETSC_TRUE;  in TSSetUp_Theta()
1090 …s_theta_theta", "Location of stage (0<Theta<=1)", "TSThetaSetTheta", th->Theta, &th->Theta, NULL));  in TSSetFromOptions_Theta()
1106     PetscCall(PetscViewerASCIIPrintf(viewer, "  Theta=%g\n", (double)th->Theta));  in TSView_Theta()
1117   *theta = th->Theta;  in TSThetaGetTheta_Theta()
1127   th->Theta = theta;  in TSThetaSetTheta_Theta()
1128   th->order = (th->Theta == 0.5) ? 2 : 1;  in TSThetaSetTheta_Theta()
1157   f   = (1.0 + (1.0 - th->Theta) * z) / (1.0 - th->Theta * z);  in TSComputeLinearStability_Theta()
1170     if (!th->endpoint && th->Theta != 1.0) *ns = 1; /* midpoint form */  in TSGetStages_Theta()
1174     if (!th->endpoint && th->Theta != 1.0) {  in TSGetStages_Theta()
1278   th->Theta       = 0.5;  in TSCreate_Theta()
1396 …PetscCheck(th->Theta == 1.0, PetscObjectComm((PetscObject)ts), PETSC_ERR_OPT_OVERWRITE, "Can not c…  in TSSetUp_BEuler()
1434 …PetscCheck(th->Theta == 0.5, PetscObjectComm((PetscObject)ts), PETSC_ERR_OPT_OVERWRITE, "Can not c…  in TSSetUp_CN()