qr() equivalent in an nxAlgebra

In order to implement alternative FIML objective functions, it would be nice to have an equivalent to R's QR decomposition.