From 5023472727ff68597033eb3a34970a3ad2f5fb5e Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Wed, 6 Nov 2013 23:20:08 +0000
Subject: [PATCH] Cleanup
---
src/ODE/ODESolvers/Rosenbrock21/Rosenbrock21.C | 4 ++--
src/ODE/ODESolvers/Rosenbrock43/Rosenbrock43.C | 2 +-
src/ODE/ODESolvers/rodas43/rodas43.C | 3 ---
src/ODE/ODESolvers/rodas43/rodas43.H | 3 +--
4 files changed, 4 insertions(+), 8 deletions(-)
diff --git a/src/ODE/ODESolvers/Rosenbrock21/Rosenbrock21.C b/src/ODE/ODESolvers/Rosenbrock21/Rosenbrock21.C
index de6b968348a..13371a8e90a 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 04b55dbbe61..061390fdbc9 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 74a3da1d80b..92b78a7ab46 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 a9b8ed7fc2c..e77c0b423b2 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,
--
GitLab