# 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)