--- rpl/lapack/lapack/zla_syrcond_c.f 2011/11/21 22:19:50 1.7
+++ rpl/lapack/lapack/zla_syrcond_c.f 2018/05/29 07:18:24 1.16
@@ -1,19 +1,19 @@
-*> \brief \b ZLA_SYRCOND_C
+*> \brief \b ZLA_SYRCOND_C computes the infinity norm condition number of op(A)*inv(diag(c)) for symmetric indefinite matrices.
*
* =========== DOCUMENTATION ===========
*
-* Online html documentation available at
-* http://www.netlib.org/lapack/explore-html/
+* Online html documentation available at
+* http://www.netlib.org/lapack/explore-html/
*
*> \htmlonly
-*> Download ZLA_SYRCOND_C + dependencies
-*>
-*> [TGZ]
-*>
-*> [ZIP]
-*>
+*> Download ZLA_SYRCOND_C + dependencies
+*>
+*> [TGZ]
+*>
+*> [ZIP]
+*>
*> [TXT]
-*> \endhtmlonly
+*> \endhtmlonly
*
* Definition:
* ===========
@@ -21,7 +21,7 @@
* DOUBLE PRECISION FUNCTION ZLA_SYRCOND_C( UPLO, N, A, LDA, AF,
* LDAF, IPIV, C, CAPPLY,
* INFO, WORK, RWORK )
-*
+*
* .. Scalar Arguments ..
* CHARACTER UPLO
* LOGICAL CAPPLY
@@ -32,7 +32,7 @@
* COMPLEX*16 A( LDA, * ), AF( LDAF, * ), WORK( * )
* DOUBLE PRECISION C( * ), RWORK( * )
* ..
-*
+*
*
*> \par Purpose:
* =============
@@ -126,12 +126,12 @@
* Authors:
* ========
*
-*> \author Univ. of Tennessee
-*> \author Univ. of California Berkeley
-*> \author Univ. of Colorado Denver
-*> \author NAG Ltd.
+*> \author Univ. of Tennessee
+*> \author Univ. of California Berkeley
+*> \author Univ. of Colorado Denver
+*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date December 2016
*
*> \ingroup complex16SYcomputational
*
@@ -140,10 +140,10 @@
$ LDAF, IPIV, C, CAPPLY,
$ INFO, WORK, RWORK )
*
-* -- LAPACK computational routine (version 3.4.0) --
+* -- LAPACK computational routine (version 3.7.0) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* December 2016
*
* .. Scalar Arguments ..
CHARACTER UPLO
@@ -162,7 +162,7 @@
INTEGER KASE
DOUBLE PRECISION AINVNM, ANORM, TMP
INTEGER I, J
- LOGICAL UP
+ LOGICAL UP, UPPER
COMPLEX*16 ZDUM
* ..
* .. Local Arrays ..
@@ -189,8 +189,15 @@
ZLA_SYRCOND_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_SYRCOND_C', -INFO )