4{
5 double y[3], x[3],
n[3];
6
7 y[0] = 0;
8 y[1] = 1;
9 y[2] = 2;
10
11 x[0] = 2;
12 x[1] = 3;
13 x[2] = 4;
14
18
19 for(
int i = 0;
i < 100000000;
i++)
20 {
21
22
23
24
25 n[0] += (y[0] - x[0]) - (y[0] - x[0]) + (y[0] - x[0]) - (y[0] - x[0])
26 + (y[0] - x[0]) - (y[0] - x[0]) + (y[0] - x[0]) - (y[0] - x[0])
27 + (y[0] - x[0]) - (y[0] - x[0]) + (y[0] - x[0]) - (y[0] - x[0])
28 + (y[0] - x[0]) - (y[0] - x[0]) + (y[0] - x[0]) - (y[0] - x[0]);
29 n[1] += (y[1] - x[1]) - (y[1] - x[1]) + (y[1] - x[1]) - (y[1] - x[1])
30 + (y[1] - x[1]) - (y[1] - x[1]) + (y[1] - x[1]) - (y[1] - x[1])
31 + (y[1] - x[1]) - (y[1] - x[1]) + (y[1] - x[1]) - (y[1] - x[1])
32 + (y[1] - x[1]) - (y[1] - x[1]) + (y[1] - x[1]) - (y[1] - x[1]);
33 n[2] += (y[2] - x[2]) - (y[2] - x[2]) + (y[2] - x[2]) - (y[2] - x[2])
34 + (y[2] - x[2]) - (y[2] - x[2]) + (y[2] - x[2]) - (y[2] - x[2])
35 + (y[2] - x[2]) - (y[2] - x[2]) + (y[2] - x[2]) - (y[2] - x[2])
36 + (y[2] - x[2]) - (y[2] - x[2]) + (y[2] - x[2]) - (y[2] - x[2]);
37
41 }
42 std::cout << y[0] << " " << y[1] << " " << y[2] << std::endl;
43}
FTensor::Index< 'i', SPACE_DIM > i
const double n
refractive index of diffusive medium