r4777 - /branches/N_state_model/maths_fns/n_state_model.py -- January 16, 2008 - 11:27

Author: bugman
Date: Wed Jan 16 11:27:02 2008
New Revision: 4777

URL: http://svn.gna.org/viewcvs/relax?rev=4777&view=rev
Log:
The vector of back-calculated reduced alignment tensors is now created.


Modified:
    branches/N_state_model/maths_fns/n_state_model.py

Modified: branches/N_state_model/maths_fns/n_state_model.py
URL: 
http://svn.gna.org/viewcvs/relax/branches/N_state_model/maths_fns/n_state_model.py?rev=4777&r1=4776&r2=4777&view=diff
==============================================================================
--- branches/N_state_model/maths_fns/n_state_model.py (original)
+++ branches/N_state_model/maths_fns/n_state_model.py Wed Jan 16 11:27:02 2008
@@ -110,5 +110,13 @@
             for i in xrange(self.num_tensors):
                 self.red_bc[i] = self.red_bc[i]  +  pc * self.RT[c] * 
self.full_tensors[i] * self.R[c]
 
+        # 5D vectorise the back-calculated tensors (create red_bc_vector 
from red_bc).
+        for i in xrange(self.num_tensors):
+            self.red_bc_vector[5*i]   = self.red_bc[i,0,0]    # Sxx.
+            self.red_bc_vector[5*i+1] = self.red_bc[i,1,1]    # Syy.
+            self.red_bc_vector[5*i+2] = self.red_bc[i,0,1]    # Sxy.
+            self.red_bc_vector[5*i+3] = self.red_bc[i,0,2]    # Sxz.
+            self.red_bc_vector[5*i+4] = self.red_bc[i,1,2]    # Syz.
+
         # Return the chi-squared value.
         return chi2(self.red_data, self.red_bc_vector, self.red_errors)