diff options
author | langou <langou@users.noreply.github.com> | 2015-06-30 10:27:18 +0000 |
---|---|---|
committer | langou <langou@users.noreply.github.com> | 2015-06-30 10:27:18 +0000 |
commit | a002ad4689cfe9647630a8eac71e24277d823938 (patch) | |
tree | 8127f396de4c38059c55eef2517a1b9fa84b281c | |
parent | ddac61b7d5e0630e5ebfdb8f311bbd111596056e (diff) | |
download | lapack-a002ad4689cfe9647630a8eac71e24277d823938.tar.gz lapack-a002ad4689cfe9647630a8eac71e24277d823938.tar.bz2 lapack-a002ad4689cfe9647630a8eac71e24277d823938.zip |
Revert commit 1474: files ?lansy.f are back as they were in revision 1473.
-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 |