r21589 - /branches/relax_disp/lib/dispersion/mq_cr72.py -- November 21, 2013 - 22:17

Author: bugman
Date: Thu Nov 21 22:17:06 2013
New Revision: 21589

URL: http://svn.gna.org/viewcvs/relax?rev=21589&view=rev
Log:
Fix for the 'MQ CR72' model for MQ-type data.

The check to prevent acos of a number less than 1 has been changed to switch 
the sign rather than to
set the back calculated R2eff to 1e99.


Modified:
    branches/relax_disp/lib/dispersion/mq_cr72.py

Modified: branches/relax_disp/lib/dispersion/mq_cr72.py
URL: 
http://svn.gna.org/viewcvs/relax/branches/relax_disp/lib/dispersion/mq_cr72.py?rev=21589&r1=21588&r2=21589&view=diff
==============================================================================
--- branches/relax_disp/lib/dispersion/mq_cr72.py (original)
+++ branches/relax_disp/lib/dispersion/mq_cr72.py Thu Nov 21 22:17:06 2013
@@ -123,8 +123,7 @@
         # Part of the equation (catch values < 1 to prevent math domain 
errors).
         part = Dpos * cosh(etapos) - Dneg * cos(etaneg)
         if part.real < 1.0:
-            back_calc[i] = 1e100
-            continue
+            part = -part
 
         # The first eigenvalue.
         lambda1 = r20_kex - cpmg_frqs[i] * arccosh(part)