Chris@16
|
1 /*
|
Chris@16
|
2 [auto_generated]
|
Chris@16
|
3 boost/numeric/odeint/integrate/integrate_adaptive.hpp
|
Chris@16
|
4
|
Chris@16
|
5 [begin_description]
|
Chris@16
|
6 Adaptive integration of ODEs.
|
Chris@16
|
7 [end_description]
|
Chris@16
|
8
|
Chris@101
|
9 Copyright 2011-2013 Karsten Ahnert
|
Chris@101
|
10 Copyright 2011-2012 Mario Mulansky
|
Chris@16
|
11
|
Chris@16
|
12 Distributed under the Boost Software License, Version 1.0.
|
Chris@16
|
13 (See accompanying file LICENSE_1_0.txt or
|
Chris@16
|
14 copy at http://www.boost.org/LICENSE_1_0.txt)
|
Chris@16
|
15 */
|
Chris@16
|
16
|
Chris@16
|
17
|
Chris@16
|
18 #ifndef BOOST_NUMERIC_ODEINT_INTEGRATE_INTEGRATE_ADAPTIVE_HPP_INCLUDED
|
Chris@16
|
19 #define BOOST_NUMERIC_ODEINT_INTEGRATE_INTEGRATE_ADAPTIVE_HPP_INCLUDED
|
Chris@16
|
20
|
Chris@16
|
21 #include <boost/type_traits/is_same.hpp>
|
Chris@16
|
22
|
Chris@16
|
23 #include <boost/numeric/odeint/stepper/stepper_categories.hpp>
|
Chris@16
|
24 #include <boost/numeric/odeint/integrate/null_observer.hpp>
|
Chris@16
|
25 #include <boost/numeric/odeint/integrate/detail/integrate_adaptive.hpp>
|
Chris@16
|
26
|
Chris@16
|
27 namespace boost {
|
Chris@16
|
28 namespace numeric {
|
Chris@16
|
29 namespace odeint {
|
Chris@16
|
30
|
Chris@16
|
31
|
Chris@16
|
32 /*
|
Chris@16
|
33 * the two overloads are needed in order to solve the forwarding problem
|
Chris@16
|
34 */
|
Chris@16
|
35 template< class Stepper , class System , class State , class Time , class Observer >
|
Chris@16
|
36 size_t integrate_adaptive(
|
Chris@16
|
37 Stepper stepper , System system , State &start_state ,
|
Chris@16
|
38 Time start_time , Time end_time , Time dt ,
|
Chris@16
|
39 Observer observer )
|
Chris@16
|
40 {
|
Chris@101
|
41 typedef typename odeint::unwrap_reference< Stepper >::type::stepper_category stepper_category;
|
Chris@16
|
42 return detail::integrate_adaptive(
|
Chris@16
|
43 stepper , system , start_state ,
|
Chris@16
|
44 start_time , end_time , dt ,
|
Chris@101
|
45 observer , stepper_category() );
|
Chris@16
|
46
|
Chris@16
|
47 /*
|
Chris@16
|
48 * Suggestion for a new extendable version:
|
Chris@16
|
49 *
|
Chris@16
|
50 * integrator_adaptive< Stepper , System, State , Time , Observer , typename Stepper::stepper_category > integrator;
|
Chris@16
|
51 * return integrator.run( stepper , system , start_state , start_time , end_time , dt , observer );
|
Chris@16
|
52 */
|
Chris@16
|
53 }
|
Chris@16
|
54
|
Chris@16
|
55 /**
|
Chris@16
|
56 * \brief Second version to solve the forwarding problem,
|
Chris@16
|
57 * can be called with Boost.Range as start_state.
|
Chris@16
|
58 */
|
Chris@16
|
59 template< class Stepper , class System , class State , class Time , class Observer >
|
Chris@16
|
60 size_t integrate_adaptive(
|
Chris@16
|
61 Stepper stepper , System system , const State &start_state ,
|
Chris@16
|
62 Time start_time , Time end_time , Time dt ,
|
Chris@16
|
63 Observer observer )
|
Chris@16
|
64 {
|
Chris@101
|
65 typedef typename odeint::unwrap_reference< Stepper >::type::stepper_category stepper_category;
|
Chris@16
|
66 return detail::integrate_adaptive(
|
Chris@16
|
67 stepper , system , start_state ,
|
Chris@16
|
68 start_time , end_time , dt ,
|
Chris@101
|
69 observer , stepper_category() );
|
Chris@16
|
70 }
|
Chris@16
|
71
|
Chris@16
|
72
|
Chris@16
|
73
|
Chris@16
|
74
|
Chris@16
|
75 /**
|
Chris@16
|
76 * \brief integrate_adaptive without an observer.
|
Chris@16
|
77 */
|
Chris@16
|
78 template< class Stepper , class System , class State , class Time >
|
Chris@16
|
79 size_t integrate_adaptive(
|
Chris@16
|
80 Stepper stepper , System system , State &start_state ,
|
Chris@16
|
81 Time start_time , Time end_time , Time dt )
|
Chris@16
|
82 {
|
Chris@16
|
83 return integrate_adaptive( stepper , system , start_state , start_time , end_time , dt , null_observer() );
|
Chris@16
|
84 }
|
Chris@16
|
85
|
Chris@16
|
86 /**
|
Chris@16
|
87 * \brief Second version to solve the forwarding problem,
|
Chris@16
|
88 * can be called with Boost.Range as start_state.
|
Chris@16
|
89 */
|
Chris@16
|
90 template< class Stepper , class System , class State , class Time >
|
Chris@16
|
91 size_t integrate_adaptive(
|
Chris@16
|
92 Stepper stepper , System system , const State &start_state ,
|
Chris@16
|
93 Time start_time , Time end_time , Time dt )
|
Chris@16
|
94 {
|
Chris@16
|
95 return integrate_adaptive( stepper , system , start_state , start_time , end_time , dt , null_observer() );
|
Chris@16
|
96 }
|
Chris@16
|
97
|
Chris@16
|
98
|
Chris@16
|
99 /************* DOXYGEN ************/
|
Chris@16
|
100
|
Chris@16
|
101 /**
|
Chris@16
|
102 * \fn integrate_adaptive( Stepper stepper , System system , State &start_state , Time start_time , Time end_time , Time dt , Observer observer )
|
Chris@16
|
103 * \brief Integrates the ODE with adaptive step size.
|
Chris@16
|
104 *
|
Chris@16
|
105 * This function integrates the ODE given by system with the given stepper.
|
Chris@16
|
106 * The observer is called after each step. If the stepper has no error
|
Chris@16
|
107 * control, the step size remains constant and the observer is called at
|
Chris@16
|
108 * equidistant time points t0+n*dt. If the stepper is a ControlledStepper,
|
Chris@16
|
109 * the step size is adjusted and the observer is called in non-equidistant
|
Chris@16
|
110 * intervals.
|
Chris@16
|
111 *
|
Chris@16
|
112 * \param stepper The stepper to be used for numerical integration.
|
Chris@16
|
113 * \param system Function/Functor defining the rhs of the ODE.
|
Chris@16
|
114 * \param start_state The initial condition x0.
|
Chris@16
|
115 * \param start_time The initial time t0.
|
Chris@16
|
116 * \param end_time The final integration time tend.
|
Chris@16
|
117 * \param dt The time step between observer calls, _not_ necessarily the
|
Chris@16
|
118 * time step of the integration.
|
Chris@16
|
119 * \param observer Function/Functor called at equidistant time intervals.
|
Chris@16
|
120 * \return The number of steps performed.
|
Chris@16
|
121 */
|
Chris@16
|
122
|
Chris@16
|
123 } // namespace odeint
|
Chris@16
|
124 } // namespace numeric
|
Chris@16
|
125 } // namespace boost
|
Chris@16
|
126
|
Chris@16
|
127
|
Chris@16
|
128
|
Chris@16
|
129 #endif // BOOST_NUMERIC_ODEINT_INTEGRATE_INTEGRATE_ADAPTIVE_HPP_INCLUDED
|