Skip to content
Snippets Groups Projects
Commit 39ee9c06 authored by Dmitri Naumov's avatar Dmitri Naumov
Browse files

[T] ODE: Add calls to ode_solver->getYDot().

parent 08caa966
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Tests/MathLib/TestCVode.cpp
+ 25
− 10
View file @ 39ee9c06
...@@ -80,9 +80,12 @@ TEST(MathLibCVodeTest, Exponential) ...@@ -80,9 +80,12 @@ TEST(MathLibCVodeTest, Exponential)
auto const y = ode_solver->getSolution(); auto const y = ode_solver->getSolution();
double time_reached = ode_solver->getTime(); double time_reached = ode_solver->getTime();
auto const y_dot = ode_solver->getYDot(time_reached, y);
std::printf("t: %14.7g, y: %14.7g, diff: %14.7g\n", time_reached, y[0], std::printf(
y[0] - exp(-15.0 * time_reached)); "t: %14.7g, y: %14.7g, diff: %14.7g, y_dot: %14.7g, diff: %14.7g\n",
time_reached, y[0], y[0] - exp(-15.0 * time_reached), y_dot[0],
y_dot[0] + 15.0 * exp(-15.0 * time_reached));
ASSERT_NEAR(time, time_reached, std::numeric_limits<double>::epsilon()); ASSERT_NEAR(time, time_reached, std::numeric_limits<double>::epsilon());
// std::printf("time: %g\n", time_reached); // std::printf("time: %g\n", time_reached);
...@@ -117,9 +120,12 @@ TEST(MathLibCVodeTest, ExponentialExtraData) ...@@ -117,9 +120,12 @@ TEST(MathLibCVodeTest, ExponentialExtraData)
auto const y = ode_solver->getSolution(); auto const y = ode_solver->getSolution();
double time_reached = ode_solver->getTime(); double time_reached = ode_solver->getTime();
auto const y_dot = ode_solver->getYDot(time_reached, y);
std::printf("t: %14.7g, y: %14.7g, diff: %14.7g\n", time_reached, y[0], std::printf(
y[0] - exp(-15.0 * time_reached)); "t: %14.7g, y: %14.7g, diff: %14.7g, y_dot: %14.7g, diff: %14.7g\n",
time_reached, y[0], y[0] - exp(-15.0 * time_reached), y_dot[0],
y_dot[0] + 15.0 * exp(-15.0 * time_reached));
ASSERT_NEAR(time, time_reached, std::numeric_limits<double>::epsilon()); ASSERT_NEAR(time, time_reached, std::numeric_limits<double>::epsilon());
// std::printf("time: %g\n", time_reached); // std::printf("time: %g\n", time_reached);
...@@ -135,9 +141,12 @@ TEST(MathLibCVodeTest, ExponentialExtraData) ...@@ -135,9 +141,12 @@ TEST(MathLibCVodeTest, ExponentialExtraData)
auto const y = ode_solver->getSolution(); auto const y = ode_solver->getSolution();
double time_reached = ode_solver->getTime(); double time_reached = ode_solver->getTime();
auto const y_dot = ode_solver->getYDot(time_reached, y);
std::printf("t: %14.7g, y: %14.7g, diff: %14.7g\n", time_reached, y[0], std::printf(
y[0] - exp(-15.0 * time_reached)); "t: %14.7g, y: %14.7g, diff: %14.7g, y_dot: %14.7g, diff: %14.7g\n",
time_reached, y[0], y[0] - exp(-15.0 * time_reached), y_dot[0],
y_dot[0] + 15.0 * exp(-15.0 * time_reached));
ASSERT_NEAR(time, time_reached, std::numeric_limits<double>::epsilon()); ASSERT_NEAR(time, time_reached, std::numeric_limits<double>::epsilon());
// std::printf("time: %g\n", time_reached); // std::printf("time: %g\n", time_reached);
...@@ -172,9 +181,12 @@ TEST(MathLibCVodeTest, ExponentialWithJacobian) ...@@ -172,9 +181,12 @@ TEST(MathLibCVodeTest, ExponentialWithJacobian)
auto const y = ode_solver->getSolution(); auto const y = ode_solver->getSolution();
double time_reached = ode_solver->getTime(); double time_reached = ode_solver->getTime();
auto const y_dot = ode_solver->getYDot(time_reached, y);
std::printf("t: %14.7g, y: %14.7g, diff: %14.7g\n", time_reached, y[0], std::printf(
y[0] - exp(-15.0 * time_reached)); "t: %14.7g, y: %14.7g, diff: %14.7g, y_dot: %14.7g, diff: %14.7g\n",
time_reached, y[0], y[0] - exp(-15.0 * time_reached), y_dot[0],
y_dot[0] + 15.0 * exp(-15.0 * time_reached));
ASSERT_NEAR(time, time_reached, std::numeric_limits<double>::epsilon()); ASSERT_NEAR(time, time_reached, std::numeric_limits<double>::epsilon());
// std::printf("time: %g\n", time_reached); // std::printf("time: %g\n", time_reached);
...@@ -212,9 +224,12 @@ TEST(MathLibCVodeTest, ExponentialWithJacobianNewton) ...@@ -212,9 +224,12 @@ TEST(MathLibCVodeTest, ExponentialWithJacobianNewton)
auto const y = ode_solver->getSolution(); auto const y = ode_solver->getSolution();
double time_reached = ode_solver->getTime(); double time_reached = ode_solver->getTime();
auto const y_dot = ode_solver->getYDot(time_reached, y);
std::printf("t: %14.7g, y: %14.7g, diff: %14.7g\n", time_reached, y[0], std::printf(
y[0] - exp(-15.0 * time_reached)); "t: %14.7g, y: %14.7g, diff: %14.7g, y_dot: %14.7g, diff: %14.7g\n",
time_reached, y[0], y[0] - exp(-15.0 * time_reached), y_dot[0],
y_dot[0] + 15.0 * exp(-15.0 * time_reached));
ASSERT_NEAR(time, time_reached, std::numeric_limits<double>::epsilon()); ASSERT_NEAR(time, time_reached, std::numeric_limits<double>::epsilon());
// std::printf("time: %g\n", time_reached); // std::printf("time: %g\n", time_reached);
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment