Subject: [Boost-commit] svn:boost r57707 - sandbox/odeint/libs/numeric/odeint/examples
From: mario.mulansky_at_[hidden]
Date: 2009-11-16 11:49:47
Author: mariomulansky
Date: 2009-11-16 11:49:46 EST (Mon, 16 Nov 2009)
New Revision: 57707
URL: http://svn.boost.org/trac/boost/changeset/57707
Log:
prepare pointer version of generic rk
Text files modified:
sandbox/odeint/libs/numeric/odeint/examples/lorenz_cmp_rk4_rk_generic.cpp | 7 +++++++
1 files changed, 7 insertions(+), 0 deletions(-)
Modified: sandbox/odeint/libs/numeric/odeint/examples/lorenz_cmp_rk4_rk_generic.cpp
==============================================================================
--- sandbox/odeint/libs/numeric/odeint/examples/lorenz_cmp_rk4_rk_generic.cpp (original)
+++ sandbox/odeint/libs/numeric/odeint/examples/lorenz_cmp_rk4_rk_generic.cpp 2009-11-16 11:49:46 EST (Mon, 16 Nov 2009)
@@ -21,6 +21,7 @@
#include <tr1/array>
#include <boost/numeric/odeint.hpp>
+#include <boost/numeric/odeint/stepper_rk_generic_test.hpp>
#define tab "\t"
@@ -73,8 +74,14 @@
c[2] = 1.0/3.0;
c[3] = 1.0/6.0;
+ const double a2[3] = {0.5, 0.5, 1.0 };
+ const double b2[6] = {0.5, 0.0, 0.5, 0.0, 0.0, 1.0};
+ const double c2[4] = {1.0/6.0, 1.0/3.0, 1.0/3.0, 1.0/6.0};
+
stepper_rk_generic< state_type > stepper_generic4(a, b, c);
+ stepper_rk_generic_test< state_type > stepper_generic4_test(a2, b2, c2, 4);
+
clock_t start , end;
double t;