Singular value decomposition of m*n matrix using boost::ublas Very naive implementation :) Seems to be working already.