Methods and systems for the design and implementation of optimal multivariable MPC controllers for fast-sampling constrained dynamic systems utilizing a primal-dual feasibility approach and/or a graph approach. The primal-dual feasibility approach can compute and store matrices defining constraints of quadratic programming problems in an off-line part in order to calculate vectors of Lagrange multipliers and an optimizer. Then primal-dual feasibility can be checked in an on-line part using the Lagrange multipliers and the optimizer can provide a unique optimal solution for the constrained dynamic system. The graph approach can compute and store the matrices and the vectors, and also prepare and store a structure of directed graph in off-line part. An optimizer for a given parameter vector can be determined in on-line part using the directed graph, the matrices and the vectors.