diff options
author | julie <julielangou@users.noreply.github.com> | 2012-04-12 00:56:00 +0000 |
---|---|---|
committer | julie <julielangou@users.noreply.github.com> | 2012-04-12 00:56:00 +0000 |
commit | d1c0c25cd21baff8edfe9273b2f975ac81da2703 (patch) | |
tree | f8e44c4112fca053cb04b9cc2d07fa905ef28c85 /SRC/zla_hercond_c.f | |
parent | 21199982e23bb604c3a4e61f8dc3b5db69009039 (diff) | |
download | lapack-d1c0c25cd21baff8edfe9273b2f975ac81da2703.tar.gz lapack-d1c0c25cd21baff8edfe9273b2f975ac81da2703.tar.bz2 lapack-d1c0c25cd21baff8edfe9273b2f975ac81da2703.zip |
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....
Diffstat (limited to 'SRC/zla_hercond_c.f')
-rw-r--r-- | SRC/zla_hercond_c.f | 11 |
1 files changed, 9 insertions, 2 deletions
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 ) |