| /third_party/boost/libs/numeric/odeint/test/ |
| D | dummy_odes.hpp | 34 void operator()( const State &x , Deriv &dxdt , const Time t ) const in operator ()()
36 dxdt[0] = 1.0; in operator ()()
43 void operator()( const State &x , Deriv &dxdt , const Time t ) const in operator ()()
45 dxdt = 1.0; in operator ()()
52 void operator()( const State &x , Deriv &dxdt , const Time t ) const in operator ()()
54 boost::fusion::at_c< 0 >( dxdt ) = boost::fusion::at_c< 0 >( x ) / Time( 1.0 ); in operator ()()
61 void operator()( const State& x , Deriv& dxdt , const Time& t ) const in operator ()()
66 dxdt[0] = sigma * ( x[1] - x[0] ); in operator ()()
67 dxdt[1] = R * x[0] - x[1] - x[0] * x[2]; in operator ()()
68 dxdt[2] = -b * x[2] + x[0] * x[1]; in operator ()()
[all …]
|
| D | stepper_with_ranges.cpp | 67 typename boost::range_iterator< Deriv >::type dxdt = boost::begin( dxdt_ ); in operator ()() local
69 dxdt[0] = x[0]; in operator ()()
70 dxdt[1] = 2.0; in operator ()()
71 dxdt[2] = 3.0; in operator ()()
78 typename boost::range_iterator< Deriv >::type dxdt = boost::begin( dxdt_ ); in operator ()() local
80 dxdt[0] = x[0]; in operator ()()
81 dxdt[1] = 2.0; in operator ()()
82 dxdt[2] = 3.0; in operator ()()
94 typename boost::range_iterator< Deriv >::type dxdt = boost::begin( dxdt_ ); in operator ()() local
96 dxdt[0] = 1.0; in operator ()()
[all …]
|
| D | stepper_with_units.cpp | 68 void oscillator( const state_type &x , deriv_type &dxdt , time_type t ) in oscillator() argument
71 fusion::at_c< 0 >( dxdt ) = fusion::at_c< 1 >( x ); in oscillator()
72 fusion::at_c< 1 >( dxdt ) = - omega * omega * fusion::at_c< 0 >( x ); in oscillator()
98 deriv_type dxdt; in check_stepper() local
99 oscillator( x , dxdt , t ); in check_stepper()
100 stepper.do_step( oscillator , x , dxdt , t , dt ); in check_stepper()
103 oscillator( x , dxdt , t ); in check_stepper()
104 stepper.do_step( oscillator , x , dxdt , t , x , dt ); in check_stepper()
130 deriv_type dxdt; in check_fsal_stepper() local
131 oscillator( x , dxdt , t ); in check_fsal_stepper()
[all …]
|
| D | bulirsch_stoer.cpp | 51 void operator()( const State &x , Deriv &dxdt , double t ) const in operator ()()
53 dxdt[0] = sigma * ( x[1] - x[0] ); in operator ()()
54 dxdt[1] = R * x[0] - x[1] - x[0] * x[2]; in operator ()()
55 dxdt[2] = -b * x[2] + x[0] * x[1]; in operator ()()
62 void operator()( const State &x , Deriv &dxdt , double t ) const in operator ()()
64 dxdt[0] = 1.0; in operator ()()
65 dxdt[1] = 1.0; in operator ()()
66 dxdt[2] = 1.0; in operator ()()
73 void operator()( const State &x , Deriv &dxdt , double t ) const in operator ()()
75 dxdt[0] = sin( x[0] ); in operator ()()
[all …]
|
| /third_party/boost/boost/numeric/odeint/stepper/detail/ |
| D | adams_moulton_call_algebra.hpp | 35 …void operator()( Algebra &algebra , const StateIn &in , StateOut &out , const DerivIn &dxdt , cons… in operator ()()
38 …algebra.for_each3( out , in , dxdt , typename Operations::template scale_sum2< value_type , Time >… in operator ()()
47 …void operator()( Algebra &algebra , const StateIn &in , StateOut &out , const DerivIn &dxdt , cons… in operator ()()
50 algebra.for_each4( out , in , dxdt , steps[0].m_v , in operator ()()
60 …void operator()( Algebra &algebra , const StateIn &in , StateOut &out , const DerivIn &dxdt , cons… in operator ()()
63 algebra.for_each5( out , in , dxdt , steps[0].m_v , steps[1].m_v , in operator ()()
73 …void operator()( Algebra &algebra , const StateIn &in , StateOut &out , const DerivIn &dxdt , cons… in operator ()()
76 algebra.for_each6( out , in , dxdt , steps[0].m_v , steps[1].m_v , steps[2].m_v , in operator ()()
87 …void operator()( Algebra &algebra , const StateIn &in , StateOut &out , const DerivIn &dxdt , cons… in operator ()()
90 … algebra.for_each7( out , in , dxdt , steps[0].m_v , steps[1].m_v , steps[2].m_v , steps[3].m_v , in operator ()()
[all …]
|
| /third_party/boost/libs/numeric/odeint/examples/mtl/ |
| D | implicit_euler_mtl.cpp | 40 void operator()( const vec_mtl4 &x , vec_mtl4 &dxdt , double t ) in operator ()()
44 dxdt[ 0 ] = -0.06*x[0]; in operator ()()
48 dxdt[ i ] = 4.2*x[i-1]-2.2*x[i]*x[i]; in operator ()()
76 void operator()( const vec_ublas &x , vec_ublas &dxdt , double t ) in operator ()()
80 dxdt[ 0 ] = -0.06*x[0]; in operator ()()
84 dxdt[ i ] = 4.2*x[i-1]-2.2*x[i]*x[i]; in operator ()()
111 void operator()( const vec_mtl4 &x , vec_mtl4 &dxdt , double t ) in operator ()()
118 dxdt[ i ] = -0.5*x[i]; in operator ()()
119 dxdt[i+1]= +25*x[i+1]*x[i+2]-740*x[i+3]*x[i+3]+4.2e-2*x[i]; in operator ()()
120 dxdt[i+2]= +25*x[i]*x[i]-740*x[i+3]*x[i+3]; in operator ()()
[all …]
|
| /third_party/boost/libs/numeric/odeint/examples/thrust/ |
| D | lorenz_parameters.cu | 70 void operator()( const State &x , Deriv &dxdt , value_type t ) const in operator ()()
78 boost::begin( dxdt ) , in operator ()()
79 boost::begin( dxdt ) + m_N , in operator ()()
80 boost::begin( dxdt ) + 2 * m_N ) ) , in operator ()()
86 boost::begin( dxdt ) + m_N , in operator ()()
87 boost::begin( dxdt ) + 2 * m_N , in operator ()()
88 boost::begin( dxdt ) + 3 * m_N ) ) , in operator ()()
124 void operator()( const State &x , Deriv &dxdt , value_type t ) const in operator ()()
137 boost::begin( dxdt ) , in operator ()()
138 boost::begin( dxdt ) + m_N , in operator ()()
[all …]
|
| /third_party/boost/libs/numeric/odeint/examples/ |
| D | harmonic_oscillator.cpp | 26 void harmonic_oscillator( const state_type &x , state_type &dxdt , const double /* t */ ) in harmonic_oscillator() argument
28 dxdt[0] = x[1]; in harmonic_oscillator()
29 dxdt[1] = -x[0] - gam*x[1]; in harmonic_oscillator()
46 void operator() ( const state_type &x , state_type &dxdt , const double /* t */ ) in operator ()() argument
48 dxdt[0] = x[1]; in operator ()()
49 dxdt[1] = -x[0] - m_gam*x[1]; in operator ()()
192 integrate_const( stepper , []( const state_type &x , state_type &dxdt , double t ) { in main() argument
193 dxdt[0] = x[1]; dxdt[1] = -x[0] - gam*x[1]; } in main()
|
| D | chaotic_system.cpp | 39 typename boost::range_iterator< Deriv >::type dxdt = boost::begin( dxdt_ ); in operator ()() local
41 dxdt[0] = sigma * ( x[1] - x[0] ); in operator ()()
42 dxdt[1] = R * x[0] - x[1] - x[0] * x[2]; in operator ()()
43 dxdt[2] = -b * x[2] + x[0] * x[1]; in operator ()()
58 void lorenz_with_lyap( const state_type &x , state_type &dxdt , double t ) in lorenz_with_lyap() argument
60 lorenz()( x , dxdt , t ); in lorenz_with_lyap()
65 double *dpert = dxdt.begin() + 3 + l * 3; in lorenz_with_lyap()
|
| D | bind_member_functions.cpp | 37 void operator()( const State &x , Deriv &dxdt , Time t ) in operator ()() argument
39 (m_obj.*m_mem)( x , dxdt , t ); in operator ()()
79 void ode( const state_type &x , state_type &dxdt , double t ) const in ode()
81 dxdt[0] = 10.0 * ( x[1] - x[0] ); in ode()
82 dxdt[1] = 28.0 * x[0] - x[1] - x[0] * x[2]; in ode()
83 dxdt[2] = -8.0 / 3.0 * x[2] + x[0] * x[1]; in ode()
|
| D | lorenz.cpp | 15 void lorenz( const state_type &x , state_type &dxdt , double t ) in lorenz() argument
17 dxdt[0] = sigma * ( x[1] - x[0] ); in lorenz()
18 dxdt[1] = R * x[0] - x[1] - x[0] * x[2]; in lorenz()
19 dxdt[2] = -b * x[2] + x[0] * x[1]; in lorenz()
|
| D | bind_member_functions_cpp11.cpp | 31 void ode( const state_type &x , state_type &dxdt , double t ) const in ode()
37 dxdt[0] = sigma * ( x[1] - x[0] ); in ode()
38 dxdt[1] = R * x[0] - x[1] - x[0] * x[2]; in ode()
39 dxdt[2] = -b * x[2] + x[0] * x[1]; in ode()
|
| /third_party/boost/boost/math/interpolators/ |
| D | vector_barycentric_rational.hpp | 39 void prime(Point& dxdt, Real t) const { in prime() argument
41 m_imp->eval_with_prime(x, dxdt, t); in prime()
50 void eval_with_prime(Point& x, Point& dxdt, Real t) const { in eval_with_prime() argument
51 m_imp->eval_with_prime(x, dxdt, t); in eval_with_prime()
57 Point dxdt; in eval_with_prime() local
58 m_imp->eval_with_prime(x, dxdt, t); in eval_with_prime()
59 return {x, dxdt}; in eval_with_prime()
|
| /third_party/boost/boost/numeric/odeint/stepper/ |
| D | extrapolation_stepper.hpp | 141 void do_step_impl( System system, const StateIn &in, const DerivIn &dxdt, in do_step_impl() argument
146 do_step_impl( system, in, dxdt, t, out, dt ); in do_step_impl()
157 void do_step_impl_io( System system, StateInOut &inout, const DerivIn &dxdt, in do_step_impl_io() argument
161 do_step_impl_io( system, inout, dxdt, t, dt ); in do_step_impl_io()
172 void do_step_impl( System system, const StateIn &in, const DerivIn &dxdt, in do_step_impl() argument
180 m_midpoint.do_step( system, in, dxdt, t, out, dt ); in do_step_impl()
184 m_midpoint.do_step( system, in, dxdt, t, m_table[k - 1].m_v, dt ); in do_step_impl()
190 void do_step_impl_io( System system, StateInOut &inout, const DerivIn &dxdt, in do_step_impl_io() argument
198 do_step_impl( system, inout, dxdt, t, m_xout.m_v, dt ); in do_step_impl_io()
203 void do_step_dxdt_impl( System system, StateInOut &x, const DerivIn &dxdt, in do_step_dxdt_impl() argument
[all …]
|
| D | runge_kutta_cash_karp54_classic.hpp | 93 …void do_step_impl( System system , const StateIn &in , const DerivIn &dxdt , time_type t , StateOu… in do_step_impl() argument
106 do_step_impl( system , in , dxdt , t , out , dt ); in do_step_impl()
109 … stepper_base_type::m_algebra.for_each6( xerr , dxdt , m_k3.m_v , m_k4.m_v , m_k5.m_v , m_k6.m_v , in do_step_impl()
117 …void do_step_impl( System system , const StateIn &in , const DerivIn &dxdt , time_type t , StateOu… in do_step_impl() argument
151 stepper_base_type::m_algebra.for_each3( m_x_tmp.m_v , in , dxdt , in do_step_impl()
156 stepper_base_type::m_algebra.for_each4( m_x_tmp.m_v , in , dxdt , m_k2.m_v , in do_step_impl()
161 stepper_base_type::m_algebra.for_each5( m_x_tmp.m_v , in , dxdt , m_k2.m_v , m_k3.m_v , in do_step_impl()
165 … stepper_base_type::m_algebra.for_each6( m_x_tmp.m_v , in , dxdt , m_k2.m_v , m_k3.m_v , m_k4.m_v , in do_step_impl()
169 …stepper_base_type::m_algebra.for_each7( m_x_tmp.m_v , in , dxdt , m_k2.m_v , m_k3.m_v , m_k4.m_v ,… in do_step_impl()
173 stepper_base_type::m_algebra.for_each6( out , in , dxdt , m_k3.m_v , m_k4.m_v , m_k6.m_v , in do_step_impl()
|
| /third_party/boost/boost/numeric/odeint/stepper/base/ |
| D | explicit_error_stepper_base.hpp | 148 do_step( System system , StateInOut &x , const DerivIn &dxdt , time_type t , time_type dt ) in do_step() argument
150 this->stepper().do_step_impl( system , x , dxdt , t , x , dt ); in do_step()
162 void do_step_dxdt_impl( System system, StateInOut &x, const DerivIn &dxdt, in do_step_dxdt_impl() argument
165 this->stepper().do_step_impl( system , x , dxdt , t , x , dt ); in do_step_dxdt_impl()
196 …do_step( System system , const StateIn &in , const DerivIn &dxdt , time_type t , StateOut &out , t… in do_step() argument
198 this->stepper().do_step_impl( system , in , dxdt , t , out , dt ); in do_step()
210 const DerivIn &dxdt, time_type t, StateOut &out, in do_step_dxdt_impl() argument
213 this->stepper().do_step_impl( system , in , dxdt , t , out , dt ); in do_step_dxdt_impl()
246 …do_step( System system , StateInOut &x , const DerivIn &dxdt , time_type t , time_type dt , Err &x… in do_step() argument
248 this->stepper().do_step_impl( system , x , dxdt , t , x , dt , xerr ); in do_step()
[all …]
|
| D | explicit_stepper_base.hpp | 135 do_step( System system , StateInOut &x , const DerivIn &dxdt , time_type t , time_type dt ) in do_step() argument
137 this->stepper().do_step_impl( system , x , dxdt , t , x , dt ); in do_step()
149 void do_step_dxdt_impl( System system, StateInOut &x, const DerivIn &dxdt, in do_step_dxdt_impl() argument
152 this->stepper().do_step_impl( system , x , dxdt , t , x , dt ); in do_step_dxdt_impl()
177 …void do_step( System system , const StateIn &in , const DerivIn &dxdt , time_type t , StateOut &ou… in do_step() argument
179 this->stepper().do_step_impl( system , in , dxdt , t , out , dt ); in do_step()
192 const DerivIn &dxdt, time_type t, StateOut &out, in do_step_dxdt_impl() argument
195 this->stepper().do_step_impl( system , in , dxdt , t , out , dt ); in do_step_dxdt_impl()
|
| /third_party/boost/libs/numeric/odeint/examples/openmp/ |
| D | phase_chain_omp_state.cpp | 37 void operator()( const state_type &x , state_type &dxdt , double /* t */ ) const in operator ()()
46 dxdt[n][m] = coupling_func( x[n][m+1] - x[n][m] ) + in operator ()()
49 dxdt[n][0] = coupling_func( x[n][1] - x[n][0] ); in operator ()()
52 dxdt[n][0] += coupling_func( x[n-1].back() - x[n].front() ); in operator ()()
54 dxdt[n][M-1] = coupling_func( x[n][M-2] - x[n][M-1] ); in operator ()()
57 dxdt[n][M-1] += coupling_func( x[n+1].front() - x[n].back() ); in operator ()()
|
| /third_party/boost/libs/numeric/odeint/performance/ |
| D | lorenz.hpp | 21 void inline operator()( const state_type &x , state_type &dxdt , const double t ) const in operator ()()
26 dxdt[0] = sigma * ( x[1] - x[0] ); in operator ()()
27 dxdt[1] = R * x[0] - x[1] - x[0] * x[2]; in operator ()()
28 dxdt[2] = x[0]*x[1] - b * x[2]; in operator ()()
|
| /third_party/boost/libs/numeric/odeint/examples/mpi/ |
| D | phase_chain.cpp | 38 void operator()( const state_type &x , state_type &dxdt , double /* t */ ) const in operator ()()
57 dxdt()[m] = coupling_func( x()[m+1] - x()[m] ) + in operator ()()
60 dxdt()[0] = coupling_func( x()[1] - x()[0] ); in operator ()()
64 dxdt()[0] += coupling_func( x_left - x().front() ); in operator ()()
66 dxdt()[M-1] = coupling_func( x()[M-2] - x()[M-1] ); in operator ()()
70 dxdt()[M-1] += coupling_func( x_right - x().back() ); in operator ()()
|
| /third_party/boost/libs/numeric/odeint/examples/ublas/ |
| D | lorenz_ublas.cpp | 18 void lorenz( const state_type &x , state_type &dxdt , const double t ) in lorenz() argument
24 dxdt[0] = sigma * ( x[1] - x[0] ); in lorenz()
25 dxdt[1] = R * x[0] - x[1] - x[0] * x[2]; in lorenz()
26 dxdt[2] = -b * x[2] + x[0] * x[1]; in lorenz()
|
| /third_party/boost/libs/numeric/odeint/test_external/gsl/ |
| D | check_gsl.cpp | 30 void lorenz( const state_type x , state_type dxdt , double t ) in lorenz() argument
32 gsl_vector_set( dxdt , 0 , sigma * ( gsl_vector_get(x , 1 ) - gsl_vector_get( x , 0 ) ) ); in lorenz()
33 …gsl_vector_set( dxdt , 1 , R * gsl_vector_get( x , 0 ) - gsl_vector_get( x , 1 ) - gsl_vector_get(… in lorenz()
34 …gsl_vector_set( dxdt , 2 , gsl_vector_get( x , 0 ) * gsl_vector_get( x , 1 ) - b * gsl_vector_get(… in lorenz()
|
| /third_party/boost/libs/numeric/odeint/examples/multiprecision/ |
| D | cmp_precision.cpp | 28 void rhs_m( const mp_50 x , mp_50 &dxdt , const mp_50 t ) in rhs_m() argument
30 dxdt = mp_50(3)/(mp_50(2)*t*t) + x/(mp_50(2)*t); in rhs_m()
33 void rhs_d( const double x , double &dxdt , const double t ) in rhs_d() argument
35 dxdt = 3.0/(2.0*t*t) + x/(2.0*t); in rhs_d()
|
| D | lorenz_mp.cpp | 33 void operator()( const state_type &x , state_type &dxdt , value_type t ) const in operator ()()
39 dxdt[0] = sigma * ( x[1] - x[0] ); in operator ()()
40 dxdt[1] = R * x[0] - x[1] - x[0] * x[2]; in operator ()()
41 dxdt[2] = -b * x[2] + x[0] * x[1]; in operator ()()
|
| /third_party/boost/libs/numeric/odeint/examples/gmpxx/ |
| D | lorenz_gmpxx.cpp | 32 void operator()( const state_type &x , state_type &dxdt , value_type t ) const in operator ()()
38 dxdt[0] = sigma * ( x[1] - x[0] ); in operator ()()
39 dxdt[1] = R * x[0] - x[1] - x[0] * x[2]; in operator ()()
40 dxdt[2] = -b * x[2] + x[0] * x[1]; in operator ()()
|