--- rpl/lapack/lapack/zla_hercond_c.f 2011/11/21 22:19:50 1.7 +++ rpl/lapack/lapack/zla_hercond_c.f 2016/08/27 15:34:54 1.13 @@ -1,4 +1,4 @@ -*> \brief \b ZLA_HERCOND_C +*> \brief \b ZLA_HERCOND_C computes the infinity norm condition number of op(A)*inv(diag(c)) for Hermitian indefinite matrices. * * =========== DOCUMENTATION =========== * @@ -131,7 +131,7 @@ *> \author Univ. of Colorado Denver *> \author NAG Ltd. * -*> \date November 2011 +*> \date September 2012 * *> \ingroup complex16HEcomputational * @@ -140,10 +140,10 @@ $ LDAF, IPIV, C, CAPPLY, $ INFO, WORK, RWORK ) * -* -- LAPACK computational routine (version 3.4.0) -- +* -- LAPACK computational routine (version 3.4.2) -- * -- LAPACK is a software package provided by Univ. of Tennessee, -- * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- -* November 2011 +* September 2012 * * .. Scalar Arguments .. CHARACTER UPLO @@ -161,7 +161,7 @@ * .. Local Scalars .. INTEGER KASE, I, J DOUBLE PRECISION AINVNM, ANORM, TMP - LOGICAL UP + LOGICAL UP, UPPER COMPLEX*16 ZDUM * .. * .. Local Arrays .. @@ -188,8 +188,15 @@ ZLA_HERCOND_C = 0.0D+0 * INFO = 0 - IF( N.LT.0 ) THEN + UPPER = LSAME( UPLO, 'U' ) + IF( .NOT.UPPER .AND. .NOT.LSAME( UPLO, 'L' ) ) THEN + INFO = -1 + ELSE IF( N.LT.0 ) THEN INFO = -2 + ELSE IF( LDA.LT.MAX( 1, N ) ) THEN + INFO = -4 + ELSE IF( LDAF.LT.MAX( 1, N ) ) THEN + INFO = -6 END IF IF( INFO.NE.0 ) THEN CALL XERBLA( 'ZLA_HERCOND_C', -INFO )