r26460 - /branches/frame_order_cleanup/lib/frame_order/pseudo_ellipse_torsionless.py -- November 09, 2014 - 11:25

Author: bugman
Date: Sun Nov  9 11:25:18 2014
New Revision: 26460

URL: http://svn.gna.org/viewcvs/relax?rev=26460&view=rev
Log:
Bug fix for the torsionless pseudo-ellipse 1st degree frame order matrix.

The 11 element was of the wrong sign.


Modified:
    branches/frame_order_cleanup/lib/frame_order/pseudo_ellipse_torsionless.py

Modified: 
branches/frame_order_cleanup/lib/frame_order/pseudo_ellipse_torsionless.py
URL: 
http://svn.gna.org/viewcvs/relax/branches/frame_order_cleanup/lib/frame_order/pseudo_ellipse_torsionless.py?rev=26460&r1=26459&r2=26460&view=diff
==============================================================================
--- 
branches/frame_order_cleanup/lib/frame_order/pseudo_ellipse_torsionless.py  
(original)
+++ 
branches/frame_order_cleanup/lib/frame_order/pseudo_ellipse_torsionless.py  
Sun Nov  9 11:25:18 2014
@@ -60,7 +60,7 @@
 
     # Numerical integration of phi of each element.
     matrix[0, 0] = fact * (2.0*pi + quad(part_int_daeg1_pseudo_ellipse_00, 
-pi, pi, args=(theta_x, theta_y), full_output=1)[0])
-    matrix[1, 1] = -fact * (2.0*pi + quad(part_int_daeg1_pseudo_ellipse_11, 
-pi, pi, args=(theta_x, theta_y), full_output=1)[0])
+    matrix[1, 1] = fact * (2.0*pi + quad(part_int_daeg1_pseudo_ellipse_11, 
-pi, pi, args=(theta_x, theta_y), full_output=1)[0])
     matrix[2, 2] = fact * quad(part_int_daeg1_pseudo_ellipse_22, -pi, pi, 
args=(theta_x, theta_y), full_output=1)[0]
 
     # Rotate and return the frame order matrix.