BGI/BASE/src/2bwt/Timing.c

   1 /*
   2
   3    Timing.c             Measuring Program running time
   4
   5    This module contains functions for measuring program running time.
   6
   7    Copyright (C) 2004, Wong Chi Kwong.
   8
   9    This program is free software; you can redistribute it and/or
  10    modify it under the terms of the GNU General Public License
  11    as published by the Free Software Foundation; either version 2
  12    of the License, or (at your option) any later version.
  13
  14    This program is distributed in the hope that it will be useful,
  15    but WITHOUT ANY WARRANTY; without even the implied warranty of
  16    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  17    GNU General Public License for more details.
  18
  19    You should have received a copy of the GNU General Public License
  20    along with this program; if not, write to the Free Software
  21    Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301, USA.
  22
  23 */
  24
  25 #include <stdio.h>
  26 #include <stdlib.h>
  27 #include <string.h>
  28
  29 #ifdef RUSAGE
  30 #include <sys/resource.h>
  31 #else
  32 #ifdef TIME_BY_CLOCK
  33 #include <time.h>
  34 #else
  35 #include <sys/time.h>
  36 #endif
  37 #endif
  38
  39 #include "Timing.h"
  40
  41
  42 double setStartTime() {
  43
  44 #ifdef RUSAGE
  45
  46         double usertime, systime;
  47         struct rusage usage;
  48
  49         getrusage(RUSAGE_SELF, &usage);
  50
  51         usertime = (double)usage.ru_utime.tv_sec + (double)usage.ru_utime.tv_usec / 1000000.0;
  52         systime = (double)usage.ru_stime.tv_sec + (double)usage.ru_stime.tv_usec / 1000000.0;
  53
  54         return(usertime + systime);
  55
  56 #else
  57 #ifdef TIME_BY_CLOCK
  58
  59         return (double)clock() / (double)CLOCKS_PER_SEC;
  60
  61 #else
  62
  63         struct timeval tp;
  64         gettimeofday(&tp, NULL);
  65         return (double)tp.tv_sec + (double)tp.tv_usec / (double)1000000;
  66
  67 #endif
  68 #endif
  69
  70 }
  71
  72 double getElapsedTime(double startTime) {
  73
  74 #ifdef RUSAGE
  75
  76         double usertime, systime;
  77         struct rusage usage;
  78
  79         getrusage(RUSAGE_SELF, &usage);
  80
  81         usertime = (double)usage.ru_utime.tv_sec + (double)usage.ru_utime.tv_usec / 1000000.0;
  82         systime = (double)usage.ru_stime.tv_sec + (double)usage.ru_stime.tv_usec / 1000000.0;
  83
  84         return (usertime + systime) - startTime;
  85
  86 #else
  87 #ifdef TIME_BY_CLOCK
  88
  89         return (double)clock() / (double)CLOCKS_PER_SEC - startTime;
  90
  91 #else
  92
  93         struct timeval tp;
  94         gettimeofday(&tp, NULL);
  95         return (double)tp.tv_sec + (double)tp.tv_usec / (double)1000000 - startTime;
  96
  97 #endif
  98 #endif
  99
 100 }
 101
 102 void printElapsedTime(FILE *file, const int printHour, const int printMin, const int printSec,
 103                                           const int secNumberOfDecimal, const double seconds) {
 104
 105         printElapsedTimeNoNewLine(file, printHour, printMin, printSec, 0, secNumberOfDecimal, seconds);
 106         fprintf(file, "\n");
 107
 108 }
 109
 110 void printElapsedTimeNoNewLine(FILE *file, const int printHour, const int printMin, const int printSec,
 111                                           const int secMinPrintLength, const int secNumberOfDecimal, const double seconds) {
 112
 113         int hour, min;
 114         double sec;
 115         char secondDisplay[8] = "%0.0f s";
 116
 117         #ifdef DEBUG
 118         if (printHour && !printMin && printSec) {
 119                 fprintf(stderr, "printElapsedTime(): Cannot skip minute only!\n");
 120                 exit(1);
 121         }
 122         if (secNumberOfDecimal > 9) {
 123                 fprintf(stderr, "printElapsedTime(): secNumberOfDecimal > 9!\n");
 124                 exit(1);
 125         }
 126         #endif
 127
 128         secondDisplay[1] = secondDisplay[1] + (char)secMinPrintLength;
 129         secondDisplay[3] = secondDisplay[3] + (char)secNumberOfDecimal;
 130
 131         sec = seconds;
 132         min = (int)(seconds / 60);
 133         if (!printSec && printMin) {
 134                 if (seconds - min * 60 >= 30) {
 135                         min++;
 136                 }
 137         }
 138         if (printMin) {
 139                 sec -= min * 60;
 140         }
 141
 142         hour = min / 60;
 143         if (!printMin) {
 144                 min = hour * 60;
 145                 if (min >= 30) {
 146                         hour++;
 147                 }
 148         }
 149         if (printHour) {
 150                 min -= hour * 60;
 151         }
 152
 153         if (printHour) {
 154         fprintf(file, "%d h ", hour);
 155         }
 156         if (printMin) {
 157                 fprintf(file, "%d m ", min);
 158         }
 159         if (printSec) {
 160                 fprintf(file, secondDisplay, sec);
 161         }
 162
 163 }
 164