mailr24332 - /branches/frame_order_cleanup/lib/frame_order/double_rotor.py


Others Months | Index by Date | Thread Index
>>   [Date Prev] [Date Next] [Thread Prev] [Thread Next]

Header


Content

Posted by edward on June 26, 2014 - 16:29:
Author: bugman
Date: Thu Jun 26 16:29:46 2014
New Revision: 24332

URL: http://svn.gna.org/viewcvs/relax?rev=24332&view=rev
Log:
Fix for the double rotor frame order matrix assembly function.

The rotations where applied in reverse, hence the 2nd degree frame order 
matrix was the Kronecker
product transpose of the correct form.


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

Modified: branches/frame_order_cleanup/lib/frame_order/double_rotor.py
URL: 
http://svn.gna.org/viewcvs/relax/branches/frame_order_cleanup/lib/frame_order/double_rotor.py?rev=24332&r1=24331&r2=24332&view=diff
==============================================================================
--- branches/frame_order_cleanup/lib/frame_order/double_rotor.py        
(original)
+++ branches/frame_order_cleanup/lib/frame_order/double_rotor.py        Thu 
Jun 26 16:29:46 2014
@@ -72,11 +72,11 @@
     matrix[8, 8] = 0.5 * sinc_2smax1p1 * sinc_2smax2p1
 
     # Off diagonal set 1.
-    matrix[4, 0] = 0.5 * sinc_2smax1n1 * sinc_2smax2n1
-    matrix[0, 8] = -sinc_2smax1n1
-    matrix[8, 0] = -0.5 * sinc_2smax1n1 * sinc_2smax2p1
-    matrix[4, 8] = -0.5 * sinc_2smax1p1 * sinc_2smax2n1
-    matrix[8, 4] = -sinc_2smax2n1
+    matrix[0, 4] = 0.5 * sinc_2smax1n1 * sinc_2smax2n1
+    matrix[0, 8] = -0.5 * sinc_2smax1n1 * sinc_2smax2p1
+    matrix[8, 0] = -sinc_2smax1n1
+    matrix[4, 8] = -sinc_2smax2n1
+    matrix[8, 4] = -0.5 * sinc_2smax1p1 * sinc_2smax2n1
 
     # Off diagonal set 2.
     matrix[2, 6] = matrix[6, 2] = sinc_smax2 * sinc_2smax1n1




Related Messages


Powered by MHonArc, Updated Thu Jun 26 16:40:02 2014