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.
_______________________________________________
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
