From a speed perspective it might be worth having different methods for extinction() for uni and bipartite network. In theory for bipartite networks a simple simplify() call will be sufficient to catch secondary extinctions - which is quicker than the current method