code simplified: redundant computations removed

Author: jip

TESTING/EIG/cget22.f

@@ -125,6 +125,7 @@
 *>          RESULT is REAL array, dimension (2)
 *>          RESULT(1) = | A E  -  E W | / ( |A| |E| ulp )
 *>          RESULT(2) = max | m-norm(E(j)) - 1 | / ( n ulp )
+*>                       j
 *> \endverbatim
 *
 *  Authors:
@@ -299,7 +300,7 @@ SUBROUTINE CGET22( TRANSA, TRANSE, TRANSW, N, A, LDA, E, LDE, W,
          RESULT( 1 ) = ( ERRNRM / ANORM ) / ULP
       ELSE
          IF( ANORM.LT.ONE ) THEN
-            RESULT( 1 ) = ( MIN( ERRNRM, ANORM ) / ANORM ) / ULP
+            RESULT( 1 ) = ONE / ULP
          ELSE
             RESULT( 1 ) = MIN( ERRNRM / ANORM, ONE ) / ULP
          END IF

TESTING/EIG/dget22.f

@@ -149,6 +149,7 @@
 *>          RESULT is DOUBLE PRECISION array, dimension (2)
 *>          RESULT(1) = | A E  -  E W | / ( |A| |E| ulp )
 *>          RESULT(2) = max | m-norm(E(j)) - 1 | / ( n ulp )
+*>                       j
 *> \endverbatim
 *
 *  Authors:
@@ -366,7 +367,7 @@ SUBROUTINE DGET22( TRANSA, TRANSE, TRANSW, N, A, LDA, E, LDE, WR,
          RESULT( 1 ) = ( ERRNRM / ANORM ) / ULP
       ELSE
          IF( ANORM.LT.ONE ) THEN
-            RESULT( 1 ) = ( MIN( ERRNRM, ANORM ) / ANORM ) / ULP
+            RESULT( 1 ) = ONE / ULP
          ELSE
             RESULT( 1 ) = MIN( ERRNRM / ANORM, ONE ) / ULP
          END IF

TESTING/EIG/sget22.f

@@ -149,6 +149,7 @@
 *>          RESULT is REAL array, dimension (2)
 *>          RESULT(1) = | A E  -  E W | / ( |A| |E| ulp )
 *>          RESULT(2) = max | m-norm(E(j)) - 1 | / ( n ulp )
+*>                       j
 *> \endverbatim
 *
 *  Authors:
@@ -366,7 +367,7 @@ SUBROUTINE SGET22( TRANSA, TRANSE, TRANSW, N, A, LDA, E, LDE, WR,
          RESULT( 1 ) = ( ERRNRM / ANORM ) / ULP
       ELSE
          IF( ANORM.LT.ONE ) THEN
-            RESULT( 1 ) = ( MIN( ERRNRM, ANORM ) / ANORM ) / ULP
+            RESULT( 1 ) = ONE / ULP
          ELSE
             RESULT( 1 ) = MIN( ERRNRM / ANORM, ONE ) / ULP
          END IF

TESTING/EIG/zget22.f

@@ -125,6 +125,7 @@
 *>          RESULT is DOUBLE PRECISION array, dimension (2)
 *>          RESULT(1) = | A E  -  E W | / ( |A| |E| ulp )
 *>          RESULT(2) = max | m-norm(E(j)) - 1 | / ( n ulp )
+*>                       j
 *> \endverbatim
 *
 *  Authors:
@@ -299,7 +300,7 @@ SUBROUTINE ZGET22( TRANSA, TRANSE, TRANSW, N, A, LDA, E, LDE, W,
          RESULT( 1 ) = ( ERRNRM / ANORM ) / ULP
       ELSE
          IF( ANORM.LT.ONE ) THEN
-            RESULT( 1 ) = ( MIN( ERRNRM, ANORM ) / ANORM ) / ULP
+            RESULT( 1 ) = ONE / ULP
          ELSE
             RESULT( 1 ) = MIN( ERRNRM / ANORM, ONE ) / ULP
          END IF