mailr21592 - /branches/relax_disp/target_functions/relax_disp.py


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

Header


Content

Posted by edward on November 21, 2013 - 22:48:
Author: bugman
Date: Thu Nov 21 22:48:28 2013
New Revision: 21592

URL: http://svn.gna.org/viewcvs/relax?rev=21592&view=rev
Log:
Fixes for the 'MQ CR72' dispersion model target function.


Modified:
    branches/relax_disp/target_functions/relax_disp.py

Modified: branches/relax_disp/target_functions/relax_disp.py
URL: 
http://svn.gna.org/viewcvs/relax/branches/relax_disp/target_functions/relax_disp.py?rev=21592&r1=21591&r2=21592&view=diff
==============================================================================
--- branches/relax_disp/target_functions/relax_disp.py (original)
+++ branches/relax_disp/target_functions/relax_disp.py Thu Nov 21 22:48:28 
2013
@@ -264,7 +264,7 @@
                     self.spin_lock_omega1_squared[exp_type_index][frq_index] 
= self.spin_lock_omega1[exp_type_index][frq_index] ** 2
 
         # The inverted relaxation delay.
-        if model in [MODEL_MMQ_2SITE, MODEL_NS_CPMG_2SITE_3D, 
MODEL_NS_CPMG_2SITE_3D_FULL, MODEL_NS_CPMG_2SITE_EXPANDED, 
MODEL_NS_CPMG_2SITE_STAR, MODEL_NS_CPMG_2SITE_STAR_FULL, 
MODEL_NS_R1RHO_2SITE]:
+        if model in [MODEL_MQ_CR72, MODEL_MMQ_2SITE, MODEL_NS_CPMG_2SITE_3D, 
MODEL_NS_CPMG_2SITE_3D_FULL, MODEL_NS_CPMG_2SITE_EXPANDED, 
MODEL_NS_CPMG_2SITE_STAR, MODEL_NS_CPMG_2SITE_STAR_FULL, 
MODEL_NS_R1RHO_2SITE]:
             self.inv_relax_times = 1.0 / relax_times
 
         # Special storage matrices for the multi-quantum CPMG 2-site 
numerical model.
@@ -1009,7 +1009,6 @@
         k_AB = pB * kex
 
         # Initialise.
-        aliased_dwH = 0.0
         chi2_sum = 0.0
 
         # Loop over the experiment types.
@@ -1026,6 +1025,7 @@
                     dwH_frq = dwH[spin_index] * 
self.frqs_H[exp_index][spin_index][frq_index]
 
                     # Alias the dw frequency combinations.
+                    aliased_dwH = 0.0
                     if self.exp_types[exp_index] == EXP_TYPE_CPMG_SQ:
                         aliased_dw = dw_frq
                     elif self.exp_types[exp_index] == 
EXP_TYPE_CPMG_PROTON_SQ:




Related Messages


Powered by MHonArc, Updated Thu Nov 21 23:00:03 2013