[gmx-developers] Inconsistent force calculating result between generic kernel and interaction-specific kernel?

[Mark Tianwu Zang]

Dear all,
I have inserted only a few lines after update() in md.c to monitor how
forces change after every step. My codes are very simple, just like:

for i=0 to md->nalloc
 fprintf f[i][0], f[i][1], f[i][2]

However, I found the results of my output become quite different after
exporting GMX_NB_GENERIC=1, which means the forces calculated by generic
kernel and interaction-specific kernel are not the same. I am a little
confused now.. It this phenomena quite normal or I ignored something
important?
Thanks a lot!
-Mark
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-developers/attachments/20130607/89cccbda/attachment.html>