# Rips
points = read_points('data/trefoil.pts')
distances = PairwiseDistances(points)
rips = Rips(distances)

f = Filtration()
rips.generate(2, 1.7, f.append)
print "Number of simplices:", len(f)

show_complex(points, f)
show_complex(points, [s for s in f if rips.eval(s) < 1.6])

f.sort(rips.cmp)

p = StaticPersistence(f)
p.pair_simplices()

dgms = init_diagrams(p, f, rips.eval)
show_diagram(dgms[:2])