/* * lorenz_gmpxx.cpp * * This example demonstrates how odeint can be used with arbitrary precision types. * * Copyright 2011-2012 Karsten Ahnert * Copyright 2011-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 #include #include #include using namespace std; using namespace boost::numeric::odeint; //[ gmpxx_lorenz typedef mpf_class value_type; typedef boost::array< value_type , 3 > state_type; struct lorenz { void operator()( const state_type &x , state_type &dxdt , value_type t ) const { const value_type sigma( 10.0 ); const value_type R( 28.0 ); const value_type b( value_type( 8.0 ) / value_type( 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]; } }; //] struct streaming_observer { std::ostream& m_out; streaming_observer( std::ostream &out ) : m_out( out ) { } template< class State , class Time > void operator()( const State &x , Time t ) const { m_out << t; for( size_t i=0 ; i() , lorenz() , x , value_type( 0.0 ) , value_type( 10.0 ) , value_type( value_type( 1.0 ) / value_type( 10.0 ) ) , streaming_observer( cout ) ); //] return 0; }