◆ zsptrs()

Solution to factorized system of linear equations AX = B for a complex symmetric matrix in packed form

Purpose: This routine solves a system of linear equations
A * X = B

with a complex symmetric matrix A stored in packed form using the factorization A = U*D*U^T or A = L*D*L^T computed by zsptrf.

Parameters

[in]	uplo	Specifies whether the details of the factorization are stored as an upper or lower triangular matrix. = 'U': Upper triangular, form is A = UDU^T. = 'L': Lower triangular, form is A = LDL^T.
[in]	n	Order of the matrix A. (n >= 0) (If n = 0, returns without computation)
[in]	nrhs	Number of right hand sides, i.e., number of columns of the matrix B. (nrhs >= 0) (if nrhs = 0, returns without computation)
[in]	ap[]	Array ap[lap] (lap >= n(n + 1)/2) The block diagonal matrix D and the multipliers used to obtain the factor U or L as computed by zsptrf, stored as a packed triangular matrix.
[in]	ipiv[]	Array ipiv[lipiv] (lipiv >= n) Details of the interchanges and the block structure of D as determined by zsptrf.
[in]	ldb	Leading dimension of the two dimensional array b[][]. (ldb >= max(1, n))
[in,out]	b[][]	Array b[lb][ldb] (lb >= nrhs) [in] Right hand side matrix B. [out] Solution matrix X.
[out]	info	= 0: Successful exit = -1: The argument uplo had an illegal value (uplo != 'U' nor 'L') = -2: The argument n had an illegal value (n < 0) = -3: The argument nrhs had an illegal value (nrhs < 0) = -6: The argument ldb had an illegal value (ldb < max(1, n))