Hi all,
I have implemented a model where I add a new property via property/atom.
It is uses the bond atom style and it creates a bondlist. Atoms can be added and removed from the bondlist at every timestep.
Once the bondlist is created it manipulates the property in the following fashion:
New_property[i]=New_property[i]/num_bond[i];
This seems to give a NAN result in parallel but correct result in serial. When I checked I found that the num_bond seems to be zero when I run it in parallel even though the 2 atoms are in the bondlist. There have been no other changes to the original source code. Does anyone have any idea why the bondlist is not updating in the case when it is run in parallel?
Thanks in advance for any help.
Liza