### Compute the centroid of a polytope

Does SAGE have a way to compute the center of mass of a polytope? I tried using polymake's centroid method, but that returned an error because the dimension of my polytope is not equal to the ambient dimension. I've read the manual and looked at the available functions but can't seem to find a SAGE command to do this. Does it exist?

EDIT: By request here's some sample _polymake_ code that fails (I'm not sure why this would be useful unless SAGE is using polymake under the ~~hood):
~~hood):

polytope > ~~$q=permutahedron(3);
~~$q=permutahedron(3);

polytope > print ~~$q->CENTROID;
~~$q->CENTROID;

polymake: WARNING: could not compute 'CENTROID' probably because of unsatisfied ~~preconditions:
~~preconditions:

precondition : DIM, AMBIENT_DIM ( CENTROID, VOLUME : VERTICES, TRIANGULATION.FACETS )

I'm guessing the error occurs because the ambient dimension of this polytope is 4 but the actual dimension of the polytope is 3; this is the same reason polymake tells me the volume of q is 0.