From a002ad4689cfe9647630a8eac71e24277d823938 Mon Sep 17 00:00:00 2001 From: langou Date: Tue, 30 Jun 2015 10:27:18 +0000 Subject: [PATCH] Revert commit 1474: files ?lansy.f are back as they were in revision 1473. --- SRC/clansy.f | 18 +++++++++--------- SRC/dlansy.f | 18 +++++++++--------- SRC/slansy.f | 18 +++++++++--------- SRC/zlansy.f | 18 +++++++++--------- 4 files changed, 36 insertions(+), 36 deletions(-) diff --git a/SRC/clansy.f b/SRC/clansy.f index aced505..86c2af7 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 d9fb2f9..bc70ab8 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 daba3e8..1640358 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 8ebdd6b..e261a3b 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 -- 2.7.4