AMBER: distcovar in ptraj

From: Zhenyu Lu <biosimu.yahoo.com>
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++){
to
< atcnt4 = (atcnt3==atcnt1 ? atcnt2 : atcnt3+1);
< for(l=(k==i ? j : k+1); l < action->state->atoms; l++){ 
---
Thanks!
Zhenyu
 
____________________________________________________________________________________
It's here! Your new message!  
Get new email alerts with the free Yahoo! Toolbar.
http://tools.search.yahoo.com/toolbar/features/mail/
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Wed Mar 28 2007 - 06:07:33 PDT
Custom Search