cannam@167: /* cannam@167: * Copyright (c) 2003, 2007-14 Matteo Frigo cannam@167: * Copyright (c) 2003, 2007-14 Massachusetts Institute of Technology cannam@167: * cannam@167: * Permission is hereby granted, free of charge, to any person obtaining cannam@167: * a copy of this software and associated documentation files (the cannam@167: * "Software"), to deal in the Software without restriction, including cannam@167: * without limitation the rights to use, copy, modify, merge, publish, cannam@167: * distribute, sublicense, and/or sell copies of the Software, and to cannam@167: * permit persons to whom the Software is furnished to do so, subject to cannam@167: * the following conditions: cannam@167: * cannam@167: * The above copyright notice and this permission notice shall be cannam@167: * included in all copies or substantial portions of the Software. cannam@167: * cannam@167: * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, cannam@167: * EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF cannam@167: * MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND cannam@167: * NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE cannam@167: * LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION cannam@167: * OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION cannam@167: * WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE. cannam@167: * cannam@167: */ cannam@167: cannam@167: cannam@167: /* machine-dependent cycle counters code. Needs to be inlined. */ cannam@167: cannam@167: /***************************************************************************/ cannam@167: /* To use the cycle counters in your code, simply #include "cycle.h" (this cannam@167: file), and then use the functions/macros: cannam@167: cannam@167: ticks getticks(void); cannam@167: cannam@167: ticks is an opaque typedef defined below, representing the current time. cannam@167: You extract the elapsed time between two calls to gettick() via: cannam@167: cannam@167: double elapsed(ticks t1, ticks t0); cannam@167: cannam@167: which returns a double-precision variable in arbitrary units. You cannam@167: are not expected to convert this into human units like seconds; it cannam@167: is intended only for *comparisons* of time intervals. cannam@167: cannam@167: (In order to use some of the OS-dependent timer routines like cannam@167: Solaris' gethrtime, you need to paste the autoconf snippet below cannam@167: into your configure.ac file and #include "config.h" before cycle.h, cannam@167: or define the relevant macros manually if you are not using autoconf.) cannam@167: */ cannam@167: cannam@167: /***************************************************************************/ cannam@167: /* This file uses macros like HAVE_GETHRTIME that are assumed to be cannam@167: defined according to whether the corresponding function/type/header cannam@167: is available on your system. The necessary macros are most cannam@167: conveniently defined if you are using GNU autoconf, via the tests: cannam@167: cannam@167: dnl --------------------------------------------------------------------- cannam@167: cannam@167: AC_C_INLINE cannam@167: AC_HEADER_TIME cannam@167: AC_CHECK_HEADERS([sys/time.h c_asm.h intrinsics.h mach/mach_time.h]) cannam@167: cannam@167: AC_CHECK_TYPE([hrtime_t],[AC_DEFINE(HAVE_HRTIME_T, 1, [Define to 1 if hrtime_t is defined in ])],,[#if HAVE_SYS_TIME_H cannam@167: #include cannam@167: #endif]) cannam@167: cannam@167: AC_CHECK_FUNCS([gethrtime read_real_time time_base_to_time clock_gettime mach_absolute_time]) cannam@167: cannam@167: dnl Cray UNICOS _rtc() (real-time clock) intrinsic cannam@167: AC_MSG_CHECKING([for _rtc intrinsic]) cannam@167: rtc_ok=yes cannam@167: AC_TRY_LINK([#ifdef HAVE_INTRINSICS_H cannam@167: #include cannam@167: #endif], [_rtc()], [AC_DEFINE(HAVE__RTC,1,[Define if you have the UNICOS _rtc() intrinsic.])], [rtc_ok=no]) cannam@167: AC_MSG_RESULT($rtc_ok) cannam@167: cannam@167: dnl --------------------------------------------------------------------- cannam@167: */ cannam@167: cannam@167: /***************************************************************************/ cannam@167: cannam@167: #if TIME_WITH_SYS_TIME cannam@167: # include cannam@167: # include cannam@167: #else cannam@167: # if HAVE_SYS_TIME_H cannam@167: # include cannam@167: # else cannam@167: # include cannam@167: # endif cannam@167: #endif cannam@167: cannam@167: #define INLINE_ELAPSED(INL) static INL double elapsed(ticks t1, ticks t0) \ cannam@167: { \ cannam@167: return (double)t1 - (double)t0; \ cannam@167: } cannam@167: cannam@167: /*----------------------------------------------------------------*/ cannam@167: /* Solaris */ cannam@167: #if defined(HAVE_GETHRTIME) && defined(HAVE_HRTIME_T) && !defined(HAVE_TICK_COUNTER) cannam@167: typedef hrtime_t ticks; cannam@167: cannam@167: #define getticks gethrtime cannam@167: cannam@167: INLINE_ELAPSED(inline) cannam@167: cannam@167: #define HAVE_TICK_COUNTER cannam@167: #endif cannam@167: cannam@167: /*----------------------------------------------------------------*/ cannam@167: /* AIX v. 4+ routines to read the real-time clock or time-base register */ cannam@167: #if defined(HAVE_READ_REAL_TIME) && defined(HAVE_TIME_BASE_TO_TIME) && !defined(HAVE_TICK_COUNTER) cannam@167: typedef timebasestruct_t ticks; cannam@167: cannam@167: static __inline ticks getticks(void) cannam@167: { cannam@167: ticks t; cannam@167: read_real_time(&t, TIMEBASE_SZ); cannam@167: return t; cannam@167: } cannam@167: cannam@167: static __inline double elapsed(ticks t1, ticks t0) /* time in nanoseconds */ cannam@167: { cannam@167: time_base_to_time(&t1, TIMEBASE_SZ); cannam@167: time_base_to_time(&t0, TIMEBASE_SZ); cannam@167: return (((double)t1.tb_high - (double)t0.tb_high) * 1.0e9 + cannam@167: ((double)t1.tb_low - (double)t0.tb_low)); cannam@167: } cannam@167: cannam@167: #define HAVE_TICK_COUNTER cannam@167: #endif cannam@167: cannam@167: /*----------------------------------------------------------------*/ cannam@167: /* cannam@167: * PowerPC ``cycle'' counter using the time base register. cannam@167: */ cannam@167: #if ((((defined(__GNUC__) && (defined(__powerpc__) || defined(__ppc__))) || (defined(__MWERKS__) && defined(macintosh)))) || (defined(__IBM_GCC_ASM) && (defined(__powerpc__) || defined(__ppc__)))) && !defined(HAVE_TICK_COUNTER) cannam@167: typedef unsigned long long ticks; cannam@167: cannam@167: static __inline__ ticks getticks(void) cannam@167: { cannam@167: unsigned int tbl, tbu0, tbu1; cannam@167: cannam@167: do { cannam@167: __asm__ __volatile__ ("mftbu %0" : "=r"(tbu0)); cannam@167: __asm__ __volatile__ ("mftb %0" : "=r"(tbl)); cannam@167: __asm__ __volatile__ ("mftbu %0" : "=r"(tbu1)); cannam@167: } while (tbu0 != tbu1); cannam@167: cannam@167: return (((unsigned long long)tbu0) << 32) | tbl; cannam@167: } cannam@167: cannam@167: INLINE_ELAPSED(__inline__) cannam@167: cannam@167: #define HAVE_TICK_COUNTER cannam@167: #endif cannam@167: cannam@167: /* MacOS/Mach (Darwin) time-base register interface (unlike UpTime, cannam@167: from Carbon, requires no additional libraries to be linked). */ cannam@167: #if defined(HAVE_MACH_ABSOLUTE_TIME) && defined(HAVE_MACH_MACH_TIME_H) && !defined(HAVE_TICK_COUNTER) cannam@167: #include cannam@167: typedef uint64_t ticks; cannam@167: #define getticks mach_absolute_time cannam@167: INLINE_ELAPSED(__inline__) cannam@167: #define HAVE_TICK_COUNTER cannam@167: #endif cannam@167: cannam@167: /*----------------------------------------------------------------*/ cannam@167: /* cannam@167: * Pentium cycle counter cannam@167: */ cannam@167: #if (defined(__GNUC__) || defined(__ICC)) && defined(__i386__) && !defined(HAVE_TICK_COUNTER) cannam@167: typedef unsigned long long ticks; cannam@167: cannam@167: static __inline__ ticks getticks(void) cannam@167: { cannam@167: ticks ret; cannam@167: cannam@167: __asm__ __volatile__("rdtsc": "=A" (ret)); cannam@167: /* no input, nothing else clobbered */ cannam@167: return ret; cannam@167: } cannam@167: cannam@167: INLINE_ELAPSED(__inline__) cannam@167: cannam@167: #define HAVE_TICK_COUNTER cannam@167: #define TIME_MIN 5000.0 /* unreliable pentium IV cycle counter */ cannam@167: #endif cannam@167: cannam@167: /* Visual C++ -- thanks to Morten Nissov for his help with this */ cannam@167: #if _MSC_VER >= 1200 && _M_IX86 >= 500 && !defined(HAVE_TICK_COUNTER) cannam@167: #include cannam@167: typedef LARGE_INTEGER ticks; cannam@167: #define RDTSC __asm __emit 0fh __asm __emit 031h /* hack for VC++ 5.0 */ cannam@167: cannam@167: static __inline ticks getticks(void) cannam@167: { cannam@167: ticks retval; cannam@167: cannam@167: __asm { cannam@167: RDTSC cannam@167: mov retval.HighPart, edx cannam@167: mov retval.LowPart, eax cannam@167: } cannam@167: return retval; cannam@167: } cannam@167: cannam@167: static __inline double elapsed(ticks t1, ticks t0) cannam@167: { cannam@167: return (double)t1.QuadPart - (double)t0.QuadPart; cannam@167: } cannam@167: cannam@167: #define HAVE_TICK_COUNTER cannam@167: #define TIME_MIN 5000.0 /* unreliable pentium IV cycle counter */ cannam@167: #endif cannam@167: cannam@167: /*----------------------------------------------------------------*/ cannam@167: /* cannam@167: * X86-64 cycle counter cannam@167: */ cannam@167: #if (defined(__GNUC__) || defined(__ICC) || defined(__SUNPRO_C)) && defined(__x86_64__) && !defined(HAVE_TICK_COUNTER) cannam@167: typedef unsigned long long ticks; cannam@167: cannam@167: static __inline__ ticks getticks(void) cannam@167: { cannam@167: unsigned a, d; cannam@167: __asm__ __volatile__ ("rdtsc" : "=a" (a), "=d" (d)); cannam@167: return ((ticks)a) | (((ticks)d) << 32); cannam@167: } cannam@167: cannam@167: INLINE_ELAPSED(__inline__) cannam@167: cannam@167: #define HAVE_TICK_COUNTER cannam@167: #define TIME_MIN 5000.0 cannam@167: #endif cannam@167: cannam@167: /* PGI compiler, courtesy Cristiano Calonaci, Andrea Tarsi, & Roberto Gori. cannam@167: NOTE: this code will fail to link unless you use the -Masmkeyword compiler cannam@167: option (grrr). */ cannam@167: #if defined(__PGI) && defined(__x86_64__) && !defined(HAVE_TICK_COUNTER) cannam@167: typedef unsigned long long ticks; cannam@167: static ticks getticks(void) cannam@167: { cannam@167: asm(" rdtsc; shl $0x20,%rdx; mov %eax,%eax; or %rdx,%rax; "); cannam@167: } cannam@167: INLINE_ELAPSED(__inline__) cannam@167: #define HAVE_TICK_COUNTER cannam@167: #define TIME_MIN 5000.0 cannam@167: #endif cannam@167: cannam@167: /* Visual C++, courtesy of Dirk Michaelis */ cannam@167: #if _MSC_VER >= 1400 && (defined(_M_AMD64) || defined(_M_X64)) && !defined(HAVE_TICK_COUNTER) cannam@167: cannam@167: #include cannam@167: #pragma intrinsic(__rdtsc) cannam@167: typedef unsigned __int64 ticks; cannam@167: #define getticks __rdtsc cannam@167: INLINE_ELAPSED(__inline) cannam@167: cannam@167: #define HAVE_TICK_COUNTER cannam@167: #define TIME_MIN 5000.0 cannam@167: #endif cannam@167: cannam@167: /*----------------------------------------------------------------*/ cannam@167: /* cannam@167: * IA64 cycle counter cannam@167: */ cannam@167: cannam@167: /* intel's icc/ecc compiler */ cannam@167: #if (defined(__EDG_VERSION) || defined(__ECC)) && defined(__ia64__) && !defined(HAVE_TICK_COUNTER) cannam@167: typedef unsigned long ticks; cannam@167: #include cannam@167: cannam@167: static __inline__ ticks getticks(void) cannam@167: { cannam@167: return __getReg(_IA64_REG_AR_ITC); cannam@167: } cannam@167: cannam@167: INLINE_ELAPSED(__inline__) cannam@167: cannam@167: #define HAVE_TICK_COUNTER cannam@167: #endif cannam@167: cannam@167: /* gcc */ cannam@167: #if defined(__GNUC__) && defined(__ia64__) && !defined(HAVE_TICK_COUNTER) cannam@167: typedef unsigned long ticks; cannam@167: cannam@167: static __inline__ ticks getticks(void) cannam@167: { cannam@167: ticks ret; cannam@167: cannam@167: __asm__ __volatile__ ("mov %0=ar.itc" : "=r"(ret)); cannam@167: return ret; cannam@167: } cannam@167: cannam@167: INLINE_ELAPSED(__inline__) cannam@167: cannam@167: #define HAVE_TICK_COUNTER cannam@167: #endif cannam@167: cannam@167: /* HP/UX IA64 compiler, courtesy Teresa L. Johnson: */ cannam@167: #if defined(__hpux) && defined(__ia64) && !defined(HAVE_TICK_COUNTER) cannam@167: #include cannam@167: typedef unsigned long ticks; cannam@167: cannam@167: static inline ticks getticks(void) cannam@167: { cannam@167: ticks ret; cannam@167: cannam@167: ret = _Asm_mov_from_ar (_AREG_ITC); cannam@167: return ret; cannam@167: } cannam@167: cannam@167: INLINE_ELAPSED(inline) cannam@167: cannam@167: #define HAVE_TICK_COUNTER cannam@167: #endif cannam@167: cannam@167: /* Microsoft Visual C++ */ cannam@167: #if defined(_MSC_VER) && defined(_M_IA64) && !defined(HAVE_TICK_COUNTER) cannam@167: typedef unsigned __int64 ticks; cannam@167: cannam@167: # ifdef __cplusplus cannam@167: extern "C" cannam@167: # endif cannam@167: ticks __getReg(int whichReg); cannam@167: #pragma intrinsic(__getReg) cannam@167: cannam@167: static __inline ticks getticks(void) cannam@167: { cannam@167: volatile ticks temp; cannam@167: temp = __getReg(3116); cannam@167: return temp; cannam@167: } cannam@167: cannam@167: INLINE_ELAPSED(inline) cannam@167: cannam@167: #define HAVE_TICK_COUNTER cannam@167: #endif cannam@167: cannam@167: /*----------------------------------------------------------------*/ cannam@167: /* cannam@167: * PA-RISC cycle counter cannam@167: */ cannam@167: #if (defined(__hppa__) || defined(__hppa)) && !defined(HAVE_TICK_COUNTER) cannam@167: typedef unsigned long ticks; cannam@167: cannam@167: # ifdef __GNUC__ cannam@167: static __inline__ ticks getticks(void) cannam@167: { cannam@167: ticks ret; cannam@167: cannam@167: __asm__ __volatile__("mfctl 16, %0": "=r" (ret)); cannam@167: /* no input, nothing else clobbered */ cannam@167: return ret; cannam@167: } cannam@167: # else cannam@167: # include cannam@167: static inline unsigned long getticks(void) cannam@167: { cannam@167: register ticks ret; cannam@167: _MFCTL(16, ret); cannam@167: return ret; cannam@167: } cannam@167: # endif cannam@167: cannam@167: INLINE_ELAPSED(inline) cannam@167: cannam@167: #define HAVE_TICK_COUNTER cannam@167: #endif cannam@167: cannam@167: /*----------------------------------------------------------------*/ cannam@167: /* S390, courtesy of James Treacy */ cannam@167: #if defined(__GNUC__) && defined(__s390__) && !defined(HAVE_TICK_COUNTER) cannam@167: typedef unsigned long long ticks; cannam@167: cannam@167: static __inline__ ticks getticks(void) cannam@167: { cannam@167: ticks cycles; cannam@167: __asm__("stck 0(%0)" : : "a" (&(cycles)) : "memory", "cc"); cannam@167: return cycles; cannam@167: } cannam@167: cannam@167: INLINE_ELAPSED(__inline__) cannam@167: cannam@167: #define HAVE_TICK_COUNTER cannam@167: #endif cannam@167: /*----------------------------------------------------------------*/ cannam@167: #if defined(__GNUC__) && defined(__alpha__) && !defined(HAVE_TICK_COUNTER) cannam@167: /* cannam@167: * The 32-bit cycle counter on alpha overflows pretty quickly, cannam@167: * unfortunately. A 1GHz machine overflows in 4 seconds. cannam@167: */ cannam@167: typedef unsigned int ticks; cannam@167: cannam@167: static __inline__ ticks getticks(void) cannam@167: { cannam@167: unsigned long cc; cannam@167: __asm__ __volatile__ ("rpcc %0" : "=r"(cc)); cannam@167: return (cc & 0xFFFFFFFF); cannam@167: } cannam@167: cannam@167: INLINE_ELAPSED(__inline__) cannam@167: cannam@167: #define HAVE_TICK_COUNTER cannam@167: #endif cannam@167: cannam@167: /*----------------------------------------------------------------*/ cannam@167: #if defined(__GNUC__) && defined(__sparc_v9__) && !defined(HAVE_TICK_COUNTER) cannam@167: typedef unsigned long ticks; cannam@167: cannam@167: static __inline__ ticks getticks(void) cannam@167: { cannam@167: ticks ret; cannam@167: __asm__ __volatile__("rd %%tick, %0" : "=r" (ret)); cannam@167: return ret; cannam@167: } cannam@167: cannam@167: INLINE_ELAPSED(__inline__) cannam@167: cannam@167: #define HAVE_TICK_COUNTER cannam@167: #endif cannam@167: cannam@167: /*----------------------------------------------------------------*/ cannam@167: #if (defined(__DECC) || defined(__DECCXX)) && defined(__alpha) && defined(HAVE_C_ASM_H) && !defined(HAVE_TICK_COUNTER) cannam@167: # include cannam@167: typedef unsigned int ticks; cannam@167: cannam@167: static __inline ticks getticks(void) cannam@167: { cannam@167: unsigned long cc; cannam@167: cc = asm("rpcc %v0"); cannam@167: return (cc & 0xFFFFFFFF); cannam@167: } cannam@167: cannam@167: INLINE_ELAPSED(__inline) cannam@167: cannam@167: #define HAVE_TICK_COUNTER cannam@167: #endif cannam@167: /*----------------------------------------------------------------*/ cannam@167: /* SGI/Irix */ cannam@167: #if defined(HAVE_CLOCK_GETTIME) && defined(CLOCK_SGI_CYCLE) && !defined(HAVE_TICK_COUNTER) && !defined(__ANDROID__) cannam@167: typedef struct timespec ticks; cannam@167: cannam@167: static inline ticks getticks(void) cannam@167: { cannam@167: struct timespec t; cannam@167: clock_gettime(CLOCK_SGI_CYCLE, &t); cannam@167: return t; cannam@167: } cannam@167: cannam@167: static inline double elapsed(ticks t1, ticks t0) cannam@167: { cannam@167: return ((double)t1.tv_sec - (double)t0.tv_sec) * 1.0E9 + cannam@167: ((double)t1.tv_nsec - (double)t0.tv_nsec); cannam@167: } cannam@167: #define HAVE_TICK_COUNTER cannam@167: #endif cannam@167: cannam@167: /*----------------------------------------------------------------*/ cannam@167: /* Cray UNICOS _rtc() intrinsic function */ cannam@167: #if defined(HAVE__RTC) && !defined(HAVE_TICK_COUNTER) cannam@167: #ifdef HAVE_INTRINSICS_H cannam@167: # include cannam@167: #endif cannam@167: cannam@167: typedef long long ticks; cannam@167: cannam@167: #define getticks _rtc cannam@167: cannam@167: INLINE_ELAPSED(inline) cannam@167: cannam@167: #define HAVE_TICK_COUNTER cannam@167: #endif cannam@167: cannam@167: /*----------------------------------------------------------------*/ cannam@167: /* MIPS ZBus */ cannam@167: #if HAVE_MIPS_ZBUS_TIMER cannam@167: #if defined(__mips__) && !defined(HAVE_TICK_COUNTER) cannam@167: #include cannam@167: #include cannam@167: #include cannam@167: cannam@167: typedef uint64_t ticks; cannam@167: cannam@167: static inline ticks getticks(void) cannam@167: { cannam@167: static uint64_t* addr = 0; cannam@167: cannam@167: if (addr == 0) cannam@167: { cannam@167: uint32_t rq_addr = 0x10030000; cannam@167: int fd; cannam@167: int pgsize; cannam@167: cannam@167: pgsize = getpagesize(); cannam@167: fd = open ("/dev/mem", O_RDONLY | O_SYNC, 0); cannam@167: if (fd < 0) { cannam@167: perror("open"); cannam@167: return NULL; cannam@167: } cannam@167: addr = mmap(0, pgsize, PROT_READ, MAP_SHARED, fd, rq_addr); cannam@167: close(fd); cannam@167: if (addr == (uint64_t *)-1) { cannam@167: perror("mmap"); cannam@167: return NULL; cannam@167: } cannam@167: } cannam@167: cannam@167: return *addr; cannam@167: } cannam@167: cannam@167: INLINE_ELAPSED(inline) cannam@167: cannam@167: #define HAVE_TICK_COUNTER cannam@167: #endif cannam@167: #endif /* HAVE_MIPS_ZBUS_TIMER */ cannam@167: cannam@167: #if defined(HAVE_ARMV7A_CNTVCT) cannam@167: typedef uint64_t ticks; cannam@167: static inline ticks getticks(void) cannam@167: { cannam@167: uint32_t Rt, Rt2 = 0; cannam@167: asm volatile("mrrc p15, 1, %0, %1, c14" : "=r"(Rt), "=r"(Rt2)); cannam@167: return ((uint64_t)Rt) | (((uint64_t)Rt2) << 32); cannam@167: } cannam@167: INLINE_ELAPSED(inline) cannam@167: #define HAVE_TICK_COUNTER cannam@167: #endif cannam@167: cannam@167: #if defined(HAVE_ARMV7A_PMCCNTR) cannam@167: typedef uint64_t ticks; cannam@167: static inline ticks getticks(void) cannam@167: { cannam@167: uint32_t r; cannam@167: asm volatile("mrc p15, 0, %0, c9, c13, 0" : "=r"(r) ); cannam@167: return r; cannam@167: } cannam@167: INLINE_ELAPSED(inline) cannam@167: #define HAVE_TICK_COUNTER cannam@167: #endif cannam@167: cannam@167: #if defined(__aarch64__) && defined(HAVE_ARMV8_CNTVCT_EL0) && !defined(HAVE_ARMV8_PMCCNTR_EL0) cannam@167: typedef uint64_t ticks; cannam@167: static inline ticks getticks(void) cannam@167: { cannam@167: uint64_t Rt; cannam@167: asm volatile("mrs %0, CNTVCT_EL0" : "=r" (Rt)); cannam@167: return Rt; cannam@167: } cannam@167: INLINE_ELAPSED(inline) cannam@167: #define HAVE_TICK_COUNTER cannam@167: #endif cannam@167: cannam@167: #if defined(__aarch64__) && defined(HAVE_ARMV8_PMCCNTR_EL0) cannam@167: typedef uint64_t ticks; cannam@167: static inline ticks getticks(void) cannam@167: { cannam@167: uint64_t cc = 0; cannam@167: asm volatile("mrs %0, PMCCNTR_EL0" : "=r"(cc)); cannam@167: return cc; cannam@167: } cannam@167: INLINE_ELAPSED(inline) cannam@167: #define HAVE_TICK_COUNTER cannam@167: #endif