diff --git a/src/ODE/ODESolvers/Rosenbrock21/Rosenbrock21.C b/src/ODE/ODESolvers/Rosenbrock21/Rosenbrock21.C index de6b968348ac0b5811d0d99bf44349e91070f4b3..13371a8e90a2bb909aa71107087a1d38606b14e0 100644 --- a/src/ODE/ODESolvers/Rosenbrock21/Rosenbrock21.C +++ b/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; } diff --git a/src/ODE/ODESolvers/Rosenbrock43/Rosenbrock43.C b/src/ODE/ODESolvers/Rosenbrock43/Rosenbrock43.C index 04b55dbbe61db7663fd3a300028e8b8a5f234be5..061390fdbc9e343251ac7d04da11001fdcdbedc4 100644 --- a/src/ODE/ODESolvers/Rosenbrock43/Rosenbrock43.C +++ b/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, diff --git a/src/ODE/ODESolvers/rodas43/rodas43.C b/src/ODE/ODESolvers/rodas43/rodas43.C index 74a3da1d80b26df6c9f2c700907c64560450dea2..92b78a7ab46aa2b21e24ac0802bcb5fce77750c2 100644 --- a/src/ODE/ODESolvers/rodas43/rodas43.C +++ b/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, diff --git a/src/ODE/ODESolvers/rodas43/rodas43.H b/src/ODE/ODESolvers/rodas43/rodas43.H index a9b8ed7fc2c82b23e82ead011c4075f17ec022c1..e77c0b423b285550779b1f19d45309cd3ed0e840 100644 --- a/src/ODE/ODESolvers/rodas43/rodas43.H +++ b/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,