mailr7183 - /branches/rdc_analysis/specific_fns/n_state_model.py


Others Months | Index by Date | Thread Index
>>   [Date Prev] [Date Next] [Thread Prev] [Thread Next]

Header


Content

Posted by edward on August 12, 2008 - 16:47:
Author: bugman
Date: Tue Aug 12 15:33:09 2008
New Revision: 7183

URL: http://svn.gna.org/viewcvs/relax?rev=7183&view=rev
Log:
Modified the N-state model methods to set the alignment tensor parameter 
rather Saupe order matrix.


Modified:
    branches/rdc_analysis/specific_fns/n_state_model.py

Modified: branches/rdc_analysis/specific_fns/n_state_model.py
URL: 
http://svn.gna.org/viewcvs/relax/branches/rdc_analysis/specific_fns/n_state_model.py?rev=7183&r1=7182&r2=7183&view=diff
==============================================================================
--- branches/rdc_analysis/specific_fns/n_state_model.py (original)
+++ branches/rdc_analysis/specific_fns/n_state_model.py Tue Aug 12 15:33:09 
2008
@@ -222,11 +222,11 @@
         if 'rdc' in data_types or 'pcs' in data_types:
             # Loop over the alignments, adding the alignment tensor 
parameters to the tensor data container.
             for i in xrange(len(cdp.align_tensors)):
-                cdp.align_tensors[i].Sxx = param_vector[5*i]
-                cdp.align_tensors[i].Syy = param_vector[5*i+1]
-                cdp.align_tensors[i].Sxy = param_vector[5*i+2]
-                cdp.align_tensors[i].Sxz = param_vector[5*i+3]
-                cdp.align_tensors[i].Syz = param_vector[5*i+4]
+                cdp.align_tensors[i].Axx = param_vector[5*i]
+                cdp.align_tensors[i].Ayy = param_vector[5*i+1]
+                cdp.align_tensors[i].Axy = param_vector[5*i+2]
+                cdp.align_tensors[i].Axz = param_vector[5*i+3]
+                cdp.align_tensors[i].Ayz = param_vector[5*i+4]
 
             # Create a new parameter vector without the tensors.
             param_vector = param_vector[5*len(cdp.align_tensors):]
@@ -593,21 +593,21 @@
             # Create a list of all the reduced alignment tensor elements and 
their errors (for the chi-squared function).
             elif tensor.red:
                 # Append the 5 unique elements.
-                red_tensor_elem.append(tensor.Sxx)
-                red_tensor_elem.append(tensor.Syy)
-                red_tensor_elem.append(tensor.Sxy)
-                red_tensor_elem.append(tensor.Sxz)
-                red_tensor_elem.append(tensor.Syz)
+                red_tensor_elem.append(tensor.Axx)
+                red_tensor_elem.append(tensor.Ayy)
+                red_tensor_elem.append(tensor.Axy)
+                red_tensor_elem.append(tensor.Axz)
+                red_tensor_elem.append(tensor.Ayz)
 
                 # Append the 5 unique error elements (if they exist).
-                if hasattr(tensor, 'Sxx_err'):
-                    red_tensor_err.append(tensor.Sxx_err)
-                    red_tensor_err.append(tensor.Syy_err)
-                    red_tensor_err.append(tensor.Sxy_err)
-                    red_tensor_err.append(tensor.Sxz_err)
-                    red_tensor_err.append(tensor.Syz_err)
-
-                # Otherwise append errors of 1.0 to convert the chi-squared 
equation to the SSE equation (for the tensors without errors).
+                if hasattr(tensor, 'Axx_err'):
+                    red_tensor_err.append(tensor.Axx_err)
+                    red_tensor_err.append(tensor.Ayy_err)
+                    red_tensor_err.append(tensor.Axy_err)
+                    red_tensor_err.append(tensor.Axz_err)
+                    red_tensor_err.append(tensor.Ayz_err)
+
+                # Otherwise append errors of 1.0 to convert the chi-squared 
equation to the ASE equation (for the tensors without errors).
                 else:
                     red_tensor_err = red_tensor_err + [1.0, 1.0, 1.0, 1.0, 
1.0]
 




Related Messages


Powered by MHonArc, Updated Tue Aug 12 17:20:18 2008