Remove COM from SETTLE
authorBerk Hess <hess@kth.se>
Fri, 5 Jun 2020 11:53:37 +0000 (11:53 +0000)
committerPaul Bauer <paul.bauer.q@gmail.com>
Fri, 5 Jun 2020 11:53:37 +0000 (11:53 +0000)
commit05556dbc8fb19d928055de6600894e719f901983
tree5f826671a8b9dd8a964f13c0b64756409a6ea38b
parentd879b6432d552604caf8c253c821ae05d26e8a69
Remove COM from SETTLE

SETTLE computed new postions of the atoms using the center of mass
of the molecule. But this adds rouding errors which lead to extra
energy drift. The use of the COM is now completely avoided, which
significantly improves energy conservation when coordinates are large
and also slighlty improves performance.
Corrected and updated the flop accounting for SETTLE.
docs/reference-manual/algorithms/constraint-algorithms.rst
docs/release-notes/2021/major/features.rst
src/gromacs/gmxlib/nrnb.cpp
src/gromacs/mdlib/settle.cpp
src/gromacs/mdlib/settle_gpu.cu