Subroutine */
| 972 | } /* zgebd2_ */ |
| 973 | |
| 974 | /* Subroutine */ int zgebrd_(integer *m, integer *n, doublecomplex *a, |
| 975 | integer *lda, doublereal *d__, doublereal *e, doublecomplex *tauq, |
| 976 | doublecomplex *taup, doublecomplex *work, integer *lwork, integer * |
| 977 | info) |
| 978 | { |
| 979 | /* System generated locals */ |
| 980 | integer a_dim1, a_offset, i__1, i__2, i__3, i__4, i__5; |
| 981 | doublereal d__1; |
| 982 | doublecomplex z__1; |
| 983 | |
| 984 | /* Local variables */ |
| 985 | static integer i__, j, nb, nx; |
| 986 | static doublereal ws; |
| 987 | static integer nbmin, iinfo, minmn; |
| 988 | extern /* Subroutine */ int zgemm_(char *, char *, integer *, integer *, |
| 989 | integer *, doublecomplex *, doublecomplex *, integer *, |
| 990 | doublecomplex *, integer *, doublecomplex *, doublecomplex *, |
| 991 | integer *), zgebd2_(integer *, integer *, |
| 992 | doublecomplex *, integer *, doublereal *, doublereal *, |
| 993 | doublecomplex *, doublecomplex *, doublecomplex *, integer *), |
| 994 | xerbla_(char *, integer *), zlabrd_(integer *, integer *, |
| 995 | integer *, doublecomplex *, integer *, doublereal *, doublereal *, |
| 996 | doublecomplex *, doublecomplex *, doublecomplex *, integer *, |
| 997 | doublecomplex *, integer *); |
| 998 | extern integer ilaenv_(integer *, char *, char *, integer *, integer *, |
| 999 | integer *, integer *, ftnlen, ftnlen); |
| 1000 | static integer ldwrkx, ldwrky, lwkopt; |
| 1001 | static logical lquery; |
| 1002 | |
| 1003 | |
| 1004 | /* |
| 1005 | -- LAPACK routine (version 3.2) -- |
| 1006 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 1007 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 1008 | November 2006 |
| 1009 | |
| 1010 | |
| 1011 | Purpose |
| 1012 | ======= |
| 1013 | |
| 1014 | ZGEBRD reduces a general complex M-by-N matrix A to upper or lower |
| 1015 | bidiagonal form B by a unitary transformation: Q**H * A * P = B. |
| 1016 | |
| 1017 | If m >= n, B is upper bidiagonal; if m < n, B is lower bidiagonal. |
| 1018 | |
| 1019 | Arguments |
| 1020 | ========= |
| 1021 | |
| 1022 | M (input) INTEGER |
| 1023 | The number of rows in the matrix A. M >= 0. |
| 1024 | |
| 1025 | N (input) INTEGER |
| 1026 | The number of columns in the matrix A. N >= 0. |
| 1027 | |
| 1028 | A (input/output) COMPLEX*16 array, dimension (LDA,N) |
| 1029 | On entry, the M-by-N general matrix to be reduced. |
| 1030 | On exit, |
| 1031 | if m >= n, the diagonal and the first superdiagonal are |