Discrete-time Contraction Constrained Nonlinear Model Predictive Control using Graph-based Geodesic Computation