author | Dmitriy Morozov <dmitriy@mrzv.org> |
Sun, 17 Jun 2012 12:14:23 -0700 | |
changeset 4 | 3315b6204a33 |
parent 2 | 4b3728f0d920 |
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)