# File lib/linalg/dmatrix/lu.rb, line 24
      def lu # :yields: p, l, u
         res = lu_private(false)
         if block_given?
            yield res
         else
            res
         end
      end