// Exemplo de arquivo de descricao de cena para POV-ray
// Last edited on 2017-05-07 20:10:09 by stolfilocal

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

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

#declare tx_plastico = 
  texture{
    pigment{ color rgb < 0.10, 0.80, 1.00 > }
    finish{ diffuse 0.8 ambient 0.1 specular 0.5 roughness 0.005 }
  }

#declare tx_fosca = 
  texture{
    pigment{ color rgb < 1.00, 0.80, 0.10 > }
    finish{ diffuse 0.9 ambient 0.1 }
  }

#declare tx_espelho = 
  texture{
    pigment{ color rgb < 1.00, 0.85, 0.30 > }
    finish{ diffuse 0.2 reflection 0.7*< 1.00, 0.85, 0.30 > ambient 0.1 }
  }

#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 }
  }


#declare tx_xadrez =
  texture{
    pigment{ checker color rgb < 0.10, 0.32, 0.60 >, color rgb < 1.00, 0.97, 0.90 > }
    finish{ diffuse 0.9 ambient 0.1 }
    scale 2.0
  }

// ======================================================================
// DESCRIÇÃO DA CENA 

#macro interpolar(t0, v0, t1, v1, tt)
	#local r = (tt - t0)/(t1 - t0);
	#local s = 1 - r;
	#local vv = s*v0 + r*v1;
	vv
#end

#macro arcobezier(t0, t1, A, B, C, D, tt)
	#local AB = interpolar(t0, A, t1, B, tt);
	#local BC = interpolar(t0, B, t1, C, tt);
	#local CD = interpolar(t0, C, t1, D, tt);

	#local ABC = interpolar(t0, AB, t1, BC, tt);
	#local BCD = interpolar(t0, BC, t1, CD, tt);

	#local ABCD = interpolar(t0, ABC, t1, BCD, tt);

	ABCD
#end

#macro suavizaZ(P, n, i, vv)
	#local j = mod(i+1, n);

	#declare P[i][2] = P[i][3] - vv;
	#declare P[j][1] = P[j][0] + vv;
#end

#macro suavizaX(P, n, i)
	#local j = mod(i+1, n);
	#local vv = (P[j][1] - P[i][2])/2;
	suavizaZ(P, n, i, vv)
#end

#macro noh(P, M)
	#local tt = 0.0;
	union {
		#local j = 0;
		#local i = 0;
		#while(i < M)
			#if ((tt >= j/4) & (tt < (j+1)/4))
				#declare centro = arcobezier(j/4, (j+1)/4, P[j][0], P[j][1], P[j][2], P[j][3], tt);

				sphere{ centro , 0.1 texture { tx_fosca }}
			#else
				#local j = j + 1;
			#end

			#local tt = tt + 1/n;
			#local i = i + 1;
		#end
	}
#end

#macro interpolaMatriz(P, Q, c0, c1, cc, R)
	#local i = 0;
	#local j = 0;
	#while (i < 4)
		#while (j < 4)
			#declare R[i][j] = interpolar(c0, P[i][j], c1, Q[i][j], cc);
			#local j = j+1;
		#end
		#local j = 0;
		#local i = i+1;
	#end
#end


#macro desenho(centro, n, cenas)
	#local i = 0;

	#declare P = array[4][4];
	#declare P1 = array[4][4];
	#declare R = array[4][4];

	#for (i, 0, 3, 1)
		#for(j, 0, 3, 1)
			#declare R[i][j] = <0,0,0>;
		#end
	#end

	// matriz P2
	#declare P[0][0] = <0,0,0>;
	#declare P[0][1] = <0.2,0,0.1>;
	#declare P[0][2] = <0.1,0,0.1>;
	#declare P[0][3] = <0.5,0,0.25>;

	#declare P[1][0] = P[0][3];
	#declare P[1][1] = <-0.5,0,-0.5>;
	#declare P[1][2] = <-0.3,0,0.5>;
	#declare P[1][3] = <0.5,0,0.4>;

	#declare P[2][0] = P[1][3];
	#declare P[2][1] = <0,0.2,-0.2>;
	#declare P[2][2] = <0,0.2,-0.2>;
	#declare P[2][3] = <0.5,-0.2,-0.4>;

	#declare P[3][0] = P[2][3];
	#declare P[3][1] = <0.5,0.2,0.5>;
	#declare P[3][2] = <0,0.2, -0.2>;
	#declare P[3][3] = P[0][0];

	// matriz P1
	#declare P1[0][0] = <0,0,0>;
	#declare P1[0][1] = <0.2,0,0.1>;
	#declare P1[0][2] = <0.1,0,0.1>;
	#declare P1[0][3] = <0.5,0,0.25>;

	#declare P1[1][0] = P1[0][3];
	#declare P1[1][1] = <-0.8,0,-0.8>;
	#declare P1[1][2] = <-0.7,0,0.7>;
	#declare P1[1][3] = <1.0,0,0.4>;

	#declare P1[2][0] = P1[1][3];
	#declare P1[2][1] = <0,0.2,-0.2>;
	#declare P1[2][2] = <0,0.2,-0.2>;
	#declare P1[2][3] = <0.5,-0.2,-0.4>;

	#declare P1[3][0] = P1[2][3];
	#declare P1[3][1] = <0.5,0.2,0.5>;
	#declare P1[3][2] = <0,0.2,-0.2>;
	#declare P1[3][3] = P1[0][0];

	#while(i < 4)
		suavizaX(P, 4, i)
		suavizaX(P1, 4, i)

		#local i = i + 1;
	#end

	#local i = 0;
	#while (i < cenas)
		interpolaMatriz(P, P1, 0, cenas, i, R)
		noh(R, n)
		#local i = i + 1;
	#end
	//noh(P, n)
#end

desenho(0, 1800, 4)

#include "camlight.inc"
#declare centro_cena = < 0, 0, 0 >;
#declare raio_cena = 2;
#declare dir_camera = < 10.00, 30.00, 5.00 >;
#declare dist_camera = raio_cena;
#declare intens_luz = 1.20;
camlight(centro_cena, raio_cena, dir_camera, dist_camera , z, intens_luz)