added to pinv - actually zero singular values if they're less than tol (when...

added to pinv - actually zero singular values if they're less than tol (when inverting them), rather than just leaving them as they are

miscmaths.cc

@@ -513,7 +513,8 @@ namespace MISCMATHS {
       float tol;
       tol = MaximumAbsoluteValue(D) * Max(mat.Nrows(),mat.Ncols()) * 1e-16;
       for (int n=1; n<=D.Nrows(); n++) {
-	if (fabs(D(n,n))>tol)  D(n,n) = 1.0/D(n,n);
+	if (fabs(D(n,n))>tol) D(n,n) = 1.0/D(n,n);
+	else D(n,n) = 0.0; // reduce the number of columns because too close to singular
       }
       Matrix pinv = V * D * U.t();
       pinv.Release();