Example #1
0
 void end_update_p(typename Hamiltonian::PointType& z,
                   Hamiltonian& hamiltonian, double epsilon,
                   callbacks::logger& logger) {
   z.p -= epsilon * hamiltonian.dphi_dq(z, logger);
 }