author | Dmitriy Morozov <dmitriy@mrzv.org> |
Sat, 16 Jun 2012 12:13:53 -0700 | |
changeset 2 | 4b3728f0d920 |
parent 0 | 05-rips.py@337c84a13184 |
permissions | -rw-r--r-- |
# Rips points = read_points('data/trefoil.pts') distances = PairwiseDistances(points) distances = ExplicitDistances(distances) 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])