Author: bugman
Date: Wed Sep 10 21:33:37 2014
New Revision: 25738

URL: http://svn.gna.org/viewcvs/relax?rev=25738&view=rev
Log:
Bug fix for the new frame_order.permute_axes user function.

The cone and torsion angles were not being correctly permuted.  Now the direct 
permutation array is
being used.  And the fact that cone_theta_x is a rotation along the y-axis and 
cone_theta_y along
the x-axis is taken into account.


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=25738&r1=25737&r2=25738&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 21:33:37 2014
@@ -188,11 +188,11 @@
     print("%-20s %-20s" % ("permutation", perm))
     print("%-20s %-20s" % ("z-axis inversion", inv))
 
-    # Permute the angles.
-    cdp.cone_theta_x = angles[perm.index(0)]
-    cdp.cone_theta_y = angles[perm.index(1)]
+    # Permute the angles (note that cone_theta_x is a rotation about the 
y-axis and cone_theta_y about the x-axis).
+    cdp.cone_theta_x = angles[perm[1]]
+    cdp.cone_theta_y = angles[perm[0]]
     if cdp.model == MODEL_PSEUDO_ELLIPSE:
-        cdp.cone_sigma_max = angles[perm.index(2)]
+        cdp.cone_sigma_max = angles[perm[2]]
 
     # Permute the axes and invert the z-axis as necessary.
     frame_new = transpose(array([inv[0]*frame[:, perm[0]], inv[1]*frame[:, 
perm[1]], inv[2]*frame[:, perm[2]]], float64))


_______________________________________________
relax (http://www.nmr-relax.com)

This is the relax-commits mailing list
[email protected]

To unsubscribe from this list, get a password
reminder, or change your subscription options,
visit the list information page at
https://mail.gna.org/listinfo/relax-commits

Reply via email to