Consider adding a test for non full rank matrix: e.g. take the svd of an outer product of a vector and see that computeSVD(1) recovers the vector