`zgetrf:`procedureThese routines compute an LU factorization of a general M-by-N matrix

*A*using partial pivoting with row interchanges. Optional argument*LDA*is the leading dimension of array*A*. The return values are:- a matrix containing the factors
*L*and*U*from the factorization*A = P*L*U*; - a vector with pivot indices: for 1 <= i <= min(M,N), row
*i*of the matrix was interchanged with row pivot(*i*)

- a matrix containing the factors