Welcome to Fractal Forums

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
}

float SAVEMANDELITER[130];

SAVEMANDELITER[0]=0;
for (int tr=1;tr<129;tr++) SAVEMANDELITER[tr]=pow(SAVEMANDELITER[tr-1],0.85)+0.1475;

bool Newton(float x, float y, float z){
	float a,b,c,d,e,f,d1,g1,e1,potential;
	float w,X1,Y1,Z1,xiter,yiter,W1,radius,radius1,angle,angle1, inside;
	char u_Counter;
	potential=0;
	if(z<=0)return (false);
	W1=0;
	a=x;
	b=y;
	inside=1000;
	radius=0;
	xiter=x;
	yiter=y;
	for(u_Counter=0;u_Counter<100;u_Counter++){
		w=W1;
		c=xiter;
		d=yiter;
		const float X2=xiter*xiter;
		const float Y2=yiter*yiter;
		const float X3=X2*xiter;
		const float Y3=Y2*yiter;
		const float X4=X2*X2;
		const float Y4=Y2*Y2;
		const float X5=X4 - 6 *X2 *Y2 + Y4 - 1;
		const float Y5=4* X3 *yiter - 4.0*xiter *Y3;
		const float X6=4 *( X3 - 3.0*xiter *Y2 );
		const float Y6= 4 *( 3 *X2 *yiter - Y3 ); 
		const float den = 1.0f / (X6*X6 + Y6*Y6);
		X1 = xiter-(X5*X6 + Y5*Y6) * den;
		Y1 = yiter-(Y5*X6 - X5*Y6) * den;
		radius1=radius;
		radius=sqrt(X1*X1+Y1*Y1);
		if(radius<inside)inside=radius;
		xiter=X1;
		yiter=Y1;
		W1=SAVEMANDELITER[u_Counter]; //W1=powf(W1,0.85)+0.1475; //if you cannot use arrays.. runs slower..
		if (fabs(radius-radius1)<0.01) break;
	}
	X1=a;
	Z1=b;
	potential=radius1;
	if (radius<10){potential=inside*5; } 
	if(W1<z)return(false);
	if (radius<10)w=inside;
	f=((z-w)/(W1-w))*5;
	angle=10-radius1;
	angle1=radius-radius1;
	c=(c/10.0)*5;
	d=(d/10.0)*5;
	radius= 7 - (cos(c + 1.61803399*f) + cos(c - 1.61803399*f) + cos(f + 1.61803399*d) + cos(f - 1.61803399*d) + cos(d - 1.61803399*c) + cos(d + 1.61803399*c));

//if(potential<0)return false; else return true;

	if (radius<4.8)return (true);
	return(false);
}

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;
}