print recursion limit

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,22 +614,23 @@ 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():
         if len(list_pr) > 1:
             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
+    return G