//gmesh18.pov
//povray 3.6 scene file originated by Greg M. Johnson

#include"transforms.inc"    
#declare RRR=seed(123);

//in your presets
#declare numsects=3;//number of defined sections as splines by you 
#declare startpts=9;//number of points defined for convenience at start

//in your constructed object
#declare numrings=58; // number of rings to be shaped along length.
#declare npts=26;//number of points evaluated around ring
#declare section=array[numsects+1]
#declare ringo=array[numrings+1]

#declare toppt=<-0,0,0>;
#declare botpt=<0,0,0> ;

#declare path=spline{ //linear_spline
natural_spline
    0,<1,0,1>,0.05,<1,0,0>,0.2,<1,4,0>,0.4,<2,8,0>,0.5,<0,9,0>,0.6,<-2,8,0>,0.8,<-1,4,0>,0.95,<-1,0,0>,1.00,<-1,0,1>}
          
    //0,<4,0,0>,0.25,<0,4,0>,0.5,<-4,0,0>,0.75,<0,-4,0>,1.0,<4,0,0>}
    //0,<0,0,4>,0.25,<0,4,0>,0.5,<0,0,-4>,0.75,<0,-4,0>,1.0,<0,0,4>}  
   
   /*#declare  n=0;
   #while(n<11)
   n/10, 3*(vnormalize(<0.5-rand(RRR),0,0.5-rand(RRR)>)+vnormalize(<0.5-rand(RRR),0.5-rand(RRR),0>)+vnormalize(<0,0.5-rand(RRR),0.5-rand(RRR)>))
   #declare n=n+1;
   #end           
   }
   */
   
#declare sectionplaces=array[numsects+2]{0,0.25,0.5,0.75,1.0}  
//#declare sectionplaces=array[numsects+2]{0,0.25,0.5,0.75,1.0}


#declare section[1]=spline{ cubic_spline
    #declare n=-1;
    #while (n<startpts+2)
    n/startpts, 0.7*<sin(n/startpts*2*pi),0,cos(n/startpts*2*pi)>
    #declare n=n+1;
    #end
     }
#declare section[2]=spline{ cubic_spline
    #declare n=-1;
    #while (n<startpts+2)
    n/startpts, 0.7*<sin(n/startpts*2*pi),0,cos(n/startpts*2*pi)>
    #declare n=n+1;
    #end
}
#declare section[3]=spline{ cubic_spline
    #declare n=-1;
    #while (n<npts+2)
    n/startpts, 0.7*<sin(n/startpts*2*pi),0,cos(n/startpts*2*pi)>
    #declare n=n+1;
    #end
}


/*#declare section[1]=spline{linear_spline   
        0,<0,0,2>,0.125,<0.5,0,0.5>, 0.25,<0.5,0,0>,0.5,<0,0,-0.5>,0.75,<-0.85,0,0>,0.875,<-0.25,0,0.25>,1.0,<0,0,2>}

#declare section[2]=spline{linear_spline   
        0,<0,0,2>,0.125,<0.5,0,0.5>, 0.25,<0.5,0,0>,0.5,<0,0,-0.5>,0.75,<-0.85,0,0>,0.875,<-0.25,0,0.25>,1.0,<0,0,2>}
                                                                                            
#declare section[3]=spline{linear_spline   
        0,<0,0,2>,0.125,<0.5,0,0.5>, 0.25,<0.5,0,0>,0.5,<0,0,-0.5>,0.75,<-0.85,0,0>,0.875,<-0.25,0,0.25>,1.0,<0,0,2>}
  */

#declare newradius=array[npts+1]


#declare n=0;
#while(n<npts)
    #declare newradius[n]=spline{
        natural_spline
        0,toppt
    #declare m=1;
    #while(m<numsects+1)
        
       
        sectionplaces[m],section[m](n/npts)
    #declare m=m+1;
    #end
}
#declare n=n+1;
#end              

#declare olda0=y;

