[gmx-users] NH2 angles wrong in calch.c

[Charles Karney]

The documentation of the hydrogen database (Gromacs Manual 3.2, p 94)
says that "type 3" hydrogens are added with bond angles of 120 (which is
right).

However the code, v 3.2.1, src/gmslib/calch.c, lines 182-187, uses the
tetrahedral angle 109.5 deg.  Surely this is wrong.

Minor gripe: why not define the tetrahedral angle exactly?

#define alfaH   (acos(-1/3.0))

-- 
Charles Karney                  Email:  ckarney at sarnoff.com
Sarnoff Corporation             URL:    http://charles.karney.info
201 Washington Rd               Phone:  +1 609 734 2312
Princeton, NJ 08543-5300        Fax:    +1 609 734 2323