r21483 - /branches/relax_disp/lib/dispersion/tap03.py -- November 17, 2013 - 13:50

Author: bugman
Date: Sun Nov 17 13:50:14 2013
New Revision: 21483

URL: http://svn.gna.org/viewcvs/relax?rev=21483&view=rev
Log:
The lib.dispersion.tap03 module can now handle negative gamma values.

This avoids fatal errors during optimisation.


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

Modified: branches/relax_disp/lib/dispersion/tap03.py
URL: 
http://svn.gna.org/viewcvs/relax/branches/relax_disp/lib/dispersion/tap03.py?rev=21483&r1=21482&r2=21483&view=diff
==============================================================================
--- branches/relax_disp/lib/dispersion/tap03.py (original)
+++ branches/relax_disp/lib/dispersion/tap03.py Sun Nov 17 13:50:14 2013
@@ -89,6 +89,11 @@
         # The gamma factor.
         gamma = 1.0 + phi_ex*(sigma2 - kex2 + spin_lock_fields2[i]) / 
(sigma2 + kex2 + spin_lock_fields2[i])
 
+        # Bad gamma.
+        if gamma < 0.0:
+            back_calc[i] = 1e100
+            continue
+
         # We assume that A resonates at 0 [s^-1], without loss of generality.
         da = Wa - offset[i]                         # Offset of spin-lock 
from A.
         db = Wb - offset[i]                         # Offset of spin-lock 
from B.