Added a SVD module:

- the decompostion code has been adfapted from JAMA
 - handles non square matrices of size MxN with M>=N
 - does not work for complex matrices
 - includes a solver where the parts corresponding to zero singular values are set to zero

Eigen/src/Core/MatrixBase.h

@@ -557,6 +557,10 @@ template<typename Derived> class MatrixBase
     EigenvaluesReturnType eigenvalues() const;
     RealScalar operatorNorm() const;
 
+/////////// SVD module ///////////
+
+    const SVD<EvalType> svd() const;
+
 /////////// Geometry module ///////////
 
     template<typename OtherDerived>