Subroutine */
| 2322 | } /* sgebd2_ */ |
| 2323 | |
| 2324 | /* Subroutine */ int sgebrd_(integer *m, integer *n, real *a, integer *lda, |
| 2325 | real *d__, real *e, real *tauq, real *taup, real *work, integer * |
| 2326 | lwork, integer *info) |
| 2327 | { |
| 2328 | /* System generated locals */ |
| 2329 | integer a_dim1, a_offset, i__1, i__2, i__3, i__4; |
| 2330 | |
| 2331 | /* Local variables */ |
| 2332 | static integer i__, j, nb, nx; |
| 2333 | static real ws; |
| 2334 | static integer nbmin, iinfo; |
| 2335 | extern /* Subroutine */ int sgemm_(char *, char *, integer *, integer *, |
| 2336 | integer *, real *, real *, integer *, real *, integer *, real *, |
| 2337 | real *, integer *); |
| 2338 | static integer minmn; |
| 2339 | extern /* Subroutine */ int sgebd2_(integer *, integer *, real *, integer |
| 2340 | *, real *, real *, real *, real *, real *, integer *), slabrd_( |
| 2341 | integer *, integer *, integer *, real *, integer *, real *, real * |
| 2342 | , real *, real *, real *, integer *, real *, integer *), xerbla_( |
| 2343 | char *, integer *); |
| 2344 | extern integer ilaenv_(integer *, char *, char *, integer *, integer *, |
| 2345 | integer *, integer *, ftnlen, ftnlen); |
| 2346 | static integer ldwrkx, ldwrky, lwkopt; |
| 2347 | static logical lquery; |
| 2348 | |
| 2349 | |
| 2350 | /* |
| 2351 | -- LAPACK routine (version 3.2) -- |
| 2352 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 2353 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 2354 | November 2006 |
| 2355 | |
| 2356 | |
| 2357 | Purpose |
| 2358 | ======= |
| 2359 | |
| 2360 | SGEBRD reduces a general real M-by-N matrix A to upper or lower |
| 2361 | bidiagonal form B by an orthogonal transformation: Q**T * A * P = B. |
| 2362 | |
| 2363 | If m >= n, B is upper bidiagonal; if m < n, B is lower bidiagonal. |
| 2364 | |
| 2365 | Arguments |
| 2366 | ========= |
| 2367 | |
| 2368 | M (input) INTEGER |
| 2369 | The number of rows in the matrix A. M >= 0. |
| 2370 | |
| 2371 | N (input) INTEGER |
| 2372 | The number of columns in the matrix A. N >= 0. |
| 2373 | |
| 2374 | A (input/output) REAL array, dimension (LDA,N) |
| 2375 | On entry, the M-by-N general matrix to be reduced. |
| 2376 | On exit, |
| 2377 | if m >= n, the diagonal and the first superdiagonal are |
| 2378 | overwritten with the upper bidiagonal matrix B; the |
| 2379 | elements below the diagonal, with the array TAUQ, represent |
| 2380 | the orthogonal matrix Q as a product of elementary |
| 2381 | reflectors, and the elements above the first superdiagonal, |