Issue with VV integrator and domain decomposition of small systems
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.