Version 0.8.0

VERSION

@@ -1 +1 @@
-0.6.0
+0.8.0

src/ecopcr.c

@@ -6,7 +6,7 @@
 #include <getopt.h>
 
 
-#define VERSION "0.6.0"
+#define VERSION "0.8.0"
 
 
 /* ----------------------------------------------- */
@@ -26,7 +26,7 @@ static void PrintHelp()
         PP      "-a    : Salt concentration in M for Tm computation (default 0.05 M)\n\n");
         PP      "-c    : Consider that the database sequences are [c]ircular\n\n");
         PP      "-d    : [D]atabase : to match the expected format, the database\n");
-        PP      "        has to be formated first by the ecoPCRFormat.py program located.\n");
+        PP      "        has to be formatted first by the ecoPCRFormat.py program located.\n");
         PP      "        in the tools directory.\n");
         PP      "        ecoPCRFormat.py creates three file types :\n");
         PP      "            .sdx : contains the sequences\n");

src/global.mk

@@ -1,14 +1,13 @@
-MACHINE=MAC_OS_X
 LIBPATH= -Llibapat -LlibecoPCR -Llibthermo
-MAKEDEPEND = gcc -D$(MACHINE) -M $(CPPFLAGS) -o $*.d $<
+MAKEDEPEND = gcc -M $(CPPFLAGS) -o $*.d $<
 
 CC=gcc
-CFLAGS= -W -Wall -O2 -g
+CFLAGS= -O3 -w
 
 default: all
 
 %.o: %.c
-	$(CC) -D$(MACHINE) $(CFLAGS) -c -o $@ $<
+	$(CC) $(CFLAGS) -c -o $@ $<
 
 %.P : %.c
 	$(MAKEDEPEND)

src/libapat/Gtypes.h

@@ -27,13 +27,9 @@
 #define PROTO   1                       /* prototypes flag              */
 #endif
 
-#ifdef MAC_OS_C
-#define Vrai    true                    /* TC boolean values            */
-#define Faux    false                   /*                              */
-#else
+
 #define Vrai    0x1                     /* bool values  = TRUE          */
 #define Faux    0x0                     /*              = FALSE         */
-#endif
 
 #define Nil     NULL                    /* nil pointer                  */
 
@@ -42,28 +38,7 @@
 #define kBigUInt16      0xffff          /* plus grand 16 bits ~signe    */
 #define kBigUInt32      0xffffffff      /* plus grand 32 bits ~signe    */
 
-#ifdef MAC_OS_C
-/* ==================================================== */
-/*  Types (for Macintosh ThinK C || MWerks)             */
-/* ==================================================== */
 
-                                        /* --- specific sizes --------- */
-typedef long            Int32;          /* Int32  = 32 bits signe       */
-typedef unsigned long   UInt32;         /* UInt32 = 32 bits ~signe      */
-typedef short           Int16;          /* Int16  = 16 bits signe       */
-typedef unsigned short  UInt16;         /* UInt32 = 16 bits ~signe      */
-typedef char            Int8;           /* Int8   = 8 bits signe        */
-typedef unsigned char   UInt8;          /* UInt8  = 8 bits ~signe       */
-
-                                        /* --- default types ---------- */
-
-typedef Boolean         Bool;           /* booleen                      */
-
-typedef long            Int;            /* 'natural' int (>= 32 bits)   */
-
-typedef void            *Ptr;           /* pointeur                     */
-
-#elif ((defined SUN) || (defined SGI) || (defined UNIX))
 /* ==================================================== */
 /*  Types (for Sun & Iris - 32 bits machines)           */
 /* ==================================================== */
@@ -84,14 +59,7 @@ typedef int             Int;            /* 'natural' int (>= 32 bits)   */
 
 typedef void            *Ptr;           /* pointeur                     */
 
-#else 
-/* ==================================================== */
-/*  Types (for undefined machines)                      */
-/* ==================================================== */
 
-#error undefined MACHINE <please edit Gmach.h>
-
-#endif
 
 /* ==================================================== */
 /*  special macro for prototypes                        */

