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
_______________________________________________
relax (http://www.nmr-relax.com)
This is the relax-devel mailing list
relax-devel@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-devel