Abstract Non-convex discrete-time optimal control problems in, e.g. , water or power systems, typically involve a large number of variables related through nonlinear equality constraints. The ideal goal is to find globally solution, and numerical experience indicates that algorithms aiming for Karush–Kuhn–Tucker points often solutions are indistinguishable from global optima. In our paper, we p...