Fixed free energy with LJ-PME
authorErik Lindahl <erik@kth.se>
Fri, 27 Jun 2014 09:04:29 +0000 (11:04 +0200)
committerGerrit Code Review <gerrit@gerrit.gromacs.org>
Sun, 29 Jun 2014 15:23:22 +0000 (17:23 +0200)
commitfc819eee87572d5608778642bcf22badaf3b4361
treebeca2948589a3dfacdc062762bfddb54fbfac022
parentec2c9197ed4b10a2ea6cb3186ffa10a3fcc3ccc3
Fixed free energy with LJ-PME

First, Only the normal LJ part of the potential was
shifted for LJ-PME with potential-shift
modifier in the free energy kernel,
and not the Ewald component. This had an effect
of roughly 0.05% on the total energy. Second,
Berk found that the verlet nbnxn SIMD LJ-PME kernel
had an issue that caused the potential and force
to be slightly off. With this patch, both
potential and force match the non-free-energy
kernel perfectly with at lambda=0.

Change-Id: I4bc53ce1bdb96ba06656b86470dbf8dbd1e81972
src/gromacs/gmxlib/ewald_util.c
src/gromacs/gmxlib/nonbonded/nb_free_energy.c
src/gromacs/mdlib/nbnxn_atomdata.c
src/gromacs/mdlib/nbnxn_search.c