remove set_slow_fast_timestep_ratio and add input file

AMReX-Codes · Jul 17, 2024 · a906cb8 · a906cb8
1 parent ffa94a0
commit a906cb8
Show file tree

Hide file tree

Showing 4 changed files with 8 additions and 48 deletions.
diff --git a/ExampleCodes/SUNDIALS/Reaction-Diffusion/Exec/inputs_forward_euler b/ExampleCodes/SUNDIALS/Reaction-Diffusion/Exec/inputs_forward_euler
diff --git a/ExampleCodes/SUNDIALS/Reaction-Diffusion/Exec/inputs_rk3 b/ExampleCodes/SUNDIALS/Reaction-Diffusion/Exec/inputs_rk3
diff --git a/...S/Reaction-Diffusion/Exec/inputs_sundials → ...action-Diffusion/Exec/inputs_sundials_mri b/...S/Reaction-Diffusion/Exec/inputs_sundials → ...action-Diffusion/Exec/inputs_sundials_mri
@@ -1,3 +1,4 @@
+##
 n_cell = 32
 max_grid_size = 16
 
@@ -6,6 +7,10 @@ plot_int = 100
 
 dt = 1.e-5
 
+# To replicate heat equation
+diffusion_coef = 1.0
+reaction_coef = 1.0
+
 # Use adaptive time stepping and set integrator relative and absolute tolerances
 # adapt_dt = true
 # reltol   = 1.0e-4
@@ -33,12 +38,9 @@ integration.type = SUNDIALS
 # documentation for the supported method names
 
 # Use the SUNDIALS default method for the chosen type (fixed or adaptive step sizes)
-integration.sundials.type = ERK
+integration.sundials.type = EX-MRI # all the types listed above work with respective default methods
 
-# Use forward Euler (fixed step sizes only)
+# The following combination of type and method works as well
 # integration.sundials.type = ERK
-# integration.sundials.method = ARKODE_FORWARD_EULER_1_1
+#integration.sundials.method = ARKODE_MRI_GARK_FORWARD_EULER #ARKODE_FORWARD_EULER_1_1
 
-# Use backward Euler (fixed step sizes only)
-# integration.sundials.type = DIRK
-# integration.sundials.method = ARKODE_BACKWARD_EULER_1_1
diff --git a/ExampleCodes/SUNDIALS/Reaction-Diffusion/Source/main.cpp b/ExampleCodes/SUNDIALS/Reaction-Diffusion/Source/main.cpp
@@ -246,10 +246,6 @@ void main_main ()
         integrator.set_time_step(dt);
     }
 
-    // This sets the ratio of slow timestep size to fast timestep size as an integer,
-    // or equivalently, the number of fast timesteps per slow timestep.
-    integrator.set_slow_fast_timestep_ratio(10);
-
     for (int step = 1; step <= nsteps; ++step)
     {
         // Set time to evolve to