# File lib/linalg/dmatrix/fit.rb, line 16
      def pseudo_inverse(epsilon = nil)
         DMatrix.fit(self, DMatrix.identity(vsize), epsilon)[0]
      end