Hi 
I am using python interface of acados. Recently, I am thinking about the following form of dynamics:
F(x_k) + G(u_k) = F(x_{k+1}). This can be interpreted as implicit variational integrators which may help convergence and reduce number of lifting states.