r23735 - /branches/disp_spin_speed/target_functions/relax_disp.py -- June 08, 2014 - 19:48

Author: tlinnet
Date: Sun Jun  8 19:48:31 2014
New Revision: 23735

URL: http://svn.gna.org/viewcvs/relax?rev=23735&view=rev
Log:
Important fix for the creation of the multi dimensional pA numpy array.

It should be created as numpy.zeros([ei][si][mi][oi]) instead of 
numpy.ones([ei][si][mi][oi]).

This allows for rapid testing of all dimensions with np.allclose(pA, 
numpy.ones(dw.shape)).
pA can have missing filled out values, when the number of dispersion points 
are different
per spectrometer frequency.

Task #7807 (https://gna.org/task/index.php?7807): Speed-up of dispersion 
models for Clustered analysis.

Modified:
    branches/disp_spin_speed/target_functions/relax_disp.py

Modified: branches/disp_spin_speed/target_functions/relax_disp.py
URL: 
http://svn.gna.org/viewcvs/relax/branches/disp_spin_speed/target_functions/relax_disp.py?rev=23735&r1=23734&r2=23735&view=diff
==============================================================================
--- branches/disp_spin_speed/target_functions/relax_disp.py     (original)
+++ branches/disp_spin_speed/target_functions/relax_disp.py     Sun Jun  8 
19:48:31 2014
@@ -411,7 +411,7 @@
             # The number of disp point can change per spectrometer, so we 
make the maximum size.
             self.R20A_a = np.ones(back_calc_shape + 
[self.max_num_disp_points])
             self.R20B_a = np.ones(back_calc_shape + 
[self.max_num_disp_points])
-            self.pA_a = np.ones(back_calc_shape + [self.max_num_disp_points])
+            self.pA_a = np.zeros(back_calc_shape + 
[self.max_num_disp_points])
             self.dw_frq_a = np.ones(back_calc_shape + 
[self.max_num_disp_points])
             self.kex_a = np.ones(back_calc_shape + 
[self.max_num_disp_points])
             self.cpmg_frqs_a = np.ones(back_calc_shape + 
[self.max_num_disp_points])