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

background{ color rgb 1.5*< 0.5, 0.5, 0.5 > }

#declare branco = 
  texture{
    pigment{ color rgb < 1, 1, 1 > }
    finish{ diffuse 0.9 ambient 0.1 }
  }
#declare preto = 
  texture{
    pigment{ color rgb < 0.5, 0.5, 0.5 > }
    finish{ diffuse 0.9 ambient 0.1 }
  }

#declare cubo =  
  box{ 
    <-1,-1,-1>,<1,1,1> texture{branco}
  }
#declare esfera =
  sphere{ 
    <0,0,0>,1 texture{branco}
  }
#declare cone_invertido =
  cone{
    <0,0,0.0>,0,<0,0,3>,1 texture{branco}
  }
#declare cone_normal=
  cone{
    <0,0,1>,1,<0,0,3>,0 texture{branco}
  }
#declare base=
  cylinder{
    <0,0,0>,<0,0,0.5>,1.5 texture{preto}
  }
#declare esfera=
  sphere{
    <0,0,0>,1 texture{branco}
  }
#declare mini_esfera=
  sphere{
    <0,0,0>,0.7 texture{branco}
  }
  
#local rinha = 0.15;
#declare esferinha=
  sphere{
    <0,0,0>,rinha texture{preto}
  }
#declare cano=
  cylinder{
    <0,0,0>,<0.5,0,0>,0.1 texture{preto}
  }
#macro tanque1()
  union{
    object{esfera} object{cano translate<1,0,0> }
  }
#end
#macro tanque2()
  union{
    object{esfera} object{cone_invertido} object{cano translate<1,0,0> } object{cano translate<-1.5,0,0>} 
  }
#end
#macro tanque4()
  union{
    object{cubo} object{cone_normal} object{cano translate<1,0,0>} object{cano translate<1,0,0> rotate 90*z} object{cano translate<1,0,0> rotate 180*z} object{cano translate<1,0,0> rotate 270*z}

  }
#end
#macro tanque3()
  union{
    object{esfera} object{mini_esfera translate<0,0,1.5>} object{cano translate<1,0,0> } object{cano translate<-1.5,0,0>} object{cano translate<0.7,0,1.5> rotate 90*z}
  }
#end
#macro interpola1(T,t0,v0,t1,v1)
  #local ss=(T-t0)/(t1-t0);
  #local vv=(1-ss)*v0+ss*v1;
  vv
#end
#macro interpola3(T,ta,tb,v0,v1,v2,v3)
  #local v01=interpola1(T,ta,v0,tb,v1);
  #local v12=interpola1(T,ta,v1,tb,v2);
  #local v23=interpola1(T,ta,v2,tb,v3);
  #local v012=interpola1(T,ta,v01,tb,v12);
  #local v123=interpola1(T,ta,v12,tb,v23);
  #local v0123=interpola1(T,ta,v012,tb,v123);
  v0123
#end
#macro testa_interpola1(obj,t0,v0,t1,v1,n)
  object{obj translate<t0,v0,0>}
  #local cont=1;
  #while(cont<n-1)
    #local position=t0+(t1-t0)*cont/(n-1);
    object{obj translate<position,interpola1(position,t0,v0,t1,v1),0>} 
    #local cont=cont+1;
  #end
  object{obj translate<t1,v1,0>}
#end
  
#include "bezpoly.inc"

#macro testa_interpola3_multi(obj,N,Pini,P1,P2,Pfim,n)
  #local t0=0;
  #local t1=N;
  object{ bezpoly_multi(N, Pini,P1,P2,Pfim, 0.9*rinha) }
  object{obj translate Pini}
  #local cont=1;
  #while(cont<n-1)
    #local position=t0+(t1-t0)*cont/(n-1);
    object{obj translate interpola3_multi(position,Pini,N,P1,P2,Pfim)} 
    #local cont=cont+1;
  #end
  object{obj translate Pfim}
#end

#macro interpola3_multi(T,Pini,n,P1,P2,Pfin)
  #local k=int(T);
  #if(k=0)
    #local Q0=Pini;
  #else
    #local Q0=(P2[k-1]+P1[k])/2;
  #end
  #local Q1=P1[k];
  #local Q2=P2[k];
  #if(k=n-1)
    #local Q3=Pfin;
  #else
    #local Q3=(P2[k]+P1[k+1])/2;
  #end
  #local Pf=interpola3(T,k,k+1,Q0,Q1,Q2,Q3);
  Pf
#end

#declare roleta=seed(34567);
#declare ponta11=<1.5,0,0>;
#declare ponta12=<-1.5,0,0>;
#declare ponta13=<0,0,3.5>;
#declare ponta21=<1.5,0,0>;
#declare ponta22=<0,1.5,0>;
#declare ponta23=<0,0,3.5>;
#declare ponta31=<1.5,0,0>;
#declare ponta32=<-1.5,0,0>;
#declare ponta33=<0,0,2.7>;

#include "eixos.inc"

#local n=3;
#local P1=array[n];
#local P2=array[n];
#local Pini=  <   0,   0,   0 >;
#local P1[0]= <   0,  10,   0 >;
#local P2[0]= <  10,  10,   0 >;

#local P1[1]= <  10,  20,   5 >;
#local P2[1]= <  20,  20,   5 >;

#local P1[2]= <  20,  10,  10 >;
#local P2[2]= <  10,  10,  10 >;

#local Pfim=  <   0,   0,   5 >;

#local n2=1;
#local P12=array[n2];
#local P22=array[n2];
#local Pini2=   <  10,  7, 10 >;
#local P12[0]=  <  11,  0, 20 >;
#local P22[0]=  <  12, 10, 20 >;
#local Pfim2=   <  13,  3, 10 >;

#include "eixos.inc"
// object{ eixos(15) }

testa_interpola3_multi(esferinha,n,Pini,P1,P2,Pfim,1000)
testa_interpola3_multi(esferinha,n2,Pini2,P12,P22,Pfim2,1000) 

#declare cmin = <  -1, -1, -1 >;
#declare cmax = < +20,+20,+22 >;

box{ cmin, <cmax.x,cmax.y,cmin.z+0.01>  texture{ pigment{ color rgb <0.80,0.70,0.60> } finish {diffuse 0.7 ambient 0.3 } } }

#include "camlight.inc"
#declare centro_cena = (cmin + cmax)/2;
#declare raio_cena = 0.55*vlength(cmin - cmax);
#declare dir_camera = < 5, 3, 4>; // <1,0,0>; //
#declare dist_camera = 5*raio_cena;
#declare intens_luz = 1.20;
camlight(centro_cena, raio_cena, dir_camera, dist_camera , z, intens_luz)