Cleanup

src/ODE/ODESolvers/Rosenbrock21/Rosenbrock21.C

@@ -42,8 +42,8 @@ const scalar
     Rosenbrock21::b2 = (1.0/2.0)/gamma,
     Rosenbrock21::e1 = b1 - 1.0/gamma,
     Rosenbrock21::e2 = b2,
-    Rosenbrock21::d1 = 1,
-    Rosenbrock21::d2 = -1;
+    Rosenbrock21::d1 = gamma,
+    Rosenbrock21::d2 = -gamma;
 }
 
 

src/ODE/ODESolvers/Rosenbrock43/Rosenbrock43.C

@@ -91,7 +91,7 @@ const scalar
     Rosenbrock43::e3 = 0,
     Rosenbrock43::e4 = 125.0/108.0,
 
-    Rosenbrock43::gamma = 0.5,
+    Rosenbrock43::gamma = 1.0/2.0,
     Rosenbrock43::c2 = 1,
     Rosenbrock43::c3  = 3.0/5.0,
 

src/ODE/ODESolvers/rodas43/rodas43.C

@@ -37,9 +37,6 @@ const scalar
     rodas43::c2 = 0.386,
     rodas43::c3 = 0.21,
     rodas43::c4 = 0.63,
-    rodas43::bet2p = 0.0317,
-    rodas43::bet3p = 0.0635,
-    rodas43::bet4p = 0.3438,
     rodas43::d1 =  0.25,
     rodas43::d2 = -0.1043,
     rodas43::d3 =  0.1035,

src/ODE/ODESolvers/rodas43/rodas43.H

@@ -25,7 +25,7 @@ Class
     Foam::rodas43
 
 Description
-    Stiffly-stable embedded Rosenbrock ODE solver of order (3)4.
+    L-stable, stiffly-accurate embedded Rosenbrock ODE solver of order (3)4.
 
     References:
     \verbatim
@@ -79,7 +79,6 @@ class rodas43
 
         static const scalar
             c2, c3, c4,
-            bet2p, bet3p, bet4p,
             d1, d2, d3, d4,
             a21, a31, a32,
             a41, a42, a43,