Purpose

LA_GESV computes the solution to a real or complex linear system of equations $AX = B$, where $A$ is a square matrix and $X$ and $B$ are rectangular matrices or vectors. Gaussian elimination with row interchanges is used to factor $A$ as $A=P \,L \,U$, where $P$ is a permutation matrix, $L$ is unit lower triangular, and $U$ is upper triangular. The factored form of $A$ is then used to solve the above system.