1*4cc4bb4cSJim Cownie /** @file kmp_stats_timing.cpp
2*4cc4bb4cSJim Cownie  * Timing functions
3*4cc4bb4cSJim Cownie  */
4*4cc4bb4cSJim Cownie 
5*4cc4bb4cSJim Cownie 
6*4cc4bb4cSJim Cownie //===----------------------------------------------------------------------===//
7*4cc4bb4cSJim Cownie //
8*4cc4bb4cSJim Cownie //                     The LLVM Compiler Infrastructure
9*4cc4bb4cSJim Cownie //
10*4cc4bb4cSJim Cownie // This file is dual licensed under the MIT and the University of Illinois Open
11*4cc4bb4cSJim Cownie // Source Licenses. See LICENSE.txt for details.
12*4cc4bb4cSJim Cownie //
13*4cc4bb4cSJim Cownie //===----------------------------------------------------------------------===//
14*4cc4bb4cSJim Cownie 
15*4cc4bb4cSJim Cownie 
16*4cc4bb4cSJim Cownie #include <stdlib.h>
17*4cc4bb4cSJim Cownie #include <unistd.h>
18*4cc4bb4cSJim Cownie 
19*4cc4bb4cSJim Cownie #include <iostream>
20*4cc4bb4cSJim Cownie #include <iomanip>
21*4cc4bb4cSJim Cownie #include <sstream>
22*4cc4bb4cSJim Cownie 
23*4cc4bb4cSJim Cownie #include "kmp_stats_timing.h"
24*4cc4bb4cSJim Cownie 
25*4cc4bb4cSJim Cownie using namespace std;
26*4cc4bb4cSJim Cownie 
27*4cc4bb4cSJim Cownie #if KMP_OS_LINUX
28*4cc4bb4cSJim Cownie # if KMP_MIC
29*4cc4bb4cSJim Cownie double tsc_tick_count::tick_time()
30*4cc4bb4cSJim Cownie {
31*4cc4bb4cSJim Cownie     // pretty bad assumption of 1GHz clock for MIC
32*4cc4bb4cSJim Cownie     return 1/((double)1000*1.e6);
33*4cc4bb4cSJim Cownie }
34*4cc4bb4cSJim Cownie # else
35*4cc4bb4cSJim Cownie #  include <string.h>
36*4cc4bb4cSJim Cownie // Extract the value from the CPUID information
37*4cc4bb4cSJim Cownie double tsc_tick_count::tick_time()
38*4cc4bb4cSJim Cownie {
39*4cc4bb4cSJim Cownie     static double result = 0.0;
40*4cc4bb4cSJim Cownie 
41*4cc4bb4cSJim Cownie     if (result == 0.0)
42*4cc4bb4cSJim Cownie     {
43*4cc4bb4cSJim Cownie         int cpuinfo[4];
44*4cc4bb4cSJim Cownie         char brand[256];
45*4cc4bb4cSJim Cownie 
46*4cc4bb4cSJim Cownie         __cpuid(cpuinfo, 0x80000000);
47*4cc4bb4cSJim Cownie         memset(brand, 0, sizeof(brand));
48*4cc4bb4cSJim Cownie         int ids = cpuinfo[0];
49*4cc4bb4cSJim Cownie 
50*4cc4bb4cSJim Cownie         for (unsigned int i=2; i<(ids^0x80000000)+2; i++)
51*4cc4bb4cSJim Cownie             __cpuid(brand+(i-2)*sizeof(cpuinfo), i | 0x80000000);
52*4cc4bb4cSJim Cownie 
53*4cc4bb4cSJim Cownie         char * start = &brand[0];
54*4cc4bb4cSJim Cownie         for (;*start == ' '; start++)
55*4cc4bb4cSJim Cownie             ;
56*4cc4bb4cSJim Cownie 
57*4cc4bb4cSJim Cownie         char * end = brand + strlen(brand) - 3;
58*4cc4bb4cSJim Cownie         uint64_t multiplier;
59*4cc4bb4cSJim Cownie 
60*4cc4bb4cSJim Cownie         if (*end == 'M') multiplier = 1000LL*1000LL;
61*4cc4bb4cSJim Cownie         else if (*end == 'G') multiplier = 1000LL*1000LL*1000LL;
62*4cc4bb4cSJim Cownie         else if (*end == 'T') multiplier = 1000LL*1000LL*1000LL*1000LL;
63*4cc4bb4cSJim Cownie         else
64*4cc4bb4cSJim Cownie         {
65*4cc4bb4cSJim Cownie             cout << "Error determining multiplier '" << *end << "'\n";
66*4cc4bb4cSJim Cownie             exit (-1);
67*4cc4bb4cSJim Cownie         }
68*4cc4bb4cSJim Cownie         *end = 0;
69*4cc4bb4cSJim Cownie         while (*end != ' ') end--;
70*4cc4bb4cSJim Cownie         end++;
71*4cc4bb4cSJim Cownie 
72*4cc4bb4cSJim Cownie         double freq = strtod(end, &start);
73*4cc4bb4cSJim Cownie         if (freq == 0.0)
74*4cc4bb4cSJim Cownie         {
75*4cc4bb4cSJim Cownie             cout << "Error calculating frequency " <<  end << "\n";
76*4cc4bb4cSJim Cownie             exit (-1);
77*4cc4bb4cSJim Cownie         }
78*4cc4bb4cSJim Cownie 
79*4cc4bb4cSJim Cownie         result = ((double)1.0)/(freq * multiplier);
80*4cc4bb4cSJim Cownie     }
81*4cc4bb4cSJim Cownie     return result;
82*4cc4bb4cSJim Cownie }
83*4cc4bb4cSJim Cownie # endif
84*4cc4bb4cSJim Cownie #endif
85*4cc4bb4cSJim Cownie 
86*4cc4bb4cSJim Cownie static bool useSI = true;
87*4cc4bb4cSJim Cownie 
88*4cc4bb4cSJim Cownie // Return a formatted string after normalising the value into
89*4cc4bb4cSJim Cownie // engineering style and using a suitable unit prefix (e.g. ms, us, ns).
90*4cc4bb4cSJim Cownie std::string formatSI(double interval, int width, char unit)
91*4cc4bb4cSJim Cownie {
92*4cc4bb4cSJim Cownie     std::stringstream os;
93*4cc4bb4cSJim Cownie 
94*4cc4bb4cSJim Cownie     if (useSI)
95*4cc4bb4cSJim Cownie     {
96*4cc4bb4cSJim Cownie         // Preserve accuracy for small numbers, since we only multiply and the positive powers
97*4cc4bb4cSJim Cownie         // of ten are precisely representable.
98*4cc4bb4cSJim Cownie         static struct { double scale; char prefix; } ranges[] = {
99*4cc4bb4cSJim Cownie             {1.e12,'f'},
100*4cc4bb4cSJim Cownie             {1.e9, 'p'},
101*4cc4bb4cSJim Cownie             {1.e6, 'n'},
102*4cc4bb4cSJim Cownie             {1.e3, 'u'},
103*4cc4bb4cSJim Cownie             {1.0,  'm'},
104*4cc4bb4cSJim Cownie             {1.e-3,' '},
105*4cc4bb4cSJim Cownie             {1.e-6,'k'},
106*4cc4bb4cSJim Cownie             {1.e-9,'M'},
107*4cc4bb4cSJim Cownie             {1.e-12,'G'},
108*4cc4bb4cSJim Cownie             {1.e-15,'T'},
109*4cc4bb4cSJim Cownie             {1.e-18,'P'},
110*4cc4bb4cSJim Cownie             {1.e-21,'E'},
111*4cc4bb4cSJim Cownie             {1.e-24,'Z'},
112*4cc4bb4cSJim Cownie             {1.e-27,'Y'}
113*4cc4bb4cSJim Cownie         };
114*4cc4bb4cSJim Cownie 
115*4cc4bb4cSJim Cownie         if (interval == 0.0)
116*4cc4bb4cSJim Cownie         {
117*4cc4bb4cSJim Cownie             os << std::setw(width-3) << std::right << "0.00" << std::setw(3) << unit;
118*4cc4bb4cSJim Cownie             return os.str();
119*4cc4bb4cSJim Cownie         }
120*4cc4bb4cSJim Cownie 
121*4cc4bb4cSJim Cownie         bool negative = false;
122*4cc4bb4cSJim Cownie         if (interval < 0.0)
123*4cc4bb4cSJim Cownie         {
124*4cc4bb4cSJim Cownie             negative = true;
125*4cc4bb4cSJim Cownie             interval = -interval;
126*4cc4bb4cSJim Cownie         }
127*4cc4bb4cSJim Cownie 
128*4cc4bb4cSJim Cownie         for (int i=0; i<(int)(sizeof(ranges)/sizeof(ranges[0])); i++)
129*4cc4bb4cSJim Cownie         {
130*4cc4bb4cSJim Cownie             if (interval*ranges[i].scale < 1.e0)
131*4cc4bb4cSJim Cownie             {
132*4cc4bb4cSJim Cownie                 interval = interval * 1000.e0 * ranges[i].scale;
133*4cc4bb4cSJim Cownie                 os << std::fixed << std::setprecision(2) << std::setw(width-3) << std::right <<
134*4cc4bb4cSJim Cownie                     (negative ? -interval : interval) << std::setw(2) << ranges[i].prefix << std::setw(1) << unit;
135*4cc4bb4cSJim Cownie 
136*4cc4bb4cSJim Cownie                 return os.str();
137*4cc4bb4cSJim Cownie             }
138*4cc4bb4cSJim Cownie         }
139*4cc4bb4cSJim Cownie     }
140*4cc4bb4cSJim Cownie     os << std::setprecision(2) << std::fixed << std::right << std::setw(width-3) << interval << std::setw(3) << unit;
141*4cc4bb4cSJim Cownie 
142*4cc4bb4cSJim Cownie     return os.str();
143*4cc4bb4cSJim Cownie }
144*4cc4bb4cSJim Cownie 
145*4cc4bb4cSJim Cownie tsc_tick_count::tsc_interval_t computeLastInLastOutInterval(timePair * times, int nTimes)
146*4cc4bb4cSJim Cownie {
147*4cc4bb4cSJim Cownie     timePair lastTimes = times[0];
148*4cc4bb4cSJim Cownie     tsc_tick_count * startp = lastTimes.get_startp();
149*4cc4bb4cSJim Cownie     tsc_tick_count * endp   = lastTimes.get_endp();
150*4cc4bb4cSJim Cownie 
151*4cc4bb4cSJim Cownie     for (int i=1; i<nTimes; i++)
152*4cc4bb4cSJim Cownie     {
153*4cc4bb4cSJim Cownie        (*startp) = startp->later(times[i].get_start());
154*4cc4bb4cSJim Cownie        (*endp)   = endp->later  (times[i].get_end());
155*4cc4bb4cSJim Cownie     }
156*4cc4bb4cSJim Cownie 
157*4cc4bb4cSJim Cownie     return lastTimes.duration();
158*4cc4bb4cSJim Cownie }
159*4cc4bb4cSJim Cownie 
160*4cc4bb4cSJim Cownie std::string timePair::format() const
161*4cc4bb4cSJim Cownie {
162*4cc4bb4cSJim Cownie     std::ostringstream oss;
163*4cc4bb4cSJim Cownie 
164*4cc4bb4cSJim Cownie     oss << start.getValue() << ":" << end.getValue() << " = " << (end-start).getValue();
165*4cc4bb4cSJim Cownie 
166*4cc4bb4cSJim Cownie     return oss.str();
167*4cc4bb4cSJim Cownie }
168