sltimer.f

Go to the documentation of this file.

      SUBROUTINE slboot()
*
*  -- ScaLAPACK tools routine (version 1.7) --
*     University of Tennessee, Knoxville, Oak Ridge National Laboratory,
*     and University of California, Berkeley.
*     May 1, 1997
*
*  Purpose
*  =======
*
*  SLBOOT (re)sets all timers to 0, and enables SLtimer.
*
*  =====================================================================
*
*     .. Parameters ..
      INTEGER            NTIMER
      parameter( ntimer = 64 )
      DOUBLE PRECISION   STARTFLAG, ZERO
      parameter( startflag = -5.0d+0, zero = 0.0d+0 )
*     ..
*     .. Local Scalars ..
      INTEGER            I
*     ..
*     .. Common Blocks ..
      LOGICAL            DISABLED
      DOUBLE PRECISION   CPUSEC( NTIMER ), CPUSTART( NTIMER ),
     $                   WALLSEC( NTIMER ), WALLSTART( NTIMER )
      COMMON /sltimer00/ cpusec, wallsec, cpustart, wallstart, disabled
*     ..
*     .. Executable Statements ..
*
      disabled = .false.
      DO 10 i = 1, ntimer
         cpusec( i )  = zero
         wallsec( i ) = zero
         cpustart( i )  = startflag
         wallstart( i ) = startflag
   10 CONTINUE
*
      RETURN
*
*     End of SLBOOT
*
      END
*
      SUBROUTINE sltimer( I )
*
*  -- ScaLAPACK tools routine (version 1.7) --
*     University of Tennessee, Knoxville, Oak Ridge National Laboratory,
*     and University of California, Berkeley.
*     May 1, 1997
*
*     .. Scalar Arguments ..
      INTEGER            I
*     ..
*
*  Purpose
*  =======
*
*  SLtimer provides a "stopwatch" functionality cpu/wall timer
*  (in seconds).  Up to 64 separate timers can be functioning at once.
*  The first call starts the timer, and the second stops it.  This
*  routine can be disenabled, so that calls to the timer are ignored.
*  This feature can be used to make sure certain sections of code do
*  not affect timings, even if they call routines which have SLtimer
*  calls in them.
*
*  Arguments
*  =========
*
*  I       (global input) INTEGER
*          The timer to stop/start.
*
*  =====================================================================
*
*     .. Parameters ..
      INTEGER            NTIMER
      parameter( ntimer = 64 )
      DOUBLE PRECISION   STARTFLAG
      parameter( startflag = -5.0d+0 )
*     ..
*     .. External Functions ..
      DOUBLE PRECISION   DCPUTIME00, DWALLTIME00
      EXTERNAL           dcputime00, dwalltime00
*     ..
*     .. Common Blocks ..
      LOGICAL            DISABLED
      DOUBLE PRECISION   CPUSEC( NTIMER ), CPUSTART( NTIMER ),
     $                   WALLSEC( NTIMER ), WALLSTART( NTIMER )
      COMMON /sltimer00/ cpusec, wallsec, cpustart, wallstart, disabled
*     ..
*     .. Executable Statements ..
*
*     If timing disabled, return
*
      IF( disabled )
     $   RETURN
*
      IF( wallstart( i ).EQ.startflag ) THEN
*
*        If timer has not been started, start it
*
         wallstart( i ) = dwalltime00()
         cpustart( i )  = dcputime00()
*
      ELSE
*
*        Stop timer and add interval to count
*
         cpusec( i ) = cpusec( i ) + dcputime00() - cpustart( i )
         wallsec( i ) = wallsec( i ) + dwalltime00() - wallstart( i )
         wallstart( i ) = startflag
*
      END IF
*
      RETURN
*
*     End of SLTIMER
*
      END
*
      SUBROUTINE slenable()
*
*  -- ScaLAPACK tools routine (version 1.7) --
*     University of Tennessee, Knoxville, Oak Ridge National Laboratory,
*     and University of California, Berkeley.
*     May 1, 1997
*
*  Purpose
*  =======
*
*  SLENABLE sets it so calls to SLtimer are not ignored.
*
*  =====================================================================
*
*     .. Parameters ..
      INTEGER            NTIMER
      parameter( ntimer = 64 )
