diff options
Diffstat (limited to 'libs/numeric/odeint/test_external/eigen/runge_kutta4.cpp')
-rw-r--r-- | libs/numeric/odeint/test_external/eigen/runge_kutta4.cpp | 94 |
1 files changed, 94 insertions, 0 deletions
diff --git a/libs/numeric/odeint/test_external/eigen/runge_kutta4.cpp b/libs/numeric/odeint/test_external/eigen/runge_kutta4.cpp new file mode 100644 index 0000000000..c4054fd15e --- /dev/null +++ b/libs/numeric/odeint/test_external/eigen/runge_kutta4.cpp @@ -0,0 +1,94 @@ +/* + [auto_generated] + libs/numeric/odeint/test_external/eigen/runge_kutta4.cpp + + [begin_description] + tba. + [end_description] + + Copyright 2013 Karsten Ahnert + Copyright 2013 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 <boost/config.hpp> +#ifdef BOOST_MSVC + #pragma warning(disable:4996) +#endif + +#define BOOST_TEST_MODULE odeint_eigen_runge_kutta4 + +#include <boost/test/unit_test.hpp> + +#include <boost/numeric/odeint/algebra/vector_space_algebra.hpp> +#include <boost/numeric/odeint/stepper/runge_kutta4.hpp> +#include <boost/numeric/odeint/external/eigen/eigen_resize.hpp> + +using namespace boost::unit_test; +using namespace boost::numeric::odeint; + +struct sys +{ + template< class State , class Deriv > + void operator()( const State &x , Deriv &dxdt , double t ) const + { + dxdt[0] = 1.0; + } +}; + + +BOOST_AUTO_TEST_SUITE( eigen_runge_kutta4 ) + +BOOST_AUTO_TEST_CASE( compile_time_matrix ) +{ + typedef Eigen::Matrix< double , 1 , 1 > state_type; + state_type x; + x[0] = 10.0; + runge_kutta4< state_type , double , state_type , double , vector_space_algebra > rk4; + rk4.do_step( sys() , x , 0.0 , 0.1 ); + BOOST_CHECK_CLOSE( x[0] , 10.1 , 1.0e-13 ); + +} + +BOOST_AUTO_TEST_CASE( runtime_matrix ) +{ + typedef Eigen::Matrix< double , Eigen::Dynamic , 1 > state_type; + state_type x( 1 ); + x[0] = 10.0; + runge_kutta4< state_type , double , state_type , double , vector_space_algebra > rk4; + rk4.do_step( sys() , x , 0.0 , 0.1 ); + BOOST_CHECK_CLOSE( x[0] , 10.1 , 1.0e-13 ); +} + + + + + + +BOOST_AUTO_TEST_CASE( compile_time_array ) +{ + typedef Eigen::Array< double , 1 , 1 > state_type; + state_type x; + x[0] = 10.0; + runge_kutta4< state_type , double , state_type , double , vector_space_algebra > rk4; + rk4.do_step( sys() , x , 0.0 , 0.1 ); + BOOST_CHECK_CLOSE( x[0] , 10.1 , 1.0e-13 ); + +} + +BOOST_AUTO_TEST_CASE( runtime_array ) +{ + typedef Eigen::Array< double , Eigen::Dynamic , 1 > state_type; + state_type x( 1 ); + x[0] = 10.0; + runge_kutta4< state_type , double , state_type , double , vector_space_algebra > rk4; + rk4.do_step( sys() , x , 0.0 , 0.1 ); + BOOST_CHECK_CLOSE( x[0] , 10.1 , 1.0e-13 ); +} + + + +BOOST_AUTO_TEST_SUITE_END() |