Buffer estimate now works with rvd-switch=0
authorBerk Hess <hess@kth.se>
Fri, 28 Feb 2014 11:37:29 +0000 (12:37 +0100)
committerGerrit Code Review <gerrit@gerrit.gromacs.org>
Wed, 5 Mar 2014 06:22:02 +0000 (07:22 +0100)
commit928cab16974d7353e9e05d32c36f0b03ca3ec297
treec216f70ccb8275def6bf9bb19c6ddcb9738d8d40
parent7fca425b52f21ec4afd1aa230db4480a05edad32
Buffer estimate now works with rvd-switch=0

With rvdw-switch=0, the Verlet buffer estimation code would divide
by 0. The estimate for vdw-switch functions now use the exact 3rd
derivative of the potential at the cut-off iso estimating the 2nd.
grompp now prints a note with rvdw-switch < 0.5*rvdw.
Updated estimate formula in the manual and added a missing factor
of 2 (the code did use the factor of 2).

Change-Id: Ia259d703515f69b0ed3390e2f7be671e0a15bd5f
manual/algorithms.tex
src/gromacs/gmxpreprocess/calc_verletbuf.c
src/gromacs/gmxpreprocess/readir.c