31 PetscErrorCode RHS_NS(TS ts, PetscReal t, Vec Q, Vec G, void *user_data) {  in RHS_NS()  argument
42   PetscCall(TSGetTimeStep(ts, &dt));  in RHS_NS()
103 PetscErrorCode IFunction_NS(TS ts, PetscReal t, Vec Q, Vec Q_dot, Vec G, void *user_data) {  in IFunction_NS()  argument
117   PetscCall(TSGetTimeStep(ts, &dt));  in IFunction_NS()
152 PetscErrorCode FormIJacobian_NS(TS ts, PetscReal t, Vec Q, Vec Q_dot, PetscReal shift, Mat J, Mat J…  in FormIJacobian_NS()  argument
241 PetscErrorCode TSMonitor_WallForce(TS ts, PetscInt step_no, PetscReal time, Vec Q, void *ctx) {  in TSMonitor_WallForce()  argument
281 // User provided TS Monitor
282 PetscErrorCode TSMonitor_NS(TS ts, PetscInt step_no, PetscReal time, Vec Q, void *ctx) {  in TSMonitor_NS()  argument
296 // TS: Create, setup, and solve
297 …ser user, AppCtx app_ctx, Physics phys, ProblemData problem, Vec *Q, PetscScalar *f_time, TS *ts) {  in TSSolve_NS()  argument
303   PetscCall(TSCreate(comm, ts));  in TSSolve_NS()
304   PetscCall(TSSetDM(*ts, dm));  in TSSolve_NS()
305   PetscCall(TSSetApplicationContext(*ts, user));  in TSSolve_NS()
307     PetscCall(TSSetType(*ts, TSBDF));  in TSSolve_NS()
309       PetscCall(TSSetIFunction(*ts, NULL, IFunction_NS, &user));  in TSSolve_NS()
311       PetscCall(TSSetRHSFunction(*ts, NULL, RHS_NS, &user));  in TSSolve_NS()
318     PetscCall(TSSetType(*ts, TSRK));  in TSSolve_NS()
319     PetscCall(TSRKSetType(*ts, TSRK5F));  in TSSolve_NS()
320     PetscCall(TSSetRHSFunction(*ts, NULL, RHS_NS, &user));  in TSSolve_NS()
322   PetscCall(TSSetMaxTime(*ts, 500. * user->units->second));  in TSSolve_NS()
323   PetscCall(TSSetExactFinalTime(*ts, TS_EXACTFINALTIME_STEPOVER));  in TSSolve_NS()
324   if (app_ctx->test_type == TESTTYPE_NONE) PetscCall(TSSetErrorIfStepFails(*ts, PETSC_FALSE));  in TSSolve_NS()
325   PetscCall(TSSetTimeStep(*ts, 1.e-2 * user->units->second));  in TSSolve_NS()
326   PetscCall(TSGetAdapt(*ts, &adapt));  in TSSolve_NS()
328   PetscCall(TSSetFromOptions(*ts));  in TSSolve_NS()
335       PetscCall(TSGetSNES(*ts, &snes));  in TSSolve_NS()
338       PetscCall(TSSetIJacobian(*ts, user->mat_ijacobian, Pmat, NULL, NULL));  in TSSolve_NS()
355     PetscCall(TSSetTime(*ts, app_ctx->cont_time * user->units->second));  in TSSolve_NS()
356     PetscCall(TSSetStepNumber(*ts, app_ctx->cont_steps));  in TSSolve_NS()
359     PetscCall(TSMonitorSet(*ts, TSMonitor_NS, user, NULL));  in TSSolve_NS()
362     PetscCall(TSMonitorSet(*ts, TSMonitor_WallForce, user, NULL));  in TSSolve_NS()
365     PetscCall(TSMonitorSet(*ts, TSMonitor_TurbulenceStatistics, user, NULL));  in TSSolve_NS()
370 …if (app_ctx->diff_filter_monitor) PetscCall(TSMonitorSet(*ts, TSMonitor_DifferentialFilter, user, …  in TSSolve_NS()
372   if (app_ctx->test_type == TESTTYPE_NONE) PetscCall(PrintRunInfo(user, user->phys, problem, *ts));  in TSSolve_NS()
376   PetscCall(TSGetTime(*ts, &start_time));  in TSSolve_NS()
377   PetscCall(TSGetStepNumber(*ts, &start_step));  in TSSolve_NS()
381   PetscCall(TSSetTime(*ts, start_time));  in TSSolve_NS()
382   PetscCall(TSSetStepNumber(*ts, start_step));  in TSSolve_NS()
390     PetscCall(TSGetSNES(*ts, &snes));  in TSSolve_NS()
393     PetscCall(TSSetSolution(*ts, Q_preload));  in TSSolve_NS()
394     PetscCall(TSStep(*ts));  in TSSolve_NS()
399     PetscCall(PetscBarrier((PetscObject)*ts));  in TSSolve_NS()
400     PetscCall(TSSolve(*ts, *Q));  in TSSolve_NS()
404   PetscCall(TSGetSolveTime(*ts, &final_time));  in TSSolve_NS()
409     PetscCall(TSGetStepNumber(*ts, &step_no));  in TSSolve_NS()