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