Hi,

There have been reports about inconsistencies between older (<= 4.0.7?) and newer versions of g_hbond, where the older seem to have been more reliable. I found and killed the bug that caused the newer versions to miscount the hbonds. Checkout release-4-5-patches to get the bugfix, or patch it yourself by commenting out line 1497, which reads "return hbNo;":

     if (bBox){
if (d>a && bMerge && (bContact || isInterchangable(hb, d, a, grpd, grpa))) { /* acceptor is also a donor a\
nd vice versa? */
            /* return hbNo; */
daSwap = TRUE; /* If so, then their history should be filed with donor and acceptor swapped. */
        }

Simple as that.

--
-----------------------------------------------
Erik Marklund, PhD student
Dept. of Cell and Molecular Biology, Uppsala University.
Husargatan 3, Box 596,    75124 Uppsala, Sweden
phone:    +46 18 471 4537        fax: +46 18 511 755
er...@xray.bmc.uu.se    http://folding.bmc.uu.se/

--
gmx-users mailing list    gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

Reply via email to