From d1c0c25cd21baff8edfe9273b2f975ac81da2703 Mon Sep 17 00:00:00 2001 From: julie Date: Thu, 12 Apr 2012 00:56:00 +0000 Subject: Follow up to rev r1130 Fix bug bug0088 reported by Mike Pont from NAG on the forum (see http://icl.cs.utk.edu/lapack-forum/viewtopic.php?f=13&t=2893) Actually there were a lot of problems regarding arguments checking. I tried to correct most of them. Apply the fix propose to all x[he/sy]rfsx.f routines - Use IGNORE_CWISE as suggested to prevent use of unitialize variable PARAMS I also appied the previous TESTING fix to all routines INFO has new value in ZHERFSX (see description of INFO between ZHESVX and ZHESVXX) This is set on line 634 (IF ( INFO .LE. N ) INFO = N + J) of zherfsx.f And this is not handled by the testing LIN/zdrvhex.f I just add .AND. INFO.LE.N at line 638 to avoid raising an error when INFO = N + J At the moment, I would recommand a further look at those routines. ZHE, ZSY led to 182 Tests failing to pass the threshold and the same for complex better than before, but still.... --- SRC/zla_hercond_c.f | 11 +++++++++-- 1 file changed, 9 insertions(+), 2 deletions(-) (limited to 'SRC/zla_hercond_c.f') diff --git a/SRC/zla_hercond_c.f b/SRC/zla_hercond_c.f index 77deaedf..586cf9b6 100644 --- a/SRC/zla_hercond_c.f +++ b/SRC/zla_hercond_c.f @@ -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 ) -- cgit v1.2.3