Subroutine */
| 6323 | } /* zgesdd_ */ |
| 6324 | |
| 6325 | /* Subroutine */ int zgesv_(integer *n, integer *nrhs, doublecomplex *a, |
| 6326 | integer *lda, integer *ipiv, doublecomplex *b, integer *ldb, integer * |
| 6327 | info) |
| 6328 | { |
| 6329 | /* System generated locals */ |
| 6330 | integer a_dim1, a_offset, b_dim1, b_offset, i__1; |
| 6331 | |
| 6332 | /* Local variables */ |
| 6333 | extern /* Subroutine */ int xerbla_(char *, integer *), zgetrf_( |
| 6334 | integer *, integer *, doublecomplex *, integer *, integer *, |
| 6335 | integer *), zgetrs_(char *, integer *, integer *, doublecomplex *, |
| 6336 | integer *, integer *, doublecomplex *, integer *, integer *); |
| 6337 | |
| 6338 | |
| 6339 | /* |
| 6340 | -- LAPACK driver routine (version 3.2) -- |
| 6341 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 6342 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 6343 | November 2006 |
| 6344 | |
| 6345 | |
| 6346 | Purpose |
| 6347 | ======= |
| 6348 | |
| 6349 | ZGESV computes the solution to a complex system of linear equations |
| 6350 | A * X = B, |
| 6351 | where A is an N-by-N matrix and X and B are N-by-NRHS matrices. |
| 6352 | |
| 6353 | The LU decomposition with partial pivoting and row interchanges is |
| 6354 | used to factor A as |
| 6355 | A = P * L * U, |
| 6356 | where P is a permutation matrix, L is unit lower triangular, and U is |
| 6357 | upper triangular. The factored form of A is then used to solve the |
| 6358 | system of equations A * X = B. |
| 6359 | |
| 6360 | Arguments |
| 6361 | ========= |
| 6362 | |
| 6363 | N (input) INTEGER |
| 6364 | The number of linear equations, i.e., the order of the |
| 6365 | matrix A. N >= 0. |
| 6366 | |
| 6367 | NRHS (input) INTEGER |
| 6368 | The number of right hand sides, i.e., the number of columns |
| 6369 | of the matrix B. NRHS >= 0. |
| 6370 | |
| 6371 | A (input/output) COMPLEX*16 array, dimension (LDA,N) |
| 6372 | On entry, the N-by-N coefficient matrix A. |
| 6373 | On exit, the factors L and U from the factorization |
| 6374 | A = P*L*U; the unit diagonal elements of L are not stored. |
| 6375 | |
| 6376 | LDA (input) INTEGER |
| 6377 | The leading dimension of the array A. LDA >= max(1,N). |
| 6378 | |
| 6379 | IPIV (output) INTEGER array, dimension (N) |
| 6380 | The pivot indices that define the permutation matrix P; |
| 6381 | row i of the matrix was interchanged with row IPIV(i). |
| 6382 |