FactorizationPLURaw

Computes an LU factorization of a general M-by-N matrix A using partial pivoting with row interchanges. The factorization has the form

   A = P * L * U

where P is a permutation matrix, L is lower triangular with unit diagonal elements (lower trapezoidal if m > n), and U is upper triangular (upper trapezoidal if m < n). LAPACK function GETRF.

Computing for type matrix<double>

bool  matrix::FactorizationPLURaw(
  matrix&       AF,           // factored matrix A
  long[]&       ipiv         // pivot indices array
  );

Computing for type matrix<float>

bool  matrixf::FactorizationPLURaw(
  matrixf&     AF,           // factored matrix A
  long[]&       ipiv         // pivot indices array
  );

Computing for type matrix<complex>

bool  matrixc::FactorizationPLURaw(
  matrixc&     AF,           // factored matrix A
  long[]&       ipiv         // pivot indices array
  );

Computing for type matrix<complexf>

bool  matrixcf::FactorizationPLURaw(
  matrixcf&     AF,           // factored matrix A
  long[]&       ipiv         // pivot indices array
  );

Parameters

AF

[out]  Factored matrix A. The factors L and U from the factorization  A = P*L*U; the unit diagonal elements of L are not stored.

ipiv

[out]  Pivot indices array of size M; row i of the matrix A was interchanged with row  ipiv[i].

 

Return Value

Return true if successful, otherwise false in case of an error.

Note

Matrix AF and pivot indices array ipiv[] are raw output of the GETRF function and can be used for further calculations with methods PLULinearEquationsSolution, PLUInverse and PLUCondNumReciprocal. In these cases original matrix A must be square.