ceph/src/boost/libs/numeric/odeint/test/regression/regression_149.cpp

   1 /*
   2
   3  [begin_description]
   4  Test case for issue 149:
   5  Error C2582 with msvc-10 when using iterator-based integration
   6  [end_description]
   7
   8  Copyright 2011-2015 Karsten Ahnert
   9  Copyright 2011-2015 Mario Mulansky
  10
  11  Distributed under the Boost Software License, Version 1.0.
  12  (See accompanying file LICENSE_1_0.txt or
  13  copy at http://www.boost.org/LICENSE_1_0.txt)
  14  */
  15
  16
  17 // disable checked iterator warning for msvc
  18
  19 #include <boost/config.hpp>
  20 #ifdef BOOST_MSVC
  21     #pragma warning(disable:4996)
  22 #endif
  23
  24 #define BOOST_TEST_MODULE odeint_regression_147
  25
  26 #include <utility>
  27 #include <iostream>
  28
  29 #include <boost/test/unit_test.hpp>
  30
  31 #include <boost/mpl/vector.hpp>
  32 #include <boost/range/algorithm/find_if.hpp>
  33
  34 #include <boost/numeric/odeint.hpp>
  35
  36 using namespace boost::unit_test;
  37 using namespace boost::numeric::odeint;
  38 namespace mpl = boost::mpl;
  39
  40 typedef std::vector<double> state_type;
  41
  42 void rhs( const state_type &x , state_type &dxdt , const double t )
  43 {
  44 }
  45
  46
  47 template<class Stepper>
  48 struct perform_test
  49 {
  50     void operator()( void )
  51     {
  52         bulirsch_stoer< state_type > stepper( 1e-9, 0.0, 0.0, 0.0 );
  53         state_type x( 3, 10.0 );
  54
  55         print( boost::find_if(
  56             make_adaptive_time_range( stepper, rhs, x, 0.0, 1.0, 0.01 ),
  57             pred ) );
  58
  59     }
  60
  61     static bool pred( const std::pair< const state_type &, double > &x )
  62     {
  63         return ( x.first[0] < 0.0 );
  64     }
  65
  66     template<class Iterator>
  67     void print( Iterator iter )
  68     {
  69         std::cout << iter->second << "\t" << iter->first[0] << "\t"
  70                   << iter->first[1] << "\t" << iter->first[2] << "\n";
  71     }
  72 };
  73
  74 typedef mpl::vector<
  75     euler< state_type > ,
  76     runge_kutta4< state_type > ,
  77     runge_kutta_cash_karp54< state_type > ,
  78     runge_kutta_dopri5< state_type > ,
  79     runge_kutta_fehlberg78< state_type > ,
  80     bulirsch_stoer< state_type >
  81     > steppers;
  82
  83
  84 BOOST_AUTO_TEST_SUITE( regression_147_test )
  85
  86 BOOST_AUTO_TEST_CASE_TEMPLATE( regression_147_test , Stepper,
  87                                steppers )
  88 {
  89     perform_test< Stepper > tester;
  90     tester();
  91 }
  92
  93 BOOST_AUTO_TEST_SUITE_END()