3 boost/numeric/odeint/stepper/rosenbrock4_dense_output.hpp
6 Dense output for Rosenbrock 4.
9 Copyright 2011-2012 Karsten Ahnert
10 Copyright 2011-2015 Mario Mulansky
11 Copyright 2012 Christoph Koke
13 Distributed under the Boost Software License, Version 1.0.
14 (See accompanying file LICENSE_1_0.txt or
15 copy at http://www.boost.org/LICENSE_1_0.txt)
19 #ifndef BOOST_NUMERIC_ODEINT_STEPPER_ROSENBROCK4_DENSE_OUTPUT_HPP_INCLUDED
20 #define BOOST_NUMERIC_ODEINT_STEPPER_ROSENBROCK4_DENSE_OUTPUT_HPP_INCLUDED
25 #include <boost/numeric/odeint/util/bind.hpp>
27 #include <boost/numeric/odeint/stepper/rosenbrock4_controller.hpp>
28 #include <boost/numeric/odeint/util/is_resizeable.hpp>
30 #include <boost/numeric/odeint/integrate/max_step_checker.hpp>
37 template< class ControlledStepper
>
38 class rosenbrock4_dense_output
43 typedef ControlledStepper controlled_stepper_type
;
44 typedef typename
controlled_stepper_type::stepper_type stepper_type
;
45 typedef typename
stepper_type::value_type value_type
;
46 typedef typename
stepper_type::state_type state_type
;
47 typedef typename
stepper_type::wrapped_state_type wrapped_state_type
;
48 typedef typename
stepper_type::time_type time_type
;
49 typedef typename
stepper_type::deriv_type deriv_type
;
50 typedef typename
stepper_type::wrapped_deriv_type wrapped_deriv_type
;
51 typedef typename
stepper_type::resizer_type resizer_type
;
52 typedef dense_output_stepper_tag stepper_category
;
54 typedef rosenbrock4_dense_output
< ControlledStepper
> dense_output_stepper_type
;
56 rosenbrock4_dense_output( const controlled_stepper_type
&stepper
= controlled_stepper_type() )
57 : m_stepper( stepper
) ,
59 m_current_state_x1( true ) ,
60 m_t() , m_t_old() , m_dt()
66 template< class StateType
>
67 void initialize( const StateType
&x0
, time_type t0
, time_type dt0
)
69 m_resizer
.adjust_size( x0
, detail::bind( &dense_output_stepper_type::template resize_impl
< StateType
> , detail::ref( *this ) , detail::_1
) );
70 get_current_state() = x0
;
75 template< class System
>
76 std::pair
< time_type
, time_type
> do_step( System system
)
78 failed_step_checker fail_checker
; // to throw a runtime_error if step size adjustment fails
79 controlled_step_result res
= fail
;
83 res
= m_stepper
.try_step( system
, get_current_state() , m_t
, get_old_state() , m_dt
);
84 fail_checker(); // check for overflow of failed steps
87 m_stepper
.stepper().prepare_dense_output();
88 this->toggle_current_state();
89 return std::make_pair( m_t_old
, m_t
);
94 * The two overloads are needed in order to solve the forwarding problem.
96 template< class StateOut
>
97 void calc_state( time_type t
, StateOut
&x
)
99 m_stepper
.stepper().calc_state( t
, x
, get_old_state() , m_t_old
, get_current_state() , m_t
);
102 template< class StateOut
>
103 void calc_state( time_type t
, const StateOut
&x
)
105 m_stepper
.stepper().calc_state( t
, x
, get_old_state() , m_t_old
, get_current_state() , m_t
);
109 template< class StateType
>
110 void adjust_size( const StateType
&x
)
112 m_stepper
.adjust_size( x
);
119 const state_type
& current_state( void ) const
121 return get_current_state();
124 time_type
current_time( void ) const
129 const state_type
& previous_state( void ) const
131 return get_old_state();
134 time_type
previous_time( void ) const
139 time_type
current_time_step( void ) const
149 state_type
& get_current_state( void )
151 return m_current_state_x1
? m_x1
.m_v
: m_x2
.m_v
;
154 const state_type
& get_current_state( void ) const
156 return m_current_state_x1
? m_x1
.m_v
: m_x2
.m_v
;
159 state_type
& get_old_state( void )
161 return m_current_state_x1
? m_x2
.m_v
: m_x1
.m_v
;
164 const state_type
& get_old_state( void ) const
166 return m_current_state_x1
? m_x2
.m_v
: m_x1
.m_v
;
169 void toggle_current_state( void )
171 m_current_state_x1
= ! m_current_state_x1
;
175 template< class StateIn
>
176 bool resize_impl( const StateIn
&x
)
178 bool resized
= false;
179 resized
|= adjust_size_by_resizeability( m_x1
, x
, typename is_resizeable
<state_type
>::type() );
180 resized
|= adjust_size_by_resizeability( m_x2
, x
, typename is_resizeable
<state_type
>::type() );
185 controlled_stepper_type m_stepper
;
186 resizer_type m_resizer
;
187 wrapped_state_type m_x1
, m_x2
;
188 bool m_current_state_x1
;
189 time_type m_t
, m_t_old
, m_dt
;
194 } // namespace odeint
195 } // namespace numeric
199 #endif // BOOST_NUMERIC_ODEINT_STEPPER_ROSENBROCK4_DENSE_OUTPUT_HPP_INCLUDED