margin and pressure in CUDA version

From: JC Gumbart (gumbart_at_physics.gatech.edu)
Date: Wed Sep 27 2017 - 11:48:52 CDT

I know this has been discussed at length previously, but I may have misunderstood something. I have a system that’s pretty small but the patch grid is still always at least 2 in any dimension. I first run a minimization/equilibration on a p100 with NAMD 2.12 with “margin 5”; the size in z stabilizes at ~48 A. I then run a continuation with no margin. The size in z quickly drops to ~42 A (within 5000 steps). This is still on the same GPU. If I turn “margin 5” on again, it stays around 48 A.

I’ve tested numerous combinations; it seems the CPU version always drops in size regardless of margin as do previous GPU versions. But what’s interesting is that in 2.12 it seems to depend on the specific GPU.

all cpu versions: 48 A -> 42 A
2.10 on k20c, 2.11 on p100: 48 -> 42
2.12 gpu on k20c: 48 -> 42
2.12 gpu on p100: 48 remains constant
2.12 gpu on k80: 48 -> 90 in 5k steps (????)

Regardless of version, if margin is less than 5, it drops. Note that for a couple hundred other runs of very similar systems before (same size within a few Angstroms in each dimension), I didn’t see this behavior. It appears to be very specific to the precise size of this one.

Any thoughts? Dave noted before that there should be other instabilities, but with the exception of the behavior on k80s, the change in size by 15% seems to be the only effect.

Thanks!
JC

This archive was generated by hypermail 2.1.6 : Mon Dec 31 2018 - 23:20:37 CST