Two improvements to Su regularization:
1) Support dynamic sized input matrices by using runtime size information
2) Test for "feedthrough" term in output equation.  Implementation assumes it's zero.