from    math        import sqrt

f = Filtration()
fill_alpha_complex(points, f)
f.sort(dim_data_cmp)

p = StaticCohomologyPersistence(f, prime = 11)
p.pair_simplices()

dgms = init_diagrams(p,f, lambda s: sqrt(s.data[0]), lambda n: n.cocycle)

while True:
    pt = show_diagram(dgms)
    if not pt: break

    #rf = Filtration((s for s in f if sqrt(s.data[0]) <= pt[0]))
    #rf = Filtration((s for s in f if sqrt(s.data[0]) < pt[1]))
    rf = Filtration((s for s in f if sqrt(s.data[0]) <= (pt[0] + pt[1])/2))
    values = circular.smooth(rf, pt[2])
    cocycle = [rf[i] for (c,i) in pt[2] if i < len(rf)]
    #print cocycle
    show_complex(points, subcomplex = cocycle)
    show_complex(points, values = values)
    #show_complex(points, rf, values = values, subcomplex = cocycle)