# File lib/linalg/dmatrix/nullspace.rb, line 65
      def rank(epsilon = (self.singleton_class.default_epsilon ||self.class.default_epsilon))
         sv = self.sv
         indexes = sv_null_indexes(sv, sv[0]*epsilon)
         indexes ? (hsize - indexes.size) : 0
      end