mirror of
https://github.com/boostorg/odeint.git
synced 2026-01-26 06:42:23 +00:00
140 lines
4.9 KiB
C++
140 lines
4.9 KiB
C++
/*
|
|
[auto_generated]
|
|
libs/numeric/odeint/test/resizing.cpp
|
|
|
|
[begin_description]
|
|
This file tests the resizing mechanism of odeint.
|
|
[end_description]
|
|
|
|
Copyright 2009-2012 Karsten Ahnert
|
|
Copyright 2009-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)
|
|
*/
|
|
|
|
// disable checked iterator warning for msvc
|
|
#include <boost/config.hpp>
|
|
#ifdef BOOST_MSVC
|
|
#pragma warning(disable:4996)
|
|
#endif
|
|
|
|
#define BOOST_TEST_MODULE odeint_resize
|
|
|
|
#include <vector>
|
|
#include <cmath>
|
|
|
|
#include <boost/array.hpp>
|
|
#include <boost/bind.hpp>
|
|
#include <boost/utility.hpp>
|
|
#include <boost/type_traits/integral_constant.hpp>
|
|
|
|
#include <boost/test/unit_test.hpp>
|
|
|
|
#include <boost/mpl/vector.hpp>
|
|
#include <boost/mpl/int.hpp>
|
|
#include <boost/mpl/at.hpp>
|
|
|
|
#include <boost/numeric/odeint/stepper/euler.hpp>
|
|
#include <boost/numeric/odeint/stepper/runge_kutta4_classic.hpp>
|
|
#include <boost/numeric/odeint/stepper/runge_kutta4.hpp>
|
|
#include <boost/numeric/odeint/algebra/vector_space_algebra.hpp>
|
|
|
|
#include <boost/numeric/odeint/util/resizer.hpp>
|
|
#include <boost/numeric/odeint/util/is_resizeable.hpp>
|
|
|
|
using namespace boost::unit_test;
|
|
using namespace boost::numeric::odeint;
|
|
|
|
namespace mpl = boost::mpl;
|
|
|
|
size_t adjust_size_count;
|
|
|
|
typedef boost::array< double , 1 > test_array_type;
|
|
|
|
namespace boost { namespace numeric { namespace odeint {
|
|
|
|
|
|
template<>
|
|
struct is_resizeable< test_array_type >
|
|
{
|
|
typedef boost::true_type type;
|
|
const static bool value = type::value;
|
|
};
|
|
|
|
template<>
|
|
struct same_size_impl< test_array_type , test_array_type >
|
|
{
|
|
static bool same_size( const test_array_type &x1 , const test_array_type &x2 )
|
|
{
|
|
return false;
|
|
}
|
|
};
|
|
|
|
template<>
|
|
struct resize_impl< test_array_type , test_array_type >
|
|
{
|
|
static void resize( test_array_type &x1 , const test_array_type &x2 )
|
|
{
|
|
adjust_size_count++;
|
|
}
|
|
};
|
|
|
|
} } }
|
|
|
|
|
|
|
|
void constant_system( const test_array_type &x , test_array_type &dxdt , double t ) { dxdt[0] = 1.0; }
|
|
|
|
|
|
BOOST_AUTO_TEST_SUITE( check_resize_test )
|
|
|
|
|
|
typedef euler< test_array_type , double , test_array_type , double , range_algebra , default_operations , never_resizer > euler_manual_type;
|
|
typedef euler< test_array_type , double , test_array_type , double , range_algebra , default_operations , initially_resizer > euler_initially_type;
|
|
typedef euler< test_array_type , double , test_array_type , double , range_algebra , default_operations , always_resizer > euler_always_type;
|
|
|
|
typedef runge_kutta4_classic< test_array_type , double , test_array_type , double , range_algebra , default_operations , never_resizer > rk4_manual_type;
|
|
typedef runge_kutta4_classic< test_array_type , double , test_array_type , double , range_algebra , default_operations , initially_resizer > rk4_initially_type;
|
|
typedef runge_kutta4_classic< test_array_type , double , test_array_type , double , range_algebra , default_operations , always_resizer > rk4_always_type;
|
|
|
|
|
|
typedef runge_kutta4< test_array_type , double , test_array_type , double , range_algebra , default_operations , never_resizer > rk4_gen_manual_type;
|
|
typedef runge_kutta4< test_array_type , double , test_array_type , double , range_algebra , default_operations , initially_resizer > rk4_gen_initially_type;
|
|
typedef runge_kutta4< test_array_type , double , test_array_type , double , range_algebra , default_operations , always_resizer > rk4_gen_always_type;
|
|
|
|
|
|
typedef mpl::vector<
|
|
mpl::vector< euler_manual_type , mpl::int_<1> , mpl::int_<0> > ,
|
|
mpl::vector< euler_initially_type , mpl::int_<1> , mpl::int_<1> > ,
|
|
mpl::vector< euler_always_type , mpl::int_<1> , mpl::int_<3> > ,
|
|
mpl::vector< rk4_manual_type , mpl::int_<5> , mpl::int_<0> > ,
|
|
mpl::vector< rk4_initially_type , mpl::int_<5> , mpl::int_<1> > ,
|
|
mpl::vector< rk4_always_type , mpl::int_<5> , mpl::int_<3> > ,
|
|
mpl::vector< rk4_gen_manual_type , mpl::int_<5> , mpl::int_<0> > ,
|
|
mpl::vector< rk4_gen_initially_type , mpl::int_<5> , mpl::int_<1> > ,
|
|
mpl::vector< rk4_gen_always_type , mpl::int_<5> , mpl::int_<3> >
|
|
>::type resize_check_types;
|
|
|
|
|
|
BOOST_AUTO_TEST_CASE_TEMPLATE( test_resize , T, resize_check_types )
|
|
{
|
|
typedef typename mpl::at< T , mpl::int_< 0 > >::type stepper_type;
|
|
const size_t resize_calls = mpl::at< T , mpl::int_< 1 > >::type::value;
|
|
const size_t multiplicity = mpl::at< T , mpl::int_< 2 > >::type::value;
|
|
adjust_size_count = 0;
|
|
|
|
stepper_type stepper;
|
|
test_array_type x;
|
|
stepper.do_step( constant_system , x , 0.0 , 0.1 );
|
|
stepper.do_step( constant_system , x , 0.0 , 0.1 );
|
|
stepper.do_step( constant_system , x , 0.0 , 0.1 );
|
|
|
|
BOOST_TEST_MESSAGE( "adjust_size_count : " << adjust_size_count );
|
|
BOOST_CHECK_MESSAGE( adjust_size_count == resize_calls * multiplicity , "adjust_size_count : " << adjust_size_count << " expected: " << resize_calls * multiplicity );
|
|
}
|
|
|
|
|
|
BOOST_AUTO_TEST_SUITE_END()
|