MODER/timing.hpp at master · jttoivon/MODER · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
/*

    MODER is a program to learn DNA binding motifs from SELEX datasets.
    Copyright (C) 2016, 2017  Jarkko Toivonen,
    Department of Computer Science, University of Helsinki

    MODER is free software; you can redistribute it and/or modify
    it under the terms of the GNU General Public License as published by
    the Free Software Foundation; either version 2 of the License, or
    (at your option) any later version.

    MODER is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU General Public License for more details.

    You should have received a copy of the GNU General Public License along
    with this program; if not, write to the Free Software Foundation, Inc.,
    51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.

*/
//
// 2003 Juha K"arkk"ainen
//
// Timer
//

#ifndef _TIMING_HPP
#define _TIMING_HPP

#include <time.h>
#include <sys/time.h>
#include <iostream>
#include <iomanip>

inline double get_cpu_time()
{
  static const double in_seconds = 1.0/static_cast<double>(CLOCKS_PER_SEC);
  return clock() * in_seconds;
}


#ifdef TIMING

#define TIME_START(t) \
  double timer_##t##_current_, \
    timer_##t##_last_ = get_cpu_time()

#define TIME_CHECK(t) \
  timer_##t##_current_ = get_cpu_time(); \
  std::cout << std::setw(10) << std::fixed << std::setprecision(2) \
            << timer_##t##_current_ - timer_##t##_last_ << std::endl; \
  timer_##t##_last_ = timer_##t##_current_

#define TIME_PRINT(format,t)		 \
  timer_##t##_current_ = get_cpu_time(); \
  printf(format, (timer_##t##_current_ - timer_##t##_last_));	\
  timer_##t##_last_ = timer_##t##_current_

#define TIME_GET(t)		 \
  ((timer_##t##_current_ = get_cpu_time()), (timer_##t##_current_ - timer_##t##_last_))


#ifdef __APPLE__ // apple machines don't define clock_gettime. Not at least the older ones.

inline
timeval
get_wall_time()
{
  struct timeval tv;
  gettimeofday(&tv, NULL);
  return tv;
}

inline
double
subtract_time(timeval& tp_new, timeval& tp_old)
{
  timeval temp = tp_new;
  temp.tv_sec -= tp_old.tv_sec;
  temp.tv_usec -= tp_old.tv_usec;
  if (temp.tv_usec < 0) {
    --temp.tv_sec;
    temp.tv_usec += 1000000;
  }
  double result = temp.tv_sec;
  result += temp.tv_usec / 1000000.0;

  return result;
}

#define WALL_TIME_START(t) \
  timeval timer_##t##_current_, \
  timer_##t##_last_ = get_wall_time()

#define WALL_TIME_PRINT(format,t)		 \
  timer_##t##_current_ = get_wall_time(); \
  printf(format, subtract_time(timer_##t##_current_, timer_##t##_last_)); \
  timer_##t##_last_ = timer_##t##_current_

#else

inline
timespec
get_wall_time()
{
  struct timespec tp;
  clock_gettime(CLOCK_REALTIME, &tp);
  return tp;
}

inline
double
subtract_time(timespec& tp_new, timespec& tp_old)
{
  timespec temp = tp_new;
  temp.tv_sec -= tp_old.tv_sec;
  temp.tv_nsec -= tp_old.tv_nsec;
  if (temp.tv_nsec < 0) {
    --temp.tv_sec;
    temp.tv_nsec += 1000000000;
  }
  double result = temp.tv_sec;
  result += temp.tv_nsec / 1000000000.0;

  return result;
}

#define WALL_TIME_START(t) \
  timespec timer_##t##_current_, \
  timer_##t##_last_ = get_wall_time()

#define WALL_TIME_PRINT(format,t)		 \
  timer_##t##_current_ = get_wall_time(); \
  printf(format, subtract_time(timer_##t##_current_, timer_##t##_last_)); \
  timer_##t##_last_ = timer_##t##_current_
#endif

#else

#define TIME_START(t)
#define TIME_CHECK(t)
#define TIME_PRINT(format,t)
#define TIME_GET(t) 0.0
#define WALL_TIME_START(t)
#define WALL_TIME_PRINT(format,t)

#endif

#endif // _TIMING_HPP