The LHS matrix of the delta-u and state constraints (W and Gamma) could be generated as a kronecker product to speed things up a bit
The LHS matrix of the delta-u and state constraints (W and Gamma) could be generated as a kronecker product to speed things up a bit