AMBER: distcovar in ptraj

From: Zhenyu Lu <>
Date: Mon, 26 Mar 2007 21:00:13 -0700 (PDT)

Hi, for "distcovar" calculations with ptraj, the code in actions.c calculates the matrix element for pairs ij, kl, with the condition k>=i, and l >=max(k+1,j). But does this skip those matrix elements satisfying j>k>i, k+1<l<j?

or should the code be changed from
> atcnt4 = (k>=j ? atcnt3+1 : atcnt2);

> for(l=(k>=j ? k+1 : j); l < action->state->atoms; l++){
< atcnt4 = (atcnt3==atcnt1 ? atcnt2 : atcnt3+1);
< for(l=(k==i ? j : k+1); l < action->state->atoms; l++){
It's here! Your new message!  
Get new email alerts with the free Yahoo! Toolbar.
The AMBER Mail Reflector
To post, send mail to
To unsubscribe, send "unsubscribe amber" to
Received on Wed Mar 28 2007 - 06:07:33 PDT
Custom Search