author | Dmitriy Morozov <dmitriy@mrzv.org> |
Tue, 12 Jun 2012 22:35:53 -0700 | |
changeset 1 | a0e045bf9248 |
parent 0 | 337c84a13184 |
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: s.data[0]) show_diagram(dgms)