Make all ARKODE implicit test functions linear.

Forcing functions don't effect stability so they should always be in the
explicit part.

diff --git a/tests/sundials/arkode_03.cc b/tests/sundials/arkode_03.cc
index 9fbdccae8a4a972e523f627078ae0f4984c24ad0..b6a54046b0e579b3376dad931ccaf4e18c9ba949 100644 (file)
--- a/tests/sundials/arkode_03.cc
+++ b/tests/sundials/arkode_03.cc
@@ -76,7 +76,7 @@ main()
     [&](double, const VectorType &y, VectorType &ydot) -> int {
     ydot[0] = 0;
     ydot[1] = 0;
-    ydot[2] = (b - y[2]) / eps;
+    ydot[2] = -y[2] / eps;
     return 0;
   };
 
@@ -85,7 +85,7 @@ main()
     [&](double, const VectorType &y, VectorType &ydot) -> int {
     ydot[0] = a - (y[2] + 1) * y[0] + y[1] * y[0] * y[0];
     ydot[1] = y[2] * y[0] - y[1] * y[0] * y[0];
-    ydot[2] = -y[2] * y[0];
+    ydot[2] = b / eps - y[2] * y[0];
     return 0;
   };
 

diff --git a/tests/sundials/arkode_04.cc b/tests/sundials/arkode_04.cc
index 4b2d0b48916737298c20cf970ddce429ada49279..1f5a284e1dd1e2b491341f540a292f94ccf1f0d6 100644 (file)
--- a/tests/sundials/arkode_04.cc
+++ b/tests/sundials/arkode_04.cc
@@ -79,7 +79,7 @@ main()
     [&](double, const VectorType &y, VectorType &ydot) -> int {
     ydot[0] = 0;
     ydot[1] = 0;
-    ydot[2] = (b - y[2]) / eps;
+    ydot[2] = -y[2] / eps;
     return 0;
   };
 
@@ -88,7 +88,7 @@ main()
     [&](double, const VectorType &y, VectorType &ydot) -> int {
     ydot[0] = a - (y[2] + 1) * y[0] + y[1] * y[0] * y[0];
     ydot[1] = y[2] * y[0] - y[1] * y[0] * y[0];
-    ydot[2] = -y[2] * y[0];
+    ydot[2] = b / eps - y[2] * y[0];
     return 0;
   };
 

diff --git a/tests/sundials/arkode_05.cc b/tests/sundials/arkode_05.cc
index 526eb7fc22b26d683417155ca822598003449c69..40d320f6d8c1af07498420889ff4ab941763b8f4 100644 (file)
--- a/tests/sundials/arkode_05.cc
+++ b/tests/sundials/arkode_05.cc
@@ -78,7 +78,7 @@ main()
     [&](double, const VectorType &y, VectorType &ydot) -> int {
     ydot[0] = 0;
     ydot[1] = 0;
-    ydot[2] = (b - y[2]) / eps;
+    ydot[2] = -y[2] / eps;
     return 0;
   };
 
@@ -87,7 +87,7 @@ main()
     [&](double, const VectorType &y, VectorType &ydot) -> int {
     ydot[0] = a - (y[2] + 1) * y[0] + y[1] * y[0] * y[0];
     ydot[1] = y[2] * y[0] - y[1] * y[0] * y[0];
-    ydot[2] = -y[2] * y[0];
+    ydot[2] = b / eps - y[2] * y[0];
     return 0;
   };
 

diff --git a/tests/sundials/arkode_06.cc b/tests/sundials/arkode_06.cc
index ba1a619afca6602b68da9337ca9414464ab0283f..35908ec72c2b8a8f9fc8a6fbfcf48a8842aa46da 100644 (file)
--- a/tests/sundials/arkode_06.cc
+++ b/tests/sundials/arkode_06.cc
@@ -85,7 +85,7 @@ main()
     [&](double, const VectorType &y, VectorType &ydot) -> int {
     ydot[0] = 0;
     ydot[1] = 0;
-    ydot[2] = (b - y[2]) / eps;
+    ydot[2] = -y[2] / eps;
     return 0;
   };
 
@@ -94,7 +94,7 @@ main()
     [&](double, const VectorType &y, VectorType &ydot) -> int {
     ydot[0] = a - (y[2] + 1) * y[0] + y[1] * y[0] * y[0];
     ydot[1] = y[2] * y[0] - y[1] * y[0] * y[0];
-    ydot[2] = -y[2] * y[0];
+    ydot[2] = b / eps - y[2] * y[0];
     return 0;
   };
 

diff --git a/tests/sundials/arkode_07.cc b/tests/sundials/arkode_07.cc
index 57779f7bab58c16d6e66984a4ba49a8e43df0854..b9cd3ad97a96864b6f51993b4462c9e8ca267abb 100644 (file)
--- a/tests/sundials/arkode_07.cc
+++ b/tests/sundials/arkode_07.cc
@@ -85,7 +85,7 @@ main()
     [&](double, const VectorType &y, VectorType &ydot) -> int {
     ydot[0] = 0;
     ydot[1] = 0;
-    ydot[2] = (b - y[2]) / eps;
+    ydot[2] = -y[2] / eps;
     return 0;
   };
 
@@ -94,7 +94,7 @@ main()
     [&](double, const VectorType &y, VectorType &ydot) -> int {
     ydot[0] = a - (y[2] + 1) * y[0] + y[1] * y[0] * y[0];
     ydot[1] = y[2] * y[0] - y[1] * y[0] * y[0];
-    ydot[2] = -y[2] * y[0];
+    ydot[2] = b / eps - y[2] * y[0];
     return 0;
   };