Welcome to Fractal Forums

// aBoxModKaliEiffie
// reference http://www.fractalforums.com/new-theories-and-research/aboxmodkali-the-2d-version/
// Created: 22/11/16
#define providesInit
#include "MathUtils.frag"
#include "de-Raytracer.frag"

// #group MixPinski
//MAIN

uniform bool swapXZ; checkbox[true]
//uniform bool bla bla; checkbox[false] // TODO later

uniform vec3 Offset1; slider[(-5.0,-5.0,-5.0),(1.0,1.0,1.0),(5,5,5)]
uniform vec3 cJ; slider[(-5.0,-5.0,-5.0),(0.0,0.0,0.0),(5,5,5)] // jula constant
Uniform vec3 c3;  slider[(-5.0,-5.0,-5.0),(1.0,1.0,1.0),(5,5,5)] // addCpixel constant multiplier
c3 = pos * c3; //  ??

uniform float minR2; slider[-5.00,0.25,5.00]
uniform float Scale1; slider[-5.00,1.50,5.00]
float mdR2 = Scale1/minR2; 

uniform int i;
	
float Dd = 1.0f; // using Dd for internal DE calc


uniform int MI; slider[0,60,250]
uniform float Bailout, slider[0,100,250]
uniform int ColorIterations; slider[0,10,250]

//Rotations
uniform int RotAStart; slider[0,1,15] 
uniform vec3 RotAVector; slider[(0,0,0),(1,1,1),(1,1,1)]
uniform float RotAAngle; slider[0.00,0,180]
uniform int RotBStart; slider[0,1,15] 
uniform vec3 RotBVector; slider[(0,0,0),(1,1,1),(1,1,1)]
uniform float RotBAngle; slider[0.00,0,180]

mat3 fracRotation1;
mat3 fracRotation2;
void init() {
	fracRotation1 = rotationMatrix3(normalize(RotAVector), RotAAngle);
	fracRotation2 = rotationMatrix3(normalize(RotBVector), RotBAngle);
}


