GPU code misses settle error check, simulation crashes with segfault without any further output
Running the 216 SPC water box on a single GPU with cut-off electrostatics segfaults without any other warning.
Find attached a stack-trace of the core dump.
- Status changed from New to Rejected
You are running NVE with a plain Coulomb cut-off at 0.4 nm, so of course the energy goes up and the system explodes.
grompp prints a note about using plain cut-off. Maybe this should be a warning?
This runs fine with reaction-field.
#2 Updated by Christian Blau 2 months ago
The issue here is that it just crashes out without any information, whereas for example in GROMACS2018.6 I get
step 30096: One or more water molecules can not be settled.
Check for bad contacts and/or reduce the timestep if appropriate.
Wrote pdb files with previous and current coordinates
and at least the pdb files with the respective coordinates
#4 Updated by Artem Zhmurov 2 months ago
Anyone knows any standard solution for checking for errors inside the CUDA kernel? I can do conditional atomicSet(..) or conditionally set the variable in host memory directly from the kernel so that it can be checked in the host code. But I don't think that this is a proper solution.