tools/ecoPCRFormat.py

@@ -1,4 +1,4 @@
-#!/usr/bin/env python
+#!/usr/bin/env python2.7
 
 import re
 import gzip
@@ -7,11 +7,8 @@ import sys
 import time
 import getopt
 
-try:
-    import psycopg2
-    _dbenable=True
-except ImportError:
-    _dbenable=False
+_dbenable=False
+
 
 #####
 #
@@ -221,56 +218,7 @@ def readTaxonomyDump(taxdir):
     
     return taxonomy,ranks,alternativeName,index
 
-def readTaxonomyDB(dbname):
-    connection = psycopg2.connect(database=dbname)
-    
-    cursor = connection.cursor()
-    cursor.execute("select numid,rank,parent from ncbi_taxonomy.taxon")
-    taxonomy=[list(x) for x in cursor]
 
-    cursor.execute("select rank_class from ncbi_taxonomy.taxon_rank_class order by rank_class")
-    ranks=cursor.fetchall()
-    ranks = dict(map(None,(x[0] for x in ranks),xrange(len(ranks))))
-    
-    print >>sys.stderr,"Sorting taxons..."
-    taxonomy.sort(taxonCmp)
-
-    print >>sys.stderr,"Indexing taxonomy..."
-    index = {}
-    for t in taxonomy:
-        index[t[0]]=bsearchTaxon(taxonomy, t[0])
-    
-    print >>sys.stderr,"Indexing parent and rank..."
-    for t in taxonomy:
-        t[1]=ranks[t[1]]
-        try:
-            t[2]=index[t[2]]
-        except KeyError,e:
-            if t[2] is None and t[0]==1:
-                t[2]=index[t[0]]
-            else:
-                raise e
-        
-    cursor.execute("select taxid,name,category from ncbi_taxonomy.name")
-    
-    alternativeName=[]
-    for taxid,name,classname in cursor:
-        alternativeName.append((name,classname,index[taxid]))
-        if classname == 'scientific name':
-            taxonomy[index[taxid]].append(name)
-    
-    cursor.execute("select old_numid,current_numid from ncbi_taxonomy.taxon_id_alias")
- 
-    print >>sys.stderr,"Adding taxid alias..."
-    for taxid,current in cursor:
-        if current is not None:
-            index[taxid]=index[current]
-        else:
-            index[taxid]=None
-            
-       
-    return taxonomy,ranks,alternativeName,index
-    
 #####
 #
 #
@@ -585,10 +533,9 @@ def ecoParseOptions(arguments):
            }
     
     o,filenames = getopt.getopt(arguments,
-                                'ht:T:n:gfe',
+                                'ht:n:gfe',
                                 ['help',
                                  'taxonomy=',
-                                 'taxonomy_db=',
                                  'name=',
                                  'genbank',
                                  'fasta',
@@ -601,9 +548,6 @@ def ecoParseOptions(arguments):
         elif name in ('-t','--taxonomy'):
             opt['taxmod']='dump'
             opt['taxdir']=value
-        elif name in ('-T','--taxonomy_db'):
-            opt['taxmod']='db'
-            opt['taxdb']=value
         elif name in ('-n','--name'):
             opt['prefix']=value
         elif name in ('-g','--genbank'):
@@ -622,6 +566,7 @@ def ecoParseOptions(arguments):
 
     return opt,filenames
 
+
 def printHelp():
     print "-----------------------------------"
     print " ecoPCRFormat.py"
@@ -641,11 +586,7 @@ if __name__ == '__main__':
     
     opt,filenames = ecoParseOptions(sys.argv[1:])
     
-    if opt['taxmod']=='dump':
-        taxonomy = readTaxonomyDump(opt['taxdir'])
-    elif opt['taxmod']=='db':
-        taxonomy = readTaxonomyDB(opt['taxdb'])
-        
-   
+    taxonomy = readTaxonomyDump(opt['taxdir'])
+    
     ecoDBWriter(opt['prefix'], taxonomy, filenames, opt['parser'])