-
template<typename T>
void gesc2(integer *n, T *a, integer *lda, T *rhs, integer *ipiv, integer *jpiv, T *scale)# GESC2 solves a system of linear equations using the LU factorization with complete pivoting computed by getc2.
Purpose:
GESC2 solves a system of linear equations A * X = scale* RHS with a general N-by-N matrix A using the LU factorization with complete pivoting computed by GETC2.
- Parameters:
N – [in]
N is INTEGER
The order of the matrix A.
A – [in]
A is REAL array, dimension (LDA,N)
On entry, the LU part of the factorization of the n-by-n matrix A computed by SGETC2: A = P * L * U * Q
LDA – [in]
LDA is INTEGER
The leading dimension of the array A. LDA >= fla_max(1, N).
RHS – [inout]
RHS is REAL array, dimension (N).
On entry, the right hand side vector b. On exit, the solution vector X.
IPIV – [in]
IPIV is INTEGER array, dimension (N).
The pivot indices; for 1 <= i <= N, row i of the matrix has been interchanged with row IPIV(i).
JPIV – [in]
JPIV is INTEGER array, dimension (N).
The pivot indices; for 1 <= j <= N, column j of the matrix has been interchanged with column JPIV(j).
SCALE – [out]
SCALE is REAL
On exit, SCALE contains the scale factor. SCALE is chosen 0 <= SCALE <= 1 to prevent overflow in the solution.