void initialize( System system , const CoorIn & qin , const VelocityIn & pin , time_type t )
{
m_resizer.adjust_size( qin ,
detail::bind( &velocity_verlet::template resize_impl< CoorIn > ,
detail::ref( *this ) , detail::_1 ) );
initialize_acc( system , qin , pin , t );
}
void
bridge_document_rep::initialize () {
int i, n= N(st);
brs= array<bridge> (n);
for (i=0; i<n; i++)
brs[i]= make_bridge (ttt, st[i], descend (ip, i));
initialize_acc ();
}
void do_step_v1( System system , StateInOut & x , time_type t , time_type dt )
{
typedef typename odeint::unwrap_reference< StateInOut >::type state_in_type;
typedef typename odeint::unwrap_reference< typename state_in_type::first_type >::type coor_in_type;
typedef typename odeint::unwrap_reference< typename state_in_type::second_type >::type momentum_in_type;
typedef typename boost::remove_reference< coor_in_type >::type xyz_type;
state_in_type & statein = x;
coor_in_type & qinout = statein.first;
momentum_in_type & pinout = statein.second;
// alloc a
if( m_resizer.adjust_size( qinout ,
detail::bind( &velocity_verlet::template resize_impl< xyz_type > ,
detail::ref( *this ) , detail::_1 ) )
|| m_first_call )
{
initialize_acc( system , qinout , pinout , t );
}
// check first
do_step( system , qinout , pinout , get_current_acc() , qinout , pinout , get_old_acc() , t , dt );
toggle_current_acc();
}
