r21195 - /branches/relax_disp/specific_analyses/relax_disp/disp_data.py -- October 21, 2013 - 09:57

Author: bugman
Date: Mon Oct 21 09:57:49 2013
New Revision: 21195

URL: http://svn.gna.org/viewcvs/relax?rev=21195&view=rev
Log:
Caught a divide by zero in the 
specific_analyses.relax_disp.disp_data.return_offset_data() function.

This was identified by turning all numpy warnings to errors.


Modified:
    branches/relax_disp/specific_analyses/relax_disp/disp_data.py

Modified: branches/relax_disp/specific_analyses/relax_disp/disp_data.py
URL: 
http://svn.gna.org/viewcvs/relax/branches/relax_disp/specific_analyses/relax_disp/disp_data.py?rev=21195&r1=21194&r2=21195&view=diff
==============================================================================
--- branches/relax_disp/specific_analyses/relax_disp/disp_data.py (original)
+++ branches/relax_disp/specific_analyses/relax_disp/disp_data.py Mon Oct 21 
09:57:49 2013
@@ -1365,7 +1365,10 @@
             # Calculate the tilt angle.
             omega1 = point * 2.0 * pi
             Delta_omega = shifts[spin_index, frq_index] - 
offsets[spin_index, frq_index, disp_pt_index]
-            theta[spin_index, frq_index, disp_pt_index] = atan(omega1 / 
Delta_omega)
+            if Delta_omega == 0.0:
+                theta[spin_index, frq_index, disp_pt_index] = pi / 2.0
+            else:
+                theta[spin_index, frq_index, disp_pt_index] = atan(omega1 / 
Delta_omega)
 
     # Return the structures.
     return shifts, offsets, theta