Author: tlinnet Date: Tue Jul 22 16:51:35 2014 New Revision: 24624 URL: http://svn.gna.org/viewcvs/relax?rev=24624&view=rev Log: Code validation of lib/dispersion/m61b.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/m61b.py Modified: branches/disp_spin_speed/lib/dispersion/m61b.py URL: http://svn.gna.org/viewcvs/relax/branches/disp_spin_speed/lib/dispersion/m61b.py?rev=24624&r1=24623&r2=24624&view=diff ============================================================================== --- branches/disp_spin_speed/lib/dispersion/m61b.py (original) +++ branches/disp_spin_speed/lib/dispersion/m61b.py Tue Jul 22 16:51:35 2014 @@ -64,6 +64,7 @@ from numpy import any, isfinite, min, sum from numpy.ma import fix_invalid, masked_where + def r1rho_M61b(r1rho_prime=None, pA=None, dw=None, kex=None, spin_lock_fields2=None, back_calc=None): """Calculate the R1rho values for the M61 skew model. @@ -108,7 +109,7 @@ denom = kex2_pA2dw2 + spin_lock_fields2 # Catch math domain error of dividing with 0. - # This is when denom=0. + # This is when denom = 0. mask_denom_zero = denom == 0.0 if any(mask_denom_zero): t_denom_zero = True