Subroutine */
| 3347 | } /* sgehd2_ */ |
| 3348 | |
| 3349 | /* Subroutine */ int sgehrd_(integer *n, integer *ilo, integer *ihi, real *a, |
| 3350 | integer *lda, real *tau, real *work, integer *lwork, integer *info) |
| 3351 | { |
| 3352 | /* System generated locals */ |
| 3353 | integer a_dim1, a_offset, i__1, i__2, i__3, i__4; |
| 3354 | |
| 3355 | /* Local variables */ |
| 3356 | static integer i__, j; |
| 3357 | static real t[4160] /* was [65][64] */; |
| 3358 | static integer ib; |
| 3359 | static real ei; |
| 3360 | static integer nb, nh, nx, iws, nbmin, iinfo; |
| 3361 | extern /* Subroutine */ int sgemm_(char *, char *, integer *, integer *, |
| 3362 | integer *, real *, real *, integer *, real *, integer *, real *, |
| 3363 | real *, integer *), strmm_(char *, char *, char *, |
| 3364 | char *, integer *, integer *, real *, real *, integer *, real *, |
| 3365 | integer *), saxpy_(integer *, |
| 3366 | real *, real *, integer *, real *, integer *), sgehd2_(integer *, |
| 3367 | integer *, integer *, real *, integer *, real *, real *, integer * |
| 3368 | ), slahr2_(integer *, integer *, integer *, real *, integer *, |
| 3369 | real *, real *, integer *, real *, integer *), slarfb_(char *, |
| 3370 | char *, char *, char *, integer *, integer *, integer *, real *, |
| 3371 | integer *, real *, integer *, real *, integer *, real *, integer * |
| 3372 | ), xerbla_(char *, integer *); |
| 3373 | extern integer ilaenv_(integer *, char *, char *, integer *, integer *, |
| 3374 | integer *, integer *, ftnlen, ftnlen); |
| 3375 | static integer ldwork, lwkopt; |
| 3376 | static logical lquery; |
| 3377 | |
| 3378 | |
| 3379 | /* |
| 3380 | -- LAPACK routine (version 3.2.1) -- |
| 3381 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 3382 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 3383 | -- April 2009 -- |
| 3384 | |
| 3385 | |
| 3386 | Purpose |
| 3387 | ======= |
| 3388 | |
| 3389 | SGEHRD reduces a real general matrix A to upper Hessenberg form H by |
| 3390 | an orthogonal similarity transformation: Q' * A * Q = H . |
| 3391 | |
| 3392 | Arguments |
| 3393 | ========= |
| 3394 | |
| 3395 | N (input) INTEGER |
| 3396 | The order of the matrix A. N >= 0. |
| 3397 | |
| 3398 | ILO (input) INTEGER |
| 3399 | IHI (input) INTEGER |
| 3400 | It is assumed that A is already upper triangular in rows |
| 3401 | and columns 1:ILO-1 and IHI+1:N. ILO and IHI are normally |
| 3402 | set by a previous call to SGEBAL; otherwise they should be |
| 3403 | set to 1 and N respectively. See Further Details. |
| 3404 | 1 <= ILO <= IHI <= N, if N > 0; ILO=1 and IHI=0, if N=0. |
| 3405 | |
| 3406 | A (input/output) REAL array, dimension (LDA,N) |