1 #define BZ_TAU_PROFILING
2
3 #include <blitz/array.h>
4
5 using namespace blitz;
6
main()7 int main()
8 {
9 TAU_PROFILE("main()", "int ()", TAU_DEFAULT);
10
11 const int N = 32;
12
13 Array<float,2> A(N,N), B(N,N), C(N,N), D(N,N), E(N,N);
14 A = 5.0;
15 B = 0.0;
16 C = 0.0;
17
18 for (int i=0; i < 20; ++i)
19 {
20 D = A + B + C;
21 D /= sum(pow2(D));
22 A = B * cos(D) + C * sin(D);
23 B += exp(-D);
24
25 float x = sum(A);
26 float y = sum(A+B);
27 float z = sum(sqr(A)+sqr(B));
28 C = x*A+y*B+z*C;
29
30 D = exp(-sqr(A)-sqr(B));
31 E = A + B + C + D;
32 float q = min(A);
33 float r = max(B);
34 }
35
36 return 0;
37 }
38
39