1c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// This file is part of Eigen, a lightweight C++ template library
2c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// for linear algebra.
3c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath//
4c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// Copyright (C) 2008-2010 Gael Guennebaud <gael.guennebaud@inria.fr>
5c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// Copyright (C) 2009 Benoit Jacob <jacob.benoit.1@gmail.com>
6c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath//
7c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// This Source Code Form is subject to the terms of the Mozilla
8c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// Public License v. 2.0. If a copy of the MPL was not distributed
9c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// with this file, You can obtain one at http://mozilla.org/MPL/2.0/.
10c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath
11c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#ifndef EIGEN_BENCH_TIMERR_H
12c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#define EIGEN_BENCH_TIMERR_H
13c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath
14c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#if defined(_WIN32) || defined(__CYGWIN__)
15c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath# ifndef NOMINMAX
16c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#   define NOMINMAX
17c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#   define EIGEN_BT_UNDEF_NOMINMAX
18c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath# endif
19c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath# ifndef WIN32_LEAN_AND_MEAN
20c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#   define WIN32_LEAN_AND_MEAN
21c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#   define EIGEN_BT_UNDEF_WIN32_LEAN_AND_MEAN
22c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath# endif
23c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath# include <windows.h>
24c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#elif defined(__APPLE__)
25c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#include <CoreServices/CoreServices.h>
26c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#include <mach/mach_time.h>
27c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#else
28c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath# include <unistd.h>
29c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#endif
30c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath
31c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#include <Eigen/Core>
32c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath
33c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathnamespace Eigen
34c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath{
35c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath
36c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathenum {
37c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  CPU_TIMER = 0,
38c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  REAL_TIMER = 1
39c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath};
40c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath
41c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath/** Elapsed time timer keeping the best try.
42c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  *
43c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  * On POSIX platforms we use clock_gettime with CLOCK_PROCESS_CPUTIME_ID.
44c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  * On Windows we use QueryPerformanceCounter
45c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  *
46c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  * Important: on linux, you must link with -lrt
47c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  */
48c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathclass BenchTimer
49c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath{
50c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathpublic:
51c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath
52c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  BenchTimer()
53c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  {
54c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#if defined(_WIN32) || defined(__CYGWIN__)
55c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    LARGE_INTEGER freq;
56c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    QueryPerformanceFrequency(&freq);
57c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    m_frequency = (double)freq.QuadPart;
58c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#endif
59c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    reset();
60c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  }
61c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath
62c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  ~BenchTimer() {}
63c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath
64c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  inline void reset()
65c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  {
66c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    m_bests.fill(1e9);
67c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    m_worsts.fill(0);
68c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    m_totals.setZero();
69c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  }
70c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  inline void start()
71c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  {
72c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    m_starts[CPU_TIMER]  = getCpuTime();
73c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    m_starts[REAL_TIMER] = getRealTime();
74c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  }
75c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  inline void stop()
76c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  {
77c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    m_times[CPU_TIMER] = getCpuTime() - m_starts[CPU_TIMER];
78c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    m_times[REAL_TIMER] = getRealTime() - m_starts[REAL_TIMER];
79c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    #if EIGEN_VERSION_AT_LEAST(2,90,0)
80c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    m_bests = m_bests.cwiseMin(m_times);
81c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    m_worsts = m_worsts.cwiseMax(m_times);
82c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    #else
83c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    m_bests(0) = std::min(m_bests(0),m_times(0));
84c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    m_bests(1) = std::min(m_bests(1),m_times(1));
85c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    m_worsts(0) = std::max(m_worsts(0),m_times(0));
86c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    m_worsts(1) = std::max(m_worsts(1),m_times(1));
87c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    #endif
88c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    m_totals += m_times;
89c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  }
90c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath
91c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  /** Return the elapsed time in seconds between the last start/stop pair
92c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    */
93c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  inline double value(int TIMER = CPU_TIMER) const
94c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  {
95c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    return m_times[TIMER];
96c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  }
97c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath
98c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  /** Return the best elapsed time in seconds
99c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    */
100c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  inline double best(int TIMER = CPU_TIMER) const
101c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  {
102c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    return m_bests[TIMER];
103c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  }
104c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath
105c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  /** Return the worst elapsed time in seconds
106c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    */
107c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  inline double worst(int TIMER = CPU_TIMER) const
108c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  {
109c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    return m_worsts[TIMER];
110c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  }
111c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath
112c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  /** Return the total elapsed time in seconds.
113c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    */
114c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  inline double total(int TIMER = CPU_TIMER) const
115c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  {
116c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    return m_totals[TIMER];
117c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  }
118c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath
119c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  inline double getCpuTime() const
120c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  {
121c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#ifdef _WIN32
122c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    LARGE_INTEGER query_ticks;
123c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    QueryPerformanceCounter(&query_ticks);
124c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    return query_ticks.QuadPart/m_frequency;
125c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#elif __APPLE__
126c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    return double(mach_absolute_time())*1e-9;
127c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#else
128c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    timespec ts;
129c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    clock_gettime(CLOCK_PROCESS_CPUTIME_ID, &ts);
130c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    return double(ts.tv_sec) + 1e-9 * double(ts.tv_nsec);
131c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#endif
132c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  }
133c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath
134c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  inline double getRealTime() const
135c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  {
136c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#ifdef _WIN32
137c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    SYSTEMTIME st;
138c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    GetSystemTime(&st);
139c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    return (double)st.wSecond + 1.e-3 * (double)st.wMilliseconds;
140c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#elif __APPLE__
141c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    return double(mach_absolute_time())*1e-9;
142c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#else
143c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    timespec ts;
144c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    clock_gettime(CLOCK_REALTIME, &ts);
145c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    return double(ts.tv_sec) + 1e-9 * double(ts.tv_nsec);
146c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#endif
147c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  }
148c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath
149c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathprotected:
150c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#if defined(_WIN32) || defined(__CYGWIN__)
151c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  double m_frequency;
152c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#endif
153c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  Vector2d m_starts;
154c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  Vector2d m_times;
155c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  Vector2d m_bests;
156c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  Vector2d m_worsts;
157c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  Vector2d m_totals;
158c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath
159c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathpublic:
160c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  EIGEN_MAKE_ALIGNED_OPERATOR_NEW
161c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath};
162c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath
163c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#define BENCH(TIMER,TRIES,REP,CODE) { \
164c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    TIMER.reset(); \
165c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    for(int uglyvarname1=0; uglyvarname1<TRIES; ++uglyvarname1){ \
166c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath      TIMER.start(); \
167c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath      for(int uglyvarname2=0; uglyvarname2<REP; ++uglyvarname2){ \
168c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath        CODE; \
169c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath      } \
170c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath      TIMER.stop(); \
171c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath    } \
172c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  }
173c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath
174c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath}
175c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath
176c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// clean #defined tokens
177c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#ifdef EIGEN_BT_UNDEF_NOMINMAX
178c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath# undef EIGEN_BT_UNDEF_NOMINMAX
179c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath# undef NOMINMAX
180c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#endif
181c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath
182c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#ifdef EIGEN_BT_UNDEF_WIN32_LEAN_AND_MEAN
183c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath# undef EIGEN_BT_UNDEF_WIN32_LEAN_AND_MEAN
184c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath# undef WIN32_LEAN_AND_MEAN
185c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#endif
186c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath
187c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#endif // EIGEN_BENCH_TIMERR_H
188