1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
|
#include <blitz/timer.h>
#include <iostream>
#include <fstream>
using namespace blitz;
using namespace std;
extern "C" {
void echo_f77tuned(int& N, int& niters, float& check, int& blockSize);
}
int main()
{
int N = 1024;
int niters = 48;
float check;
double Mflops = niters * 9;
Timer timer;
ofstream ofs("echotune.log");
cout << "This program decides on the best block size for a typical 2D "
<< endl << "stencil operation. Pick the block size which has the "
<< endl << "maximum Mflops/s." << endl << endl;
cout << "Block size\tMflops/s" << endl;
int blockSize;
for (blockSize=1; blockSize < 32; ++blockSize)
{
timer.start();
echo_f77tuned(N, niters, check, blockSize);
timer.stop();
cout << blockSize << "\t" << (Mflops/timer.elapsedSeconds()) << endl;
ofs << blockSize << "\t" << (Mflops/timer.elapsedSeconds()) << endl;
}
for (; blockSize < 1024; blockSize += 32)
{
timer.start();
echo_f77tuned(N, niters, check, blockSize);
timer.stop();
cout << blockSize << "\t" << (Mflops/timer.elapsedSeconds()) << endl;
ofs << blockSize << "\t" << (Mflops/timer.elapsedSeconds()) << endl;
}
return 0;
}
|
