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