[gmx-developers] Problem with dihedral restraint in gromacs 4.6

[Landry Charlier]

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.
>
> Cheers,
> Bogdan
> --
> gmx-developers mailing list
> gmx-developers at gromacs.org
> http://lists.gromacs.org/mailman/listinfo/gmx-developers
> Please don't post (un)subscribe requests to the list. Use the
> www interface or send it to gmx-developers-request at gromacs.org.
>