Fix Bug 0073: [Re] Add dummy routine in timer to avoid compiler problem

Bug reported by John Tellefson to LAPACK Mailing List on February 7th.
Dummy routine was removed in LAPACK 3.3.0 because it seems it was working on all machines.
[Re] Incorporation of the dummy routine and incorporate John's suggested cosmectic changes to allow
an easier change in number of operation to perform.

"modified secondtst.f from lapack-3.2.2 with the old
code that makes it reliably report timing loop results. I also added a
variable TOTAL to simplify changes in setting array size and iterations
which affect total flops and made few other changes.

diff --git a/INSTALL/dsecndtst.f b/INSTALL/dsecndtst.f
index f8765d8eeaaf2a65f678ab1d683d8ecb623a1137..f272d37f4fb8b86272affad140e75b62af7a2e28 100644 (file)
--- a/INSTALL/dsecndtst.f
+++ b/INSTALL/dsecndtst.f
@@ -1,16 +1,16 @@
       PROGRAM TEST5
 *
-*  -- LAPACK test routine (version 3.2) --
+*  -- LAPACK test routine (version 3.3.1) --
 *     Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd..
-*     November 2006
+*     February 2010
 *
 *     .. Parameters ..
       INTEGER            NMAX, ITS
-      PARAMETER          ( NMAX = 100, ITS = 5000 )
+      PARAMETER          ( NMAX = 1000, ITS = 50000 )
 *     ..
 *     .. Local Scalars ..
       INTEGER            I, J
-      DOUBLE PRECISION   ALPHA, AVG, T1, T2, TNOSEC
+      DOUBLE PRECISION   ALPHA, AVG, T1, T2, TNOSEC, TOTAL
 *     ..
 *     .. Local Arrays ..
       DOUBLE PRECISION   X( NMAX ), Y( NMAX )
@@ -24,6 +24,8 @@
 *     ..
 *     .. Executable Statements ..
 *
+*    .. Figure TOTAL flops ..
+      TOTAL = DBLE(NMAX) * DBLE(ITS) * 2.0
 *
 *     Initialize X and Y
 *
@@ -33,7 +35,7 @@
    10 CONTINUE
       ALPHA = 0.315D0
 *
-*     Time 1,000,000 DAXPY operations
+*     Time TOTAL SAXPY operations
 *
       T1 = DSECND( )
       DO 30 J = 1, ITS
@@ -43,15 +45,15 @@
          ALPHA = -ALPHA
    30 CONTINUE
       T2 = DSECND( )
-      WRITE( 6, 9999 )T2 - T1
-      IF( T2-T1.GT.0.0D0 ) THEN
-         WRITE( 6, 9998 )1.0D0 / ( T2-T1 )
+      TNOSEC = T2 - T1
+      WRITE( 6, 9999 )TOTAL, TNOSEC
+      IF( TNOSEC.GT.0.0 ) THEN
+         WRITE( 6, 9998 )(TOTAL/1.0D6)/TNOSEC
       ELSE
          WRITE( 6, 9994 )
       END IF
-      TNOSEC = T2 - T1
 *
-*     Time 1,000,000 DAXPY operations with DSECND in the outer loop
+*     Time TOTAL DAXPY operations with DSECND in the outer loop
 *
       T1 = DSECND( )
       DO 50 J = 1, ITS
@@ -62,26 +64,28 @@
          T2 = DSECND( )
    50 CONTINUE
 *
-*     Compute the time in milliseconds used by an average call
+*     Compute the time used in milliseconds used by an average call
 *     to DSECND.
 *
       WRITE( 6, 9997 )T2 - T1
-      AVG = ( ( T2-T1 )-TNOSEC )*1000.D0 / DBLE( ITS )
-      WRITE( 6, 9996 )AVG
+      AVG = ( ( T2-T1 ) - TNOSEC ) * 1000.0D+00/DBLE( ITS )
+      IF( AVG.GT.0.0)
+     $   WRITE( 6, 9996 )AVG
 *
 *     Compute the equivalent number of floating point operations used
 *     by an average call to DSECND.
 *
-      IF( TNOSEC.GT.0.0D0 )
-     $   WRITE( 6, 9995 )1000.D0*AVG / TNOSEC
+      IF(( AVG.GT.0.0 ).AND.( TNOSEC.GT.0.0 ))
+     $   WRITE( 6, 9995 )(AVG/1000) * TOTAL / TNOSEC
 *
- 9999 FORMAT( ' Time for 1,000,000 DAXPY ops  = ', G10.3, ' seconds' )
+ 9999 FORMAT( ' Time for ', G10.3,' DAXPY ops = ', G10.3, ' seconds' )
  9998 FORMAT( ' DAXPY performance rate        = ', G10.3, ' mflops ' )
  9997 FORMAT( ' Including DSECND, time        = ', G10.3, ' seconds' )
  9996 FORMAT( ' Average time for DSECND       = ', G10.3,
      $      ' milliseconds' )
  9995 FORMAT( ' Equivalent floating point ops = ', G10.3, ' ops' )
