r10110 - /1.3/maths_fns/rotation_matrix.py -- January 05, 2010 - 19:14

Author: bugman
Date: Tue Jan  5 19:14:09 2010
New Revision: 10110

URL: http://svn.gna.org/viewcvs/relax?rev=10110&view=rev
Log:
Bug fix for the euler_xyx_to_R() function.

The R[y, x] element was incorrect!!!


Modified:
    1.3/maths_fns/rotation_matrix.py

Modified: 1.3/maths_fns/rotation_matrix.py
URL: 
http://svn.gna.org/viewcvs/relax/1.3/maths_fns/rotation_matrix.py?rev=10110&r1=10109&r2=10110&view=diff
==============================================================================
--- 1.3/maths_fns/rotation_matrix.py (original)
+++ 1.3/maths_fns/rotation_matrix.py Tue Jan  5 19:14:09 2010
@@ -458,7 +458,7 @@
     According to wikipedia 
(http://en.wikipedia.org/wiki/Euler_angles#Table_of_matrices), the rotation 
matrix for the xyx convention is::
 
               |  cb                  sa*sb               ca*sb            |
-        R  =  |  sa*sb               ca*cg - sa*cb*sg   -sa*cg - ca*cb*sg |,
+        R  =  |  sb*sg               ca*cg - sa*cb*sg   -sa*cg - ca*cb*sg |,
               | -sb*cg               ca*sg + sa*cb*cg   -sa*sg + ca*cb*cg |
 
     where::
@@ -491,7 +491,7 @@
 
     # The unit mux vector component of the rotation matrix.
     R[0, 0] =  cos_b
-    R[1, 0] =  sin_a * sin_b
+    R[1, 0] =  sin_b * sin_g
     R[2, 0] = -sin_b * cos_g
 
     # The unit muy vector component of the rotation matrix.