mailr21684 - /trunk/specific_analyses/relax_disp/optimisation.py


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

Header


Content

Posted by edward on November 27, 2013 - 17:31:
Author: bugman
Date: Wed Nov 27 17:31:13 2013
New Revision: 21684

URL: http://svn.gna.org/viewcvs/relax?rev=21684&view=rev
Log:
Fix for the specific_analyses.relax_disp.optimisation.back_calc_r2eff() 
function.

The R2eff error data structure when the cpmg_frqs or spin_lock_nu1 argument 
is supplied was all
zeros, whereas it should all be ones.  This was causing many divide by zero 
numpy warnings to appear
on certain operating systems (Mac OS X).


Modified:
    trunk/specific_analyses/relax_disp/optimisation.py

Modified: trunk/specific_analyses/relax_disp/optimisation.py
URL: 
http://svn.gna.org/viewcvs/relax/trunk/specific_analyses/relax_disp/optimisation.py?rev=21684&r1=21683&r2=21684&view=diff
==============================================================================
--- trunk/specific_analyses/relax_disp/optimisation.py (original)
+++ trunk/specific_analyses/relax_disp/optimisation.py Wed Nov 27 17:31:13 
2013
@@ -25,7 +25,7 @@
 # Python module imports.
 from minfx.generic import generic_minimise
 from minfx.grid import grid
-from numpy import dot, float64, int32, zeros
+from numpy import dot, float64, int32, ones, zeros
 from numpy.linalg import inv
 from re import search
 import sys
@@ -110,7 +110,7 @@
                     else:
                         num = len(spin_lock_nu1[exp_type_index][frq_index])
                     values[exp_type_index][0].append(zeros(num, float64))
-                    errors[exp_type_index][0].append(zeros(num, float64))
+                    errors[exp_type_index][0].append(ones(num, float64))
                     missing[exp_type_index][0].append(zeros(num, int32))
 
     # Initialise the relaxation dispersion fit functions.




Related Messages


Powered by MHonArc, Updated Wed Nov 27 18:20:02 2013