mailr23646 - /branches/frame_order_cleanup/test_suite/shared_data/frame_order/cam/generate_base.py


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

Header


Content

Posted by edward on June 05, 2014 - 10:50:
Author: bugman
Date: Thu Jun  5 10:50:12 2014
New Revision: 23646

URL: http://svn.gna.org/viewcvs/relax?rev=23646&view=rev
Log:
Updated the CaM frame order data generation base script to print out more 
information.

This is for the first axis system so that the same amount of information as 
the second system is
printed.


Modified:
    
branches/frame_order_cleanup/test_suite/shared_data/frame_order/cam/generate_base.py

Modified: 
branches/frame_order_cleanup/test_suite/shared_data/frame_order/cam/generate_base.py
URL: 
http://svn.gna.org/viewcvs/relax/branches/frame_order_cleanup/test_suite/shared_data/frame_order/cam/generate_base.py?rev=23646&r1=23645&r2=23646&view=diff
==============================================================================
--- 
branches/frame_order_cleanup/test_suite/shared_data/frame_order/cam/generate_base.py
        (original)
+++ 
branches/frame_order_cleanup/test_suite/shared_data/frame_order/cam/generate_base.py
        Thu Jun  5 10:50:12 2014
@@ -602,7 +602,7 @@
         self.axes = dot(R, axes)
         alpha, beta, gamma = R_to_euler_zyz(self.axes)
 
-        # Print out.
+        # Printout.
         print("Tilt axis: %s, norm = %s" % (repr(tilt_axis), 
norm(tilt_axis)))
         print("CoM-pivot axis: %s, norm = %s" % (repr(axis_z), norm(axis_z)))
         print("Rotation axis: %s, norm = %s" % (repr(self.axes[:, 2]), 
norm(self.axes[:, 2])))
@@ -627,3 +627,10 @@
 
         # The eigenframe.
         self.axes = transpose(array([axis_x, axis_y, axis_z]))
+        alpha, beta, gamma = R_to_euler_zyz(self.axes)
+
+        # Printout.
+        print("CoM-pivot axis: %s, norm = %s" % (repr(axis_z), norm(axis_z)))
+        print("Rotation axis: %s, norm = %s" % (repr(self.axes[:, 2]), 
norm(self.axes[:, 2])))
+        print("Full axis system:\n%s" % self.axes)
+        print("Full axis system Euler angles:\n\talpha: %s\n\tbeta: 
%s\n\tgamma: %s" % (repr(alpha), repr(beta), repr(gamma)))




Related Messages


Powered by MHonArc, Updated Thu Jun 05 11:20:02 2014