diff options
Diffstat (limited to 'libs/numeric/odeint/performance/generic_odeint_rk4_lorenz.cpp')
-rw-r--r-- | libs/numeric/odeint/performance/generic_odeint_rk4_lorenz.cpp | 84 |
1 files changed, 0 insertions, 84 deletions
diff --git a/libs/numeric/odeint/performance/generic_odeint_rk4_lorenz.cpp b/libs/numeric/odeint/performance/generic_odeint_rk4_lorenz.cpp deleted file mode 100644 index 0d0975646..000000000 --- a/libs/numeric/odeint/performance/generic_odeint_rk4_lorenz.cpp +++ /dev/null @@ -1,84 +0,0 @@ -/* - * Copyright 2011 Mario Mulansky - * Copyright 2012 Karsten Ahnert - * - * Distributed under the Boost Software License, Version 1.0. - * (See accompanying file LICENSE_1_0.txt or - * copy at http://www.boost.org/LICENSE_1_0.txt) - */ - - -#include <boost/array.hpp> - -//#include <boost/numeric/odeint/stepper/explicit_generic_rk.hpp> -#include <boost/numeric/odeint/stepper/runge_kutta4.hpp> -#include <boost/numeric/odeint/algebra/array_algebra.hpp> - -#include "rk_performance_test_case.hpp" - -#include "lorenz.hpp" - -using namespace boost::numeric::odeint; - -typedef boost::array< double , 3 > state_type; - -/* -typedef explicit_generic_rk< 4 , 4 , state_type , double , state_type , double , array_algebra > rk4_type; - -typedef rk4_type::coef_a_type coef_a_type; -typedef rk4_type::coef_b_type coef_b_type; -typedef rk4_type::coef_c_type coef_c_type; - -const boost::array< double , 1 > a1 = {{ 0.5 }}; -const boost::array< double , 2 > a2 = {{ 0.0 , 0.5 }}; -const boost::array< double , 3 > a3 = {{ 0.0 , 0.0 , 1.0 }}; - -const coef_a_type a = fusion::make_vector( a1 , a2 , a3 ); -const coef_b_type b = {{ 1.0/6 , 1.0/3 , 1.0/3 , 1.0/6 }}; -const coef_c_type c = {{ 0.0 , 0.5 , 0.5 , 1.0 }}; -*/ - -typedef runge_kutta4< state_type , double , state_type , double , array_algebra > rk4_type; - - -class rk4_wrapper -{ - -public: - - rk4_wrapper() - // : m_stepper( a , b , c ) - {} - - void reset_init_cond() - { - m_x[0] = 10.0 * rand() / RAND_MAX; - m_x[1] = 10.0 * rand() / RAND_MAX; - m_x[2] = 10.0 * rand() / RAND_MAX; - m_t = 0.0; - } - - inline void do_step( const double dt ) - { - m_stepper.do_step( lorenz(), m_x , m_t , dt ); - } - - double state( const size_t i ) const - { return m_x[i]; } - -private: - state_type m_x; - double m_t; - rk4_type m_stepper; -}; - - - -int main() -{ - srand( 12312354 ); - - rk4_wrapper stepper; - - run( stepper ); -} |