bg = GRIMMReader.get_breakpoint_graph(source) bg_graphs[file_name.split(".")[0]] = bg bg = BreakpointGraph() for br_gr in bg_graphs.values(): bg.update(br_gr, merge_edges=True) target_multicolor = Multicolor("Anguilla_japonica") print("Breakpoint graph stats:") print( "\t", "non-infinity nodes count:", len( list(node for node in bg.nodes() if not BGVertex.is_infinity_vertex(node)))) normal_edges, infinity_edges = [], [] for edge in bg.edges(): if edge.is_infinity_edge: infinity_edges.append(edge) else: normal_edges.append(edge) print("\t", "non-infinity edges count:", len(normal_edges)) print("\t", "infinity edges count:", len(infinity_edges)) print("\t", "connected component count:", len(list(bg.connected_components_subgraphs()))) print() print() cc_cnt_with_exactly_two_targeted_infinity_edges = 0 for cnt, cc in enumerate(bg.connected_components_subgraphs()): infinity_edges = defaultdict(int)