Example #1
0
 void update_q(typename Hamiltonian::PointType& z,
               Hamiltonian& hamiltonian, double epsilon,
               callbacks::logger& logger) {
   z.q += epsilon * hamiltonian.dtau_dp(z);
   hamiltonian.update_potential_gradient(z, logger);
 }