- 9994 FORMAT( ' *** Error:  Time for operations was zero' )
+ 9994 FORMAT( ' *** Warning:  Time for operations was less or equal',
+     $        ' than zero => timing in TESTING might be dubious' )
       CALL MYSUB(NMAX,X,Y)
       END
       SUBROUTINE MYSUB(N,X,Y)

diff --git a/INSTALL/secondtst.f b/INSTALL/secondtst.f
index ebd09a77046114a5fb32e7f146a8e5aba8ccf413..693924080349c9a50c3414318865eb4fb857ce3a 100644 (file)
--- a/INSTALL/secondtst.f
+++ b/INSTALL/secondtst.f
@@ -1,16 +1,16 @@
       PROGRAM TEST4
 *
-*  -- LAPACK test routine (version 3.2.2) --
+*  -- LAPACK test routine (version 3.3.1) --
 *     Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd..
-*     June 2010
+*     February 2010
 *
 *     .. Parameters ..
       INTEGER            NMAX, ITS
-      PARAMETER          ( NMAX = 1000, ITS = 5000 )
+      PARAMETER          ( NMAX = 1000, ITS = 50000 )
 *     ..
 *     .. Local Scalars ..
       INTEGER            I, J
-      REAL               ALPHA, AVG, T1, T2, TNOSEC
+      REAL               ALPHA, AVG, T1, T2, TNOSEC, TOTAL
 *     ..
 *     .. Local Arrays ..
       REAL               X( NMAX ), Y( NMAX )
@@ -24,6 +24,8 @@
 *     ..
 *     .. Executable Statements ..
 *
+*    .. Figure TOTAL flops ..
+      TOTAL = REAL(NMAX) * REAL(ITS) * 2.0
 *
 *     Initialize X and Y
 *
@@ -33,7 +35,7 @@
    10 CONTINUE
       ALPHA = 0.315
 *
-*     Time 10,000,000 SAXPY operations
+*     Time TOTAL SAXPY operations
 *
       T1 = SECOND( )
       DO 30 J = 1, ITS
@@ -43,15 +45,15 @@
          ALPHA = -ALPHA
    30 CONTINUE
       T2 = SECOND( )
-      WRITE( 6, 9999 )T2 - T1
-      IF( T2-T1.GT.0.0 ) THEN
-         WRITE( 6, 9998 )10.0 / ( T2-T1 )
+      TNOSEC = T2 - T1
+      WRITE( 6, 9999 )TOTAL, TNOSEC
+      IF( TNOSEC.GT.0.0 ) THEN
+         WRITE( 6, 9998 )(TOTAL/1.0E6)/TNOSEC
       ELSE
          WRITE( 6, 9994 )
       END IF
-      TNOSEC = T2 - T1
 *
-*     Time 1,000,000 SAXPY operations with SECOND in the outer loop
+*     Time TOTAL SAXPY operations with SECOND in the outer loop
 *
       T1 = SECOND( )
       DO 50 J = 1, ITS
@@ -74,15 +76,20 @@
 *     by an average call to SECOND.
 *
       IF(( AVG.GT.0.0 ).AND.( TNOSEC.GT.0.0 ))
-     $   WRITE( 6, 9995 )1000.*AVG / TNOSEC
+     $   WRITE( 6, 9995 )(AVG/1000) * TOTAL / TNOSEC
 *
- 9999 FORMAT( ' Time for 10,000,000 SAXPY ops  = ', G10.3, ' seconds' )
- 9998 FORMAT( ' SAXPY performance rate         = ', G10.3, ' mflops ' )
- 9997 FORMAT( ' Including SECOND, time         = ', G10.3, ' seconds' )
- 9996 FORMAT( ' Average time for SECOND        = ', G10.3,
+ 9999 FORMAT( ' Time for ', G10.3,' SAXPY ops = ', G10.3, ' seconds' )
+ 9998 FORMAT( ' SAXPY performance rate        = ', G10.3, ' mflops ' )
+ 9997 FORMAT( ' Including SECOND, time        = ', G10.3, ' seconds' )
+ 9996 FORMAT( ' Average time for SECOND       = ', G10.3,
      $      ' milliseconds' )
- 9995 FORMAT( ' Equivalent floating point ops  = ', G10.3, ' ops' )
+ 9995 FORMAT( ' Equivalent floating point ops = ', G10.3, ' ops' )
  9994 FORMAT( ' *** Warning:  Time for operations was less or equal',
      $        ' than zero => timing in TESTING might be dubious' )
+      CALL MYSUB(NMAX,X,Y)
+      END
+      SUBROUTINE MYSUB(N,X,Y)
+      INTEGER N
+      REAL X(N), Y(N)
+      RETURN
       END
-