// Dragons sans utiliser la récursivité

#include "colors.inc"
#include "metals.inc" 

//-------------- n ITERATIONS---------------


camera  { location <-0.5,-0,-2> look_at<-0.5,-0.2,0>}
light_source { <-3,3,-100> color White } 
light_source { <2,3,-50> color Yellow} 

background { rgb<1,0.95,0.85>*1 } //<0.,0.5,0.95>*1.0

//---------------- macros -------------------
#macro sim(a,b,e,f,g) // pour similitude appliquée à <u,v>
   <a*e-b*f,a*f+b*e,g>
  #end

//----------- MACRO DRAGON------------------------
#macro dragon(s1x,s1y,s2x,s2y,transx,transy,n)
   
   #local tab = array [pow(2,n)]  ; 
   #declare tab[0]=<0,0,0>;
   #declare tab[1]=<1,0,0>; 
   #local trx=transx;
   #local try=transy;
   #local be=1;

   #while(be<n)                   // boucle (1) be=1..n-1
      #declare centre= tab[pow(2,be)-1];
      #declare LL=vlength(tab[pow(2,be)-1]-tab[0]);  // on ramène ainsi le diamètre à environ 2
      #declare bi=0;                                 // donc caméra, lumières fixes pour tout n
      #while (bi<pow(2,be))              // boucle (2) bi=0..2^be-1
          #declare temp=tab[bi]-centre; 
          // similitude s2 puis translation trans2 à la branche d'ordre be
          #declare tab[pow(2,be+1)-bi-1]=(centre +sim(s2x,s2y,temp.x,temp.y,temp.z))/LL;// 
          #declare bi=bi+1;
      #end                              // fin boucle (2)
   
      #declare bi=0;
      #while (bi<pow(2,be))               // fin boucle (3) bi=0..2^be-1
          // similitude s1 à la branche  d'ordre be
          #declare temp=tab[bi];
          #declare tab[bi]=(sim(s1x,s1y,temp.x,temp.y,temp.z)+<trx,try,0>)/LL;//
          #declare bi=bi+1;
      #end 
      #declare trx=trx/LL;
      #declare try=try/LL;                               // fin boucle (3)
      #declare be=be+1;
   #end
   //-------------------------------------------   
   
   #local r=0.25/pow(1.4,n) ;
   #declare cp=0; 
   union {
        #while(cp<pow(2,n)-1 )
          #if ( vlength(tab[cp]-tab[cp+1])!= 0 )
             sphere {tab[cp] ,r}
             cylinder { tab[cp],tab[cp+1] ,r}
          #end
        #declare cp=cp+1;
        #end
        sphere {tab[pow(2,n)-1] ,r} 
        texture {T_Copper_3B}
        }

    //sphere {0 ,r*1.01  pigment {Red}  }

#end // macro DRAGON

//--------------INITIALISATIONS---------------


  // dragon classique
#declare n=10;
object{ dragon(1, 0, 0, 1,0,0,10) rotate <0,0,(n+3)*45>
        }
         
/* */ 
/* // dragon enflammé !
#declare n=10;
object{ dragon(1, 0, 0, 1,cos(pi/4),sin(pi/4),n)
        rotate <0,0,(n-4)*45-22>
        translate <0.1,0.,0>
        finish{ ambient 0.6}
        }
*/      

/* dragon sans points multiples
#declare n=12;
object{ dragon(1, 0, 0, 1,0,0.8,n)
        rotate <0,0,(n+3)*45>
        translate <0.1,0.,0>
        scale 1
        }