changes for overcomplete matrices

miscmaths.cc

@@ -512,10 +512,13 @@ namespace MISCMATHS {
       }
     }
 
-  ReturnMatrix pinv(const Matrix& mat)
+  ReturnMatrix pinv(const Matrix& mat2)
     {
       // calculates the psuedo-inverse using SVD
       // note that the right-pinv(x') = pinv(x).t()
+      Matrix mat(mat2);
+      if ( mat2.Ncols() > mat2.Nrows() )
+	mat=mat.t();
       Tracer tr("pinv");
       DiagonalMatrix D;
       Matrix U, V;
@@ -527,6 +530,8 @@ namespace MISCMATHS {
 	else D(n,n) = 0.0; // reduce the number of columns because too close to singular
       }
       Matrix pinv = V * D * U.t();
+      if ( mat2.Ncols() > mat2.Nrows() )
+	pinv=pinv.t();
       pinv.Release();
       return pinv;
     }