[gmx-developers] Problem with dihedral restraint in gromacs 4.6
landry.charlier at univ-montp1.fr
Mon Feb 11 14:49:51 CET 2013
So it's not a bug but a mistake in the manual.
Thank you for your help Bodgan.
Quoting Bogdan Costescu <bcostescu at gmail.com>:
> On Mon, Feb 11, 2013 at 11:51 AM, Landry Charlier
> <landry.charlier at univ-montp1.fr> wrote:
>> Incorrect number of parameters - found 2, expected 3 or 6 for Dih. Rest..
>> grompp waits a third parameter but the manual p 133 says that only 2
>> parameters (phi & dphi) are required.
> The code expects indeed 3 (or 6 for free energy calculations)
> parameters. They are phi, dphi and kfac. You can look at the code in
> src/gmxlib/bondfree.c::dihres() for how they are used.
> gmx-developers mailing list
> gmx-developers at gromacs.org
> Please don't post (un)subscribe requests to the list. Use the
> www interface or send it to gmx-developers-request at gromacs.org.
More information about the gromacs.org_gmx-developers