Author: bugman Date: Wed Aug 4 17:17:01 2010 New Revision: 11406 URL: http://svn.gna.org/viewcvs/relax?rev=11406&view=rev Log: Modified rotate_daeg() to operate on the Kronecker product form of the frame order matrix. Modified: 1.3/maths_fns/frame_order_matrix_ops.py Modified: 1.3/maths_fns/frame_order_matrix_ops.py URL: http://svn.gna.org/viewcvs/relax/1.3/maths_fns/frame_order_matrix_ops.py?rev=11406&r1=11405&r2=11406&view=diff ============================================================================== --- 1.3/maths_fns/frame_order_matrix_ops.py (original) +++ 1.3/maths_fns/frame_order_matrix_ops.py Wed Aug 4 17:17:01 2010 @@ -1252,6 +1252,9 @@ def rotate_daeg(matrix, R): """Rotate the given frame order matrix. + It is assumed that the frame order matrix is in the Kronecker product form. + + @param matrix: The Frame Order matrix, 2nd degree to be populated. @type matrix: numpy 9D, rank-2 array @param R: The rotation matrix to be populated. @@ -1261,14 +1264,8 @@ # The outer product of R. R_kron = kron_prod(R, R) - # Perform the T23 transpose to obtain the Kronecker product matrix! - transpose_23(matrix) - # Rotate. matrix_rot = dot(R_kron, dot(matrix, transpose(R_kron))) - - # Perform T23 again to return back. - transpose_23(matrix_rot) # Return the matrix. return matrix_rot