Just a quick follow-up for the posterity. As for v2.5.0, the bounds on the slack variable \epsilon, the input increments \mathbf{\Delta U} and the terminal states \mathbf{\hat{x}}(k+H_p) are now treated as box constraints. And pretty-printing an MPC object now shows the dimensions of the optimization problem, which is just clearer for the user.