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]