diff options
| -rw-r--r-- | SRC/clansy.f | 18 | ||||
| -rw-r--r-- | SRC/dlansy.f | 18 | ||||
| -rw-r--r-- | SRC/slansy.f | 18 | ||||
| -rw-r--r-- | SRC/zlansy.f | 18 |
4 files changed, 36 insertions, 36 deletions
diff --git a/SRC/clansy.f b/SRC/clansy.f index aced505b..86c2af74 100644 --- a/SRC/clansy.f +++ b/SRC/clansy.f @@ -187,24 +187,24 @@ * Find normI(A) ( = norm1(A), since A is symmetric). * VALUE = ZERO - DO 50 I = 1, N - WORK( I ) = ZERO - 50 CONTINUE IF( LSAME( UPLO, 'U' ) ) THEN - DO 70 J = 1, N + DO 60 J = 1, N SUM = ZERO - DO 60 I = 1, J - 1 + DO 50 I = 1, J - 1 ABSA = ABS( A( I, J ) ) SUM = SUM + ABSA WORK( I ) = WORK( I ) + ABSA - 60 CONTINUE + 50 CONTINUE WORK( J ) = SUM + ABS( A( J, J ) ) - 70 CONTINUE - DO 80 I = 1, N + 60 CONTINUE + DO 70 I = 1, N SUM = WORK( I ) IF( VALUE .LT. SUM .OR. SISNAN( SUM ) ) VALUE = SUM - 80 CONTINUE + 70 CONTINUE ELSE + DO 80 I = 1, N + WORK( I ) = ZERO + 80 CONTINUE DO 100 J = 1, N SUM = WORK( J ) + ABS( A( J, J ) ) DO 90 I = J + 1, N diff --git a/SRC/dlansy.f b/SRC/dlansy.f index d9fb2f98..bc70ab8e 100644 --- a/SRC/dlansy.f +++ b/SRC/dlansy.f @@ -185,24 +185,24 @@ * Find normI(A) ( = norm1(A), since A is symmetric). * VALUE = ZERO - DO 50 I = 1, N - WORK( I ) = ZERO - 50 CONTINUE IF( LSAME( UPLO, 'U' ) ) THEN - DO 70 J = 1, N + DO 60 J = 1, N SUM = ZERO - DO 60 I = 1, J - 1 + DO 50 I = 1, J - 1 ABSA = ABS( A( I, J ) ) SUM = SUM + ABSA WORK( I ) = WORK( I ) + ABSA - 60 CONTINUE + 50 CONTINUE WORK( J ) = SUM + ABS( A( J, J ) ) - 70 CONTINUE - DO 80 I = 1, N + 60 CONTINUE + DO 70 I = 1, N SUM = WORK( I ) IF( VALUE .LT. SUM .OR. DISNAN( SUM ) ) VALUE = SUM - 80 CONTINUE + 70 CONTINUE ELSE + DO 80 I = 1, N + WORK( I ) = ZERO + 80 CONTINUE DO 100 J = 1, N SUM = WORK( J ) + ABS( A( J, J ) ) DO 90 I = J + 1, N diff --git a/SRC/slansy.f b/SRC/slansy.f index daba3e87..1640358c 100644 --- a/SRC/slansy.f +++ b/SRC/slansy.f @@ -185,24 +185,24 @@ * Find normI(A) ( = norm1(A), since A is symmetric). * VALUE = ZERO - DO 50 I = 1, N - WORK( I ) = ZERO - 50 CONTINUE IF( LSAME( UPLO, 'U' ) ) THEN - DO 70 J = 1, N + DO 60 J = 1, N SUM = ZERO - DO 60 I = 1, J - 1 + DO 50 I = 1, J - 1 ABSA = ABS( A( I, J ) ) SUM = SUM + ABSA WORK( I ) = WORK( I ) + ABSA - 60 CONTINUE + 50 CONTINUE WORK( J ) = SUM + ABS( A( J, J ) ) - 70 CONTINUE - DO 80 I = 1, N + 60 CONTINUE + DO 70 I = 1, N SUM = WORK( I ) IF( VALUE .LT. SUM .OR. SISNAN( SUM ) ) VALUE = SUM - 80 CONTINUE + 70 CONTINUE ELSE + DO 80 I = 1, N + WORK( I ) = ZERO + 80 CONTINUE DO 100 J = 1, N SUM = WORK( J ) + ABS( A( J, J ) ) DO 90 I = J + 1, N diff --git a/SRC/zlansy.f b/SRC/zlansy.f index 8ebdd6b3..e261a3b9 100644 --- a/SRC/zlansy.f +++ b/SRC/zlansy.f @@ -187,24 +187,24 @@ * Find normI(A) ( = norm1(A), since A is symmetric). * VALUE = ZERO - DO 50 I = 1, N - WORK( I ) = ZERO - 50 CONTINUE IF( LSAME( UPLO, 'U' ) ) THEN - DO 70 J = 1, N + DO 60 J = 1, N SUM = ZERO - DO 60 I = 1, J - 1 + DO 50 I = 1, J - 1 ABSA = ABS( A( I, J ) ) SUM = SUM + ABSA WORK( I ) = WORK( I ) + ABSA - 60 CONTINUE + 50 CONTINUE WORK( J ) = SUM + ABS( A( J, J ) ) - 70 CONTINUE - DO 80 I = 1, N + 60 CONTINUE + DO 70 I = 1, N SUM = WORK( I ) IF( VALUE .LT. SUM .OR. DISNAN( SUM ) ) VALUE = SUM - 80 CONTINUE + 70 CONTINUE ELSE + DO 80 I = 1, N + WORK( I ) = ZERO + 80 CONTINUE DO 100 J = 1, N SUM = WORK( J ) + ABS( A( J, J ) ) DO 90 I = J + 1, N |