1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556 |
- /*
- [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 <array>
- #include <type_traits>
- #include <boost/numeric/odeint.hpp>
- namespace odeint = boost::numeric::odeint;
- typedef std::array< double , 3 > state_type;
- struct lorenz
- {
- void ode( const state_type &x , state_type &dxdt , double t ) const
- {
- const double sigma = 10.0;
- const double R = 28.0;
- const double b = 8.0 / 3.0;
- dxdt[0] = sigma * ( x[1] - x[0] );
- dxdt[1] = R * x[0] - x[1] - x[0] * x[2];
- dxdt[2] = -b * x[2] + x[0] * x[1];
- }
- };
- int main( int argc , char *argv[] )
- {
- using namespace boost::numeric::odeint;
- //[ bind_member_function_cpp11
- namespace pl = std::placeholders;
- state_type x = {{ 10.0 , 10.0 , 10.0 }};
- integrate_const( runge_kutta4< state_type >() ,
- std::bind( &lorenz::ode , lorenz() , pl::_1 , pl::_2 , pl::_3 ) ,
- x , 0.0 , 10.0 , 0.01 );
- //]
- return 0;
- }
|