From 5a29a32fcd08eacabb255e9082e5bf246ba77b5e Mon Sep 17 00:00:00 2001 From: Sina Majidian Date: Sat, 19 Oct 2024 11:11:50 -0400 Subject: [PATCH] print recursion limit --- FastOMA/_utils_roothog.py | 12 +++++++----- 1 file changed, 7 insertions(+), 5 deletions(-) diff --git a/FastOMA/_utils_roothog.py b/FastOMA/_utils_roothog.py index 1c8c114..87d179c 100644 --- a/FastOMA/_utils_roothog.py +++ b/FastOMA/_utils_roothog.py @@ -5,6 +5,7 @@ import pickle from os import listdir import os +import sys from .zoo.unionfind import UnionFind from ._wrappers import logger @@ -613,14 +614,15 @@ def merge_rhogs2(hogmaps, rhogs_prots, conf_infer_roothogs): cluster_rhogs_list = cluster_rhogs(candidates_pair) print("There are " + str(len(cluster_rhogs_list)) + " clusters.") + logger.debug("There are " + str(len(cluster_rhogs_list)) + " clusters.") + print("** the recursion limit is "+str( sys.getrecursionlimit())) + logger.debug("** the recursion limit is "+str( sys.getrecursionlimit())) cluster_rhogs_list = cluster_rhogs_nx(cluster_rhogs_list, candidates_pair) print("There are " + str(len(cluster_rhogs_list)) + " selected clusters.") - #cluster_rhogs_list = cluster_rhogs(candidates_pair) - - + # cluster_rhogs_list = cluster_rhogs(candidates_pair) num_rhog_g1 = 0 for rhg, list_pr in rhogs_prots.items(): @@ -628,7 +630,7 @@ def merge_rhogs2(hogmaps, rhogs_prots, conf_infer_roothogs): num_rhog_g1 += 1 logger.debug("There are " + str(num_rhog_g1) + " rhogs (size>1) before merging.") - print(len(rhogs_prots)) + print("len(rhogs_prots) is ", len(rhogs_prots)) file_out_merge = open("merging_rhogs.txt","w") file_out_merge.write("#first column is the host hog, the rest will be merged here.\n") @@ -1205,4 +1207,4 @@ def HCS(G): G = nx.compose(H, _H) - return G \ No newline at end of file + return G