summaryrefslogtreecommitdiff
path: root/inference-engine/thirdparty/clDNN/common/boost/1.64.0/include/boost-1_64/boost/numeric/odeint/stepper/implicit_euler.hpp
diff options
context:
space:
mode:
Diffstat (limited to 'inference-engine/thirdparty/clDNN/common/boost/1.64.0/include/boost-1_64/boost/numeric/odeint/stepper/implicit_euler.hpp')
-rw-r--r--inference-engine/thirdparty/clDNN/common/boost/1.64.0/include/boost-1_64/boost/numeric/odeint/stepper/implicit_euler.hpp170
1 files changed, 0 insertions, 170 deletions
diff --git a/inference-engine/thirdparty/clDNN/common/boost/1.64.0/include/boost-1_64/boost/numeric/odeint/stepper/implicit_euler.hpp b/inference-engine/thirdparty/clDNN/common/boost/1.64.0/include/boost-1_64/boost/numeric/odeint/stepper/implicit_euler.hpp
deleted file mode 100644
index e1c64164b..000000000
--- a/inference-engine/thirdparty/clDNN/common/boost/1.64.0/include/boost-1_64/boost/numeric/odeint/stepper/implicit_euler.hpp
+++ /dev/null
@@ -1,170 +0,0 @@
-/*
- [auto_generated]
- boost/numeric/odeint/stepper/implicit_euler.hpp
-
- [begin_description]
- Impementation of the implicit Euler method. Works with ublas::vector as state type.
- [end_description]
-
- Copyright 2010-2012 Mario Mulansky
- Copyright 2010-2012 Karsten Ahnert
- Copyright 2012 Christoph Koke
-
- 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_STEPPER_IMPLICIT_EULER_HPP_INCLUDED
-#define BOOST_NUMERIC_ODEINT_STEPPER_IMPLICIT_EULER_HPP_INCLUDED
-
-
-#include <utility>
-
-#include <boost/numeric/odeint/util/bind.hpp>
-#include <boost/numeric/odeint/util/unwrap_reference.hpp>
-#include <boost/numeric/odeint/stepper/stepper_categories.hpp>
-
-#include <boost/numeric/odeint/util/ublas_wrapper.hpp>
-#include <boost/numeric/odeint/util/is_resizeable.hpp>
-#include <boost/numeric/odeint/util/resizer.hpp>
-
-#include <boost/numeric/ublas/vector.hpp>
-#include <boost/numeric/ublas/matrix.hpp>
-#include <boost/numeric/ublas/lu.hpp>
-
-namespace boost {
-namespace numeric {
-namespace odeint {
-
-
-
-
-
-
-
-
-template< class ValueType , class Resizer = initially_resizer >
-class implicit_euler
-{
-
-public:
-
- typedef ValueType value_type;
- typedef value_type time_type;
- typedef boost::numeric::ublas::vector< value_type > state_type;
- typedef state_wrapper< state_type > wrapped_state_type;
- typedef state_type deriv_type;
- typedef state_wrapper< deriv_type > wrapped_deriv_type;
- typedef boost::numeric::ublas::matrix< value_type > matrix_type;
- typedef state_wrapper< matrix_type > wrapped_matrix_type;
- typedef boost::numeric::ublas::permutation_matrix< size_t > pmatrix_type;
- typedef state_wrapper< pmatrix_type > wrapped_pmatrix_type;
- typedef Resizer resizer_type;
- typedef stepper_tag stepper_category;
- typedef implicit_euler< ValueType , Resizer > stepper_type;
-
- implicit_euler( value_type epsilon = 1E-6 )
- : m_epsilon( epsilon )
- { }
-
-
- template< class System >
- void do_step( System system , state_type &x , time_type t , time_type dt )
- {
- typedef typename odeint::unwrap_reference< System >::type system_type;
- typedef typename odeint::unwrap_reference< typename system_type::first_type >::type deriv_func_type;
- typedef typename odeint::unwrap_reference< typename system_type::second_type >::type jacobi_func_type;
- system_type &sys = system;
- deriv_func_type &deriv_func = sys.first;
- jacobi_func_type &jacobi_func = sys.second;
-
- m_resizer.adjust_size( x , detail::bind( &stepper_type::template resize_impl<state_type> , detail::ref( *this ) , detail::_1 ) );
-
- for( size_t i=0 ; i<x.size() ; ++i )
- m_pm.m_v[i] = i;
-
- t += dt;
-
- // apply first Newton step
- deriv_func( x , m_dxdt.m_v , t );
-
- m_b.m_v = dt * m_dxdt.m_v;
-
- jacobi_func( x , m_jacobi.m_v , t );
- m_jacobi.m_v *= dt;
- m_jacobi.m_v -= boost::numeric::ublas::identity_matrix< value_type >( x.size() );
-
- solve( m_b.m_v , m_jacobi.m_v );
-
- m_x.m_v = x - m_b.m_v;
-
- // iterate Newton until some precision is reached
- // ToDo: maybe we should apply only one Newton step -> linear implicit one-step scheme
- while( boost::numeric::ublas::norm_2( m_b.m_v ) > m_epsilon )
- {
- deriv_func( m_x.m_v , m_dxdt.m_v , t );
- m_b.m_v = x - m_x.m_v + dt*m_dxdt.m_v;
-
- // simplified version, only the first Jacobian is used
- // jacobi( m_x , m_jacobi , t );
- // m_jacobi *= dt;
- // m_jacobi -= boost::numeric::ublas::identity_matrix< value_type >( x.size() );
-
- solve( m_b.m_v , m_jacobi.m_v );
-
- m_x.m_v -= m_b.m_v;
- }
- x = m_x.m_v;
- }
-
- template< class StateType >
- void adjust_size( const StateType &x )
- {
- resize_impl( x );
- }
-
-
-private:
-
- template< class StateIn >
- bool resize_impl( const StateIn &x )
- {
- bool resized = false;
- resized |= adjust_size_by_resizeability( m_dxdt , x , typename is_resizeable<deriv_type>::type() );
- resized |= adjust_size_by_resizeability( m_x , x , typename is_resizeable<state_type>::type() );
- resized |= adjust_size_by_resizeability( m_b , x , typename is_resizeable<deriv_type>::type() );
- resized |= adjust_size_by_resizeability( m_jacobi , x , typename is_resizeable<matrix_type>::type() );
- resized |= adjust_size_by_resizeability( m_pm , x , typename is_resizeable<pmatrix_type>::type() );
- return resized;
- }
-
-
- void solve( state_type &x , matrix_type &m )
- {
- int res = boost::numeric::ublas::lu_factorize( m , m_pm.m_v );
- if( res != 0 ) std::exit(0);
- boost::numeric::ublas::lu_substitute( m , m_pm.m_v , x );
- }
-
-private:
-
- value_type m_epsilon;
- resizer_type m_resizer;
- wrapped_deriv_type m_dxdt;
- wrapped_state_type m_x;
- wrapped_deriv_type m_b;
- wrapped_matrix_type m_jacobi;
- wrapped_pmatrix_type m_pm;
-
-
-};
-
-
-} // odeint
-} // numeric
-} // boost
-
-
-#endif // BOOST_NUMERIC_ODEINT_STEPPER_IMPLICIT_EULER_HPP_INCLUDED