source: GTP/trunk/Lib/Vis/shared/EvalStats/EvalStats.cpp @ 744

Revision 744, 7.0 KB checked in by mattausch, 19 years ago (diff)
Line 
1// EvalStats.cpp : Defines the entry point for the console application.
2//
3#include "stdafx.h"
4
5#include <string>
6#include <vector>
7#include <istream>
8#include <iostream>
9#include <fstream>
10#include <math.h>
11#include <algorithm>
12
13using namespace std;
14
15/** This is a small function which takes two log files and computes the
16        difference in percent of the weighted render cost.
17*/
18
19
20struct RenderStats
21{
22        float mRenderCost;
23        float mAvgRenderCost;
24};
25
26
27typedef vector<RenderStats> StatsContainer;
28
29
30enum {NO_TAG, RENDER_COST, AVG_RENDER_COST};
31
32
33bool extractRenderStats(ifstream &file, StatsContainer &renderStats)
34{
35        if (!file.is_open())
36                return false;
37               
38        string buf;
39
40        int currentTag = NO_TAG;
41
42        while (!(getline(file, buf)).eof())
43        {
44                //cout << "buf: " << buf << endl;
45
46                if (buf[0] == '#')
47                {
48                        char entry[50];
49
50                        sscanf(buf.c_str(), "#%s", entry);
51                       
52                        // new entry
53                        if (strcmp(entry, "Pass") == 0)
54                        {
55                                //cout << "\nnew render stat " << endl;
56                                renderStats.push_back(RenderStats());
57                        }
58                        else if (strcmp(entry, "TotalRenderCost") == 0)
59                        {
60                                currentTag = RENDER_COST;
61                        }
62                        else if (strcmp(entry, "AvgRenderCost") == 0)
63                        {
64                                currentTag = AVG_RENDER_COST;
65                        }
66                }
67                else
68                {
69                        float val;
70                        sscanf(buf.c_str(), "%f", &val);
71                       
72                        switch (currentTag)
73                        {
74                        case RENDER_COST:
75                                //cout << "render cost: " << val << endl;
76                renderStats.back().mRenderCost = val;
77                                break;
78                        case AVG_RENDER_COST:
79                                //cout << "avg render cost: " << val << endl;
80                                renderStats.back().mAvgRenderCost = val;
81                                break;
82                        default:
83                                break;
84                        }
85                       
86                        currentTag = NO_TAG;
87                }
88        }
89
90        return true;
91}
92
93
94
95void ComputeStats(ofstream &statsOut,
96                                  const RenderStats &firstStats,
97                                  const RenderStats &currentStats,
98                                  const int index)
99{
100        float denom = firstStats.mRenderCost ? firstStats.mRenderCost : 0.000001f;
101        float costRatio = currentStats.mRenderCost / denom;
102                       
103        denom = firstStats.mAvgRenderCost ? firstStats.mAvgRenderCost : 0.000001f;
104        float avgCostRatio = currentStats.mAvgRenderCost / denom;
105
106        statsOut << "#ViewCells\n" << index + 1 << endl
107                         //<< "#TotalRenderCostGain\n" << 100.0f - costRatio * 100.0f << endl
108                         //<< "#AvgRenderCostGain\n" << 100.0f - avgCostRatio * 100.0f << endl << endl;
109                         << "#TotalRenderCostRatio\n" << costRatio << endl
110                         << "#AvgRenderCostRatio\n" << avgCostRatio << endl << endl;
111}
112
113
114inline bool vlt(const RenderStats &c1, const RenderStats &c2)
115{
116        return c1.mRenderCost > c2.mRenderCost;
117}
118
119
120// evaluate number of view cells needed for same rendercost
121void EvalNumViewCells(ofstream &outstream,
122                                          const StatsContainer &firstStats,
123                                          const StatsContainer &currentStats)
124{
125       
126        const int n = min((int)firstStats.size(), (int)currentStats.size());
127        //const int n = (int)firstStats.size();
128       
129        StatsContainer::const_iterator it, it_end = currentStats.end();
130
131        int i = 0;
132
133        // currentStats would be the other (necessarily worse?) method, firstStats would be our method
134        // loop through view cells of currentStats
135        // Compare render cost with render cost of other view cell
136        // compute ratio of view cells
137        for (it = currentStats.begin(); it != it_end; ++ it, ++ i)
138        {
139               
140                const float renderCost = (*it).mRenderCost;
141   
142                // find equivalent render cost in the other stats container.
143                // don't stop until cursor is placed one element behind
144                // or end of vector is reached
145                StatsContainer::const_iterator equalCostIt = std::upper_bound(firstStats.begin(), firstStats.end(), *it, vlt);
146
147                int j = (int)(equalCostIt - firstStats.begin());
148
149                float val;
150
151                // soecial cases
152                if (j == 0)
153                {
154                        val = (float)j;
155                }
156                else if (j >= n)
157                {
158                        val = (float)(j - 1);
159                }
160                else
161                // interpolate linearly. NOTE: probably big error because of steep curve
162                // => intepolate logarithmically
163                {
164#if 1
165                        const float rcu = log(firstStats[j - 1].mRenderCost);
166                        const float rcl = log(firstStats[j].mRenderCost);
167                        const float rc = log(renderCost);
168#else
169                        const float rcu = firstStats[j - 1].mRenderCost;
170                        const float rcl = firstStats[j].mRenderCost;
171                        const float rc = renderCost;
172#endif
173                        const float factor = (rcu - rcl != 0) ?
174                                (rcu - rc) / (rcu - rcl) : 1;
175
176                        // view cells needed for same render cost
177                        val = (float)j - 1 + factor;
178                }
179
180                // lower bound
181                //int j = max (0, i - 1);
182                //cout << "i: " << i << " j: " << j << endl;
183                float ratio = (i && val) ? (float)i / val : 1;
184
185                //cout << "ratio: " << ratio << endl;
186                outstream << "#Pass\n" << i << endl;
187                outstream << "#RenderCost\n" << renderCost << endl;
188                outstream << "#ViewCellsRatio\n" << ratio << endl << endl;     
189        }
190}
191
192
193int _tmain(int argc, _TCHAR* argv[])
194{
195        vector<StatsContainer> renderStats;
196       
197        if (argc < 3)
198        {
199                cerr << "arguments missing. Usage: input1 input2 ... intputn" << endl;
200                exit(1);
201        }
202
203        vector<string> outFilenames;
204        vector<string> outFilenames2;
205
206        // read input files from command line
207        for (int i = 1; i < argc; ++ i)
208        {
209                StatsContainer currentStats;
210     
211                ifstream file;
212                file.open(argv[i]);
213
214                // extract the render cost
215                cout << "extracting render cost of file " << argv[i] << endl;
216       
217                if (extractRenderStats(file, currentStats))
218                {
219                        renderStats.push_back(currentStats);
220
221                        // create output file name
222                        string fn = argv[i];
223               
224                        string::size_type pos = fn.find(".log", 0);
225                        fn.erase(pos, 4);
226                        //sscanf(argv[i], "%s.log", fn);
227                        string outFilename = string(fn) + string("-ratio.log");
228                        string outFilename2 = string(fn) + string("-reverse.log");
229
230                        outFilenames.push_back(outFilename);
231                        outFilenames2.push_back(outFilename2);
232                        cout << "new filen: " << outFilename2 << endl;
233                }
234                else
235                {
236                        cout << "could not open file!" << endl;
237                }
238        }
239
240        //statsOut << "rc1 size: " << (int)renderCosts1.size() << " " << (int)renderCosts2.size() << endl << endl;
241        vector<StatsContainer>::const_iterator it = renderStats.begin(), it_end = renderStats.end();
242
243        vector<string>::const_iterator sit, sit_end = outFilenames.end();
244        vector<string>::const_iterator sit2, sit2_end = outFilenames2.end();
245
246        sit = outFilenames.begin();
247        sit2 = outFilenames2.begin();
248
249        // compare all values with this vector
250        StatsContainer &firstStats = renderStats[0];
251
252        // don't compare with itself
253        //++ it;
254        //++ sit;
255
256        for (it; it != it_end; ++ it, ++ sit, ++ sit2)
257        {
258                cout << "now writing output to file " << *sit << endl;
259                ofstream statsOut;
260                statsOut.open((*sit).c_str());
261                ofstream statsOut2;
262                statsOut2.open((*sit2).c_str());
263
264                cout << "opening new file: " << (*sit2).c_str() << endl;
265                // compute size of output vector
266                const int n = min((int)firstStats.size(), (int)(*it).size());
267       
268                for (int i = 0; i < n; ++ i)
269                {
270                        //cout << "rc: " << (*it)[i].mRenderCost << endl;
271                        ComputeStats(statsOut, firstStats[i], (*it)[i], i);
272                }
273
274                // evaluate number of view cells needed for same rendercost
275                EvalNumViewCells(statsOut2, firstStats, (*it));
276
277                statsOut.close();
278                statsOut2.close();
279        }
280
281        return 0;
282}
283
Note: See TracBrowser for help on using the repository browser.