END OF AN ERA, FRACTALFORUMS.COM IS CONTINUED ON FRACTALFORUMS.ORG

tester00 {
; ----------------

init:
  z = #pixel
  float x = real(z)
  float y = imag(z)
  ;float x2 = x
  float zz = @zinit
  float rad = 0
  float rr = 0
  int i = @itercnt
loop:

  ; here is the main fmla
  ; -------------------
  while i>1
  x = x-@p
  y = y-@q
  zz = zz-@r
  rad=(x*x+y*y+zz*zz)
  x=x/rad
  y=-y/rad
  zz = (zz+@r)/rad
  i = i-1
  endwhile
  ; -------------------
  ;while i>1
  rr=sqrt(@q*@q+@r*@r)
  y=y+2*@r/rr
  zz=zz-2*@q/rr
  rad=(x*x+y*y+zz*zz)
  x=4*rr*x/rad
  y=4*rr*(y/rad)-@r
  zz=4*rr*(zz/rad)+@q
  ;i = i-1
  ;endwhile
  ; -------------------
  z = x + flip(y)
  ;ends here
bailout:
   false
   ;|z| > @bailout
default:
  title = "Test function"
  center = (0,0)
  float param p
    caption = "p"
    default = 0
  endparam
  float param q
    caption = "q"
    default = 1.9
  endparam
  float param r
    caption = "r"
    default = .04
  endparam
  float param zinit
    caption = "z init"
    default = 0
  endparam
  int param itercnt
    caption = "Iter count"
    default = 30
  endparam
  float param bailout
    caption = "Bailout"
    hint = "-disabled-"
    default = 1024
  endparam
}

bool Mandel3TriZ(double x, double y, double z){
	//space coordinates transformations..
	//Turns the 3D space in an infinite amount 
	//of 2D slices 
	
	const double w=1/(x*x+y*y+z*z);
	const double dem=1-(z)*w; 
	const double ctx =((x)*w)/dem;
	const double cty =((-y)*w)/dem;
	double Xp=0,Yp=0;
	double radius;

	//put your favorite fractal iteration here
	//A simple Mandelbrot is show..but you can 
	//try others..
	for(unsigned char Counter=0;Counter<20;Counter++){
	const double xs=Xp*Xp-Yp*Yp+ctx;
	const double ys=2*Xp*Yp+cty;
	Xp =xs;
	Yp = ys;
	radius=Xp*Xp+Yp*Yp;
	if (radius>100) break;
	const double w1=1/(xs*xs+ys*ys+z*z);
	const double dem1=1-(z)*w; 
	Xp =((xs)*w1)/dem1;
	Yp =((-ys)*w1)/dem1;

	}
	//false, the iteration took the blue pill
	//true, the iteration took the red pill
	if(radius<100)
	return false;
	else
	return true;
}

//Based on the MandelBox..4D
bool CheeseBox(double x, double y, double z){
	
  const double ctx =(x   + y   + z);
  const double cty = (-x -y + z);
  const double ctz =(-x + y   -z);
  const double ctw = (x   -y -z);
	float w=0;//cheesebox
	//float w=pow(x*x+y*y+z*z,0.333); //cybertronian box
	for(unsigned char Counter=0;Counter<15;Counter++){
	if(x>1)x=2-x;
	else 
	if(x<-1)x=-2-x;
	if(y>1)y=2-y;
	else 
	if(y<-1)y=-2-y;
	if(z>1)z=2-z;
	else 
	if(z<-1)z=-2-z;
	if(w>1)w=2-w;
	else 
	if(w<-1)w=-2-w;
const float vm=sqrt(x*x+y*y+z*z+w*w);

if(vm<0.5){
x=x*4;
y=y*4;
z=z*4;
w=w*4;
} else 
if(vm<1){
const float vsq=vm*vm;
x/=vsq;
y/=vsq;
z/=vsq;
w/=vsq;
}
x=x*2+ctx;
y=y*2+cty;
z=z*2+ctz;
w=w*2+ctw;

radius=sqrt(x*x+y*y+z*z+w*w);

	if (radius>60) break;
	
	}

	
	if(radius<60)
	return true;
	else
	return false;
}