In this technical note we introduce a multivariable extension of a recently proposed model reduction algorithm for linear time-invariant systems. Reduced models are found by solving a convex. optimization problem with Linear Matrix Inequality constraints given a state space model or frequency samples of a model. In order to illustrate the method we apply it to a large-scale model of a deformable mirror.