xref: /aosp_15_r20/external/eigen/bench/btl/data/regularize.cxx (revision bf2c37156dfe67e5dfebd6d394bad8b2ab5804d4)
1*bf2c3715SXin Li //=====================================================
2*bf2c3715SXin Li // File   :  regularize.cxx
3*bf2c3715SXin Li // Author :  L. Plagne <[email protected])>
4*bf2c3715SXin Li // Copyright (C) EDF R&D,  lun sep 30 14:23:15 CEST 2002
5*bf2c3715SXin Li //=====================================================
6*bf2c3715SXin Li //
7*bf2c3715SXin Li // This program is free software; you can redistribute it and/or
8*bf2c3715SXin Li // modify it under the terms of the GNU General Public License
9*bf2c3715SXin Li // as published by the Free Software Foundation; either version 2
10*bf2c3715SXin Li // of the License, or (at your option) any later version.
11*bf2c3715SXin Li //
12*bf2c3715SXin Li // This program is distributed in the hope that it will be useful,
13*bf2c3715SXin Li // but WITHOUT ANY WARRANTY; without even the implied warranty of
14*bf2c3715SXin Li // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
15*bf2c3715SXin Li // GNU General Public License for more details.
16*bf2c3715SXin Li // You should have received a copy of the GNU General Public License
17*bf2c3715SXin Li // along with this program; if not, write to the Free Software
18*bf2c3715SXin Li // Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.
19*bf2c3715SXin Li //
20*bf2c3715SXin Li #include "utilities.h"
21*bf2c3715SXin Li #include <vector>
22*bf2c3715SXin Li #include <string>
23*bf2c3715SXin Li #include <iostream>
24*bf2c3715SXin Li #include <fstream>
25*bf2c3715SXin Li #include "bench_parameter.hh"
26*bf2c3715SXin Li #include <set>
27*bf2c3715SXin Li 
28*bf2c3715SXin Li using namespace std;
29*bf2c3715SXin Li 
30*bf2c3715SXin Li void read_xy_file(const string & filename, vector<int> & tab_sizes, vector<double> & tab_mflops);
31*bf2c3715SXin Li void regularize_curve(const string & filename,
32*bf2c3715SXin Li 		      const vector<double> & tab_mflops,
33*bf2c3715SXin Li 		      const vector<int> & tab_sizes,
34*bf2c3715SXin Li 		      int start_cut_size, int stop_cut_size);
35*bf2c3715SXin Li /////////////////////////////////////////////////////////////////////////////////////////////////
36*bf2c3715SXin Li 
main(int argc,char * argv[])37*bf2c3715SXin Li int main( int argc , char *argv[] )
38*bf2c3715SXin Li {
39*bf2c3715SXin Li 
40*bf2c3715SXin Li   // input data
41*bf2c3715SXin Li 
42*bf2c3715SXin Li   if (argc<4){
43*bf2c3715SXin Li     INFOS("!!! Error ... usage : main filename start_cut_size stop_cut_size regularize_filename");
44*bf2c3715SXin Li     exit(0);
45*bf2c3715SXin Li   }
46*bf2c3715SXin Li   INFOS(argc);
47*bf2c3715SXin Li 
48*bf2c3715SXin Li   int start_cut_size=atoi(argv[2]);
49*bf2c3715SXin Li   int stop_cut_size=atoi(argv[3]);
50*bf2c3715SXin Li 
51*bf2c3715SXin Li   string filename=argv[1];
52*bf2c3715SXin Li   string regularize_filename=argv[4];
53*bf2c3715SXin Li 
54*bf2c3715SXin Li   INFOS(filename);
55*bf2c3715SXin Li   INFOS("start_cut_size="<<start_cut_size);
56*bf2c3715SXin Li 
57*bf2c3715SXin Li   vector<int> tab_sizes;
58*bf2c3715SXin Li   vector<double> tab_mflops;
59*bf2c3715SXin Li 
60*bf2c3715SXin Li   read_xy_file(filename,tab_sizes,tab_mflops);
61*bf2c3715SXin Li 
62*bf2c3715SXin Li   // regularizeing
63*bf2c3715SXin Li 
64*bf2c3715SXin Li   regularize_curve(regularize_filename,tab_mflops,tab_sizes,start_cut_size,stop_cut_size);
65*bf2c3715SXin Li 
66*bf2c3715SXin Li 
67*bf2c3715SXin Li }
68*bf2c3715SXin Li 
69*bf2c3715SXin Li //////////////////////////////////////////////////////////////////////////////////////
70*bf2c3715SXin Li 
regularize_curve(const string & filename,const vector<double> & tab_mflops,const vector<int> & tab_sizes,int start_cut_size,int stop_cut_size)71*bf2c3715SXin Li void regularize_curve(const string & filename,
72*bf2c3715SXin Li 		      const vector<double> & tab_mflops,
73*bf2c3715SXin Li 		      const vector<int> & tab_sizes,
74*bf2c3715SXin Li 		      int start_cut_size, int stop_cut_size)
75*bf2c3715SXin Li {
76*bf2c3715SXin Li   int size=tab_mflops.size();
77*bf2c3715SXin Li   ofstream output_file (filename.c_str(),ios::out) ;
78*bf2c3715SXin Li 
79*bf2c3715SXin Li   int i=0;
80*bf2c3715SXin Li 
81*bf2c3715SXin Li   while(tab_sizes[i]<start_cut_size){
82*bf2c3715SXin Li 
83*bf2c3715SXin Li     output_file << tab_sizes[i] << " " <<  tab_mflops[i] << endl ;
84*bf2c3715SXin Li     i++;
85*bf2c3715SXin Li 
86*bf2c3715SXin Li   }
87*bf2c3715SXin Li 
88*bf2c3715SXin Li   output_file << endl ;
89*bf2c3715SXin Li 
90*bf2c3715SXin Li   while(tab_sizes[i]<stop_cut_size){
91*bf2c3715SXin Li 
92*bf2c3715SXin Li     i++;
93*bf2c3715SXin Li 
94*bf2c3715SXin Li   }
95*bf2c3715SXin Li 
96*bf2c3715SXin Li   while(i<size){
97*bf2c3715SXin Li 
98*bf2c3715SXin Li     output_file << tab_sizes[i] << " " <<  tab_mflops[i] << endl ;
99*bf2c3715SXin Li     i++;
100*bf2c3715SXin Li 
101*bf2c3715SXin Li   }
102*bf2c3715SXin Li 
103*bf2c3715SXin Li   output_file.close();
104*bf2c3715SXin Li 
105*bf2c3715SXin Li }
106*bf2c3715SXin Li 
107*bf2c3715SXin Li ///////////////////////////////////////////////////////////////////////////////////////////////////////////////////
108*bf2c3715SXin Li 
read_xy_file(const string & filename,vector<int> & tab_sizes,vector<double> & tab_mflops)109*bf2c3715SXin Li void read_xy_file(const string & filename, vector<int> & tab_sizes, vector<double> & tab_mflops){
110*bf2c3715SXin Li 
111*bf2c3715SXin Li   ifstream input_file (filename.c_str(),ios::in) ;
112*bf2c3715SXin Li 
113*bf2c3715SXin Li   if (!input_file){
114*bf2c3715SXin Li     INFOS("!!! Error opening "<<filename);
115*bf2c3715SXin Li     exit(0);
116*bf2c3715SXin Li   }
117*bf2c3715SXin Li 
118*bf2c3715SXin Li   int nb_point=0;
119*bf2c3715SXin Li   int size=0;
120*bf2c3715SXin Li   double mflops=0;
121*bf2c3715SXin Li 
122*bf2c3715SXin Li   while (input_file >> size >> mflops ){
123*bf2c3715SXin Li     nb_point++;
124*bf2c3715SXin Li     tab_sizes.push_back(size);
125*bf2c3715SXin Li     tab_mflops.push_back(mflops);
126*bf2c3715SXin Li   }
127*bf2c3715SXin Li   SCRUTE(nb_point);
128*bf2c3715SXin Li 
129*bf2c3715SXin Li   input_file.close();
130*bf2c3715SXin Li }
131*bf2c3715SXin Li 
132