14cc4bb4cSJim Cownie /** @file kmp_stats_timing.cpp
24cc4bb4cSJim Cownie  * Timing functions
34cc4bb4cSJim Cownie  */
44cc4bb4cSJim Cownie 
54cc4bb4cSJim Cownie 
64cc4bb4cSJim Cownie //===----------------------------------------------------------------------===//
74cc4bb4cSJim Cownie //
84cc4bb4cSJim Cownie //                     The LLVM Compiler Infrastructure
94cc4bb4cSJim Cownie //
104cc4bb4cSJim Cownie // This file is dual licensed under the MIT and the University of Illinois Open
114cc4bb4cSJim Cownie // Source Licenses. See LICENSE.txt for details.
124cc4bb4cSJim Cownie //
134cc4bb4cSJim Cownie //===----------------------------------------------------------------------===//
144cc4bb4cSJim Cownie 
154cc4bb4cSJim Cownie 
164cc4bb4cSJim Cownie #include <stdlib.h>
174cc4bb4cSJim Cownie #include <unistd.h>
184cc4bb4cSJim Cownie 
194cc4bb4cSJim Cownie #include <iostream>
204cc4bb4cSJim Cownie #include <iomanip>
214cc4bb4cSJim Cownie #include <sstream>
224cc4bb4cSJim Cownie 
234cc4bb4cSJim Cownie #include "kmp_stats_timing.h"
244cc4bb4cSJim Cownie 
254cc4bb4cSJim Cownie using namespace std;
264cc4bb4cSJim Cownie 
274cc4bb4cSJim Cownie #if KMP_OS_LINUX
284cc4bb4cSJim Cownie # if KMP_MIC
294cc4bb4cSJim Cownie double tsc_tick_count::tick_time()
304cc4bb4cSJim Cownie {
314cc4bb4cSJim Cownie     // pretty bad assumption of 1GHz clock for MIC
324cc4bb4cSJim Cownie     return 1/((double)1000*1.e6);
334cc4bb4cSJim Cownie }
344cc4bb4cSJim Cownie # else
354cc4bb4cSJim Cownie #  include <string.h>
364cc4bb4cSJim Cownie // Extract the value from the CPUID information
374cc4bb4cSJim Cownie double tsc_tick_count::tick_time()
384cc4bb4cSJim Cownie {
394cc4bb4cSJim Cownie     static double result = 0.0;
404cc4bb4cSJim Cownie 
414cc4bb4cSJim Cownie     if (result == 0.0)
424cc4bb4cSJim Cownie     {
434cc4bb4cSJim Cownie         int cpuinfo[4];
444cc4bb4cSJim Cownie         char brand[256];
454cc4bb4cSJim Cownie 
464cc4bb4cSJim Cownie         __cpuid(cpuinfo, 0x80000000);
474cc4bb4cSJim Cownie         memset(brand, 0, sizeof(brand));
484cc4bb4cSJim Cownie         int ids = cpuinfo[0];
494cc4bb4cSJim Cownie 
504cc4bb4cSJim Cownie         for (unsigned int i=2; i<(ids^0x80000000)+2; i++)
514cc4bb4cSJim Cownie             __cpuid(brand+(i-2)*sizeof(cpuinfo), i | 0x80000000);
524cc4bb4cSJim Cownie 
534cc4bb4cSJim Cownie         char * start = &brand[0];
544cc4bb4cSJim Cownie         for (;*start == ' '; start++)
554cc4bb4cSJim Cownie             ;
564cc4bb4cSJim Cownie 
57*74bf17b8SAndrey Churbanov         char * end = brand + KMP_STRLEN(brand) - 3;
584cc4bb4cSJim Cownie         uint64_t multiplier;
594cc4bb4cSJim Cownie 
604cc4bb4cSJim Cownie         if (*end == 'M') multiplier = 1000LL*1000LL;
614cc4bb4cSJim Cownie         else if (*end == 'G') multiplier = 1000LL*1000LL*1000LL;
624cc4bb4cSJim Cownie         else if (*end == 'T') multiplier = 1000LL*1000LL*1000LL*1000LL;
634cc4bb4cSJim Cownie         else
644cc4bb4cSJim Cownie         {
654cc4bb4cSJim Cownie             cout << "Error determining multiplier '" << *end << "'\n";
664cc4bb4cSJim Cownie             exit (-1);
674cc4bb4cSJim Cownie         }
684cc4bb4cSJim Cownie         *end = 0;
694cc4bb4cSJim Cownie         while (*end != ' ') end--;
704cc4bb4cSJim Cownie         end++;
714cc4bb4cSJim Cownie 
724cc4bb4cSJim Cownie         double freq = strtod(end, &start);
734cc4bb4cSJim Cownie         if (freq == 0.0)
744cc4bb4cSJim Cownie         {
754cc4bb4cSJim Cownie             cout << "Error calculating frequency " <<  end << "\n";
764cc4bb4cSJim Cownie             exit (-1);
774cc4bb4cSJim Cownie         }
784cc4bb4cSJim Cownie 
794cc4bb4cSJim Cownie         result = ((double)1.0)/(freq * multiplier);
804cc4bb4cSJim Cownie     }
814cc4bb4cSJim Cownie     return result;
824cc4bb4cSJim Cownie }
834cc4bb4cSJim Cownie # endif
844cc4bb4cSJim Cownie #endif
854cc4bb4cSJim Cownie 
864cc4bb4cSJim Cownie static bool useSI = true;
874cc4bb4cSJim Cownie 
884cc4bb4cSJim Cownie // Return a formatted string after normalising the value into
894cc4bb4cSJim Cownie // engineering style and using a suitable unit prefix (e.g. ms, us, ns).
904cc4bb4cSJim Cownie std::string formatSI(double interval, int width, char unit)
914cc4bb4cSJim Cownie {
924cc4bb4cSJim Cownie     std::stringstream os;
934cc4bb4cSJim Cownie 
944cc4bb4cSJim Cownie     if (useSI)
954cc4bb4cSJim Cownie     {
964cc4bb4cSJim Cownie         // Preserve accuracy for small numbers, since we only multiply and the positive powers
974cc4bb4cSJim Cownie         // of ten are precisely representable.
984cc4bb4cSJim Cownie         static struct { double scale; char prefix; } ranges[] = {
994cc4bb4cSJim Cownie             {1.e12,'f'},
1004cc4bb4cSJim Cownie             {1.e9, 'p'},
1014cc4bb4cSJim Cownie             {1.e6, 'n'},
1024cc4bb4cSJim Cownie             {1.e3, 'u'},
1034cc4bb4cSJim Cownie             {1.0,  'm'},
1044cc4bb4cSJim Cownie             {1.e-3,' '},
1054cc4bb4cSJim Cownie             {1.e-6,'k'},
1064cc4bb4cSJim Cownie             {1.e-9,'M'},
1074cc4bb4cSJim Cownie             {1.e-12,'G'},
1084cc4bb4cSJim Cownie             {1.e-15,'T'},
1094cc4bb4cSJim Cownie             {1.e-18,'P'},
1104cc4bb4cSJim Cownie             {1.e-21,'E'},
1114cc4bb4cSJim Cownie             {1.e-24,'Z'},
1124cc4bb4cSJim Cownie             {1.e-27,'Y'}
1134cc4bb4cSJim Cownie         };
1144cc4bb4cSJim Cownie 
1154cc4bb4cSJim Cownie         if (interval == 0.0)
1164cc4bb4cSJim Cownie         {
1174cc4bb4cSJim Cownie             os << std::setw(width-3) << std::right << "0.00" << std::setw(3) << unit;
1184cc4bb4cSJim Cownie             return os.str();
1194cc4bb4cSJim Cownie         }
1204cc4bb4cSJim Cownie 
1214cc4bb4cSJim Cownie         bool negative = false;
1224cc4bb4cSJim Cownie         if (interval < 0.0)
1234cc4bb4cSJim Cownie         {
1244cc4bb4cSJim Cownie             negative = true;
1254cc4bb4cSJim Cownie             interval = -interval;
1264cc4bb4cSJim Cownie         }
1274cc4bb4cSJim Cownie 
1284cc4bb4cSJim Cownie         for (int i=0; i<(int)(sizeof(ranges)/sizeof(ranges[0])); i++)
1294cc4bb4cSJim Cownie         {
1304cc4bb4cSJim Cownie             if (interval*ranges[i].scale < 1.e0)
1314cc4bb4cSJim Cownie             {
1324cc4bb4cSJim Cownie                 interval = interval * 1000.e0 * ranges[i].scale;
1334cc4bb4cSJim Cownie                 os << std::fixed << std::setprecision(2) << std::setw(width-3) << std::right <<
1344cc4bb4cSJim Cownie                     (negative ? -interval : interval) << std::setw(2) << ranges[i].prefix << std::setw(1) << unit;
1354cc4bb4cSJim Cownie 
1364cc4bb4cSJim Cownie                 return os.str();
1374cc4bb4cSJim Cownie             }
1384cc4bb4cSJim Cownie         }
1394cc4bb4cSJim Cownie     }
1404cc4bb4cSJim Cownie     os << std::setprecision(2) << std::fixed << std::right << std::setw(width-3) << interval << std::setw(3) << unit;
1414cc4bb4cSJim Cownie 
1424cc4bb4cSJim Cownie     return os.str();
1434cc4bb4cSJim Cownie }
1444cc4bb4cSJim Cownie 
1454cc4bb4cSJim Cownie tsc_tick_count::tsc_interval_t computeLastInLastOutInterval(timePair * times, int nTimes)
1464cc4bb4cSJim Cownie {
1474cc4bb4cSJim Cownie     timePair lastTimes = times[0];
1484cc4bb4cSJim Cownie     tsc_tick_count * startp = lastTimes.get_startp();
1494cc4bb4cSJim Cownie     tsc_tick_count * endp   = lastTimes.get_endp();
1504cc4bb4cSJim Cownie 
1514cc4bb4cSJim Cownie     for (int i=1; i<nTimes; i++)
1524cc4bb4cSJim Cownie     {
1534cc4bb4cSJim Cownie        (*startp) = startp->later(times[i].get_start());
1544cc4bb4cSJim Cownie        (*endp)   = endp->later  (times[i].get_end());
1554cc4bb4cSJim Cownie     }
1564cc4bb4cSJim Cownie 
1574cc4bb4cSJim Cownie     return lastTimes.duration();
1584cc4bb4cSJim Cownie }
1594cc4bb4cSJim Cownie 
1604cc4bb4cSJim Cownie std::string timePair::format() const
1614cc4bb4cSJim Cownie {
1624cc4bb4cSJim Cownie     std::ostringstream oss;
1634cc4bb4cSJim Cownie 
1644cc4bb4cSJim Cownie     oss << start.getValue() << ":" << end.getValue() << " = " << (end-start).getValue();
1654cc4bb4cSJim Cownie 
1664cc4bb4cSJim Cownie     return oss.str();
1674cc4bb4cSJim Cownie }
168