remove Ros2 (Verwer 1999)

Author: BenWibking

integration/Rosenbrock/_parameters

@@ -8,8 +8,7 @@ X_reject_buffer              real         1.0
 #   0 = Rodas5P 8-stage method from DifferentialEquations.jl
 #   1 = Rodas4P 6-stage method from DifferentialEquations.jl
 #   2 = Rodas3P 5-stage method from DifferentialEquations.jl
-#   3 = Hamkar et al. (2012) 3-stage ROS2S tableau
-#   4 = Verwer et al. (1999) 2-stage ROS2 tableau
+#   3 = ROS2S   3-stage method from Hamkar et al. (2012)
 rosenbrock_tableau           int          0
 
 # H211b timestep controller parameters
@@ -18,4 +17,3 @@ h211b_k                      real         2.5
 h211b_fac_min                real         0.2
 h211b_fac_max                real         6.0
 h211b_reduction_fac          real         0.5
-

integration/Rosenbrock/rosenbrock_integrator.H

@@ -597,9 +597,6 @@ int rosenbrock_integrator (BurnT& state, rosenbrock_t<integrator_neqs<BurnT>()>&
     if (integrator_rp::rosenbrock_tableau == 3) {
         return rosenbrock_integrator<BurnT, rosenbrock::ros2s_tableau>(state, rstate);
     }
-    if (integrator_rp::rosenbrock_tableau == 4) {
-        return rosenbrock_integrator<BurnT, rosenbrock::ros2_tableau>(state, rstate);
-    }
 
     return rosenbrock_integrator<BurnT, rosenbrock::rodas5p_tableau>(state, rstate);
 }

integration/Rosenbrock/rosenbrock_tableau.H

@@ -66,47 +66,6 @@ struct ros2s_tableau {
     }
 };
 
-struct ros2_tableau {
-    static constexpr int stages = 2;
-    static constexpr amrex::Real gamma = 1.7071067811865475_rt;
-
-    AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE
-    static constexpr amrex::Real ct (const int i) {
-        (void) i;
-        return 1.0_rt;
-    }
-
-    AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE
-    static constexpr amrex::Real a (const int i, const int j) {
-        if (i == 2 && j == 1) {
-            return 0.585786437626905_rt;
-        }
-        return 0.0_rt;
-    }
-
-    AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE
-    static constexpr amrex::Real c (const int i, const int j) {
-        if (i == 2 && j == 1) {
-            return -1.1715728752538102_rt;
-        }
-        return 0.0_rt;
-    }
-
-    AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE
-    static constexpr amrex::Real b (const int i) {
-        if (i == 1) {
-            return 0.8786796564403575_rt;
-        }
-        return 0.2928932188134525_rt;
-    }
-
-    AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE
-    static constexpr amrex::Real e (const int i) {
-        (void) i;
-        return 0.2928932188134525_rt;
-    }
-};
-
 struct rodas3p_tableau {
     static constexpr int stages = 5;
     static constexpr amrex::Real gamma = 1.0_rt / 3.0_rt;