Welcome to Fractal Forums

  //combine with tubular
void polygonator (inout vec3 z) {   
	float rCyz=(atan(z.z,z.y));
	float i=1.;
	while (rCyz>pi/sides && i<sides) {rCyz-=2.*pi/sides; i++;}   //define pi and sides outside.
	 while (rCyz<-pi/sides && i<sides) {rCyz+=2.*pi/sides;i++;}  // sides = 3 for.. triangle.. :D
	z.yz*=cos(rCyz);  
 // if we really want to do this right, we need to pick a part of the fractal with reflective symmetry
 // to "divide" from... the above while loops just split it at pi (180 degrees).. which is basically fine
}
void Tubular (inout vec3 z) {   //tubular actually transforms a square to a circle,  but since I did it on 1 axis of a 3d system... 
float rCyz= (z.y*z.y)/(z.z*z.z);
if (rCyz<1.) {rCyz=sqrt(rCyz+1.);} else {rCyz=sqrt(1./rCyz+1.);}
z.yz*=rCyz;  
}

Description:
Warps a cube to a sphere; transform made by M.Benesi, optimized by
Luca.


// spheric code
Cvector3 temp = z;
z *= z;
double rCyz = z.y / z.z;

double rCxy z= (z.y + z.z ) / z.x ;


if (rCyz < 1.0) rCyz = sqrt(rCyz + 1.0);
else rCyz=sqrt( 1.0 / rCyz + 1.0);

if (rCxyz < 1.0) rCxyz =s qrt(rCxyz + 1.0);
else rCxyz=sqrt(1.0/ rCxyz + 1.0);

z.yz *=rCyz; 


z *= rCxyz/sr32;

void fullpolytest (inout vec3 z) {  
	

	float  rCyz=(atan(z.z,z.y));
	float i=1.;
	if (samesides) {   //lock number of sides between this and the other section below- so both use sides2 variable....
		while (rCyz>pi/sides2 && i<sides2) {rCyz-=pi2/sides2; i++;}
	 	while (rCyz<-pi/sides2 && i<sides2) {rCyz+=pi2/sides2;i++;}
	} else {
		while (rCyz>pi/sides && i<sides) {rCyz-=pi2/sides; i++;}
	 	while (rCyz<-pi/sides && i<sides) {rCyz+=pi2/sides;i++;}
	}
	z.yz*=cos(rCyz);
	i=1.;	
	if (sides2==3) {		// for the tetrahedron... it was skewed, need to make it regular  
		z.x/=sqrt(2);     // square root of 2 appears correct... and it makes a bit of sense, but I'd like confirmation 
				// of the math... 
	}
	float  rCxyz= atan(sqrt(z.y*z.y+z.z*z.z),z.x);
	i=1.;
	
	while (rCxyz>pi/sides2 && i<sides2) {rCxyz-=pi2/sides2; i++;}
	while (rCxyz<-pi/sides2 && i<sides2) {rCxyz+=pi2/sides2;i++;}
	z*=cos(rCxyz);
//  this is the spheric part... need it because the above transform is for sphere to polyhedra
	rCyz= (z.y*z.y)/(z.z*z.z);
	rCxyz= (z.y*z.y+z.z*z.z)/(z.x*z.x);
	if (rCyz<1.) {rCyz=sqrt(rCyz+1.);} else {rCyz=sqrt(1./rCyz+1.);}
	if (rCxyz<1.) {rCxyz=sqrt(rCxyz+1.);} else {rCxyz=sqrt(1./rCxyz+1.);}

	z.yz*=rCyz;  
	z*=rCxyz;  
}