knowing which atoms trigger the "Non-numeric coordinates" issue

Hi all,

Just a brief question on working with this issue. Is there a 'simple'
way to know which atom (i.e. its coordinates, ID etc) triggers the NaN
coordinates in the calculation? I have been able to make due with
dumping at every single time step and iterate through the output
one-by-one around the time of crash, which is notoriously inefficient
just from the sound of it. Is there a more convenient way to achieve
this?

Thanks,
QT

Hi all,

Just a brief question on working with this issue. Is there a 'simple'
way to know which atom (i.e. its coordinates, ID etc) triggers the NaN
coordinates in the calculation? I have been able to make due with
dumping at every single time step and iterate through the output
one-by-one around the time of crash, which is notoriously inefficient
just from the sound of it. Is there a more convenient way to achieve
this?

how about this?

<define a dump assuming the label is '1'>

variable fsum atom fx+fy+fz
dump_modify 1 thresh v_fsum > 1.0e100

the non-numeric number is larger than the largest numeric value, thus
the dump will *only* output atoms that have NaN as forces.

axel.