author | Dmitriy Morozov <dmitriy@mrzv.org> |
Mon, 18 Jun 2012 10:04:58 -0700 | |
changeset 5 | 9acb920a945c |
parent 4 | 3315b6204a33 |
permissions | -rw-r--r-- |
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)