# File lib/linalg/dmatrix/nullspace.rb, line 97
      def regular?(epsilon = (self.singleton_class.default_epsilon ||self.class.default_epsilon))
         square? and rank(epsilon) == hsize
      end