//boid04 was ideal in smooth collision avoidance swarming.
//bodi05 attempts to make them go somewhere 
//boid08b attempts to make them be "near center, not at center, with x*(1/x)          

//trail03 made a growing mess of string.  trail 04 will attempt to make something that looks like worms.

// by Greg M. Johnson   http://members.xoom.com/gregjohn/animation.html.

#include "colors.inc"
#include "textures.inc"
#include "stones.inc"
#include "metals.inc" 


#macro Reorient(Axis1,Axis2)
                #local vX1=vnormalize(Axis1);                   //john vansickle
                #local vX2=vnormalize(Axis2);
                #local vY=vnormalize(vcross(vX1,vX2));
                #local vZ1=vnormalize(vcross(vX1,vY));
                #local vZ2=vnormalize(vcross(vX2,vY));
                matrix < vX1.x, vY.x,vZ1.x, vX1.y,vY.y,vZ1.y, vX1.z,vY.z, vZ1.z, 0,0,0 >
                matrix < vX2.x,vX2.y,vX2.z,  vY.x,vY.y, vY.z, vZ2.x,vZ2.y,vZ2.z, 0,0,0 >
              #end

background {White}



#declare num=100 ;                      //number of actors
#declare framenum=375;
               
#declare RRR=seed(2022) ;        
#declare UMAX=.007 ;                    //maximum initial postion       
#declare VMAX=.07 ;                     //maximum initial velocity

#declare colorp= array[num+1]          

#declare actorp= array[num+1]           
#declare actorv= array[num+1]
#declare actoravoid= array[num+1]
                             

#fopen MyFile2 "c:\pov31\include\trailchase04.inc" write

#write(MyFile2,"#declare position=array[", framenum, "][", num, "][6] {\n")
 

                             
 
//initialize    
#declare n=0    ;
#while (n<num)

#declare colorp[n]=<rand(RRR),rand(RRR),rand(RRR)>;                                 //color assignment  of each actor


#declare actorp[n]=UMAX*<rand(RRR)-0.5,rand(RRR)-0.5,rand(RRR)-0.5>;               //this is position of each actor.     
#declare actorv[n]=VMAX*vnormalize(<rand(RRR)-0.5,rand(RRR)-0.5,rand(RRR)-0.5>);   // this is velocity of each actor.

//#declare tee=tee+1;
//#end

#declare actoravoid[n]=<0,0,0>;                                                    // this is the vector to be avoided. 
#declare n=n+1   ;
#end                                    
                            
//#declare actave=<0,0,0>; 


#declare tt=1;//



#while(tt<framenum)//*clock+17)
         

#write(MyFile2,"{ ")



//#debug "tt", str(tt,2,2)                     

#declare n=0    ;
#while (n<num)
#declare actoravoid[n]=<0,0,0>;
#declare n=n+1   ;
#end                                      
                    

                                 
#declare i=0    ;
#while (i<num)


#declare j=0    ;
#while (j<num)

#if (i!=j)
#declare actoravoid[j]=actoravoid[j]+vnormalize(actorp[i]-actorp[j])/pow(vlength(actorp[j]-actorp[i]),3);

#else
#end 
#declare j=j+1;
#end  

//#declare actave=actave+actorp[i];
#declare i=i+1;
#end     


//#declare actave=<0,0,0>//2*cos(tt/375*3*pi),0,2*sin(tt/375*3*pi)>;
       
       

#declare n=0    ;
#while (n<num)
  
//#debug str(n,2,2)
//#debug "\n"
  
#declare chasee=mod(n+1,num);//50)+int(n/50.5)*50;                      actor N chases actor N+1

#declare actave=actorp[chasee];                                 //actave is the target that each actor will seek!
                          
#declare actorp[n]=actorp[n]+.2*actorv[n];     //was 0.4*
#declare actorv[n]=vnormalize(.8*actorv[n]+3.9*(-actoravoid[n])+0.05*vnormalize(-actorp[n]+actave)*(vlength(-actorp[n]+actave)));//+vlength(1/(-actorp[n]+actave);
//                                         6.9*(-actoravoid    +0.02  was here!                               
                                         
//other users: Be careful how you mess with this equation!


//#write(MyFile,"{ ")                          
                          
//#write (MyFile2,tt,",",n,",",actorp[n][0].x,",",actorp[n][0].y,",",actorp[n][0].z,",",actorv[n][0].x,",",actorv[n][0].y,",",actorv[n][0].z,",","\n")    
//#write (MyFile2,"{ ",actorp[n][0],actorv[n][0],",","} ","\n")    



#write (MyFile2,"{",actorp[n].x,",",actorp[n].y,",",actorp[n].z,",",actorv[n].x,",",actorv[n].y,",",actorv[n].z,"}",",","\n")                          

// 350 clusters of 34 rows of 6 variables  
  
                      
/*                      
 #if (iI>0) #write(MyFile,",") #end
  #write(MyFile,"{ ")
  #local iJ=0; #while(iJ<COLUMNS)
    #if (iJ>0) #write(MyFile,",") #end
    #write(MyFile,MyDataArray[iI][iJ])
  #local iJ=iJ+1; #end
  #write(MyFile, " }\n" )
#local iI=iI+1; #end
#write(MyFile, "}\n" ) 
 */
 
//#debug str(n,2,2)
 
#declare n=n+1   ;
#end                                    

#write(MyFile2,"}, ")                          


//#declare actorv[0]=vnormalize(<sin(tt/20),0,cos(tt/20)>);
  
  
  


#declare tt=tt+1   ;



#debug "            "
#debug str(tt,2,2) 
#debug "\n"

#end                              

#write(MyFile2,"} ")                          

#fclose Myfile2 

 
 
/*


//display      
       
#declare n=0    ;
#while (n<num)
        
#declare tee=0;        
#while (tee<16)        
 //cone {<0,0,0>,.6,<2,0,0>,0
 //sphere {<0,0,0>,.25         
 cylinder {<0,0,0>,<0.25,0,0>,0.2//5
     Reorient(x,actorv[n][tee])
     translate actorp[n][tee]
     //pigment{Coral}
     pigment {color colorp[n]} 
     finish{ambient 0.3}
     //pigment {rgb 1/pow(tt,0.15)}
     }
     
#declare tee=tee+1;
#end     
      
#declare n=n+1   ;
#end          
                     

*/
                            
camera {               
  //location <0, 100, 0>   
  location <0, -5,-35>*(1+tt/350/6)
  look_at  <0, 0, 0>
  angle 30//5         
  rotate <0,60*tt/350,0>
//  Reorient(-z,actorv[0])
 // translate actorp[0]
}                               

light_source{<10,200,-50>rgb 1.2}
light_source{<10,-200,50>rgb .2}
    
fog { fog_type 1 distance 2000 rgb <.2,.3,.7>}