r18036 - /branches/frame_order_testing/specific_fns/frame_order.py -- November 27, 2012 - 15:41

Author: bugman
Date: Tue Nov 27 15:41:36 2012
New Revision: 18036

URL: http://svn.gna.org/viewcvs/relax?rev=18036&view=rev
Log:
The angles of the cone axis for the frame order eigenframe are now wrapped 
between 0 and 2pi.

This is needed for the numerical stability of the trigonometric functions, 
and for better
visualisation.


Modified:
    branches/frame_order_testing/specific_fns/frame_order.py

Modified: branches/frame_order_testing/specific_fns/frame_order.py
URL: 
http://svn.gna.org/viewcvs/relax/branches/frame_order_testing/specific_fns/frame_order.py?rev=18036&r1=18035&r2=18036&view=diff
==============================================================================
--- branches/frame_order_testing/specific_fns/frame_order.py (original)
+++ branches/frame_order_testing/specific_fns/frame_order.py Tue Nov 27 
15:41:36 2012
@@ -1451,9 +1451,9 @@
             if eigen_gamma != None:
                 cdp.eigen_gamma_sim[sim_index] = wrap_angles(eigen_gamma, 
cdp.eigen_gamma-pi, cdp.eigen_gamma+pi)
             if axis_theta != None:
-                cdp.axis_theta_sim[sim_index] = axis_theta
+                cdp.axis_theta_sim[sim_index] = wrap_angles(axis_theta, 
cdp.axis_theta-pi, cdp.axis_theta+pi)
             if axis_phi != None:
-                cdp.axis_phi_sim[sim_index] = axis_phi
+                cdp.axis_phi_sim[sim_index] = wrap_angles(axis_phi, 
cdp.axis_phi-pi, cdp.axis_phi+pi)
 
             # Cone parameters.
             if cone_theta != None:
@@ -1494,9 +1494,9 @@
             if eigen_gamma != None:
                 cdp.eigen_gamma = wrap_angles(eigen_gamma, 0.0, 2.0*pi)
             if axis_theta != None:
-                cdp.axis_theta = axis_theta
+                cdp.axis_theta = wrap_angles(axis_theta, 0.0, 2.0*pi)
             if axis_phi != None:
-                cdp.axis_phi = axis_phi
+                cdp.axis_phi = wrap_angles(axis_phi, 0.0, 2.0*pi)
 
             # Cone parameters.
             if cone_theta != None: