Skip to content
Snippets Groups Projects
Commit b57b0bd3 authored by nina.marthe_ird.fr's avatar nina.marthe_ird.fr
Browse files

commented debug print lines

parent 8efb565b
No related branches found
No related tags found
No related merge requests found
......@@ -8,20 +8,15 @@ from .Functions_output import *
from .argparser import *
#from inference import *
from pathlib import Path
import psutil
#import psutil
def main():
#import tracemalloc
#tracemalloc.start()
args=arg()
read_args(args)
print(psutil.Process().memory_info().rss)
#print(tracemalloc.get_traced_memory())
#tracemalloc.reset_peak()
#print(psutil.Process().memory_info().rss)
# intersect is in the output directory
intersect=Path(args.outdir.joinpath("intersect")).resolve()
......@@ -29,7 +24,7 @@ def main():
load_intersect(intersect.as_posix(),args.verbose)
#print(psutil.Process().memory_info().rss / 1024**2) # B -> MiB
print(psutil.Process().memory_info().rss)
#print(psutil.Process().memory_info().rss)
segments=args.segment_coordinates_path.joinpath("segments.txt")
......@@ -62,7 +57,7 @@ def main():
if args.verbose:
print(f" Genomes found : {args.target}")
print(psutil.Process().memory_info().rss) # B -> MiB
#print(psutil.Process().memory_info().rss)
# get haplotypes for each target genome
if args.haplotype:
......@@ -86,7 +81,7 @@ def main():
# build a dictionnary with the segment sizes to compute the coverage and id
if not args.graph_gaf:
seg_size=get_segments_length(segments,args.verbose)
print(psutil.Process().memory_info().rss) # B -> MiB
#print(psutil.Process().memory_info().rss)
genome_index=0
for target_genome in args.target:
print(f'\n{target_genome} transfer :')
......@@ -101,7 +96,7 @@ def main():
print(f' Loading the walks for the genome {target_genome}')
walks_path=args.segment_coordinates_path.joinpath("walks.txt").as_posix()
target_genome_paths=get_paths(walks_path,target_genome,args.haplotype)
print(psutil.Process().memory_info().rss) # B -> MiB
#print(psutil.Process().memory_info().rss)
if not args.verbose:
print(" Loading the segments coordinates")
for file in segment_coord_files:
......@@ -118,9 +113,9 @@ def main():
list_feat_absent=[]
# do the annotation transfer (or var/aln)
print(psutil.Process().memory_info().rss) # B -> MiB
#print(psutil.Process().memory_info().rss)
transfer_on_target(segments,out_target_gff,out_target_var,out_target_aln,target_genome,target_genome_paths,list_feat_absent,seg_size,args)
print(psutil.Process().memory_info().rss) # B -> MiB
#print(psutil.Process().memory_info().rss)
# if pav matrix is asked, add the information of this transfer on the matrix
if args.pav_matrix:
for feat in list_feat_absent:
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment