code fixed: take absolute value to make comparison correct

jip · jip · commit 4bcc555d318d · 2020-06-02T22:26:49.000+10:00
diff --git a/TESTING/EIG/cget52.f b/TESTING/EIG/cget52.f
@@ -277,7 +277,7 @@ SUBROUTINE CGET52( LEFT, N, A, LDA, B, LDB, E, LDE, ALPHA, BETA,
          DO 20 J = 1, N
             TEMP1 = MAX( TEMP1, ABS1( E( J, JVEC ) ) )
    20    CONTINUE
-         ENRMER = MAX( ENRMER, TEMP1-ONE )
+         ENRMER = MAX( ENRMER, ABS( TEMP1-ONE ) )
    30 CONTINUE
 *
 *     Compute RESULT(2) : the normalization error in E.
diff --git a/TESTING/EIG/dget52.f b/TESTING/EIG/dget52.f
@@ -360,14 +360,14 @@ SUBROUTINE DGET52( LEFT, N, A, LDA, B, LDB, E, LDE, ALPHAR,
                DO 20 J = 1, N
                   TEMP1 = MAX( TEMP1, ABS( E( J, JVEC ) ) )
    20          CONTINUE
-               ENRMER = MAX( ENRMER, TEMP1-ONE )
+               ENRMER = MAX( ENRMER, ABS( TEMP1-ONE ) )
             ELSE
                ILCPLX = .TRUE.
                DO 30 J = 1, N
                   TEMP1 = MAX( TEMP1, ABS( E( J, JVEC ) )+
      $                    ABS( E( J, JVEC+1 ) ) )
    30          CONTINUE
-               ENRMER = MAX( ENRMER, TEMP1-ONE )
+               ENRMER = MAX( ENRMER, ABS( TEMP1-ONE ) )
             END IF
          END IF
    40 CONTINUE
diff --git a/TESTING/EIG/sget52.f b/TESTING/EIG/sget52.f
@@ -360,14 +360,14 @@ SUBROUTINE SGET52( LEFT, N, A, LDA, B, LDB, E, LDE, ALPHAR,
                DO 20 J = 1, N
                   TEMP1 = MAX( TEMP1, ABS( E( J, JVEC ) ) )
    20          CONTINUE
-               ENRMER = MAX( ENRMER, TEMP1-ONE )
+               ENRMER = MAX( ENRMER, ABS( TEMP1-ONE ) )
             ELSE
                ILCPLX = .TRUE.
                DO 30 J = 1, N
                   TEMP1 = MAX( TEMP1, ABS( E( J, JVEC ) )+
      $                    ABS( E( J, JVEC+1 ) ) )
    30          CONTINUE
-               ENRMER = MAX( ENRMER, TEMP1-ONE )
+               ENRMER = MAX( ENRMER, ABS( TEMP1-ONE ) )
             END IF
          END IF
    40 CONTINUE
diff --git a/TESTING/EIG/zget52.f b/TESTING/EIG/zget52.f
@@ -278,7 +278,7 @@ SUBROUTINE ZGET52( LEFT, N, A, LDA, B, LDB, E, LDE, ALPHA, BETA,
          DO 20 J = 1, N
             TEMP1 = MAX( TEMP1, ABS1( E( J, JVEC ) ) )
    20    CONTINUE
-         ENRMER = MAX( ENRMER, TEMP1-ONE )
+         ENRMER = MAX( ENRMER, ABS( TEMP1-ONE ) )
    30 CONTINUE
 *
 *     Compute RESULT(2) : the normalization error in E.

Original file line number	Diff line number	Diff line change
`@@ -277,7 +277,7 @@ SUBROUTINE CGET52( LEFT, N, A, LDA, B, LDB, E, LDE, ALPHA, BETA,`
`277`	`277`	`DO 20 J = 1, N`
`278`	`278`	`TEMP1 = MAX( TEMP1, ABS1( E( J, JVEC ) ) )`
`279`	`279`	`20 CONTINUE`
`280`		`- ENRMER = MAX( ENRMER, TEMP1-ONE )`
	`280`	`+ ENRMER = MAX( ENRMER, ABS( TEMP1-ONE ) )`
`281`	`281`	`30 CONTINUE`
`282`	`282`	`*`
`283`	`283`	`* Compute RESULT(2) : the normalization error in E.`
Original file line number	Diff line number	Diff line change
`@@ -278,7 +278,7 @@ SUBROUTINE ZGET52( LEFT, N, A, LDA, B, LDB, E, LDE, ALPHA, BETA,`
`278`	`278`	`DO 20 J = 1, N`
`279`	`279`	`TEMP1 = MAX( TEMP1, ABS1( E( J, JVEC ) ) )`
`280`	`280`	`20 CONTINUE`
`281`		`- ENRMER = MAX( ENRMER, TEMP1-ONE )`
	`281`	`+ ENRMER = MAX( ENRMER, ABS( TEMP1-ONE ) )`
`282`	`282`	`30 CONTINUE`
`283`	`283`	`*`
`284`	`284`	`* Compute RESULT(2) : the normalization error in E.`