Get the QR factorization of a given NumPy array