mailRe: r24628 - /branches/disp_spin_speed/lib/dispersion/ns_cpmg_2site_expanded.py


Others Months | Index by Date | Thread Index
>>   [Date Prev] [Date Next] [Thread Prev] [Thread Next]

Header


Content

Posted by Edward d'Auvergne on July 22, 2014 - 17:37:
Hi Troels,

For this change, it might be better to revert it.  This changes the
original Maple script from Nikolai Skrynnikov which I simply copied
and pasted unmodified.  It's not very important, but it might be
useful to preserve the original form.

Cheers,

Edward



On 22 July 2014 16:51,  <tlinnet@xxxxxxxxxxxxx> wrote:
Author: tlinnet
Date: Tue Jul 22 16:51:42 2014
New Revision: 24628

URL: http://svn.gna.org/viewcvs/relax?rev=24628&view=rev
Log:
Code validation of lib/dispersion/ns_cpmg_2site_expanded.py.

Task #7807 (https://gna.org/task/index.php?7807): Speed-up of dispersion 
models for Clustered analysis.

Modified:
    branches/disp_spin_speed/lib/dispersion/ns_cpmg_2site_expanded.py

Modified: branches/disp_spin_speed/lib/dispersion/ns_cpmg_2site_expanded.py
URL: 
http://svn.gna.org/viewcvs/relax/branches/disp_spin_speed/lib/dispersion/ns_cpmg_2site_expanded.py?rev=24628&r1=24627&r2=24628&view=diff
==============================================================================
--- branches/disp_spin_speed/lib/dispersion/ns_cpmg_2site_expanded.py   
(original)
+++ branches/disp_spin_speed/lib/dispersion/ns_cpmg_2site_expanded.py   Tue 
Jul 22 16:51:42 2014
@@ -74,9 +74,9 @@
     #tcp:=0.040/N;

     Ksym:=sqrt(Ka*Kb);
-    #dX:=(Ka-Kb+I*dw)/2;       # Ra=Rb
+    #dX:=(Ka-Kb+I*dw)/2;    # Ra=Rb
     dX:=((Ra-Rb)+(Ka-Kb)+I*dw)/2;
-
+
     L:=([[-dX, Ksym], [Ksym, dX]]);

     # in the end everything is multiplied by 
exp(-0.5*(Ra+Rb+Ka+Kb)*(Tc+2*tpalmer))
@@ -117,7 +117,7 @@
     cGG:=evalm(GG2&*Pspalmer&*GG1);

     #s0:=array([Kb, Ka]);
-    s0:=array([sqrt(Kb),sqrt(Ka)]);    # accounts for exchange 
symmetrization
+    s0:=array([sqrt(Kb),sqrt(Ka)]); # accounts for exchange symmetrization
     st:=evalm(cGG&*s0);
     #obs:=(1/(Ka+Kb))*st[1];
     obs:=(sqrt(Kb)/(Ka+Kb))*st[1];  # accounts for exchange symmetrization


_______________________________________________
relax (http://www.nmr-relax.com)

This is the relax-commits mailing list
relax-commits@xxxxxxx

To unsubscribe from this list, get a password
reminder, or change your subscription options,
visit the list information page at
https://mail.gna.org/listinfo/relax-commits



Related Messages


Powered by MHonArc, Updated Tue Jul 22 18:00:12 2014