[gmx-users] Velocity Verlet integrator

[Michael Shirts]

Putting both velocity Verlet and leapfrog Verlet both in Gromacs turns
out to be non-trivial for the bookkeeping.  The easiest way to do this
was split the velocity Verlet updates.

Also, the additional computational cost of two half steps for
velocities is trivial compared to the cost of the forces for almost
all systems.

On Sun, Mar 23, 2014 at 8:39 AM, Mario Fernández Pendás
<mariofp77 at gmail.com> wrote:
> Dear all,
>
> In terms of computational efficiency, why the velocity Verlet integrator is
> implemented in GROMACS in one full step for positions and two half steps
> for velocities?
>
> Would it be more efficient to merge the second halft step for velocities
> with the first half step of the following scheme, ie, integrating in one
> full step for velocities?
>
> Thank you very much,
> Mario Fernánez-Pendás
> --
> Gromacs Users mailing list
>
> * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!
>
> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>
> * For (un)subscribe requests visit
> https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-request at gromacs.org.