--- rpl/lapack/lapack/zla_porcond_c.f 2011/11/21 20:43:14 1.6
+++ rpl/lapack/lapack/zla_porcond_c.f 2018/05/29 07:18:23 1.16
@@ -1,27 +1,27 @@
-*> \brief \b ZLA_PORCOND_C
+*> \brief \b ZLA_PORCOND_C computes the infinity norm condition number of op(A)*inv(diag(c)) for Hermitian positive-definite 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_PORCOND_C + dependencies
-*>
-*> [TGZ]
-*>
-*> [ZIP]
-*>
+*> Download ZLA_PORCOND_C + dependencies
+*>
+*> [TGZ]
+*>
+*> [ZIP]
+*>
*> [TXT]
-*> \endhtmlonly
+*> \endhtmlonly
*
* Definition:
* ===========
*
-* DOUBLE PRECISION FUNCTION ZLA_PORCOND_C( UPLO, N, A, LDA, AF,
+* DOUBLE PRECISION FUNCTION ZLA_PORCOND_C( UPLO, N, A, LDA, AF,
* LDAF, C, CAPPLY, INFO,
* WORK, RWORK )
-*
+*
* .. Scalar Arguments ..
* CHARACTER UPLO
* LOGICAL CAPPLY
@@ -31,7 +31,7 @@
* COMPLEX*16 A( LDA, * ), AF( LDAF, * ), WORK( * )
* DOUBLE PRECISION C( * ), RWORK( * )
* ..
-*
+*
*
*> \par Purpose:
* =============
@@ -118,24 +118,24 @@
* 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 complex16POcomputational
*
* =====================================================================
- DOUBLE PRECISION FUNCTION ZLA_PORCOND_C( UPLO, N, A, LDA, AF,
+ DOUBLE PRECISION FUNCTION ZLA_PORCOND_C( UPLO, N, A, LDA, AF,
$ LDAF, 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
@@ -153,7 +153,7 @@
INTEGER KASE
DOUBLE PRECISION AINVNM, ANORM, TMP
INTEGER I, J
- LOGICAL UP
+ LOGICAL UP, UPPER
COMPLEX*16 ZDUM
* ..
* .. Local Arrays ..
@@ -180,8 +180,15 @@
ZLA_PORCOND_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_PORCOND_C', -INFO )