5 static PetscErrorCode RHSFunction(TS ts, PetscReal t, Vec X, Vec F, PetscCtx ctx)  in RHSFunction()  argument
24 PetscErrorCode TestCheckStage(TSAdapt adapt, TS ts, PetscReal t, Vec X, PetscBool *accept)  in TestCheckStage()  argument
29   PetscCall(TSGetStepNumber(ts, &step));  in TestCheckStage()
34 static PetscErrorCode TestExplicitTS(TS ts, PetscInt order, const char subtype[])  in TestExplicitTS()  argument
51   PetscCall(TSGetType(ts, &type));  in TestExplicitTS()
52   PetscCall(TSGetSolution(ts, &U));  in TestExplicitTS()
54   PetscCall(TSSetStepNumber(ts, 0));  in TestExplicitTS()
55   PetscCall(TSSetTime(ts, 0));  in TestExplicitTS()
56   PetscCall(TSSetTimeStep(ts, dt));  in TestExplicitTS()
57   PetscCall(TSSetMaxTime(ts, Tf));  in TestExplicitTS()
58   PetscCall(TSSetExactFinalTime(ts, TS_EXACTFINALTIME_MATCHSTEP));  in TestExplicitTS()
59   PetscCall(TSSolve(ts, NULL));  in TestExplicitTS()
60   PetscCall(TSRollBack(ts));  in TestExplicitTS()
61   PetscCall(TSSolve(ts, NULL));  in TestExplicitTS()
62   PetscCall(TSGetTime(ts, &t));  in TestExplicitTS()
64   PetscCall(TSGetSolution(ts, &U));  in TestExplicitTS()
66   PetscCall(TSEvaluateStep(ts, order, X, NULL));  in TestExplicitTS()
75   PetscCall(TSGetSolution(ts, &U));  in TestExplicitTS()
77   PetscCall(TSEvaluateStep(ts, order - 1, Y, &done));  in TestExplicitTS()
86   PetscCall(TSGetAdapt(ts, &adapt));  in TestExplicitTS()
88   PetscCall(TSSetErrorIfStepFails(ts, PETSC_FALSE));  in TestExplicitTS()
89   PetscCall(TSSetStepNumber(ts, 0));  in TestExplicitTS()
90   PetscCall(TSSetTime(ts, 0));  in TestExplicitTS()
91   PetscCall(TSSetTimeStep(ts, dt));  in TestExplicitTS()
92   PetscCall(TSSolve(ts, NULL));  in TestExplicitTS()
94   PetscCall(TSSetErrorIfStepFails(ts, PETSC_TRUE));  in TestExplicitTS()
95   PetscCall(TSGetStepNumber(ts, &step));  in TestExplicitTS()
96   PetscCall(TSGetConvergedReason(ts, &reason));  in TestExplicitTS()
102 static PetscErrorCode TestTSRK(TS ts, TSRKType type)  in TestTSRK()  argument
111   PetscCall(TSRKSetType(ts, type));  in TestTSRK()
112   PetscCall(TSRKGetType(ts, &type));  in TestTSRK()
113   PetscCall(TSRKGetOrder(ts, &order));  in TestTSRK()
115   PetscCall(TSGetAdapt(ts, &adapt));  in TestTSRK()
123   PetscCall(TestExplicitTS(ts, order, type));  in TestTSRK()
125   PetscCall(TSGetAdapt(ts, &adapt));  in TestTSRK()
132   TS       ts;  in main()  local
139   PetscCall(TSCreate(PETSC_COMM_SELF, &ts));  in main()
140   PetscCall(TSSetType(ts, TSRK));  in main()
141   PetscCall(TSSetRHSFunction(ts, NULL, RHSFunction, NULL));  in main()
143   PetscCall(TSSetSolution(ts, X));  in main()
145   PetscCall(TSSetFromOptions(ts));  in main()
147   PetscCall(TestTSRK(ts, TSRK1FE));  in main()
148   PetscCall(TestTSRK(ts, TSRK2A));  in main()
149   PetscCall(TestTSRK(ts, TSRK3));  in main()
150   PetscCall(TestTSRK(ts, TSRK3BS));  in main()
151   PetscCall(TestTSRK(ts, TSRK4));  in main()
152   PetscCall(TestTSRK(ts, TSRK5F));  in main()
153   PetscCall(TestTSRK(ts, TSRK5DP));  in main()
154   PetscCall(TestTSRK(ts, TSRK5BS));  in main()
155   PetscCall(TestTSRK(ts, TSRK6VR));  in main()
156   PetscCall(TestTSRK(ts, TSRK7VR));  in main()
157   PetscCall(TestTSRK(ts, TSRK8VR));  in main()
159   PetscCall(TSRollBack(ts));  in main()
160   PetscCall(TSDestroy(&ts));  in main()