p = StaticPersistence(f) p.pair_simplices() dgms = init_diagrams(p, f) for i, dgm in enumerate(dgms): print "Dimension:", i print dgm