mailr25475 - /trunk/specific_analyses/relax_disp/estimate_r2eff.py


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

Header


Content

Posted by tlinnet on August 31, 2014 - 08:49:
Author: tlinnet
Date: Sun Aug 31 08:49:57 2014
New Revision: 25475

URL: http://svn.gna.org/viewcvs/relax?rev=25475&view=rev
Log:
In module for estimating R2eff errors, removed " values, errors" to be send 
to function for gradient, since they are not used.

task #7822(https://gna.org/task/index.php?7822): Implement user function to 
estimate R2eff and associated errors for exponential curve fitting.

Modified:
    trunk/specific_analyses/relax_disp/estimate_r2eff.py

Modified: trunk/specific_analyses/relax_disp/estimate_r2eff.py
URL: 
http://svn.gna.org/viewcvs/relax/trunk/specific_analyses/relax_disp/estimate_r2eff.py?rev=25475&r1=25474&r2=25475&view=diff
==============================================================================
--- trunk/specific_analyses/relax_disp/estimate_r2eff.py        (original)
+++ trunk/specific_analyses/relax_disp/estimate_r2eff.py        Sun Aug 31 
08:49:57 2014
@@ -414,17 +414,13 @@
         return Kw
 
 
-    def func_exp_grad(self, params=None, times=None, values=None, 
errors=None):
+    def func_exp_grad(self, params=None, times=None):
         """The gradient (Jacobian matrix) of func_exp for Co-variance 
calculation.
 
         @param params:  The vector of parameter values.
         @type params:   numpy rank-1 float array
         @keyword times: The time points.
         @type times:    numpy array
-        @param values:  The measured values.
-        @type values:   numpy array
-        @param errors:  The standard deviation of the measured intensity 
values per time point.
-        @type errors:   numpy array
         @return:        The Jacobian matrix with 'm' rows of function 
derivatives per 'n' columns of parameters.
         @rtype:         numpy array
         """
@@ -490,7 +486,7 @@
         back_calc = self.func_exp(params=params, times=times)
 
         # Get the Jacobian, with partial derivative, with respect to r2eff 
and i0.
-        exp_grad = self.func_exp_grad(params=params, times=times, 
values=values, errors=errors)
+        exp_grad = self.func_exp_grad(params=params, times=times)
 
         # Transpose back, to get rows.
         exp_grad_t = transpose(exp_grad)
@@ -890,7 +886,7 @@
 
         else:
             # Use the direct Jacobian from python.
-            jacobian_matrix_exp = E.func_exp_grad(params=param_vector, 
times=E.times, values=E.values, errors=E.errors)
+            jacobian_matrix_exp = E.func_exp_grad(params=param_vector, 
times=E.times)
             weights = 1. / E.errors**2
 
     pcov = multifit_covar(J=jacobian_matrix_exp, weights=weights)




Related Messages


Powered by MHonArc, Updated Sun Aug 31 09:00:06 2014