| |
| /* |
| [auto_generated] |
| boost/numeric/odeint/iterator/adaptive_iterator.hpp |
| |
| [begin_description] |
| Iterator for iterating throught the solution of an ODE with adaptive step size. The dereferenced types containes also the time. |
| [end_description] |
| |
| Copyright 2012-2013 Karsten Ahnert |
| Copyright 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) |
| */ |
| |
| |
| #ifndef BOOST_NUMERIC_ODEINT_ITERATOR_ADAPTIVE_ITERATOR_HPP_INCLUDED |
| #define BOOST_NUMERIC_ODEINT_ITERATOR_ADAPTIVE_ITERATOR_HPP_INCLUDED |
| |
| #include <boost/numeric/odeint/util/stepper_traits.hpp> |
| #include <boost/numeric/odeint/util/unit_helper.hpp> |
| #include <boost/numeric/odeint/stepper/stepper_categories.hpp> |
| #include <boost/numeric/odeint/stepper/controlled_step_result.hpp> |
| #include <boost/numeric/odeint/iterator/impl/adaptive_iterator_impl.hpp> |
| |
| namespace boost { |
| namespace numeric { |
| namespace odeint { |
| |
| /* use the adaptive_iterator_impl with the right tags */ |
| template< class Stepper , class System , class State |
| #ifndef DOXYGEN_SKIP |
| , class StepperTag = typename base_tag< typename traits::stepper_category< Stepper >::type >::type |
| #endif |
| > |
| class adaptive_iterator : public adaptive_iterator_impl< |
| adaptive_iterator< Stepper , System , State , StepperTag > , |
| Stepper , System , State , detail::ode_state_iterator_tag , StepperTag |
| > |
| { |
| typedef typename traits::time_type< Stepper >::type time_type; |
| typedef adaptive_iterator< Stepper , System , State , StepperTag > iterator_type; |
| |
| public: |
| adaptive_iterator( Stepper stepper , System sys , State &s , time_type t_start , time_type t_end , time_type dt ) |
| : adaptive_iterator_impl< iterator_type , Stepper , System , State , detail::ode_state_iterator_tag , StepperTag >( stepper , sys , s , t_start , t_end , dt ) |
| {} |
| |
| adaptive_iterator( Stepper stepper , System sys , State &s ) |
| : adaptive_iterator_impl< iterator_type , Stepper , System , State , detail::ode_state_iterator_tag , StepperTag >( stepper , sys , s ) |
| {} |
| }; |
| |
| |
| |
| |
| template< class Stepper , class System , class State > |
| adaptive_iterator< Stepper , System , State > make_adaptive_iterator_begin( |
| Stepper stepper , |
| System system , |
| State &x , |
| typename traits::time_type< Stepper >::type t_start , |
| typename traits::time_type< Stepper >::type t_end , |
| typename traits::time_type< Stepper >::type dt ) |
| { |
| return adaptive_iterator< Stepper , System , State >( stepper , system , x , t_start , t_end , dt ); |
| } |
| |
| |
| template< class Stepper , class System , class State > |
| adaptive_iterator< Stepper , System , State > make_adaptive_iterator_end( |
| Stepper stepper , |
| System system , |
| State &x ) |
| { |
| return adaptive_iterator< Stepper , System , State >( stepper , system , x ); |
| } |
| |
| |
| template< class Stepper , class System , class State > |
| std::pair< adaptive_iterator< Stepper , System , State > , adaptive_iterator< Stepper , System , State > > |
| make_adaptive_range( |
| Stepper stepper , |
| System system , |
| State &x , |
| typename traits::time_type< Stepper >::type t_start , |
| typename traits::time_type< Stepper >::type t_end , |
| typename traits::time_type< Stepper >::type dt ) |
| { |
| return std::make_pair( |
| adaptive_iterator< Stepper , System , State >( stepper , system , x , t_start , t_end , dt ) , |
| adaptive_iterator< Stepper , System , State >( stepper , system , x ) |
| ); |
| } |
| |
| /** |
| * \class adaptive_iterator |
| * |
| * \brief ODE Iterator with adaptive step size. The value type of this iterator is the state type of the stepper. |
| * |
| * Implements an iterator representing the solution of an ODE from t_start |
| * to t_end evaluated at steps with an adaptive step size dt. |
| * After each iteration the iterator dereferences to the state x at the next |
| * time t+dt where dt is controlled by the stepper. |
| * This iterator can be used with ControlledSteppers and |
| * DenseOutputSteppers and it always makes use of the all the given steppers |
| * capabilities. A for_each over such an iterator range behaves similar to |
| * the integrate_adaptive routine. |
| * |
| * adaptive_iterator is a model of single-pass iterator. |
| * |
| * The value type of this iterator is the state type of the stepper. Hence one can only access the state and not the current time. |
| * |
| * \tparam Stepper The stepper type which should be used during the iteration. |
| * \tparam System The type of the system function (ODE) which should be solved. |
| * \tparam State The state type of the ODE. |
| */ |
| |
| |
| |
| |
| |
| |
| /** |
| * \fn make_adaptive_iterator_begin( Stepper stepper , System system , State &x , |
| typename traits::time_type< Stepper >::type t_start , |
| typename traits::time_type< Stepper >::type t_end , |
| typename traits::time_type< Stepper >::type dt ) |
| * |
| * \brief Factory function for adaptive_iterator. Constructs a begin iterator. |
| * |
| * \param stepper The stepper to use during the iteration. |
| * \param system The system function (ODE) to solve. |
| * \param x The initial state. |
| * \param t_start The initial time. |
| * \param t_end The end time, at which the iteration should stop. |
| * \param dt The initial time step. |
| * \returns The adaptive iterator. |
| */ |
| |
| |
| /** |
| * \fn make_adaptive_iterator_end( Stepper stepper , System system , State &x ) |
| * \brief Factory function for adaptive_iterator. Constructs a end iterator. |
| * |
| * \param stepper The stepper to use during the iteration. |
| * \param system The system function (ODE) to solve. |
| * \param x The initial state. |
| * \returns The adaptive iterator. |
| */ |
| |
| |
| /** |
| * \fn make_adaptive_range( Stepper stepper , System system , State &x , |
| typename traits::time_type< Stepper >::type t_start , |
| typename traits::time_type< Stepper >::type t_end , |
| typename traits::time_type< Stepper >::type dt ) |
| * |
| * \brief Factory function to construct a single pass range of adaptive iterators. A range is here a pair of adaptive_iterator. |
| * |
| * \param stepper The stepper to use during the iteration. |
| * \param system The system function (ODE) to solve. |
| * \param x The initial state. |
| * \param t_start The initial time. |
| * \param t_end The end time, at which the iteration should stop. |
| * \param dt The initial time step. |
| * \returns The adaptive range. |
| */ |
| |
| |
| |
| |
| |
| |
| } // namespace odeint |
| } // namespace numeric |
| } // namespace boost |
| |
| |
| |
| #endif // BOOST_NUMERIC_ODEINT_ITERATOR_ADAPTIVE_ITERATOR_HPP_INCLUDED |