[gmx-users] CG Force and Potential through do_force()

[Sikandar Mashayak]

Hi

I want to do Coarse-Grained simulation in Gromacs. I have one CG force
formula, which I plan to implement by adding a short code in do_force()
routine in Gromacs 4.0.7 source code.

After browsing through the do_force() code, I am now familiar with data
structure relevant to add the new force. I mean, I can specify new force by
adding to f[indx][XX] and so on.

But I am not very sure about potential, whether do I have to explicitly
specify the potential on atoms due to this force or gromacs automatically
takes care of it by knowing the forces on each atoms?

thanks
sikandar
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20100324/77887aed/attachment.html>