author | Dmitriy Morozov <dmitriy@mrzv.org> |
Mon, 18 Jun 2012 10:04:58 -0700 | |
changeset 5 | 9acb920a945c |
parent 2 | 4b3728f0d920 |
permissions | -rw-r--r-- |
# Alpha shapes points = read_points('data/trefoil.pts') f = Filtration() fill_alpha_complex(points, f) for s in f: print s, s.data from math import sqrt show_complex(points, [s for s in f if sqrt(s.data[0]) < .5]) show_complex(points, [s for s in f if sqrt(s.data[0]) < .8]) f.sort(dim_data_cmp) p = StaticPersistence(f) p.pair_simplices() dgms = init_diagrams(p, f, lambda s: sqrt(s.data[0])) show_diagram(dgms)