Towards a theory of non-commutative optimization: geodesic 1st and 2nd order methods for moment maps and polytopes