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
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
|
#include <blitz/tinyvec.h>
#include <blitz/tinymat.h>
#include <blitz/timer.h>
using namespace blitz;
template<class T>
void optimizationSink(T&);
int main()
{
TinyMatrix<float,3,3> A1, A2;
TinyVector<float,3> b1, b2, c1, c2;
Timer timer;
const int iters = 9000000;
for (int i=0; i < 3; ++i)
{
for (int j=0; j < 3; ++j)
A1(i,j) = 1.0;
b1(i) = 1.0;
b2(i) = 1.0;
}
optimizationSink(A1);
optimizationSink(A2);
optimizationSink(b1);
optimizationSink(b2);
timer.start();
for (int i=0; i < iters; ++i)
{
c1 = product(A1,b1);
c2 = product(A1,b2);
b1 = product(A1,c1);
c2 = product(A1,c2);
c1 = product(A1,b1);
c2 = product(A1,b2);
b1 = product(A1,c1);
c2 = product(A1,c2);
c1 = product(A1,b1);
c2 = product(A1,b2);
b1 = product(A1,c1);
c2 = product(A1,c2);
c1 = product(A1,b1);
c2 = product(A1,b2);
b1 = product(A1,c1);
c2 = product(A1,c2);
c1 = product(A1,b1);
c2 = product(A1,b2);
b1 = product(A1,c1);
c2 = product(A1,c2);
c1 = product(A1,b1);
c2 = product(A1,b2);
b1 = product(A1,c1);
c2 = product(A1,c2);
c1 = product(A1,b1);
c2 = product(A1,b2);
b1 = product(A1,c1);
c2 = product(A1,c2);
c1 = product(A1,b1);
c2 = product(A1,b2);
b1 = product(A1,c1);
c2 = product(A1,c2);
c1 = product(A1,b1);
c2 = product(A1,b2);
b1 = product(A1,c1);
c2 = product(A1,c2);
c1 = product(A1,b1);
c2 = product(A1,b2);
b1 = product(A1,c1);
c2 = product(A1,c2);
c1 = product(A1,b1);
c2 = product(A1,b2);
b1 = product(A1,c1);
c2 = product(A1,c2);
c1 = product(A1,b1);
c2 = product(A1,b2);
b1 = product(A1,c1);
c2 = product(A1,c2);
c1 = product(A1,b1);
c2 = product(A1,b2);
b1 = product(A1,c1);
c2 = product(A1,c2);
c1 = product(A1,b1);
c2 = product(A1,b2);
b1 = product(A1,c1);
c2 = product(A1,c2);
c1 = product(A1,b1);
c2 = product(A1,b2);
b1 = product(A1,c1);
c2 = product(A1,c2);
c1 = product(A1,b1);
c2 = product(A1,b2);
b1 = product(A1,c1);
c2 = product(A1,c2);
}
timer.stop();
double ops = 64.0 * iters;
double flops = ops * 15;
double seconds = timer.elapsedSeconds();
double timePerOp = seconds / ops;
double cycles = timePerOp * (100.0 * 1e+6);
cout << "ops = " << ops << endl
<< "seconds = " << seconds << endl;
cout << "timePerOp = " << timePerOp << endl
<< "cycles = " << cycles << endl;
double Mflops = flops / seconds / 1.0e+6;
cout << "Mflops = " << Mflops << endl;
optimizationSink(c1);
optimizationSink(c2);
return 0;
}
template<class T>
void optimizationSink(T&)
{
}
|
