// Last edited on 2021-01-12 22:57:07 by jstolfi
//  Argumentos da otimização:

//  Deve carregar antes "direcao_do_sol.inc", "dimensoes.inc", e "parametros.inc".

#declare edir = array[nx][ny][nz]  // Normais dos espelhos.

//Direcao do sol ao meio-dia no equinócio:
#local dir_sol = direcao_do_sol(lat,0,12.00);

#local ix = 0;
#while (ix < nx)
  #local iy = 0;
  #while (iy < ny)
    #local iz = 0;
    #while (iz < nz)
    
      #local CE = epos[ix][iy][iz];
      #local dir_cal = vnormalize(CC - CE);
      
      // Normal do espelho = reflete sol na direção da caldeira:
      #local DE = vnormalize(dir_sol + dir_cal);
      #declare edir[ix][iy][iz] = DE;

      #local iz = iz + 1;
    #end
    #local iy = iy + 1;
  #end
  #local ix = ix + 1;
#end