Makes NaN an illegal value for ANORM on input of GECON. Let NaNs prop…

…agate in the computation of RCOND
Reference-LAPACK · Nov 8, 2023 · 49bfa63 · 49bfa63
1 parent 09abbd0
commit 49bfa63
Show file tree

Hide file tree

Showing 4 changed files with 20 additions and 16 deletions.
diff --git a/SRC/cgecon.f b/SRC/cgecon.f
@@ -105,7 +105,9 @@
 *> \verbatim
 *>          INFO is INTEGER
 *>          = 0:  successful exit
-*>          < 0:  if INFO = -i, the i-th argument had an illegal value
+*>          < 0:  if INFO = -i, the i-th argument had an illegal value.
+*>                NaNs are illegal inputs for ANORM.
+*>                Infinity values are valid inputs for ANORM.
 *>          = 1:  RCOND = NaN.
 *> \endverbatim
 *
@@ -201,7 +203,7 @@ SUBROUTINE CGECON( NORM, N, A, LDA, ANORM, RCOND, WORK, RWORK,
          RETURN
       ELSE IF( SISNAN( ANORM ) ) THEN
          RCOND = ANORM
-         INFO = 1
+         INFO = -6
          RETURN
       END IF
 *
@@ -260,8 +262,7 @@ SUBROUTINE CGECON( NORM, N, A, LDA, ANORM, RCOND, WORK, RWORK,
 *
 *     Compute the estimate of the reciprocal condition number.
 *
-      IF( AINVNM.NE.ZERO )
-     $   RCOND = ( ONE / AINVNM ) / ANORM
+      RCOND = ( ONE / AINVNM ) / ANORM
 *
 *     Check for NaNs
 *

diff --git a/SRC/dgecon.f b/SRC/dgecon.f
@@ -105,7 +105,9 @@
 *> \verbatim
 *>          INFO is INTEGER
 *>          = 0:  successful exit
-*>          < 0:  if INFO = -i, the i-th argument had an illegal value
+*>          < 0:  if INFO = -i, the i-th argument had an illegal value.
+*>                NaNs are illegal inputs for ANORM.
+*>                Infinity values are valid inputs for ANORM.
 *>          = 1:  RCOND = NaN.
 *> \endverbatim
 *
@@ -194,7 +196,7 @@ SUBROUTINE DGECON( NORM, N, A, LDA, ANORM, RCOND, WORK, IWORK,
          RETURN
       ELSE IF( DISNAN( ANORM ) ) THEN
          RCOND = ANORM
-         INFO = 1
+         INFO = -6
          RETURN
       END IF
 *
@@ -252,8 +254,7 @@ SUBROUTINE DGECON( NORM, N, A, LDA, ANORM, RCOND, WORK, IWORK,
 *
 *     Compute the estimate of the reciprocal condition number.
 *
-      IF( AINVNM.NE.ZERO )
-     $   RCOND = ( ONE / AINVNM ) / ANORM
+      RCOND = ( ONE / AINVNM ) / ANORM
 *
 *     Check for NaNs
 *

diff --git a/SRC/sgecon.f b/SRC/sgecon.f
@@ -105,7 +105,9 @@
 *> \verbatim
 *>          INFO is INTEGER
 *>          = 0:  successful exit
-*>          < 0:  if INFO = -i, the i-th argument had an illegal value
+*>          < 0:  if INFO = -i, the i-th argument had an illegal value.
+*>                NaNs are illegal inputs for ANORM.
+*>                Infinity values are valid inputs for ANORM.
 *>          = 1:  RCOND = NaN.
 *> \endverbatim
 *
@@ -194,7 +196,7 @@ SUBROUTINE SGECON( NORM, N, A, LDA, ANORM, RCOND, WORK, IWORK,
          RETURN
       ELSE IF( SISNAN( ANORM ) ) THEN
          RCOND = ANORM
-         INFO = 1
+         INFO = -6
          RETURN
       END IF
 *
@@ -252,8 +254,7 @@ SUBROUTINE SGECON( NORM, N, A, LDA, ANORM, RCOND, WORK, IWORK,
 *
 *     Compute the estimate of the reciprocal condition number.
 *
-      IF( AINVNM.NE.ZERO )
-     $   RCOND = ( ONE / AINVNM ) / ANORM
+      RCOND = ( ONE / AINVNM ) / ANORM
 *
 *     Check for NaNs
 *

diff --git a/SRC/zgecon.f b/SRC/zgecon.f
@@ -105,7 +105,9 @@
 *> \verbatim
 *>          INFO is INTEGER
 *>          = 0:  successful exit
-*>          < 0:  if INFO = -i, the i-th argument had an illegal value
+*>          < 0:  if INFO = -i, the i-th argument had an illegal value.
+*>                NaNs are illegal inputs for ANORM.
+*>                Infinity values are valid inputs for ANORM.
 *>          = 1:  RCOND = NaN.
 *> \endverbatim
 *
@@ -201,7 +203,7 @@ SUBROUTINE ZGECON( NORM, N, A, LDA, ANORM, RCOND, WORK, RWORK,
          RETURN
       ELSE IF( DISNAN( ANORM ) ) THEN
          RCOND = ANORM
-         INFO = 1
+         INFO = -6
          RETURN
       END IF
 *
@@ -260,8 +262,7 @@ SUBROUTINE ZGECON( NORM, N, A, LDA, ANORM, RCOND, WORK, RWORK,
 *
 *     Compute the estimate of the reciprocal condition number.
 *
-      IF( AINVNM.NE.ZERO )
-     $   RCOND = ( ONE / AINVNM ) / ANORM
+      RCOND = ( ONE / AINVNM ) / ANORM
 *
 *     Check for NaNs
 *