diff options
Diffstat (limited to 'tests')
| -rw-r--r-- | tests/Makefile.am | 2 | ||||
| -rw-r--r-- | tests/Vector.cpp | 40 |
2 files changed, 42 insertions, 0 deletions
diff --git a/tests/Makefile.am b/tests/Makefile.am index 1cd184d..745562c 100644 --- a/tests/Makefile.am +++ b/tests/Makefile.am @@ -28,6 +28,7 @@ check_PROGRAMS = Euler-test \ RKF45-test \ CK45-test \ DP45-test \ + Vector-test \ EquationSystem-test TESTS = $(check_PROGRAMS) @@ -40,4 +41,5 @@ BS23_test_SOURCES = BS23.cpp RKF45_test_SOURCES = RKF45.cpp CK45_test_SOURCES = CK45.cpp DP45_test_SOURCES = DP45.cpp +Vector_test_SOURCES = Vector.cpp EquationSystem_test_SOURCES = EquationSystem.cpp diff --git a/tests/Vector.cpp b/tests/Vector.cpp new file mode 100644 index 0000000..7f53352 --- /dev/null +++ b/tests/Vector.cpp @@ -0,0 +1,40 @@ +#include "vZ.hpp" +#include <cmath> +#include <cstdlib> +#include <iostream> +#include <iomanip> + +// y' = -y (y == C*exp(-t)) +vZ::Vector<3> +f(double t, vZ::Vector<3> y) +{ + return -y; +} + +int +main() +{ + vZ::GenericDP45Integrator<vZ::Vector<3> > integrator(f); + integrator.tol(1e-6).y(vZ::Vector<3>(1.0, 1.0, 1.0)).x(0.0).h(0.06); + + integrator.integrate(2.0); + + double actual = integrator.y().x(); + double expected = std::exp(-2.0); + + std::cout << std::setprecision(10) + << "Numerical: " << actual << std::endl + << "Expected: " << expected << std::endl + << "h: " << integrator.h() << std::endl + << "Iterations: " << integrator.iterations() << std::endl + << "Rejections: " << integrator.rejections() << std::endl; + + double error = std::fabs(expected - actual)/expected; + if (error > 1.5e-6) { + std::cerr << "Error: " << 100.0*error << "%" << std::endl; + return EXIT_FAILURE; + } else { + std::cout << "Error: " << 100.0*error << "%" << std::endl; + return EXIT_SUCCESS; + } +} |