load("Diagonal.sage")

# The following routine considers the polytope P3, its preferred colouring P3_colouring() and state P3_state(), 
# and classifies all the ascending and descending links. It groups these links into isomorphism classes, and for each class
# it computes a presentation of the fundamental group (that turns out to be always trivial) and its Betti numbers
# It makes a final check: every link contributes to the final Euler characteristic of the manifold.

get_links_topology(P3(), P3_colouring(), P3_state())

# The same routine for P4, ..., P8 is obtained by substituting "P3" with "P4", ... , "P8"