Author: bugman Date: Wed Sep 10 09:22:20 2014 New Revision: 25727 URL: http://svn.gna.org/viewcvs/relax?rev=25727&view=rev Log: Added extensive printouts to the frame_order.permute_axes user function. Modified: branches/frame_order_cleanup/specific_analyses/frame_order/uf.py Modified: branches/frame_order_cleanup/specific_analyses/frame_order/uf.py URL: http://svn.gna.org/viewcvs/relax/branches/frame_order_cleanup/specific_analyses/frame_order/uf.py?rev=25727&r1=25726&r2=25727&view=diff ============================================================================== --- branches/frame_order_cleanup/specific_analyses/frame_order/uf.py (original) +++ branches/frame_order_cleanup/specific_analyses/frame_order/uf.py Wed Sep 10 09:22:20 2014 @@ -122,6 +122,16 @@ frame = zeros((3, 3), float64) euler_to_R_zyz(cdp.eigen_alpha, cdp.eigen_beta, cdp.eigen_gamma, frame) + # Start printout. + print("\nOriginal parameters:") + print("%-20s %20.10f" % ("eigen_alpha", cdp.eigen_alpha)) + print("%-20s %20.10f" % ("eigen_beta", cdp.eigen_beta)) + print("%-20s %20.10f" % ("eigen_gamma", cdp.eigen_gamma)) + print("%-20s\n%s" % ("eigenframe", frame)) + print("%-20s %20.10f" % ("cone_theta_x", cdp.cone_theta_x)) + print("%-20s %20.10f" % ("cone_theta_y", cdp.cone_theta_y)) + print("%-20s %20.10f" % ("cone_sigma_max", cone_sigma_max)) + # The permutation with the condition that cone_theta_x <= cone_theta_y. if angles[1] <= angles[2]: perm = [2, 0, 1] @@ -139,6 +149,17 @@ # Convert the permuted frame to Euler angles and store them. cdp.eigen_alpha, cdp.eigen_beta, cdp.eigen_gamma = R_to_euler_zyz(frame_new) + + # End printout. + print("\nPermuted parameters:") + print("%-20s %20s" % ("permutation", perm)) + print("%-20s %20.10f" % ("eigen_alpha", cdp.eigen_alpha)) + print("%-20s %20.10f" % ("eigen_beta", cdp.eigen_beta)) + print("%-20s %20.10f" % ("eigen_gamma", cdp.eigen_gamma)) + print("%-20s\n%s" % ("eigenframe", frame_new)) + print("%-20s %20.10f" % ("cone_theta_x", cdp.cone_theta_x)) + print("%-20s %20.10f" % ("cone_theta_y", cdp.cone_theta_y)) + print("%-20s %20.10f" % ("cone_sigma_max", cdp.cone_sigma_max)) def pivot(pivot=None, order=1, fix=False):