version 1.1, 2015/11/26 11:44:16
|
version 1.7, 2023/08/07 08:38:51
|
Line 211
|
Line 211
|
*> \author Univ. of Colorado Denver |
*> \author Univ. of Colorado Denver |
*> \author NAG Ltd. |
*> \author NAG Ltd. |
* |
* |
*> \date January 2015 |
|
* |
|
*> \ingroup doubleOTHERcomputational |
*> \ingroup doubleOTHERcomputational |
* |
* |
*> \par Further Details: |
*> \par Further Details: |
Line 220
|
Line 218
|
*> |
*> |
*> \verbatim |
*> \verbatim |
*> |
*> |
*> This routine reduces A to Hessenberg form and maintains B in |
*> This routine reduces A to Hessenberg form and maintains B in triangular form |
*> using a blocked variant of Moler and Stewart's original algorithm, |
*> using a blocked variant of Moler and Stewart's original algorithm, |
*> as described by Kagstrom, Kressner, Quintana-Orti, and Quintana-Orti |
*> as described by Kagstrom, Kressner, Quintana-Orti, and Quintana-Orti |
*> (BIT 2008). |
*> (BIT 2008). |
Line 230
|
Line 228
|
SUBROUTINE DGGHD3( COMPQ, COMPZ, N, ILO, IHI, A, LDA, B, LDB, Q, |
SUBROUTINE DGGHD3( COMPQ, COMPZ, N, ILO, IHI, A, LDA, B, LDB, Q, |
$ LDQ, Z, LDZ, WORK, LWORK, INFO ) |
$ LDQ, Z, LDZ, WORK, LWORK, INFO ) |
* |
* |
* -- LAPACK computational routine (version 3.6.0) -- |
* -- LAPACK computational routine -- |
* -- LAPACK is a software package provided by Univ. of Tennessee, -- |
* -- LAPACK is a software package provided by Univ. of Tennessee, -- |
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
* January 2015 |
|
* |
* |
IMPLICIT NONE |
IMPLICIT NONE |
* |
* |
Line 266
|
Line 263
|
EXTERNAL ILAENV, LSAME |
EXTERNAL ILAENV, LSAME |
* .. |
* .. |
* .. External Subroutines .. |
* .. External Subroutines .. |
EXTERNAL DGGHRD, DLARTG, DLASET, DORM22, DROT, XERBLA |
EXTERNAL DGGHRD, DLARTG, DLASET, DORM22, DROT, DGEMM, |
|
$ DGEMV, DTRMV, DLACPY, XERBLA |
* .. |
* .. |
* .. Intrinsic Functions .. |
* .. Intrinsic Functions .. |
INTRINSIC DBLE, MAX |
INTRINSIC DBLE, MAX |
Line 277
|
Line 275
|
* |
* |
INFO = 0 |
INFO = 0 |
NB = ILAENV( 1, 'DGGHD3', ' ', N, ILO, IHI, -1 ) |
NB = ILAENV( 1, 'DGGHD3', ' ', N, ILO, IHI, -1 ) |
LWKOPT = 6*N*NB |
LWKOPT = MAX( 6*N*NB, 1 ) |
WORK( 1 ) = DBLE( LWKOPT ) |
WORK( 1 ) = DBLE( LWKOPT ) |
INITQ = LSAME( COMPQ, 'I' ) |
INITQ = LSAME( COMPQ, 'I' ) |
WANTQ = INITQ .OR. LSAME( COMPQ, 'V' ) |
WANTQ = INITQ .OR. LSAME( COMPQ, 'V' ) |