Author: tlinnet Date: Fri Aug 1 18:09:34 2014 New Revision: 24908 URL: http://svn.gna.org/viewcvs/relax?rev=24908&view=rev Log: Changed to the matrix_exponential function for NS R1rho 2site. Modified: trunk/lib/dispersion/ns_r1rho_2site.py Modified: trunk/lib/dispersion/ns_r1rho_2site.py URL: http://svn.gna.org/viewcvs/relax/trunk/lib/dispersion/ns_r1rho_2site.py?rev=24908&r1=24907&r2=24908&view=diff ============================================================================== --- trunk/lib/dispersion/ns_r1rho_2site.py (original) +++ trunk/lib/dispersion/ns_r1rho_2site.py Fri Aug 1 18:09:34 2014 @@ -54,7 +54,7 @@ from numpy.ma import fix_invalid, masked_less # relax module imports. -from lib.dispersion.matrix_exponential import matrix_exponential_rank_NE_NS_NM_NO_ND_x_x +from lib.dispersion.matrix_exponential import matrix_exponential # Repetitive calculations (to speed up calculations). m_r1rho_prime = array([ @@ -235,7 +235,7 @@ R_mat = rr1rho_3d_2site_rankN(R1=r1, r1rho_prime=r1rho_prime, dw=dw, omega=omega, offset=offset, w1=spin_lock_fields, k_AB=k_AB, k_BA=k_BA, relax_time=relax_time) # This matrix is a propagator that will evolve the magnetization with the matrix R. - Rexpo_mat = matrix_exponential_rank_NE_NS_NM_NO_ND_x_x(R_mat) + Rexpo_mat = matrix_exponential(R_mat) # Magnetization evolution. Rexpo_M0_mat = einsum('...ij, ...jk', Rexpo_mat, M0)