[gmx-developers] Re: g_dist and vsite

Ilya Chorny ichorny at gmail.com
Wed Mar 10 18:51:21 CET 2010


I have been having problems calculating the distance between CA and CB atoms
in ALA using g_dist. The trajectory was run using vsite so the mass of CB in
ALA is set to 0. Since g_dist calculates distances between the center of
mass of the groups if one of the groups has a mass of zero there will be an
error.

I am going to edit gmx_dist as follows. Does anyone see any glaring issues
with doing this?

Thanks,

Ilya

previous:

com[g][d] /= mass[g];

new:

if mass[g] > 0.0
        com[g][d] /= mass[g];


-- 
Ilya Chorny Ph.D.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-developers/attachments/20100310/bff36ccf/attachment.html>


More information about the gromacs.org_gmx-developers mailing list