File indexing completed on 2025-01-18 09:42:55
0001
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014
0015
0016
0017
0018 #ifndef BOOST_NUMERIC_ODEINT_STEPPER_GENERATION_GENERATION_DENSE_OUTPUT_RUNGE_KUTTA_HPP_INCLUDED
0019 #define BOOST_NUMERIC_ODEINT_STEPPER_GENERATION_GENERATION_DENSE_OUTPUT_RUNGE_KUTTA_HPP_INCLUDED
0020
0021 #include <boost/numeric/odeint/stepper/controlled_runge_kutta.hpp>
0022 #include <boost/numeric/odeint/stepper/dense_output_runge_kutta.hpp>
0023 #include <boost/numeric/odeint/stepper/generation/make_dense_output.hpp>
0024
0025 namespace boost {
0026 namespace numeric {
0027 namespace odeint {
0028
0029
0030 template< class Stepper >
0031 struct dense_output_factory< Stepper , dense_output_runge_kutta< controlled_runge_kutta< Stepper > > >
0032 {
0033 typedef Stepper stepper_type;
0034 typedef controlled_runge_kutta< stepper_type > controller_type;
0035 typedef typename controller_type::error_checker_type error_checker_type;
0036 typedef typename controller_type::step_adjuster_type step_adjuster_type;
0037 typedef typename stepper_type::value_type value_type;
0038 typedef typename stepper_type::time_type time_type;
0039 typedef dense_output_runge_kutta< controller_type > dense_output_type;
0040
0041 dense_output_type operator()( value_type abs_error , value_type rel_error , const stepper_type &stepper )
0042 {
0043 return dense_output_type( controller_type( error_checker_type( abs_error , rel_error ) ,
0044 step_adjuster_type() , stepper ) );
0045 }
0046
0047 dense_output_type operator()( value_type abs_error , value_type rel_error ,
0048 time_type max_dt , const stepper_type &stepper )
0049 {
0050 return dense_output_type(
0051 controller_type( error_checker_type( abs_error , rel_error) ,
0052 step_adjuster_type( max_dt ) , stepper ) );
0053 }
0054 };
0055
0056
0057
0058
0059
0060 }
0061 }
0062 }
0063
0064
0065 #endif