Skip to content

Refactor time_stepper #1012

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
sbryngelson merged 11 commits into from
Oct 13, 2025
Merged

Refactor time_stepper #1012

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
11 commits
Select commit Hold shift + click to select a range
7699ee3
refactor time_stepper routines
Oct 6, 2025
9bcad63
format
Oct 6, 2025
d349c91
Update src/simulation/m_time_steppers.fpp
hyeoksu-lee Oct 6, 2025
bdd0382
Update src/simulation/m_time_steppers.fpp
hyeoksu-lee Oct 6, 2025
1cf0390
fix precision-related bug
Oct 10, 2025
5ab9193
format
Oct 10, 2025
34f2307
remove unused var
Oct 10, 2025
37460c7
fix cantera version compatibility for MacOS runner
Oct 10, 2025
1d236ce
add gpu_update for rk_coef
Oct 10, 2025
1787d0c
format
Oct 10, 2025
d06b2ac
Merge branch 'master' into dts
wilfonba Oct 12, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
10 changes: 2 additions & 8 deletions src/simulation/m_start_up.fpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1149,14 +1149,8 @@ contains
mytime = mytime + dt

! Total-variation-diminishing (TVD) Runge-Kutta (RK) time-steppers
if (time_stepper == 1) then
call s_1st_order_tvd_rk(t_step, time_avg)
elseif (time_stepper == 2) then
call s_2nd_order_tvd_rk(t_step, time_avg)
elseif (time_stepper == 3 .and. (.not. adap_dt)) then
call s_3rd_order_tvd_rk(t_step, time_avg)
elseif (time_stepper == 3 .and. adap_dt) then
call s_strang_splitting(t_step, time_avg)
if (any(time_stepper == (/1, 2, 3/))) then
call s_tvd_rk(t_step, time_avg, time_stepper)
end if

if (relax) call s_infinite_relaxation_k(q_cons_ts(1)%vf)
Expand Down
Loading
Loading