Camilo Calderon wrote:
Hi all!
I am new to the GROMACS-world, and I am mystified.
How do I implement the 3 body potentials advertised in sec. 4.2.8 and
sec. 4.2.9 of the GROMACS 3.3 manual?
Right now I have the following:
[ cross_cross_bond ]
;ai aj funct r1e r2e krr
<numbers here>
[ cross_bond_angle ]
;ai aj ak funct r1e r2e r3e krt
<more numbers here>
when I grompp the topol.top, what I get is a Error, with a complaint
"invalid directive", and then a Fatal Error, with the complaint
"Incorrect number of parameters - found 5, expected 2 or 4." What is
really tough about this is that I've been poring over the source code,
trying to decypher what input format ought to be - there is nothing on
this either in the Manual, or in the various lists!
The contents of either gmxlib/ifunc.c or kernel/ifunc.c (sorry can't
remember which it is) may help you work out how many parameters it needs
- text search for CROSS, and then find the definition of the macro
that is being used there.
Can you give us more detail from the .top file and the error messages
grompp returns?
Mark
_______________________________________________
gmx-users mailing list gmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php