Zgbcon Zgbtrf Zgbtrs Zgecon Zgetrf Zgetri Zgetrs Zgtcon Zgttrf Zgttrs Zspcon Zsptrf Zsptri Zsptrs Zsycon Zsytrf Zsytri Zsytrs

Zsptrs() Sub Zsptrs ( Uplo As  String, N As  Long, Ap() As  Complex, IPiv() As  Long, B() As  Complex, Info As  Long, Optional Nrhs As  Long = 1  ) 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 = U*D*U^T. = "L": Lower triangular, form is A = L*D*L^T. [in] N Order of the matrix A. (N >= 0) (If N = 0, returns without computation) [in] Ap() Array Ap(LAp - 1) (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 - 1) (LIPiv >= N) Details of the interchanges and the block structure of D as determined by Zsptrf. [in,out] B() Array B(LB1 - 1, LB2 - 1) (LB1 >= max(1, N), LB2 >= Nrhs) (2D array) or B(LB - 1) (LB >= max(1, N), Nrhs = 1) (1D array) [in] N x Nrhs matrix of right hand side matrix B. [out] If Info = 0, the N x Nrhs 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 Ap() is invalid. = -4: The argument IPiv() is invalid. = -5: The argument B() is invalid. = -7: The argument Nrhs had an illegal value. (Nrhs < 0) [in] Nrhs (Optional) Number of right hand sides, i.e., number of columns of the matrix B. (Nrhs >= 0) (If Nrhs = 0, returns without computation) (default = 1) Reference LAPACK Example Program See example of Zsptrf.