float DE(vec3 z) {
	
	float r2=0.;

	int i=0;
	for(; i<MI; i++) {


		
z = abs(z + Offset1) - abs(z - Offset1) - z;// tglad fold

if (swapXZ)//  default ON
{
	z = vec3(z.z, z.y, z.x); // swizzle swap
}
if(i>=RotAStart) z.xyz*=fracRotation1;

if (z.y > z.x) z = vec3(z.y, z.x, z.z); // conditional swizzle

// Abox type spherical fold, NO fixed radius
// no DE adjustment as we are modifying the scale value instead		
float r2 = dot(z,z);		
float m = Scale1;
		
if (r2 < minR2) m =  mdR2;

else
    
if (r2 < 1.0f) m =  m / r2;



z =  z * m + c3 + cJ;//   z * scale + addCpixel + addJuliaConstant  

if(i>=RotBStart) z.xyz*=fracRotation2;   

//part of distance estimation calculation

Dd = Dd * abs(m) + 1.0f;

//for Bailout
r2 = dot(z,z);

if(r2>Bailout)
 		
{
	float r = sqrt(r2);
	return = r / abs(Dd);
}

#include "DE-Raytracer.frag"

uniform int MI; slider[1,5,100]
uniform int bailout; slider[1,16,1024]
uniform float Scale; slider[-5,2,5]
uniform float Min_R; slider[0,0.5,5]
uniform float X_add; slider[-5,0,5]
uniform float Y_add; slider[-5,0,5]
uniform float Z_add; slider[-5,0,5]
uniform float Julia_Cx; slider[-5,0,5]
uniform float Julia_Cy; slider[-5,0,5]
uniform float Julia_Cz; slider[-5,0,5]

float DE(vec3 p) {
	vec3 z=vec3(p);
	float m;
	float rr=0.; 
	float Dd = 1.0; 

	for(int i=0; i<MI && rr<bailout; i++){	
	z.x = X_add - abs(z.x);
	z.y = Y_add - abs(z.y);
	z.z = Z_add - abs(z.z);
	rr = z.x*z.x + z.y*z.y + z.z*z.z;
	if (rr < sqrt(Min_R)) {m = Scale/sqrt(Min_R); };
	if (rr < 1) {m = Scale/rr;} else { m = Scale; };
	z.x = z.x * m + Julia_Cx;
	z.y = z.y * m + Julia_Cy;
	z.z = z.z * m + Julia_Cz;
	Dd *= Scale; // for DE calc		
	}
float r = sqrt(rr);
return (r) / abs(Dd);
}


#preset default
FOV = 0.4
Eye = 0.2248197,-0.1026156,-1.696023
Target = -2.466143,0.075118,7.93347
Up = -0.0682873,0.9969614,-0.0374841
AutoFocus = false
Gamma = 2
ToneMapping = 4
Exposure = 1
Brightness = 1
Contrast = 1
Saturation = 1
GaussianWeight = 1
AntiAliasScale = 2
DepthToAlpha = false
Detail = -3.2932
FudgeFactor = 0.1545559
Dither = 0.5
NormalBackStep = 1
CamLight = 1,0.9647059,0.9215686,1
BaseColor = 1,0.9803922,0.9411765
OrbitStrength = 0.8399072
X = 0.5,0.6,0.6,0.4040632
Y = 1,0.6,0,0.4
Z = 0.8,0.78,1,0.5
R = 0.4,0.7,1,0.12
BackgroundColor = 0.6,0.6,0.45
GradientBackground = 0.3
CycleColors = false
Cycles = 1.1
EnableFloor = false
FloorNormal = 0,0,1
FloorHeight = 0
FloorColor = 1,1,1
EquiRectangular = false
FocalPlane = 1
Aperture = 0
ShowDepth = false
DepthMagnitude = 1
DetailAO = -1.166667
MaxDistance = 100
MaxRaySteps = 1669
AO = 0,0,0,0.928489
Specular = 0.4
SpecularExp = 16
SpecularMax = 10
SpotLight = 1,0.9607843,0.9215686,0.5639188
SpotLightDir = 1,0.1
CamLightMin = 0
Glow = 1,0.8862745,0.8078431,0.0849825
GlowMax = 222
Fog = 0.1160714
HardShadow = 0.7913247 NotLocked
ShadowSoft = 19.76471
QualityShadows = false
Reflection = 0 NotLocked
DebugSun = false NotLocked
MI = 16
Scale = 2.152391
Min_R = 0.5
X_add = 0.0223714
Y_add = 0
Z_add = 0.3131991
bailout = 1024
Julia_Cx = 0.4072398
Julia_Cy = 1.289593
Julia_Cz = 1.237288
#endpreset

#preset 1
FOV = 0.4
Eye = -2.107985,0.0854123,-0.2484874
Target = -11.64699,1.814691,2.204474
Up = 0.1201045,0.9689693,-0.216041
AutoFocus = false
Gamma = 2
ToneMapping = 4
Exposure = 1
Brightness = 1
Contrast = 1
Saturation = 1
GaussianWeight = 1
AntiAliasScale = 2
DepthToAlpha = false
Detail = -3.06795
FudgeFactor = 0.1242673
Dither = 0.4180791
NormalBackStep = 0
CamLight = 1,0.9647059,0.9215686,1
BaseColor = 1,0.9803922,0.9411765
OrbitStrength = 0.8399072
X = 0.5,0.6,0.6,0.4040632
Y = 1,0.6,0,0.4
Z = 0.8,0.78,1,0.5
R = 0.4,0.7,1,0.12
BackgroundColor = 0.6,0.6,0.45
GradientBackground = 0.3
CycleColors = false
Cycles = 1.1
EnableFloor = false
FloorNormal = 0,0,1
FloorHeight = 0
FloorColor = 1,1,1
EquiRectangular = false
FocalPlane = 1
Aperture = 0
ShowDepth = false
DepthMagnitude = 1
DetailAO = -1.166667
MaxDistance = 100
MaxRaySteps = 502
AO = 0,0,0,0.928489
Specular = 0.4
SpecularExp = 16
SpecularMax = 10
SpotLight = 1,0.9607843,0.9215686,0.5639188
SpotLightDir = 1,0.1
CamLightMin = 0
Glow = 1,0.8862745,0.8078431,0.0849825
GlowMax = 222
Fog = 0.1160714
HardShadow = 0.7913247 NotLocked
ShadowSoft = 19.76471
QualityShadows = false
Reflection = 0 NotLocked
DebugSun = false NotLocked
MI = 16
bailout = 16
Scale = 1.915254
Min_R = 0.5
X_add = 0
Y_add = 0
Z_add = 0
Julia_Cx = 2.678161
Julia_Cy = 0.6206897
Julia_Cz = 0.6716418
#endpreset

// aBoxModKaliEiffie
// reference http://www.fractalforums.com/new-theories-and-research/aboxmodkali-the-2d-version/
// Created: 22/11/16
// INITIAL CONDITIONS

int swapXZ =  consts->fractal.custom[14]; 
float3 Offset1  = (float3){consts->fractal.custom[6], consts->fractal.custom[7], consts->fractal.custom[8]};
float minR2 =  consts->fractal.custom[9]; //default 0.25f (was 0.5)	

float Scale1 = consts->fractal.mandelbox.scale;
 //default 1.5f

float mdR2 = Scale1/minR2;

// addC, c3 default (0,0,0), cJ default (0,0,0)
float3 c3  = (float3){consts->fractal.custom[0], consts->fractal.custom[1], consts->fractal.custom[2]};

c3 = point * c3; // applying a scale factor to original position c

float3 cJ  = (float3){consts->fractal.custom[3], consts->fractal.custom[4], consts->fractal.custom[5]};// julia_C

float minimumR = 100.0f; // for coloring calculation	

float DE = 1.0f;

// aBoxModKaliEiffie
// reference http://www.fractalforums.com/new-theories-and-research/aboxmodkali-the-2d-version/
// Created: 22/11/16
// ITERATION LOOP

//a  box fold
z = fabs(z + Offset1) - fabs(z - Offset1) - z;// tglad fold

// the swizzling
if (swapXZ == 1)//  default ON
{
	z = (float3){z.z, z.y, z.x}; // swizzle swap
}

if (z.y > z.x) z = (float3){z.y, z.x, z.z}; // conditional swizzle

// a randomly placed 3D rotation
z = Matrix33MulFloat3(consts->fractal.mandelbox.mainRot, z);

// Abox type spherical fold, NO fixed radius
// no DE adjustment as we are modifying the scale value instead		
float r2 = dot(z,z);		
float m = Scale1;	
if (r2 < minR2) m =  mdR2;
else if (r2 < 1.0f) m =  m / r2;

z =  z * m + c3 + cJ;//   z * scale + addCpixel + addJuliaConstant  
 
//part of distance estimation calculation
DE = DE * fabs(m) + 1.0f;

//for Bailout
r2 = dot(z,z);
// for coloring
if (r2 < minimumR) minimumR = r2; // to be checked calibrated


//needed for calculation of fake lights based on orbit traps
#if _orbitTrapsEnabled
if (i >= consts->fractal.fakeLightsMinIter && i <= consts->fractal.fakeLightsMaxIter) distFromOrbitTrap += (1.0f/(r2));
#endif


if(r2>100.0f|| any(isinf(z)))	
{
	float r = sqrt(r2);
	dist = r / fabs(DE);

	out.colourIndex =   minimumR * 5000.0f;
	break;
		}

if(absC_sym	)
	{
		if (absC_ZX) c3.x = abs(c3.x);
		if (absC_ZY) c3.y = abs(c3.y);
		if (absC_ZZ) c3.z = abs(c3.z);

		c3.x = sign(p3.x) * c3.x;
		c3.y = sign(p3.y) * c3.y;
		c3.z = sign(p3.z) * c3.z;
	}	

uniform int MI; slider[1,30,100]
uniform int ColorIterations;  slider[0,10,100]
uniform int bailout; slider[1,1024,1024]
uniform float Scale; slider[-5,1.5,5]
uniform float MinR; slider[0,0.6,3]
uniform vec3 xyzAdd; slider[(-5,-5,-5),(0.5,0.5,0.5),(5,5,5)]
uniform bool Julia;  checkbox[false]
uniform vec3 Julia_C; slider[(-5,-5,-5),(0,0,0),(5,5,5)]

float DE(vec3 p) {
	vec3 z=vec3(p);
	float m;
	float mf=m/MinR;
	float rr=0; 
	float Dd = 1.0f; 
		for(int i=0; i<MI && rr<bailout; i++){
		m=Scale;
		z = xyzAdd - abs(z);
		rr = z.x*z.x + z.y*z.y + z.z*z.z;
		if (rr < MinR) {m = mf; } else 
			if (rr < 1) {m/=rr;};
		z*=m;
		if (Julia) {z+=Julia_C; }; 
		Dd = Dd * abs(m) + 1.0f;
		if (i < ColorIterations) orbitTrap = min(orbitTrap, abs(vec4(p,rr)));	
	}
float r = sqrt(rr);
return (r) / abs(Dd);
}