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