where A is a known m -by- n matrix of rank r and b is a known m -vector. We note that this set of constraint equations encompass a very general variety of holonomic and non-holonomic equality constraints. For example, holonomic constraints of the form