Bug #1669
Issue with VV integrator and domain decomposition of small systems
Description
At the first MD step the virial computation for the VV integrator requires a backup of the velocities. The backup array was allocated using the global instead of the local atom count. For small systems the local atom count, including communicated and constraint atoms, can be larger than the global atom count. This can cause crashes due to a buffer out of range.
Associated revisions
History
#1 Updated by Gerrit Code Review Bot over 4 years ago
Gerrit received a related patchset '1' for Issue #1669.
Uploader: Berk Hess (hess@kth.se)
Change-Id: I64947405c138f601db7daa4f9628a04cff9fa8bb
Gerrit URL: https://gerrit.gromacs.org/4359
#2 Updated by Mark Abraham about 4 years ago
- Status changed from New to Closed
Fix allocation issue with VV integrator
Allocation of temp buffers for VV and VVAK now uses state->natoms
instead of top_global->natoms.
Fixes #1669.
Change-Id: I64947405c138f601db7daa4f9628a04cff9fa8bb