123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126 |
- /*
- [auto_generated]
- libs/numeric/odeint/examples/bind_member_functions.hpp
- [begin_description]
- tba.
- [end_description]
- Copyright 2012 Karsten Ahnert
- Copyright 2012 Mario Mulansky
- 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 <iostream>
- #include <boost/numeric/odeint.hpp>
- namespace odeint = boost::numeric::odeint;
- typedef boost::array< double , 3 > state_type;
- //[ ode_wrapper
- template< class Obj , class Mem >
- class ode_wrapper
- {
- Obj m_obj;
- Mem m_mem;
- public:
- ode_wrapper( Obj obj , Mem mem ) : m_obj( obj ) , m_mem( mem ) { }
- template< class State , class Deriv , class Time >
- void operator()( const State &x , Deriv &dxdt , Time t )
- {
- (m_obj.*m_mem)( x , dxdt , t );
- }
- };
- template< class Obj , class Mem >
- ode_wrapper< Obj , Mem > make_ode_wrapper( Obj obj , Mem mem )
- {
- return ode_wrapper< Obj , Mem >( obj , mem );
- }
- //]
- template< class Obj , class Mem >
- class observer_wrapper
- {
- Obj m_obj;
- Mem m_mem;
- public:
- observer_wrapper( Obj obj , Mem mem ) : m_obj( obj ) , m_mem( mem ) { }
- template< class State , class Time >
- void operator()( const State &x , Time t )
- {
- (m_obj.*m_mem)( x , t );
- }
- };
- template< class Obj , class Mem >
- observer_wrapper< Obj , Mem > make_observer_wrapper( Obj obj , Mem mem )
- {
- return observer_wrapper< Obj , Mem >( obj , mem );
- }
- //[ bind_member_function
- struct lorenz
- {
- void ode( const state_type &x , state_type &dxdt , double t ) const
- {
- dxdt[0] = 10.0 * ( x[1] - x[0] );
- dxdt[1] = 28.0 * x[0] - x[1] - x[0] * x[2];
- dxdt[2] = -8.0 / 3.0 * x[2] + x[0] * x[1];
- }
- };
- int main( int argc , char *argv[] )
- {
- using namespace boost::numeric::odeint;
- state_type x = {{ 10.0 , 10.0 , 10.0 }};
- integrate_const( runge_kutta4< state_type >() , make_ode_wrapper( lorenz() , &lorenz::ode ) ,
- x , 0.0 , 10.0 , 0.01 );
- return 0;
- }
- //]
- /*
- struct lorenz
- {
- void ode( const state_type &x , state_type &dxdt , double t ) const
- {
- dxdt[0] = 10.0 * ( x[1] - x[0] );
- dxdt[1] = 28.0 * x[0] - x[1] - x[0] * x[2];
- dxdt[2] = -8.0 / 3.0 * x[2] + x[0] * x[1];
- }
- void obs( const state_type &x , double t ) const
- {
- std::cout << t << " " << x[0] << " " << x[1] << " " << x[2] << "\n";
- }
- };
- int main( int argc , char *argv[] )
- {
- using namespace boost::numeric::odeint;
- state_type x = {{ 10.0 , 10.0 , 10.0 }};
- integrate_const( runge_kutta4< state_type >() ,
- make_ode_wrapper( lorenz() , &lorenz::ode ) ,
- x , 0.0 , 10.0 , 0.01 ,
- make_observer_wrapper( lorenz() , &lorenz::obs ) );
- return 0;
- }
- */
|