r22170 - /branches/double_rotor/test_suite/shared_data/frame_order/cam/double_rotor/eigen_system.py -- February 12, 2014 - 09:13

Author: bugman
Date: Wed Feb 12 09:13:08 2014
New Revision: 22170

URL: http://svn.gna.org/viewcvs/relax?rev=22170&view=rev
Log:
Added a simple script for analysing the Eigensystem of the CaM frame order 
double rotor test model.


Added:
    
branches/double_rotor/test_suite/shared_data/frame_order/cam/double_rotor/eigen_system.py

Added: 
branches/double_rotor/test_suite/shared_data/frame_order/cam/double_rotor/eigen_system.py
URL: 
http://svn.gna.org/viewcvs/relax/branches/double_rotor/test_suite/shared_data/frame_order/cam/double_rotor/eigen_system.py?rev=22170&view=auto
==============================================================================
--- 
branches/double_rotor/test_suite/shared_data/frame_order/cam/double_rotor/eigen_system.py
 (added)
+++ 
branches/double_rotor/test_suite/shared_data/frame_order/cam/double_rotor/eigen_system.py
 Wed Feb 12 09:13:08 2014
@@ -1,0 +1,33 @@
+# Script for determining the full eigensystem of the double rotor system.
+
+# Python module imports
+from numpy import array, cross, dot, float64, transpose
+from numpy.linalg import inv, norm
+
+
+# The axes of the double rotor.
+axis1 = array([ -7.778375610280605e-01,   6.284649244351433e-01, 
-7.532653237683726e-04], float64)
+axis2 = array([-0.487095774865268, -0.60362450312215 , -0.63116968030708 ], 
float64)
+
+# The cross product.
+axis3 = cross(axis1, axis2)
+axis3 = axis3 / norm(axis3)
+
+# Printout.
+print("\n")
+print("Axis 1:       %s" % axis1)
+print("Axis 2:       %s" % axis2)
+print("Axis 3:       %s" % axis3)
+print("Axis 1 norm:  %s" % norm(axis1))
+print("Axis 2 norm:  %s" % norm(axis2))
+print("Axis 3 norm:  %s" % norm(axis3))
+
+# Eigenframe.
+frame = transpose(array([axis1, axis2, axis3]))
+print("Eigenframe:\n%s" % frame)
+
+# Dot products.
+print("Dot 1 to 2:   %s" % dot(axis1, axis2))
+print("Dot 1 to 3:   %s" % dot(axis1, axis3))
+print("Dot 2 to 3:   %s" % dot(axis2, axis3))
+print inv(frame) - transpose(frame)