diff options
| author | Tavian Barnes <tavianator@gmail.com> | 2010-10-06 21:52:59 -0400 |
|---|---|---|
| committer | Tavian Barnes <tavianator@gmail.com> | 2010-10-06 21:52:59 -0400 |
| commit | 25cc2bf981d52f0c7688b1e380b795f0cda8a852 (patch) | |
| tree | e88a72cd50f93cbf472ad83d449a26a1d5a342db | |
| parent | f1983adb487bb531c1c68226596b26eebf99876a (diff) | |
| download | vz-25cc2bf981d52f0c7688b1e380b795f0cda8a852.tar.xz | |
Add AdaptiveIntegrator::tol() method to adjust both tolerances at once.
| -rw-r--r-- | src/vZ/Adaptive.hpp | 2 | ||||
| -rw-r--r-- | tests/HE12.cpp | 3 |
2 files changed, 3 insertions, 2 deletions
diff --git a/src/vZ/Adaptive.hpp b/src/vZ/Adaptive.hpp index 864b038..74de0db 100644 --- a/src/vZ/Adaptive.hpp +++ b/src/vZ/Adaptive.hpp @@ -33,6 +33,8 @@ namespace vZ typedef typename GenericRKIntegrator<Y>::Scalar Scalar; typedef typename GenericRKIntegrator<Y>::Function Function; + GenericAdaptiveIntegrator& tol(Scalar tol) + { m_atol = tol; m_rtol = tol; return *this; } GenericAdaptiveIntegrator& atol(Scalar tol) { m_atol = tol; return *this; } GenericAdaptiveIntegrator& rtol(Scalar tol) { m_rtol = tol; return *this; } diff --git a/tests/HE12.cpp b/tests/HE12.cpp index 9a556ef..44d6943 100644 --- a/tests/HE12.cpp +++ b/tests/HE12.cpp @@ -14,9 +14,8 @@ f(double t, double y) int main() { - double tol = 1e-6; vZ::HE12Integrator integrator(f); - integrator.atol(tol).rtol(tol).y(1.0).x(0.0).h(0.02); + integrator.tol(1e-6).y(1.0).x(0.0).h(0.02); integrator.integrate(2.0); |