This function computes QR factorization of matrix \(A\)
\[A = Q R\]
where \(A\) is a dense matrix of size \(m \times n\), \(Q\) is a \(m \times n\) matrix with orthonormal columns, and \(R\) is an
upper triangular matrix.
The maximum matrix size supported in FPGA is templated by NRMAX and NCMAX.