| NFFT 3.5.0

- Main Page
- Modules
- Data Structures
- Files

- File List
- Globals

cycle.h
1 /*
2 * Copyright (c) 2003, 2007-14 Matteo Frigo
3 * Copyright (c) 2003, 2007-14 Massachusetts Institute of Technology
4 *
5 * Permission is hereby granted, free of charge, to any person obtaining
6 * a copy of this software and associated documentation files (the
7 * "Software"), to deal in the Software without restriction, including
8 * without limitation the rights to use, copy, modify, merge, publish,
9 * distribute, sublicense, and/or sell copies of the Software, and to
10 * permit persons to whom the Software is furnished to do so, subject to
11 * the following conditions:
12 *
13 * The above copyright notice and this permission notice shall be
14 * included in all copies or substantial portions of the Software.
15 *
16 * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
17 * EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
18 * MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
19 * NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE
20 * LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
21 * OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION
22 * WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
23 *
24 */
25
26
27 /* machine-dependent cycle counters code. Needs to be inlined. */
28
29 /***************************************************************************/
30 /* To use the cycle counters in your code, simply #include "cycle.h" (this
31 file), and then use the functions/macros:
32
33 ticks getticks(void);
34
35 ticks is an opaque typedef defined below, representing the current time.
36 You extract the elapsed time between two calls to gettick() via:
37
38 double elapsed(ticks t1, ticks t0);
39
40 which returns a double-precision variable in arbitrary units. You
41 are not expected to convert this into human units like seconds; it
42 is intended only for *comparisons* of time intervals.
43
44 (In order to use some of the OS-dependent timer routines like
45 Solaris' gethrtime, you need to paste the autoconf snippet below
46 into your configure.ac file and #include "config.h" before cycle.h,
47 or define the relevant macros manually if you are not using autoconf.)
48 */
49
50 /***************************************************************************/
51 /* This file uses macros like HAVE_GETHRTIME that are assumed to be
52 defined according to whether the corresponding function/type/header
53 is available on your system. The necessary macros are most
54 conveniently defined if you are using GNU autoconf, via the tests:
55
56 dnl ---------------------------------------------------------------------
57
58 AC_C_INLINE
59 AC_HEADER_TIME
60 AC_CHECK_HEADERS([sys/time.h c_asm.h intrinsics.h mach/mach_time.h])
61
62 AC_CHECK_TYPE([hrtime_t],[AC_DEFINE(HAVE_HRTIME_T, 1, [Define to 1 if hrtime_t is defined in <sys/time.h>])],,[#if HAVE_SYS_TIME_H
63 #include <sys/time.h>
64 #endif])
65
66 AC_CHECK_FUNCS([gethrtime read_real_time time_base_to_time clock_gettime mach_absolute_time])
67
68 dnl Cray UNICOS _rtc() (real-time clock) intrinsic
69 AC_MSG_CHECKING([for _rtc intrinsic])
70 rtc_ok=yes
71 AC_TRY_LINK([#ifdef HAVE_INTRINSICS_H
72 #include <intrinsics.h>
73 #endif], [_rtc()], [AC_DEFINE(HAVE__RTC,1,[Define if you have the UNICOS _rtc() intrinsic.])], [rtc_ok=no])
74 AC_MSG_RESULT($rtc_ok)
75
76 dnl ---------------------------------------------------------------------
77 */
78
79 /***************************************************************************/
80
81 #if TIME_WITH_SYS_TIME
82 # include <sys/time.h>
83 # include <time.h>
84 #else
85 # if HAVE_SYS_TIME_H
86 # include <sys/time.h>
87 # else
88 # include <time.h>
89 # endif
90 #endif
91
92 #define INLINE_ELAPSED(INL) static INL double elapsed(ticks t1, ticks t0) \
93 { \
94 return (double)t1 - (double)t0; \
95 }
96
97 /*----------------------------------------------------------------*/
98 /* Solaris */
99 #if defined(HAVE_GETHRTIME) && defined(HAVE_HRTIME_T) && !defined(HAVE_TICK_COUNTER)
100 typedef hrtime_t ticks;
101
102 #define getticks gethrtime
103
104 INLINE_ELAPSED(inline)
105
106 #define HAVE_TICK_COUNTER
107 #endif
108
109 /*----------------------------------------------------------------*/
110 /* AIX v. 4+ routines to read the real-time clock or time-base register */
111 #if defined(HAVE_READ_REAL_TIME) && defined(HAVE_TIME_BASE_TO_TIME) && !defined(HAVE_TICK_COUNTER)
112 typedef timebasestruct_t ticks;
113
114 static __inline ticks getticks(void)
115 {
116 ticks t;
117 read_real_time(&t, TIMEBASE_SZ);
118 return t;
119 }
120
121 static __inline double elapsed(ticks t1, ticks t0) /* time in nanoseconds */
122 {
123 time_base_to_time(&t1, TIMEBASE_SZ);
124 time_base_to_time(&t0, TIMEBASE_SZ);
125 return (((double)t1.tb_high - (double)t0.tb_high) * 1.0e9 +
126 ((double)t1.tb_low - (double)t0.tb_low));
127 }
128
129 #define HAVE_TICK_COUNTER
130 #endif
131
132 /*----------------------------------------------------------------*/
133 /*
134 * PowerPC ``cycle'' counter using the time base register.
135 */
136 #if ((((defined(__GNUC__) && (defined(__powerpc__) || defined(__ppc__))) || (defined(__MWERKS__) && defined(macintosh)))) || (defined(__IBM_GCC_ASM) && (defined(__powerpc__) || defined(__ppc__)))) && !defined(HAVE_TICK_COUNTER)
137 typedef unsigned long long ticks;
138
139 static __inline__ ticks getticks(void)
140 {
141 unsigned int tbl, tbu0, tbu1;
142
143 do {
144 __asm__ __volatile__ ("mftbu %0" : "=r"(tbu0));
145 __asm__ __volatile__ ("mftb %0" : "=r"(tbl));
146 __asm__ __volatile__ ("mftbu %0" : "=r"(tbu1));
147 } while (tbu0 != tbu1);
148
149 return (((unsigned long long)tbu0) << 32) | tbl;
150 }
151
152 INLINE_ELAPSED(__inline__)
153
154 #define HAVE_TICK_COUNTER
155 #endif
156
157 /* MacOS/Mach (Darwin) time-base register interface (unlike UpTime,
158 from Carbon, requires no additional libraries to be linked). */
159 #if defined(HAVE_MACH_ABSOLUTE_TIME) && defined(HAVE_MACH_MACH_TIME_H) && !defined(HAVE_TICK_COUNTER)
160 #include <mach/mach_time.h>
161 typedef uint64_t ticks;
162 #define getticks mach_absolute_time
163 INLINE_ELAPSED(__inline__)
164 #define HAVE_TICK_COUNTER
165 #endif
166
167 /*----------------------------------------------------------------*/
168 /*
169 * Pentium cycle counter
170 */
171 #if (defined(__GNUC__) || defined(__ICC)) && defined(__i386__) && !defined(HAVE_TICK_COUNTER)
172 typedef unsigned long long ticks;
173
174 static __inline__ ticks getticks(void)
175 {
176 ticks ret;
177
178 __asm__ __volatile__("rdtsc": "=A" (ret));
179 /* no input, nothing else clobbered */
180 return ret;
181 }
182
183 INLINE_ELAPSED(__inline__)
184
185 #define HAVE_TICK_COUNTER
186 #define TIME_MIN 5000.0 /* unreliable pentium IV cycle counter */
187 #endif
188
189 /* Visual C++ -- thanks to Morten Nissov for his help with this */
190 #if defined(_MSC_VER) && _MSC_VER >= 1200 && _M_IX86 >= 500 && !defined(HAVE_TICK_COUNTER)
191 #include <windows.h>
192 typedef LARGE_INTEGER ticks;
193 #define RDTSC __asm __emit 0fh __asm __emit 031h /* hack for VC++ 5.0 */
194
195 static __inline ticks getticks(void)
196 {
197 ticks retval;
198
199 __asm {
200 RDTSC
201 mov retval.HighPart, edx
202 mov retval.LowPart, eax
203 }
204 return retval;
205 }
206
207 static __inline double elapsed(ticks t1, ticks t0)
208 {
209 return (double)t1.QuadPart - (double)t0.QuadPart;
210 }
211
212 #define HAVE_TICK_COUNTER
213 #define TIME_MIN 5000.0 /* unreliable pentium IV cycle counter */
214 #endif
215
216 /*----------------------------------------------------------------*/
217 /*
218 * X86-64 cycle counter
219 */
220 #if (defined(__GNUC__) || defined(__ICC) || defined(__SUNPRO_C)) && defined(__x86_64__) && !defined(HAVE_TICK_COUNTER)
221 typedef unsigned long long ticks;
222
223 static __inline__ ticks getticks(void)
224 {
225 unsigned a, d;
226 asm volatile("rdtsc" : "=a" (a), "=d" (d));
227 return ((ticks)a) | (((ticks)d) << 32);
228 }
229
230 INLINE_ELAPSED(__inline__)
231
232 #define HAVE_TICK_COUNTER
233 #define TIME_MIN 5000.0
234 #endif
235
236 /* PGI compiler, courtesy Cristiano Calonaci, Andrea Tarsi, & Roberto Gori.
237 NOTE: this code will fail to link unless you use the -Masmkeyword compiler
238 option (grrr). */
239 #if defined(__PGI) && defined(__x86_64__) && !defined(HAVE_TICK_COUNTER)
240 typedef unsigned long long ticks;
241 static ticks getticks(void)
242 {
243 asm(" rdtsc; shl $0x20,%rdx; mov %eax,%eax; or %rdx,%rax; ");
244 }
245 INLINE_ELAPSED(__inline__)
246 #define HAVE_TICK_COUNTER
247 #define TIME_MIN 5000.0
248 #endif
249
250 /* Visual C++, courtesy of Dirk Michaelis */
251 #if defined(_MSC_VER) && _MSC_VER >= 1400 && (defined(_M_AMD64) || defined(_M_X64)) && !defined(HAVE_TICK_COUNTER)
252
253 #include <intrin.h>
254 #pragma intrinsic(__rdtsc)
255 typedef unsigned __int64 ticks;
256 #define getticks __rdtsc
257 INLINE_ELAPSED(__inline)
258
259 #define HAVE_TICK_COUNTER
260 #define TIME_MIN 5000.0
261 #endif
262
263 /*----------------------------------------------------------------*/
264 /*
265 * IA64 cycle counter
266 */
267
268 /* intel's icc/ecc compiler */
269 #if (defined(__EDG_VERSION) || defined(__ECC)) && defined(__ia64__) && !defined(HAVE_TICK_COUNTER)
270 typedef unsigned long ticks;
271 #include <ia64intrin.h>
272
273 static __inline__ ticks getticks(void)
274 {
275 return __getReg(_IA64_REG_AR_ITC);
276 }
277
278 INLINE_ELAPSED(__inline__)
279
280 #define HAVE_TICK_COUNTER
281 #endif
282
283 /* gcc */
284 #if defined(__GNUC__) && defined(__ia64__) && !defined(HAVE_TICK_COUNTER)
285 typedef unsigned long ticks;
286
287 static __inline__ ticks getticks(void)
288 {
289 ticks ret;
290
291 __asm__ __volatile__ ("mov %0=ar.itc" : "=r"(ret));
292 return ret;
293 }
294
295 INLINE_ELAPSED(__inline__)
296
297 #define HAVE_TICK_COUNTER
298 #endif
299
300 /* HP/UX IA64 compiler, courtesy Teresa L. Johnson: */
301 #if defined(__hpux) && defined(__ia64) && !defined(HAVE_TICK_COUNTER)
302 #include <machine/sys/inline.h>
303 typedef unsigned long ticks;
304
305 static inline ticks getticks(void)
306 {
307 ticks ret;
308
309 ret = _Asm_mov_from_ar (_AREG_ITC);
310 return ret;
311 }
312
313 INLINE_ELAPSED(inline)
314
315 #define HAVE_TICK_COUNTER
316 #endif
317
318 /* Microsoft Visual C++ */
319 #if defined(_MSC_VER) && defined(_M_IA64) && !defined(HAVE_TICK_COUNTER)
320 typedef unsigned __int64 ticks;
321
322 # ifdef __cplusplus
323 extern "C"
324 # endif
325 ticks __getReg(int whichReg);
326 #pragma intrinsic(__getReg)
327
328 static __inline ticks getticks(void)
329 {
330 volatile ticks temp;
331 temp = __getReg(3116);
332 return temp;
333 }
334
335 INLINE_ELAPSED(inline)
336
337 #define HAVE_TICK_COUNTER
338 #endif
339
340 /*----------------------------------------------------------------*/
341 /*
342 * PA-RISC cycle counter
343 */
344 #if (defined(__hppa__) || defined(__hppa)) && !defined(HAVE_TICK_COUNTER)
345 typedef unsigned long ticks;
346
347 # ifdef __GNUC__
348 static __inline__ ticks getticks(void)
349 {
350 ticks ret;
351
352 __asm__ __volatile__("mfctl 16, %0": "=r" (ret));
353 /* no input, nothing else clobbered */
354 return ret;
355 }
356 # else
357 # include <machine/inline.h>
358 static inline unsigned long getticks(void)
359 {
360 register ticks ret;
361 _MFCTL(16, ret);
362 return ret;
363 }
364 # endif
365
366 INLINE_ELAPSED(inline)
367
368 #define HAVE_TICK_COUNTER
369 #endif
370
371 /*----------------------------------------------------------------*/
372 /* S390, courtesy of James Treacy */
373 #if defined(__GNUC__) && defined(__s390__) && !defined(HAVE_TICK_COUNTER)
374 typedef unsigned long long ticks;
375
376 static __inline__ ticks getticks(void)
377 {
378 ticks cycles;
379 __asm__("stck 0(%0)" : : "a" (&(cycles)) : "memory", "cc");
380 return cycles;
381 }
382
383 INLINE_ELAPSED(__inline__)
384
385 #define HAVE_TICK_COUNTER
386 #endif
387 /*----------------------------------------------------------------*/
388 #if defined(__GNUC__) && defined(__alpha__) && !defined(HAVE_TICK_COUNTER)
389 /*
390 * The 32-bit cycle counter on alpha overflows pretty quickly,
391 * unfortunately. A 1GHz machine overflows in 4 seconds.
392 */
393 typedef unsigned int ticks;
394
395 static __inline__ ticks getticks(void)
396 {
397 unsigned long cc;
398 __asm__ __volatile__ ("rpcc %0" : "=r"(cc));
399 return (cc & 0xFFFFFFFF);
400 }
401
402 INLINE_ELAPSED(__inline__)
403
404 #define HAVE_TICK_COUNTER
405 #endif
406
407 /*----------------------------------------------------------------*/
408 #if defined(__GNUC__) && defined(__sparc_v9__) && !defined(HAVE_TICK_COUNTER)
409 typedef unsigned long ticks;
410
411 static __inline__ ticks getticks(void)
412 {
413 ticks ret;
414 __asm__ __volatile__("rd %%tick, %0" : "=r" (ret));
415 return ret;
416 }
417
418 INLINE_ELAPSED(__inline__)
419
420 #define HAVE_TICK_COUNTER
421 #endif
422
423 /*----------------------------------------------------------------*/
424 #if (defined(__DECC) || defined(__DECCXX)) && defined(__alpha) && defined(HAVE_C_ASM_H) && !defined(HAVE_TICK_COUNTER)
425 # include <c_asm.h>
426 typedef unsigned int ticks;
427
428 static __inline ticks getticks(void)
429 {
430 unsigned long cc;
431 cc = asm("rpcc %v0");
432 return (cc & 0xFFFFFFFF);
433 }
434
435 INLINE_ELAPSED(__inline)
436
437 #define HAVE_TICK_COUNTER
438 #endif
439 /*----------------------------------------------------------------*/
440 /* SGI/Irix */
441 #if defined(HAVE_CLOCK_GETTIME) && defined(CLOCK_SGI_CYCLE) && !defined(HAVE_TICK_COUNTER)
442 typedef struct timespec ticks;
443
444 static inline ticks getticks(void)
445 {
446 struct timespec t;
447 clock_gettime(CLOCK_SGI_CYCLE, &t);
448 return t;
449 }
450
451 static inline double elapsed(ticks t1, ticks t0)
452 {
453 return ((double)t1.tv_sec - (double)t0.tv_sec) * 1.0E9 +
454 ((double)t1.tv_nsec - (double)t0.tv_nsec);
455 }
456 #define HAVE_TICK_COUNTER
457 #endif
458
459 /*----------------------------------------------------------------*/
460 /* Cray UNICOS _rtc() intrinsic function */
461 #if defined(HAVE__RTC) && !defined(HAVE_TICK_COUNTER)
462 #ifdef HAVE_INTRINSICS_H
463 # include <intrinsics.h>
464 #endif
465
466 typedef long long ticks;
467
468 #define getticks _rtc
469
470 INLINE_ELAPSED(inline)
471
472 #define HAVE_TICK_COUNTER
473 #endif
474
475 /*----------------------------------------------------------------*/
476 /* MIPS ZBus */
477 #if defined(HAVE_MIPS_ZBUS_TIMER) && HAVE_MIPS_ZBUS_TIMER
478 #if defined(__mips__) && !defined(HAVE_TICK_COUNTER)
479 #include <sys/mman.h>
480 #include <unistd.h>
481 #include <fcntl.h>
482
483 typedef uint64_t ticks;
484
485 static inline ticks getticks(void)
486 {
487 static uint64_t* addr = 0;
488
489 if (addr == 0)
490 {
491 uint32_t rq_addr = 0x10030000;
492 int fd;
493 int pgsize;
494
495 pgsize = getpagesize();
496 fd = open ("/dev/mem", O_RDONLY | O_SYNC, 0);
497 if (fd < 0) {
498 perror("open");
499 return NULL;
500 }
501 addr = mmap(0, pgsize, PROT_READ, MAP_SHARED, fd, rq_addr);
502 close(fd);
503 if (addr == (uint64_t *)-1) {
504 perror("mmap");
505 return NULL;
506 }
507 }
508
509 return *addr;
510 }
511
512 INLINE_ELAPSED(inline)
513
514 #define HAVE_TICK_COUNTER
515 #endif
516 #endif /* HAVE_MIPS_ZBUS_TIMER */
517