version 1.8, 2011/11/21 20:42:53
|
version 1.14, 2016/08/27 15:34:24
|
Line 108
|
Line 108
|
*> \verbatim |
*> \verbatim |
*> H is DOUBLE PRECISION array, dimension (LDH,N) |
*> H is DOUBLE PRECISION array, dimension (LDH,N) |
*> The upper Hessenberg matrix H. |
*> The upper Hessenberg matrix H. |
|
*> If a NaN is detected in H, the routine will return with INFO=-6. |
*> \endverbatim |
*> \endverbatim |
*> |
*> |
*> \param[in] LDH |
*> \param[in] LDH |
Line 243
|
Line 244
|
*> \author Univ. of Colorado Denver |
*> \author Univ. of Colorado Denver |
*> \author NAG Ltd. |
*> \author NAG Ltd. |
* |
* |
*> \date November 2011 |
*> \date November 2013 |
* |
* |
*> \ingroup doubleOTHERcomputational |
*> \ingroup doubleOTHERcomputational |
* |
* |
Line 262
|
Line 263
|
$ VL, LDVL, VR, LDVR, MM, M, WORK, IFAILL, |
$ VL, LDVL, VR, LDVR, MM, M, WORK, IFAILL, |
$ IFAILR, INFO ) |
$ IFAILR, INFO ) |
* |
* |
* -- LAPACK computational routine (version 3.4.0) -- |
* -- LAPACK computational routine (version 3.5.0) -- |
* -- 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..-- |
* November 2011 |
* November 2013 |
* |
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
CHARACTER EIGSRC, INITV, SIDE |
CHARACTER EIGSRC, INITV, SIDE |
Line 291
|
Line 292
|
$ WKR |
$ WKR |
* .. |
* .. |
* .. External Functions .. |
* .. External Functions .. |
LOGICAL LSAME |
LOGICAL LSAME, DISNAN |
DOUBLE PRECISION DLAMCH, DLANHS |
DOUBLE PRECISION DLAMCH, DLANHS |
EXTERNAL LSAME, DLAMCH, DLANHS |
EXTERNAL LSAME, DLAMCH, DLANHS, DISNAN |
* .. |
* .. |
* .. External Subroutines .. |
* .. External Subroutines .. |
EXTERNAL DLAEIN, XERBLA |
EXTERNAL DLAEIN, XERBLA |
Line 423
|
Line 424
|
* has not ben computed before. |
* has not ben computed before. |
* |
* |
HNORM = DLANHS( 'I', KR-KL+1, H( KL, KL ), LDH, WORK ) |
HNORM = DLANHS( 'I', KR-KL+1, H( KL, KL ), LDH, WORK ) |
IF( HNORM.GT.ZERO ) THEN |
IF( DISNAN( HNORM ) ) THEN |
|
INFO = -6 |
|
RETURN |
|
ELSE IF( HNORM.GT.ZERO ) THEN |
EPS3 = HNORM*ULP |
EPS3 = HNORM*ULP |
ELSE |
ELSE |
EPS3 = SMLNUM |
EPS3 = SMLNUM |