while True:
    pt = show_diagram(dgms)
    if not pt: break
    print pt
    i = pt[2]
    smap = persistence.make_simplex_map(f)
    chain = [smap[ii] for ii in i.chain]
    pair_cycle = [smap[ii] for ii in i.pair().cycle]
    pair_chain = [smap[ii] for ii in i.pair().chain]
    show_complex(elephant_points, subcomplex = chain)
    if pt[1] != float('inf'):
        show_complex(elephant_points, subcomplex = pair_cycle + pair_chain)
    #show_complex(elephant_points, [s for s in f if s.dimension() != 1], subcomplex = cycle)