mailr7604 - /1.3/generic_fns/dasha.py


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

Header


Content

Posted by edward on October 11, 2008 - 21:25:
Author: bugman
Date: Sat Oct 11 21:25:49 2008
New Revision: 7604

URL: http://svn.gna.org/viewcvs/relax?rev=7604&view=rev
Log:
Converted the extract() function to the new design.


Modified:
    1.3/generic_fns/dasha.py

Modified: 1.3/generic_fns/dasha.py
URL: 
http://svn.gna.org/viewcvs/relax/1.3/generic_fns/dasha.py?rev=7604&r1=7603&r2=7604&view=diff
==============================================================================
--- 1.3/generic_fns/dasha.py (original)
+++ 1.3/generic_fns/dasha.py Sat Oct 11 21:25:49 2008
@@ -32,7 +32,7 @@
 
 # relax module imports.
 from data import Relax_data_store; ds = Relax_data_store()
-from generic_fns import angles, diffusion_tensor, pipes
+from generic_fns import angles, diffusion_tensor, pipes, value
 from generic_fns.mol_res_spin import exists_mol_res_spin_data, 
first_residue_num, last_residue_num, residue_loop, spin_loop
 from relax_errors import RelaxDirError, RelaxError, RelaxFileError, 
RelaxNoPdbError, RelaxNoSequenceError, RelaxNoTensorError, RelaxNucleusError
 from relax_io import mkdir_nofail, open_write_file, test_binary
@@ -415,7 +415,11 @@
 
 
 def extract(dir):
-    """Function for extracting the Dasha results out of the 'dasha_results' 
file."""
+    """Extract the Dasha results out of the 'dasha_results' file.
+
+    @param dir:     The optional directory where the results file is located.
+    @type dir:      str or None
+    """
 
     # Test if sequence data is loaded.
     if not exists_mol_res_spin_data():
@@ -423,7 +427,7 @@
 
     # The directory.
     if dir == None:
-        dir = run
+        dir = pipes.cdp_name()
     if not access(dir, F_OK):
         raise RelaxDirError, ('Dasha', dir)
 
@@ -440,25 +444,25 @@
         if param in ['te', 'tf', 'ts']:
             scaling = 1e-9
         elif param == 'Rex':
-            scaling = 1.0 / (2.0 * pi * ds.frq[self.run][0]) ** 2
+            scaling = 1.0 / (2.0 * pi * cdp.frq[0]) ** 2
         else:
             scaling = 1.0
 
         # Set the values.
-        self.relax.generic.value.read(self.run, param=param, 
scaling=scaling, file=file_name, num_col=0, name_col=None, data_col=1, 
error_col=2)
+        value.read(param=param, scaling=scaling, file=file_name, 
res_num_col=0, res_name_col=None, data_col=1, error_col=2)
 
         # Clean up of non-existant parameters (set the parameter to None!).
-        for i in xrange(len(ds.res[self.run])):
-            # Skip deselected residues.
-            if not ds.res[self.run][i].select:
+        for spin in spin_loop():
+            # Skip deselected spins.
+            if not spin.select:
                 continue
 
             # Skip the residue (don't set the parameter to None) if the 
parameter exists in the model.
-            if param in ds.res[self.run][i].params:
+            if param in spin.params:
                 continue
 
             # Set the parameter to None.
-            setattr(ds.res[self.run][i], lower(param), None)
+            setattr(spin, lower(param), None)
 
     # Extract the chi-squared values.
     file_name = dir + '/chi2.out'
@@ -468,4 +472,4 @@
         raise RelaxFileError, ('Dasha', file_name)
 
     # Set the values.
-    self.relax.generic.value.read(self.run, param='chi2', file=file_name, 
num_col=0, name_col=None, data_col=1, error_col=2)
+    value.read(param='chi2', file=file_name, res_num_col=0, 
res_name_col=None, data_col=1, error_col=2)




Related Messages


Powered by MHonArc, Updated Sat Oct 11 22:20:03 2008