version 1.11, 2012/12/14 14:22:48
|
version 1.17, 2018/05/29 07:18:22
|
Line 2
|
Line 2
|
* |
* |
* =========== DOCUMENTATION =========== |
* =========== DOCUMENTATION =========== |
* |
* |
* Online html documentation available at |
* Online html documentation available at |
* http://www.netlib.org/lapack/explore-html/ |
* http://www.netlib.org/lapack/explore-html/ |
* |
* |
*> \htmlonly |
*> \htmlonly |
*> Download ZHSEIN + dependencies |
*> Download ZHSEIN + dependencies |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/zhsein.f"> |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/zhsein.f"> |
*> [TGZ]</a> |
*> [TGZ]</a> |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/zhsein.f"> |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/zhsein.f"> |
*> [ZIP]</a> |
*> [ZIP]</a> |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/zhsein.f"> |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/zhsein.f"> |
*> [TXT]</a> |
*> [TXT]</a> |
*> \endhtmlonly |
*> \endhtmlonly |
* |
* |
* Definition: |
* Definition: |
* =========== |
* =========== |
Line 21
|
Line 21
|
* SUBROUTINE ZHSEIN( SIDE, EIGSRC, INITV, SELECT, N, H, LDH, W, VL, |
* SUBROUTINE ZHSEIN( SIDE, EIGSRC, INITV, SELECT, N, H, LDH, W, VL, |
* LDVL, VR, LDVR, MM, M, WORK, RWORK, IFAILL, |
* LDVL, VR, LDVR, MM, M, WORK, RWORK, IFAILL, |
* IFAILR, INFO ) |
* IFAILR, INFO ) |
* |
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
* CHARACTER EIGSRC, INITV, SIDE |
* CHARACTER EIGSRC, INITV, SIDE |
* INTEGER INFO, LDH, LDVL, LDVR, M, MM, N |
* INTEGER INFO, LDH, LDVL, LDVR, M, MM, N |
Line 33
|
Line 33
|
* COMPLEX*16 H( LDH, * ), VL( LDVL, * ), VR( LDVR, * ), |
* COMPLEX*16 H( LDH, * ), VL( LDVL, * ), VR( LDVR, * ), |
* $ W( * ), WORK( * ) |
* $ W( * ), WORK( * ) |
* .. |
* .. |
* |
* |
* |
* |
*> \par Purpose: |
*> \par Purpose: |
* ============= |
* ============= |
Line 104
|
Line 104
|
*> \verbatim |
*> \verbatim |
*> H is COMPLEX*16 array, dimension (LDH,N) |
*> H is COMPLEX*16 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 220
|
Line 221
|
* Authors: |
* Authors: |
* ======== |
* ======== |
* |
* |
*> \author Univ. of Tennessee |
*> \author Univ. of Tennessee |
*> \author Univ. of California Berkeley |
*> \author Univ. of California Berkeley |
*> \author Univ. of Colorado Denver |
*> \author Univ. of Colorado Denver |
*> \author NAG Ltd. |
*> \author NAG Ltd. |
* |
* |
*> \date November 2011 |
*> \date December 2016 |
* |
* |
*> \ingroup complex16OTHERcomputational |
*> \ingroup complex16OTHERcomputational |
* |
* |
Line 244
|
Line 245
|
$ LDVL, VR, LDVR, MM, M, WORK, RWORK, IFAILL, |
$ LDVL, VR, LDVR, MM, M, WORK, RWORK, IFAILL, |
$ IFAILR, INFO ) |
$ IFAILR, INFO ) |
* |
* |
* -- LAPACK computational routine (version 3.4.0) -- |
* -- LAPACK computational routine (version 3.7.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 |
* December 2016 |
* |
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
CHARACTER EIGSRC, INITV, SIDE |
CHARACTER EIGSRC, INITV, SIDE |
Line 276
|
Line 277
|
COMPLEX*16 CDUM, WK |
COMPLEX*16 CDUM, WK |
* .. |
* .. |
* .. External Functions .. |
* .. External Functions .. |
LOGICAL LSAME |
LOGICAL LSAME, DISNAN |
DOUBLE PRECISION DLAMCH, ZLANHS |
DOUBLE PRECISION DLAMCH, ZLANHS |
EXTERNAL LSAME, DLAMCH, ZLANHS |
EXTERNAL LSAME, DLAMCH, ZLANHS, DISNAN |
* .. |
* .. |
* .. External Subroutines .. |
* .. External Subroutines .. |
EXTERNAL XERBLA, ZLAEIN |
EXTERNAL XERBLA, ZLAEIN |
Line 399
|
Line 400
|
* has not ben computed before. |
* has not ben computed before. |
* |
* |
HNORM = ZLANHS( 'I', KR-KL+1, H( KL, KL ), LDH, RWORK ) |
HNORM = ZLANHS( 'I', KR-KL+1, H( KL, KL ), LDH, RWORK ) |
IF( HNORM.GT.RZERO ) THEN |
IF( DISNAN( HNORM ) ) THEN |
|
INFO = -6 |
|
RETURN |
|
ELSE IF( HNORM.GT.RZERO ) THEN |
EPS3 = HNORM*ULP |
EPS3 = HNORM*ULP |
ELSE |
ELSE |
EPS3 = SMLNUM |
EPS3 = SMLNUM |