r25543 - /branches/est_par_error/specific_analyses/relax_disp/estimate_r2eff.py -- September 02, 2014 - 14:03

Author: tlinnet
Date: Tue Sep  2 14:03:37 2014
New Revision: 25543

URL: http://svn.gna.org/viewcvs/relax?rev=25543&view=rev
Log:
Fix for dimension extractions, when estimating the parameter errors.

task #7824(https://gna.org/task/index.php?7824): Model parameter ERROR 
estimation from Jacobian and Co-variance matrix of dispersion models.

Modified:
    branches/est_par_error/specific_analyses/relax_disp/estimate_r2eff.py

Modified: 
branches/est_par_error/specific_analyses/relax_disp/estimate_r2eff.py
URL: 
http://svn.gna.org/viewcvs/relax/branches/est_par_error/specific_analyses/relax_disp/estimate_r2eff.py?rev=25543&r1=25542&r2=25543&view=diff
==============================================================================
--- branches/est_par_error/specific_analyses/relax_disp/estimate_r2eff.py     
  (original)
+++ branches/est_par_error/specific_analyses/relax_disp/estimate_r2eff.py     
  Tue Sep  2 14:03:37 2014
@@ -250,7 +250,7 @@
         # NM: Number of spectrometer frequencies.
         # NO: Maximum number of offsets.
         # ND: Number of dispersion(data) points.
-        NJ, NE, NS, NM, NO, ND = func_jacobian.shape
+        ND, NE, NS, NM, NO, NJ = func_jacobian.shape
 
         if NS != 1:
             raise RelaxError("The number of spins does not fit.")
@@ -268,7 +268,7 @@
             cur_weights = weights[ei, si, mi, oi]
 
             # Extract every column/row from the first to last columns. Is 
this correct?
-            cur_jacobian = func_jacobian[0:NJ:1, ei, si, mi, oi]
+            cur_jacobian = func_jacobian[0:ND:1, ei, si, mi, oi]
 
             # Get the co-variance
             pcov = multifit_covar(J=cur_jacobian, weights=cur_weights)