r14361 - /1.3/data/align_tensor.py -- August 11, 2011 - 16:38

Author: bugman
Date: Thu Aug 11 16:38:36 2011
New Revision: 14361

URL: http://svn.gna.org/viewcvs/relax?rev=14361&view=rev
Log:
Bug fix for the data.align_tensor.calc_eta() function - divide by zeros are 
now avoided.

When Azz is zero, eta is now set to NaN.


Modified:
    1.3/data/align_tensor.py

Modified: 1.3/data/align_tensor.py
URL: 
http://svn.gna.org/viewcvs/relax/1.3/data/align_tensor.py?rev=14361&r1=14360&r2=14361&view=diff
==============================================================================
--- 1.3/data/align_tensor.py (original)
+++ 1.3/data/align_tensor.py Thu Aug 11 16:38:36 2011
@@ -1,6 +1,6 @@
 
###############################################################################
 #                                                                            
 #
-# Copyright (C) 2003-2004, 2006-2010 Edward d'Auvergne                       
 #
+# Copyright (C) 2003-2011 Edward d'Auvergne                                  
 #
 #                                                                            
 #
 # This file is part of the program relax.                                    
 #
 #                                                                            
 #
@@ -29,6 +29,7 @@
 
 # relax module imports.
 from data_classes import Element
+from float import nan
 from maths_fns.rotation_matrix import R_to_euler_zyz
 from relax_errors import RelaxError
 from relax_xml import fill_object_contents, xml_to_object
@@ -265,6 +266,10 @@
     @rtype:         float
     """
 
+    # Zero Azz value, so return NaN.
+    if A_diag[2, 2] == 0:
+        return nan
+
     # Return eta.
     return (A_diag[0, 0] - A_diag[1, 1]) / A_diag[2, 2]