diff options
author | Tavian Barnes <tavianator@gmail.com> | 2010-10-07 01:54:17 -0400 |
---|---|---|
committer | Tavian Barnes <tavianator@gmail.com> | 2010-10-07 01:54:17 -0400 |
commit | 180d804d9199b79bcad951638105250c77a002bc (patch) | |
tree | c8da91cfcf86d72b6fd3e983e231cf1ff3d20ba8 /tests/CK45.cpp | |
parent | 1f25636c64f0a7c07fa268c0d66c084948b2dd15 (diff) | |
download | vz-180d804d9199b79bcad951638105250c77a002bc.tar.xz |
Add CK45 method.
Diffstat (limited to 'tests/CK45.cpp')
-rw-r--r-- | tests/CK45.cpp | 40 |
1 files changed, 40 insertions, 0 deletions
diff --git a/tests/CK45.cpp b/tests/CK45.cpp new file mode 100644 index 0000000..e726bc6 --- /dev/null +++ b/tests/CK45.cpp @@ -0,0 +1,40 @@ +#include "vZ.hpp" +#include <cmath> +#include <cstdlib> +#include <iostream> +#include <iomanip> + +// y' = y (y == C*exp(t)) +double +f(double t, double y) +{ + return y; +} + +int +main() +{ + vZ::CK45Integrator integrator(f); + integrator.tol(1e-6).y(1.0).x(0.0).h(0.06); + + integrator.integrate(2.0); + + double actual = integrator.y(); + 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 > 2.7e-6) { + std::cerr << "Error: " << 100.0*error << "%" << std::endl; + return EXIT_FAILURE; + } else { + std::cout << "Error: " << 100.0*error << "%" << std::endl; + return EXIT_SUCCESS; + } +} |