mailr23048 - /trunk/lib/dispersion/b14.py


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

Header


Content

Posted by tlinnet on May 07, 2014 - 15:54:
Author: tlinnet
Date: Wed May  7 15:54:49 2014
New Revision: 23048

URL: http://svn.gna.org/viewcvs/relax?rev=23048&view=rev
Log:
Replaced "v2pPdN" with v5, to follow paper and manual.

sr #3154: (https://gna.org/support/?3154) Implementation of Baldwin (2014) 
B14 model - 2-site exact solution model for all time scales.

This follows the tutorial for adding relaxation dispersion models at:
http://wiki.nmr-relax.com/Tutorial_for_adding_relaxation_dispersion_models_to_relax#Debugging

Modified:
    trunk/lib/dispersion/b14.py

Modified: trunk/lib/dispersion/b14.py
URL: 
http://svn.gna.org/viewcvs/relax/trunk/lib/dispersion/b14.py?rev=23048&r1=23047&r2=23048&view=diff
==============================================================================
--- trunk/lib/dispersion/b14.py (original)
+++ trunk/lib/dispersion/b14.py Wed May  7 15:54:49 2014
@@ -222,9 +222,9 @@
     y = power( (v1c - v3) / (v1c + v3), ncyc)
 
     # Off diagonal common factor. sinh fuctions.
-    v2pPdN = (-deltaR2 + kex + dw*1j) * v1s - 2. * (oGt2 + k_AB * t1pt2) * 
ex1c
-
-    Tog = 0.5 * (1. + y) + (1. - y) * v2pPdN / (2. * v3 * N )
+    v5 = (-deltaR2 + kex + dw*1j) * v1s - 2. * (oGt2 + k_AB * t1pt2) * ex1c
+
+    Tog = 0.5 * (1. + y) + (1. - y) * v5 / (2. * v3 * N )
 
     # Estimate R2eff. relax_time = Trel = 1/inv_tcpmg.
     Minty = Rpre - ncyc * inv_tcpmg * arccosh(v1c.real) - inv_tcpmg * 
log(Tog.real)




Related Messages


Powered by MHonArc, Updated Wed May 07 16:00:02 2014