mailr8117 - /1.3/generic_fns/spectrum.py


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

Header


Content

Posted by sebastien . morin . 1 on December 04, 2008 - 02:09:
Author: semor
Date: Thu Dec  4 02:09:50 2008
New Revision: 8117

URL: http://svn.gna.org/viewcvs/relax?rev=8117&view=rev
Log:
Reverted revisions 8096 and 8097.

The command used was:
svn merge -r8097:8095 .

This follows a discussion started by Edward d'Auvergne at:
https://mail.gna.org/public/relax-users/2008-12/msg00000.html
(Message ID: 7f080ed10812030146w7e73ce24r75d32121e5f1c05f@xxxxxxxxxxxxxx)


Modified:
    1.3/generic_fns/spectrum.py

Modified: 1.3/generic_fns/spectrum.py
URL: 
http://svn.gna.org/viewcvs/relax/1.3/generic_fns/spectrum.py?rev=8117&r1=8116&r2=8117&view=diff
==============================================================================
--- 1.3/generic_fns/spectrum.py (original)
+++ 1.3/generic_fns/spectrum.py Thu Dec  4 02:09:50 2008
@@ -38,7 +38,6 @@
 from relax_errors import RelaxError, RelaxArgNotInListError, 
RelaxImplementError, RelaxNoSequenceError
 from relax_io import extract_data, strip
 from relax_warnings import RelaxWarning, RelaxNoSpinWarning
-from specific_fns.relax_fit import Relax_fit
 
 
 def __errors_height_no_repl():
@@ -411,7 +410,7 @@
             __errors_repl()
 
 
-def intensity_generic(line, int_col, file_data_header):
+def intensity_generic(line, int_col):
     """Function for returning relevant data from the generic peak intensity 
line.
 
     The residue number, heteronucleus and proton names, and peak intensity 
will be returned.
@@ -425,46 +424,8 @@
     @raises RelaxError: When the expected peak intensity is not a float.
     """
 
-    # Determine the number of delays (and associated intensities).
-    i = 5
-    while 1:
-        i = i + 1
-        try:
-            current_field = line[i-1]
-        except:
-            num_delays = int(i - 6)
-            if num_delays == 0:
-                raise RelaxError, "Generic file with no associated delays 
(and intensities)."
-            break
-
-    # The residue number.
-    res_num = ''
-    try:
-        res_num = int(line[1])
-    except:
-        raise RelaxError, "Improperly formatted generic file."
-
-    # Nuclei names.
-    x_name = ''
-    x_name = line[4]
-    #h_name = ''
-    #h_name = line[5]
-
-    # Extract both delays and associated intensities.
-    i = 0
-    while i < num_delays:
-        i = i + 1
-        delay = file_data_header[0][i + 4]
-        intensity = line[i + 4]
-        Relax_fit.relax_time(time=file_data_header[0][i + 6], 
spectrum_id=line[i + 6])
-
-    ## Intensity.
-    #try:
-    #    intensity = float(line[int_col])
-    #except ValueError:
-    #    raise RelaxError, "The peak intensity value " + `intensity` + " 
from the line " + `line` + " is invalid."
-    ## Return the data.
-    #return res_num, h_name, x_name, intensity
+
+    # Not implemented yet...
 
 
 def intensity_nmrview(line, int_col):
@@ -782,11 +743,6 @@
     num = number_of_header_lines(file_data, format, int_col, intensity_fn)
     print "Number of header lines found: " + `num`
 
-    # Store the header (if using the generic file type).
-    if format == 'generic':
-        file_data_header = file_data[:num]
-        file_data_header = strip(file_data_header)
-
     # Remove the header.
     file_data = file_data[num:]
 
@@ -808,10 +764,7 @@
     # Loop over the peak intensity data.
     for i in xrange(len(file_data)):
         # Extract the data.
-        if format == 'generic':
-            res_num, H_name, X_name, intensity = intensity_fn(file_data[i], 
int_col, file_data_header)
-        else:
-            res_num, H_name, X_name, intensity = intensity_fn(file_data[i], 
int_col)
+        res_num, H_name, X_name, intensity = intensity_fn(file_data[i], 
int_col)
 
         # Skip data.
         if X_name != heteronuc or H_name != proton:




Related Messages


Powered by MHonArc, Updated Thu Dec 04 02:20:02 2008