I have written some code which runs k rounds of Sherali-Adams relaxation hierarchy on a system Ax >= b. It was enough of a pain that I would like to potentially save others from it in the future :-)So I am wondering if I should contribute it here. I thin I could easily adapt it so that it could take a sage Polyhedron object and return the resulting extended Polyhedron. Does this seem like something Sage would want?