*     ..
*     .. Common Blocks ..
      LOGICAL            DISABLED
      DOUBLE PRECISION   CPUSEC( NTIMER ), CPUSTART( NTIMER ),
     $                   WALLSEC( NTIMER ), WALLSTART( NTIMER )
      COMMON /sltimer00/ cpusec, wallsec, cpustart, wallstart, disabled
*     ..
*     .. Executable Statements ..
*
      disabled = .false.
*
      RETURN
*
      END
*
      SUBROUTINE sldisable()
*
*  -- ScaLAPACK tools routine (version 1.7) --
*     University of Tennessee, Knoxville, Oak Ridge National Laboratory,
*     and University of California, Berkeley.
*     May 1, 1997
*
*  Purpose
*  =======
*
*  SLDISABLE sets it so calls to SLTIMER are ignored.
*
*  =====================================================================
*
*     .. Parameters ..
      INTEGER            NTIMER
      parameter( ntimer = 64 )
*     ..
*     .. Common Blocks ..
      LOGICAL            DISABLED
      DOUBLE PRECISION   CPUSEC( NTIMER ), CPUSTART( NTIMER ),
     $                   WALLSEC( NTIMER ), WALLSTART( NTIMER )
      COMMON /sltimer00/ cpusec, wallsec, cpustart, wallstart, disabled
*     ..
*     .. Executable Statements ..
*
      disabled = .true.
*
      RETURN
*
*     End of SLDISABLE
*
      END
*
      DOUBLE PRECISION FUNCTION slinquire( TIMETYPE, I )
*
*  -- ScaLAPACK tools routine (version 1.7) --
*     University of Tennessee, Knoxville, Oak Ridge National Laboratory,
*     and University of California, Berkeley.
*     May 1, 1997
*
*     .. Scalar Arguments ..
      CHARACTER*1        timetype
      INTEGER            i
*     ..
*
*  Purpose
*  =======
*
*  SLINQUIRE returns wall or cpu time that has accumulated in timer I.
*
*  Arguments
*  =========
*
*  TIMETYPE (global input) CHARACTER
*           Controls what time will be returned:
*           = 'W': wall clock time is returned,
*           = 'C': CPU time is returned (default).
*
*  I        (global input) INTEGER
*           The timer to return.
*
*  =====================================================================
*
*     .. Parameters ..
      INTEGER            ntimer
      parameter( ntimer = 64 )
      DOUBLE PRECISION   errflag
      parameter( errflag = -1.0d+0 )
*     ..
*     .. Local Scalars ..
      DOUBLE PRECISION   time
*     ..
*     .. External Functions ..
      LOGICAL            lsame
      DOUBLE PRECISION   dcputime00, dwalltime00
      EXTERNAL           dcputime00, dwalltime00, lsame
*     ..
*     .. Common Blocks ..
      LOGICAL            disabled
      DOUBLE PRECISION   cpusec( ntimer ), cpustart( ntimer ),
     $                   wallsec( ntimer ), wallstart( ntimer )
      COMMON /sltimer00/ cpusec, wallsec, cpustart, wallstart, disabled
*     ..
*     .. Executable Statements ..
*
      IF( lsame( timetype, 'W' ) ) THEN
*
*        If walltime not available on this machine, return -1 flag
*
         IF( dwalltime00().EQ.errflag ) THEN
            time = errflag
         ELSE
            time = wallsec( i )
         END IF
      ELSE
         IF( dcputime00().EQ.errflag ) THEN
            time = errflag
         ELSE
            time = cpusec( i )
         END IF
      END IF
*
      slinquire = time
*
      RETURN
*
*     End of SLINQUIRE
*
      END
*
      SUBROUTINE slcombine( ICTXT, SCOPE, OP, TIMETYPE, N, IBEG,
     $                      TIMES )
*
*  -- ScaLAPACK tools routine (version 1.7) --
*     University of Tennessee, Knoxville, Oak Ridge National Laboratory,
*     and University of California, Berkeley.
*     May 1, 1997
*
*     .. Scalar Arguments ..
      CHARACTER          OP, SCOPE, TIMETYPE
      INTEGER            IBEG, ICTXT, N