#declare n=0;
#while(n<numrings)
    #declare ringo[n]=spline{
        linear_spline
    #declare o=0;
    #while(o<npts)                      
    
        #declare a01= path(n/numrings+0.001)-path(n/numrings-0.01);
        #declare a02= newradius[o](n/numrings);    
        
        //#declare a03= vtransform(a02,Reorient_Trans(a01.y,x+0.00000001));
        #declare a04=               vtransform(a02,Reorient_Trans(y,a01)); 
        
        #declare v01=18/pi*atan2(a01.x,a01.y+0.00001);//z
        #declare v02=18/pi*atan2(-a01.z,a01.y+0.00001);//
                                                    
        
        //#debug str (a01.y,8,4) #debug str(v01,8,3) #debug "\n" 
                                                              
                                                                               
        //#declare a04=   vaxis_rotate(a02,-z,v01);                               
        //#declare a04=   vaxis_rotate(a04,-x,v02);
         
        //o/npts,  vtransform(a02,Reorient_Trans(y,a01))+path(n/numrings)                    // all -y messed up.
        
        o/npts,  a04+path(n/numrings) //lateset   //pinch off at +/- transform.
        
        //o/npts,  vtransform(a02,Point_At_Trans(a01))+path(n/numrings)     //frequently messed up                          
        
    #declare o=o+1;
    #end
}
#declare n=n+1;
#end





  
  
  
//#declare botpt=vtransform(<0,-6,0>,chin_transform);


#declare themesh=
mesh2{
    vertex_vectors{
        npts*numrings+2
        toppt,
        #declare nspacer=0.5;
     
        #declare nring=0;
        #while(nring<numrings)
            #declare nspacer=mod(nspacer+0.5,1);
            #declare n=0;
            #while(n<npts)
                ringo[nring](mod((nspacer+n)/npts,1)),

            //#debug str(ring[ringo](mod(nspacer+n,1)/npts).x,6,3)
            //#debug str(ring[ringo](mod(nspacer+n,1)/npts).y,6,3)
            //#debug str(ring[ringo](mod(nspacer+n,1)/npts).z,6,3)
            //#debug "\n"

            #declare n=n+1;
            #end
        #declare nring=nring+1;
        #end
        botpt
        }

    face_indices{
        npts*2+(numrings-1)*2*npts
        #declare nspacer=0.5;
        #declare n=0;
        #while(n<npts)
            <0,n+1,1+mod(n+1,npts)>,
        #declare n=n+1;
        #end          
      
        #declare a=1;
        #declare anen=0;
        #while(anen<numrings-1)
            #if(nspacer>0.25)
                #declare np=0;
                #while(np<npts-1)
                    <a+np,a+np+1,a+np+npts>,
                  <a+1+np,a+np+npts,a+np+1+npts>,
                #declare np=np+1;
                #end
                      <a,a+npts-1,a+npts+npts-1>,
                       //<a,a+npts,a+npts+npts>,   //the problem??
                       <a,a+npts,a+npts+npts-1>,   //the problem??
            #else
                //#//if(nspacer=0)
                #declare np=0;
                #while(np<npts-1)
                   <a+np,a+np+1,a+np+npts+1>,
                  <a+np,a+np+npts,a+np+npts+1>,
                #declare np=np+1;
                #end                                   
                   //<a,npts-1,npts>,
                    //<1,2,0>//<npts-1,npts+npts-2,npts>,
                  // <a+17,a+8,a+9>,        //9: <a+17,a+8,a+9>
                   <a+npts+npts-1,a+npts-1,a+npts>,        //9: <a+17,a+8,a+9>
                 
                   //<a+2,a+npts-1,a+npts>,  //
                   <a,a+npts-1,a+npts>,  //<a+1,a+9,a+8+2>,
                 
            #end
            
            
            #declare nspacer=mod(nspacer+0.5,1);
            
            
            
        #declare a=a+npts;
        #declare anen=anen+1;
        #end

        #declare n=0;
        #while(n<npts)
            <numrings*npts+1,  numrings*npts+1-npts+n,  numrings*npts+1-npts+mod(n+1,npts)>,
        #declare n=n+1;
        #end
    }
       pigment {rgb 1}
}

object{themesh rotate 160*y translate -8*x}
object{themesh rotate -90*x translate 13*y+8*z}
object{themesh rotate -104*y translate 11*x+.5*z}

light_source{<0,200,-50> color rgb 1}
light_source{<0,200,0> color rgb x rotate -80*x rotate 45*y}
light_source{<0,200,0> color rgb z rotate -80*x rotate -45*y} 
light_source{<0,0,0> color rgb 1 translate 10*x }
background{rgb 1}
camera{location <2,8,-22> look_at <0,4,0>}