def METHOD_NAME(cg, edges: Iterable, time_stamp: Optional[float] = None):
    from ...utils.general import in2d
    coords0 = chunk_utils.get_chunk_coordinates_multiple(cg.meta, edges[:, 0])
    coords1 = chunk_utils.get_chunk_coordinates_multiple(cg.meta, edges[:, 1])
    coords = np.concatenate([np.array(coords0), np.array(coords1)])
    bbox = [np.min(coords, axis=0), np.max(coords, axis=0)]
    bbox[1] += 1
    root_ids = set(
        cg.get_roots(edges.ravel(), assert_roots=True, time_stamp=time_stamp)
    )
    sg_edges = cg.get_subgraph(
        root_ids,
        bbox=bbox,
        bbox_is_coordinate=False,
        edges_only=True,
    )
    existence_status = in2d(edges, sg_edges)
    edge_layers = cg.get_cross_chunk_edges_layer(edges)
    active_status = []
    for layer in np.unique(edge_layers):
        layer_edges = edges[edge_layers == layer]
        edges_parents = cg.get_roots(
            layer_edges.ravel(), time_stamp=time_stamp, stop_layer=layer + 1
        ).reshape(-1, 2)
        mask = edges_parents[:, 0] == edges_parents[:, 1]
        active_status.extend(mask)
    active_status = np.array(active_status, dtype=bool)
    return existence_status, active_status