1//===================================================== 2// File : portable_perf_analyzer.hh 3// Author : L. Plagne <laurent.plagne@edf.fr)> 4// Copyright (C) EDF R&D, mar d�c 3 18:59:35 CET 2002 5//===================================================== 6// 7// This program is free software; you can redistribute it and/or 8// modify it under the terms of the GNU General Public License 9// as published by the Free Software Foundation; either version 2 10// of the License, or (at your option) any later version. 11// 12// This program is distributed in the hope that it will be useful, 13// but WITHOUT ANY WARRANTY; without even the implied warranty of 14// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 15// GNU General Public License for more details. 16// You should have received a copy of the GNU General Public License 17// along with this program; if not, write to the Free Software 18// Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. 19// 20#ifndef _PORTABLE_PERF_ANALYZER_HH 21#define _PORTABLE_PERF_ANALYZER_HH 22 23#include "utilities.h" 24#include "timers/portable_timer.hh" 25 26template <class Action> 27class Portable_Perf_Analyzer{ 28public: 29 Portable_Perf_Analyzer( void ):_nb_calc(1),_nb_init(1),_chronos(){ 30 MESSAGE("Portable_Perf_Analyzer Ctor"); 31 }; 32 Portable_Perf_Analyzer( const Portable_Perf_Analyzer & ){ 33 INFOS("Copy Ctor not implemented"); 34 exit(0); 35 }; 36 ~Portable_Perf_Analyzer( void ){ 37 MESSAGE("Portable_Perf_Analyzer Dtor"); 38 }; 39 40 41 42 inline double eval_mflops(int size) 43 { 44 45 Action action(size); 46 47// double time_baseline = time_init(action); 48// while (time_baseline < MIN_TIME_INIT) 49// { 50// _nb_init *= 2; 51// time_baseline = time_init(action); 52// } 53// 54// // optimize 55// for (int i=1; i<NB_TRIES; ++i) 56// time_baseline = std::min(time_baseline, time_init(action)); 57// 58// time_baseline = time_baseline/(double(_nb_init)); 59 60 double time_action = time_calculate(action); 61 while (time_action < MIN_TIME) 62 { 63 _nb_calc *= 2; 64 time_action = time_calculate(action); 65 } 66 67 // optimize 68 for (int i=1; i<NB_TRIES; ++i) 69 time_action = std::min(time_action, time_calculate(action)); 70 71// INFOS("size="<<size); 72// INFOS("_nb_init="<<_nb_init); 73// INFOS("_nb_calc="<<_nb_calc); 74 75 time_action = time_action / (double(_nb_calc)); 76 77 action.check_result(); 78 79 80 double time_baseline = time_init(action); 81 for (int i=1; i<NB_TRIES; ++i) 82 time_baseline = std::min(time_baseline, time_init(action)); 83 time_baseline = time_baseline/(double(_nb_init)); 84 85 86 87// INFOS("time_baseline="<<time_baseline); 88// INFOS("time_action="<<time_action); 89 90 time_action = time_action - time_baseline; 91 92// INFOS("time_corrected="<<time_action); 93 94 return action.nb_op_base()/(time_action*1000000.0); 95 } 96 97 inline double time_init(Action & action) 98 { 99 // time measurement 100 _chronos.start(); 101 for (int ii=0; ii<_nb_init; ii++) 102 action.initialize(); 103 _chronos.stop(); 104 return _chronos.user_time(); 105 } 106 107 108 inline double time_calculate(Action & action) 109 { 110 // time measurement 111 _chronos.start(); 112 for (int ii=0;ii<_nb_calc;ii++) 113 { 114 action.initialize(); 115 action.calculate(); 116 } 117 _chronos.stop(); 118 return _chronos.user_time(); 119 } 120 121 unsigned long long get_nb_calc( void ) 122 { 123 return _nb_calc; 124 } 125 126 127private: 128 unsigned long long _nb_calc; 129 unsigned long long _nb_init; 130 Portable_Timer _chronos; 131 132}; 133 134#endif //_PORTABLE_PERF_ANALYZER_HH 135