r18461 - /branches/frame_order_testing/specific_fns/frame_order.py -- February 13, 2013 - 15:45

Author: bugman
Date: Wed Feb 13 15:45:08 2013
New Revision: 18461

URL: http://svn.gna.org/viewcvs/relax?rev=18461&view=rev
Log:
Bug fix for the RDCs used in the frame order analysis when multiple models 
are present.

The unit vectors used in the RDC calculation were not of length 1!


Modified:
    branches/frame_order_testing/specific_fns/frame_order.py

Modified: branches/frame_order_testing/specific_fns/frame_order.py
URL: 
http://svn.gna.org/viewcvs/relax/branches/frame_order_testing/specific_fns/frame_order.py?rev=18461&r1=18460&r2=18461&view=diff
==============================================================================
--- branches/frame_order_testing/specific_fns/frame_order.py (original)
+++ branches/frame_order_testing/specific_fns/frame_order.py Wed Feb 13 
15:45:08 2013
@@ -28,7 +28,7 @@
 from minfx.generic import generic_minimise
 from minfx.grid import grid_point_array
 from numpy import arccos, array, dot, eye, float64, identity, ones, 
transpose, zeros
-from numpy.linalg import inv
+from numpy.linalg import inv, norm
 from re import search
 import sys
 from warnings import warn
@@ -830,6 +830,9 @@
 
                 # Store.
                 unit_vect.append(ave_vector)
+
+            # Normalise (to be sure).
+            unit_vect[-1] = unit_vect[-1] / norm(unit_vect[-1])
 
             # Gyromagnetic ratios.
             g1 = return_gyromagnetic_ratio(spin1.isotope)