r26474 - /branches/frame_order_cleanup/lib/frame_order/pseudo_ellipse_torsionless.py -- November 09, 2014 - 22:23

Author: bugman
Date: Sun Nov  9 22:23:54 2014
New Revision: 26474

URL: http://svn.gna.org/viewcvs/relax?rev=26474&view=rev
Log:
Fix for r26460 - the second factor also needs to be inverted.


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=26474&r1=26473&r2=26474&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 22:23:54 2014
@@ -97,8 +97,10 @@
     # Invert.
     if fact == 0.0:
         fact = 1e100
+        fact2 = 1e100
     else:
         fact = 1.0 / fact
+        fact2 = 1.0 / fact2
 
     # Diagonal.
     matrix[0, 0] = fact2 * (3.0*pi + 
quad(part_int_daeg2_pseudo_ellipse_torsionless_00, -pi, pi, args=(theta_x, 
theta_y), full_output=1)[0])