gkowal/amun-code
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

EVOLUTION: Do not increase CFL coefficient in some SSPRK methods.

Browse Source 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
Grzegorz Kowal 2024-07-10 20:54:55 -03:00
parent 0e703cf9b5
commit 94ecb53bc4
1 changed files with 0 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
sources/evolution.F90
View File

@ -231,7 +231,6 @@ module evolution
evolve => evolve_ssprk34
order = 3
registers = 2
cfl = 2.0d+00 * cfl
name_int = "3ʳᵈ order 4-step SSPRK(4,3)"
case ("ssprk(5,3)", "SSPRK(5,3)", "rk3.5", "RK3.5")
@ -239,7 +238,6 @@ module evolution
evolve => evolve_ssprk35
order = 3
registers = 3
cfl = 2.65062919143939d+00 * cfl
name_int = "3ʳᵈ order 5-step SSPRK(5,3)"
case ("ssprk(5,4)", "SSPRK(5,4)", "rk4.5", "RK4.5")
@ -247,7 +245,6 @@ module evolution
evolve => evolve_ssprk45
order = 4
registers = 3
cfl = 1.50818d+00 * cfl
name_int = "4th order 5-step SSPRK(5,4)"
case ("ssprk(10,4)", "SSPRK(10,4)", "rk4.10", "RK4.10")
@ -255,7 +252,6 @@ module evolution
evolve => evolve_ssprk410
order = 4
registers = 3
cfl = 6.0d+00 * cfl
name_int = "4th order 10-step SSPRK(10,4)"
case ("ssprk(m,2)", "SSPRK(m,2)", "ssprk2(1)m", "SSPRK2(1)M")
@ -265,7 +261,6 @@ module evolution
order = 2
registers = 3
stages = max(2, min(9, stages))
cfl = (stages - 1) * cfl
write(name_int, "('2ⁿᵈ-order ',i0,'-step embedded " // &
"SSPRK(',i0,',2) method')") stages, stages
@ -298,7 +293,6 @@ module evolution
order = 3
registers = 3
stages = 4
cfl = 2.0d+00 * cfl
name_int = "3ʳᵈ-order 4-step embedded SSP3(2)4 method"
betas(:) = [ 5.5d-01, -2.7d-01, 5.0d-02 ]
@ -321,7 +315,6 @@ module evolution
end do
n = n - 1
stages = max(4, n**2)
cfl = (n - 1) * n * cfl
write(name_int, "('3ʳᵈ-order order ',i0,'-step embedded SSPRK3(2)n² method')") stages
betas(:) = [ 5.5d-01, -2.7d-01, 5.0d-02 ]