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

kl = new V2D(1.95, 0.03);
b_a = kl.y / kl.x;
box = new V2D(1.0);
extraIter = 5;
balls = true;
fourGen = false;
clamp = new V2D(0.3, 2.0);
doInv = false;
invCenter = new V3D(1, 0.96, 0);
invR = 1;
recenter = new V3D(0);

void DE() {
V3D p = previous.ergebnis;
V3D q = p.copy();
V3D lz = q.getAdd(V3D.V1);
V3D llz = q.getSub(V3D.V1);
DF = 1.0;
double DE = 1e10;
double f = Math.signum(kl.y);
for(int i = 0 ; i < maxIterM ; i++){
q.x = q.x + b_a * q.y;
if(fourGen){
q.x = wrap(q.x, 2.0 * box.x, -box.x);
q.x = wrap(q.y, kl.x, 0);
q.x = wrap(q.z, 2.0 * box.y, -box.y);
}else{
q.x = wrap(q.x, 2.0 * box.x, -box.x);
q.x = wrap(q.z, 2.0 * box.y, -box.y);
}
q.x = q.x - b_a * q.y;
if(q.y >= kl.x * (0.5 + f * (2.0 * kl.x - 1.95) * 0.25 * Math.signum(q.x + kl.y * 0.5) * (1.0 - Math
.exp(-(7.2 - (1.95 - kl.x) * 15.0) * Math.abs(q.x + kl.y * 0.5))))){
q.x = -kl.y - q.x;
q.y = kl.x - q.y;
q.z = -q.z;
}
q = TransA(q);
V3D tmp = q.getSub(llz);
if(tmp.dot(tmp) < 1e-5){
break;
}
llz = lz;
lz = q;
}
double y = balls ? Math.min(q.y, kl.x - q.y) : q.y;
DE = Math.min(DE, Math.min(y, clamp.x) / Math.max(DF, clamp.y));
dist = DE;
hitMat = mat;
}

private V3D TransA(V3D q) {
double r = 1.0 / q.dot(q);
q.mult(-r);
q.x = -kl.y - q.x;
q.y = kl.x + q.y;
DF *= r;
return q;
}

private double wrap(double x, double a, double s) {
x -= s;
return x - a * Math.floor(x / a) + s;
}