summaryrefslogtreecommitdiff
path: root/SRC/clansb.f
diff options
context:
space:
mode:
authorjames <james@8a072113-8704-0410-8d35-dd094bca7971>2012-07-16 14:21:07 +0000
committerjames <james@8a072113-8704-0410-8d35-dd094bca7971>2012-07-16 14:21:07 +0000
commitb0ccd072d7381538a5947bc93908abeabaaf83d8 (patch)
treefc3e26f33d620d6968d56441ca92a256a071c3e6 /SRC/clansb.f
parent8989738d77f02c7b50839e69ab8a2ed8aac75024 (diff)
downloadlapack-b0ccd072d7381538a5947bc93908abeabaaf83d8.tar.gz
lapack-b0ccd072d7381538a5947bc93908abeabaaf83d8.tar.bz2
lapack-b0ccd072d7381538a5947bc93908abeabaaf83d8.zip
matrix norms now consistently return nan when input matrix contains a nan;
previously, this was (usually) only the case for the 'F' norm
Diffstat (limited to 'SRC/clansb.f')
-rw-r--r--SRC/clansb.f15
1 files changed, 9 insertions, 6 deletions
diff --git a/SRC/clansb.f b/SRC/clansb.f
index efc46dae..4bb46103 100644
--- a/SRC/clansb.f
+++ b/SRC/clansb.f
@@ -155,8 +155,8 @@
REAL ABSA, SCALE, SUM, VALUE
* ..
* .. External Functions ..
- LOGICAL LSAME
- EXTERNAL LSAME
+ LOGICAL LSAME, SISNAN
+ EXTERNAL LSAME, SISNAN
* ..
* .. External Subroutines ..
EXTERNAL CLASSQ
@@ -176,13 +176,15 @@
IF( LSAME( UPLO, 'U' ) ) THEN
DO 20 J = 1, N
DO 10 I = MAX( K+2-J, 1 ), K + 1
- VALUE = MAX( VALUE, ABS( AB( I, J ) ) )
+ SUM = ABS( AB( I, J ) )
+ IF( VALUE .LT. SUM .OR. SISNAN( SUM ) ) VALUE = SUM
10 CONTINUE
20 CONTINUE
ELSE
DO 40 J = 1, N
DO 30 I = 1, MIN( N+1-J, K+1 )
- VALUE = MAX( VALUE, ABS( AB( I, J ) ) )
+ SUM = ABS( AB( I, J ) )
+ IF( VALUE .LT. SUM .OR. SISNAN( SUM ) ) VALUE = SUM
30 CONTINUE
40 CONTINUE
END IF
@@ -204,7 +206,8 @@
WORK( J ) = SUM + ABS( AB( K+1, J ) )
60 CONTINUE
DO 70 I = 1, N
- VALUE = MAX( VALUE, WORK( I ) )
+ SUM = WORK( I )
+ IF( VALUE .LT. SUM .OR. SISNAN( SUM ) ) VALUE = SUM
70 CONTINUE
ELSE
DO 80 I = 1, N
@@ -218,7 +221,7 @@
SUM = SUM + ABSA
WORK( I ) = WORK( I ) + ABSA
90 CONTINUE
- VALUE = MAX( VALUE, SUM )
+ IF( VALUE .LT. SUM .OR. SISNAN( SUM ) ) VALUE = SUM
100 CONTINUE
END IF
ELSE IF( ( LSAME( NORM, 'F' ) ) .OR. ( LSAME( NORM, 'E' ) ) ) THEN