// Exemplo de arquivo de descricao de cena para POV-ray
// Last edited on 2010-03-04 15:44:01 by stolfi

// ======================================================================
// CORES E TEXTURAS

background{ color rgb < 0.75, 0.80, 0.85 > }

#declare tx_tronco = 
  texture{
    pigment{ color rgb < 0.65, 0.16, 0.16 > }
    finish{ diffuse 0.9 ambient 0.1 }
  }

#declare tx_chao = 
  texture{
    pigment{ color rgb < 0.15, 0.16, 0.16 > }
    finish{ diffuse 0.9 ambient 0.1 }
  }

#declare tx_branco = 
  texture{
    pigment{ color rgb < 0.8, 0.8, 0.8 > }
    finish{ diffuse 0.9 ambient 0.8 }
  }

#declare tx_vermelha = 
  texture{
    pigment{ color rgb < 0.89, 0.15, 0.13 > }
    finish{ diffuse 0.9 ambient 0.8 }
  }

#declare tx_xadrez =
  texture{
    pigment{ checker color rgb < 0.2, 0.3, 0.6 >, color rgb < 1, 1, 1 >}
    finish{ diffuse 0.9 ambient 0.1 }
    scale 2.0
  }
#declare tx_vidro = 
  texture{
    pigment{ color rgb < 0.85, 0.95, 1.00 > filter 0.70 }
    finish{ diffuse 0.03 reflection 0.25 ambient 0.02 specular 0.25 roughness 0.005 }
  }


// ======================================================================
// DESCRICAO DA CENA

// Partes da cena:

#include "eixos.inc"
#include "retalho.inc"

#declare chao = 
  box{ <-20,-20,-1>, <+20,+20,0> }

#macro interpola(f0, v0, f1, v1, f)
  #local ss = (f - f0) / (f1 - f0);
  #local rr = 1 - ss;
  (rr * v0 + ss * v1)
#end

#macro interpola4(t0, v0, v1, v2, v3, t1, tt)
	#local v01 = interpola(t0, v0, t1, v1, tt);
	#local v12 = interpola(t0, v1, t1, v2, tt);
	#local v23 = interpola(t0, v2, t1, v3, tt);

	#local v012 = interpola(t0, v01, t1, v12, tt);
	#local v123 = interpola(t0, v12, t1, v23, tt);

	#local vv = interpola(t0, v012, t1, v123, tt);
	vv
#end

#macro pingo(vv)
	union {
		object {
			sphere {
				vv
				0.2
			}
		}
	}
#end

#macro grafico(x0, y0, y1, y2, y3, x3, N)
	union {
		#local i = 0;
		#local p = (x3 - x0) / N;
		#while(i < N)
			#local xx = x0 + (p * i);
			#local vv = interpola4(x0, y0, y1, y2, y3, x3, xx);
			pingo(vv)
			#local i = i + 1;
		#end
	}
#end

#macro supergrafico(pts_por_curva)
	union {
		#local x0 = 0;
		#local x1 = 10;
		#local x2 = 30;
		#local x3 = 10;
		#local x4 = 0;
		#local y0 = <1, 3, 1>;
		#local y1 = <3, 8, 16>;
		#local y2 = <15, 2, 6>;
		#local y3 = <-4, -6, 12>;
		#local y4 = <5, 6, 4>;
		#local y5 = <-1, -4, 0>;
		#local y6 = <12, 4, 2>;
		#local y7 = <-8, -12, -18>;
		#local y8 = <-5, -8, -10>;
		#local y9 = <1, 3, 1>;
		#local b1 = (y2+y3)/2;
		#local b2 = (y4+y5)/2;
		#local b3 = (y6+y7)/2;

		grafico(x0, y0, y1, y2, b1, x1, pts_por_curva)
		grafico(x1, b1, y3, y4, b2, x2, pts_por_curva)
		grafico(x2, b2, y5, y6, b3, x3, pts_por_curva)
		grafico(x3, b3, y7, y8, y9, x4, pts_por_curva)
	}
#end

#declare p21 = <1, 8, 0>;
#declare p22 = <3, 10, 0>;
#declare p23 = <5, 11, 0>;
#declare p24 = <7, 12, 0>;
#declare p31 = <0, 4, 0>;
#declare p32 = <2, 6, 0>;
#declare p33 = <4, 5, 0>;
#declare p34 = <6, 8, 0>;
#declare p41 = <1, 0, 0>;
#declare p42 = <3, 2, 0>;
#declare p43 = <5, 4, 0>;
#declare p44 = <7, 2, 0>;

#declare q11 = <0, 24, 0>;
#declare q12 = <2, 28, 0>;
#declare q13 = <4, 26, 0>;
#declare q14 = <6, 27, 0>;
#declare q21 = <2, 18, 0>;
#declare q22 = <3, 20, 0>;
#declare q23 = <5, 21, 0>;
#declare q24 = <6, 22, 0>;
#declare q31 = <1, 14, 0>;
#declare q32 = <3, 16, 0>;
#declare q33 = <4, 15, 0>;
#declare q34 = <7, 18, 0>;

#declare p11 = (p21 + q31) / 2;
#declare q41 = p11;
#declare p12 = (p22 + q32) / 2;
#declare q42 = p12;
#declare p13 = (p23 + q33) / 2;
#declare q43 = p13;
#declare p14 = (p24 + q34) / 2;
#declare q44 = p14;
 
#macro superretalho()
	union {
		retalho(p11, p12, p13, p14, p21, p22, p23, p24, p31, p32, p33, p34, p41, p42, p43, p44, 0.1, tx_branco, tx_vermelha)
		retalho(q11, q12, q13, q14, q21, q22, q23, q24, q31, q32, q33, q34, q41, q42, q43, q44, 0.1, tx_branco, tx_vermelha)
	}
#end

union{
	//object { eixos(10.0) }
	object { supergrafico(50) translate<-12, 0, 0> texture{ tx_vidro } }
	object { superretalho() translate<8, -14, 0> }
}

#include "camlight.inc"
#declare centro_cena = < 0.00, 1.00, 0.00 >;
#declare raio_cena = 30.0;
#declare dir_camera = < 0.00, 1.0, 10.00 >;
#declare dist_camera = 30.0;
#declare intens_luz = 1.00;
camlight(centro_cena, raio_cena, dir_camera, dist_camera , z, intens_luz)