Fixed small bug (--sl merging)

felipelouza · Sep 16, 2021 · b190421 · b190421
1 parent 2da5546
commit b190421
Showing 1 changed file with 8 additions and 3 deletions.
diff --git a/eGap b/eGap
@@ -203,9 +203,6 @@ def check_input(args):
     if (args.sa):
       print("SA construction not supported for merging BWT files!")
       sys.exit(1)
-    if (args.sl):
-      print("SL construction not supported for merging BWT files!")
-      sys.exit(1)
     args.basename = args.out
   # ---- if the input are concatenated texts there is a single file
   else:
@@ -283,6 +280,14 @@ def phase1(args,logfile, logfile_name):
           filename = os.path.splitext(name)[0]+".{byts}.da".format(byts=args.dbytes)
           with open(filename,'rb') as fd:
             shutil.copyfileobj(fd, dablfile, 1024*1024*10) # 10 MBs buffer
+    if(args.sl): 
+      # concatenate .sl files in a single .sl_bl file 
+      print("==== creating .sl_bl file")
+      with open(args.basename+".{byts}.sl_bl".format(byts=args.slbytes),"wb") as slblfile:
+        for name in args.input:
+          filename = os.path.splitext(name)[0]+".{byts}.sl".format(byts=args.slbytes)
+          with open(filename,'rb') as fd:
+            shutil.copyfileobj(fd, slblfile, 1024*1024*10) # 10 MBs buffer
     return True # everything fine
   else:        
     # ---- gSACAK