*     ..
*     .. Array Arguments ..
      DOUBLE PRECISION   TIMES( N )
*     ..
*
*  Purpose
*  =======
*
*  SLCOMBINE takes the timing information stored on a scope of processes
*  and combines them into the user's TIMES array.
*
*  Arguments
*  =========
*
*  ICTXT    (local input) INTEGER
*           The BLACS context handle.
*
*  SCOPE    (global input) CHARACTER
*           Controls what processes in grid participate in combine.
*           Options are 'Rowwise', 'Columnwise', or 'All'.
*
*  OP       (global input) CHARACTER
*           Controls what combine should be done:
*           = '>': get maximal time on any process (default),
*           = '<': get minimal time on any process,
*           = '+': get sum of times across processes.
*
*  TIMETYPE (global input) CHARACTER
*           Controls what time will be returned in TIMES:
*           = 'W': wall clock time,
*           = 'C': CPU time (default).
*
*  N        (global input) INTEGER
*           The number of timers to combine.
*
*  IBEG     (global input) INTEGER
*           The first timer to be combined.
*
*  TIMES    (global output) DOUBLE PRECISION array, dimension (N)
*           The requested timing information is returned in this array.
*
*  =====================================================================
*
*     .. Parameters ..
      INTEGER            NTIMER
      parameter( ntimer = 64 )
      DOUBLE PRECISION   ERRFLAG
      parameter( errflag = -1.0d+0 )
*     ..
*     .. Local Scalars ..
      LOGICAL            TMPDIS
      INTEGER            I
*     ..
*     .. External Subroutines ..
      EXTERNAL           dgamx2d, dgamn2d, dgsum2d
*     ..
*     .. External Functions ..
      LOGICAL            LSAME
      DOUBLE PRECISION   DCPUTIME00, DWALLTIME00
      EXTERNAL           dcputime00, dwalltime00, lsame
*     ..
*     .. Common Blocks ..
      LOGICAL            DISABLED
      DOUBLE PRECISION   CPUSEC( NTIMER ), CPUSTART( NTIMER ),
     $                   wallsec( ntimer ), wallstart( ntimer )
      COMMON /sltimer00/ cpusec, wallsec, cpustart, wallstart, disabled
*     ..
*     .. Executable Statements ..
*
*     Disable timer for combine operation
*
      tmpdis = disabled
      disabled = .true.
*
*     Copy timer information into user's times array
*
      IF( lsame( timetype, 'W' ) ) THEN
*
*        If walltime not available on this machine, fill in times
*        with -1 flag, and return
*
         IF( dwalltime00().EQ.errflag ) THEN
            DO 10 i = 1, n
               times( i ) = errflag
   10       CONTINUE
            RETURN
         ELSE
            DO 20 i = 1, n
               times( i ) = wallsec( ibeg + i - 1 )
   20       CONTINUE
         END IF
      ELSE
         IF( dcputime00().EQ.errflag ) THEN
            DO 30 i = 1, n
               times( i ) = errflag
   30       CONTINUE
            RETURN
         ELSE
            DO 40 i = 1, n
               times( i ) = cpusec( ibeg + i - 1 )
   40       CONTINUE
         END IF
      ENDIF
*
*     Combine all nodes' information, restore disabled, and return
*
      IF( op.EQ.'>' ) THEN
         CALL dgamx2d( ictxt, scope, ' ', n, 1, times, n, -1, -1,
     $                 -1, -1, 0 )
      ELSE IF( op.EQ.'<' ) THEN
         CALL dgamn2d( ictxt, scope, ' ', n, 1, times, n, -1, -1,
     $                 -1, -1, 0 )
      ELSE IF( op.EQ.'+' ) THEN
         CALL dgsum2d( ictxt, scope, ' ', n, 1, times, n, -1, 0 )
      ELSE
         CALL dgamx2d( ictxt, scope, ' ', n, 1, times, n, -1, -1,
     $                 -1, -1, 0 )
      END IF
*
      disabled = tmpdis
*
      RETURN
*
*     End of SLCOMBINE
*
      END