r17291 - /branches/frame_order_testing/maths_fns/frame_order.py -- July 20, 2012 - 16:58

Author: bugman
Date: Fri Jul 20 16:58:06 2012
New Revision: 17291

URL: http://svn.gna.org/viewcvs/relax?rev=17291&view=rev
Log:
Changed the structures of r17286 to be of numpy.float32 types instead of 
float16.

This is to speed up the calculation due to faster type conversions.


Modified:
    branches/frame_order_testing/maths_fns/frame_order.py

Modified: branches/frame_order_testing/maths_fns/frame_order.py
URL: 
http://svn.gna.org/viewcvs/relax/branches/frame_order_testing/maths_fns/frame_order.py?rev=17291&r1=17290&r2=17291&view=diff
==============================================================================
--- branches/frame_order_testing/maths_fns/frame_order.py (original)
+++ branches/frame_order_testing/maths_fns/frame_order.py Fri Jul 20 16:58:06 
2012
@@ -258,13 +258,13 @@
             self.paramag_unit_vect = zeros(atomic_pos.shape, float64)
             self.paramag_dist = zeros(self.num_spins, float64)
             self.pcs_const = zeros(self.num_align, float64)
-            self.r_pivot_atom = zeros((3, self.num_spins), float16)
-            self.r_pivot_atom_rev = zeros((3, self.num_spins), float16)
-            self.r_ln_pivot = zeros((3, self.num_spins), float16)
+            self.r_pivot_atom = zeros((3, self.num_spins), float32)
+            self.r_pivot_atom_rev = zeros((3, self.num_spins), float32)
+            self.r_ln_pivot = zeros((3, self.num_spins), float32)
             for j in xrange(self.num_spins):
                 self.r_ln_pivot[:, j] = pivot - self.paramag_centre
             if self.paramag_centre == None:
-                self.paramag_centre = zeros(3, float16)
+                self.paramag_centre = zeros(3, float32)
 
             # Set up the paramagnetic constant (without the interatomic 
distance and in Angstrom units).
             for align_index in range(self.num_align):