r24632 - /branches/disp_spin_speed/lib/dispersion/ns_r1rho_2site.py -- July 22, 2014 - 16:51

Author: tlinnet
Date: Tue Jul 22 16:51:50 2014
New Revision: 24632

URL: http://svn.gna.org/viewcvs/relax?rev=24632&view=rev
Log:
Code validation of lib/dispersion/ns_r1rho_2site.py.

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

Modified:
    branches/disp_spin_speed/lib/dispersion/ns_r1rho_2site.py

Modified: branches/disp_spin_speed/lib/dispersion/ns_r1rho_2site.py
URL: 
http://svn.gna.org/viewcvs/relax/branches/disp_spin_speed/lib/dispersion/ns_r1rho_2site.py?rev=24632&r1=24631&r2=24632&view=diff
==============================================================================
--- branches/disp_spin_speed/lib/dispersion/ns_r1rho_2site.py   (original)
+++ branches/disp_spin_speed/lib/dispersion/ns_r1rho_2site.py   Tue Jul 22 
16:51:50 2014
@@ -162,8 +162,8 @@
     #d = W - offset
 
     # Alias to original parameter name.
-    wA=dA
-    wB=dB
+    wA = dA
+    wB = dB
 
     # Multiply and expand.
     mat_r1rho_prime = multiply.outer( r1rho_prime * relax_time, 
m_r1rho_prime )
@@ -238,10 +238,10 @@
     Rexpo_mat = matrix_exponential_rank_NE_NS_NM_NO_ND_x_x(R_mat)
 
     # Magnetization evolution.
-    Rexpo_M0_mat = einsum('...ij,...jk', Rexpo_mat, M0)
+    Rexpo_M0_mat = einsum('...ij, ...jk', Rexpo_mat, M0)
 
     # Magnetization evolution, which include all dimensions.
-    MA_mat = einsum('...ij,...jk', M0_T, Rexpo_M0_mat)[:, :, :, :, :, 0, 0]
+    MA_mat = einsum('...ij, ...jk', M0_T, Rexpo_M0_mat)[:, :, :, :, :, 0, 0]
 
     # Insert safe checks.
     if min(MA_mat) < 0.0: