Author: bugman Date: Wed Aug 4 14:30:13 2010 New Revision: 11402 URL: http://svn.gna.org/viewcvs/relax?rev=11402&view=rev Log: Removal of the reduce_alignment_tensor_reverse() frame order function. This is not being used, and is of no use! 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=11402&r1=11401&r2=11402&view=diff ============================================================================== --- 1.3/maths_fns/frame_order_matrix_ops.py (original) +++ 1.3/maths_fns/frame_order_matrix_ops.py Wed Aug 4 14:30:13 2010 @@ -1227,56 +1227,6 @@ red_tensor[4] = red_tensor[4] + (D[4, 8] + D[5, 7])*A[4] -def reduce_alignment_tensor_reverse(D, A, red_tensor): - """Calculate the reduction in the alignment tensor caused by the Frame Order matrix. - - This is the reverse rotation notation. - - @param D: The Frame Order matrix, 2nd degree to be populated. - @type D: numpy 9D, rank-2 array - @param A: The full alignment tensor in {Axx, Ayy, Axy, Axz, Ayz} notation. - @type A: numpy 5D, rank-1 array - @param red_tensor: The structure in {Axx, Ayy, Axy, Axz, Ayz} notation to place the reduced - alignment tensor. - @type red_tensor: numpy 5D, rank-1 array - """ - - # The reduced tensor element A0. - red_tensor[0] = (D[0, 0] - D[6, 6])*A[0] - red_tensor[0] = red_tensor[0] + (D[3, 3] - D[6, 6])*A[1] - red_tensor[0] = red_tensor[0] + 2.0*D[0, 3]*A[2] - red_tensor[0] = red_tensor[0] + 2.0*D[0, 6]*A[3] - red_tensor[0] = red_tensor[0] + 2.0*D[3, 6]*A[4] - - # The reduced tensor element A1. - red_tensor[1] = (D[1, 1] - D[7, 7])*A[0] - red_tensor[1] = red_tensor[1] + (D[4, 4] - D[7, 7])*A[1] - red_tensor[1] = red_tensor[1] + 2.0*D[1, 4]*A[2] - red_tensor[1] = red_tensor[1] + 2.0*D[1, 7]*A[3] - red_tensor[1] = red_tensor[1] + 2.0*D[4, 7]*A[4] - - # The reduced tensor element A2. - red_tensor[2] = (D[0, 1] - D[6, 7])*A[0] - red_tensor[2] = red_tensor[2] + (D[3, 4] - D[6, 7])*A[1] - red_tensor[2] = red_tensor[2] + (D[0, 4] + D[1, 3])*A[2] - red_tensor[2] = red_tensor[2] + (D[0, 7] + D[1, 6])*A[3] - red_tensor[2] = red_tensor[2] + (D[3, 7] + D[4, 6])*A[4] - - # The reduced tensor element A3. - red_tensor[3] = (D[0, 2] - D[6, 8])*A[0] - red_tensor[3] = red_tensor[3] + (D[3, 5] - D[6, 8])*A[1] - red_tensor[3] = red_tensor[3] + (D[0, 5] + D[2, 3])*A[2] - red_tensor[3] = red_tensor[3] + (D[0, 8] + D[2, 6])*A[3] - red_tensor[3] = red_tensor[3] + (D[3, 8] + D[5, 6])*A[4] - - # The reduced tensor element A4. - red_tensor[4] = (D[1, 2] - D[7, 8])*A[0] - red_tensor[4] = red_tensor[4] + (D[4, 5] - D[7, 8])*A[1] - red_tensor[4] = red_tensor[4] + (D[1, 5] + D[2, 4])*A[2] - red_tensor[4] = red_tensor[4] + (D[1, 8] + D[2, 7])*A[3] - red_tensor[4] = red_tensor[4] + (D[4, 8] + D[5, 7])*A[4] - - def rotate_daeg(matrix, R): """Rotate the given frame order matrix.