Lower-Upper Decomposition (GETRF)

This function computes the LU decomposition (with partial pivoting) of matrix $A$

$$A = L U$$

where $A$ is a dense matrix of size $n \times n$, $L$ is a lower triangular matrix with unit diagonal, and $U$ is a upper triangular matrix. This function implement partial pivoting. The maximum matrix size supported in FPGA is templated by NMAX.