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])
_______________________________________________
relax (http://www.nmr-relax.com)
This is the relax-commits mailing list
[email protected]
To unsubscribe from this list, get a password
reminder, or change your subscription options,
visit the list information page at
https://mail.gna.org/listinfo/relax-commits