Deltographos::Blog

float i,j,d=.1,e; vec3 p;
for(i=0.;i++<64.&&d>1e-4;){
    p=vec3((FC.xy-.5*r)/r.y,e*d)*4.;
    p*=rotate3D(t,vec3(2,1,-2));
    for(j=0.;j++<12.;)
        p=abs(p*p)*rotate3D(j/i,vec3(1)),
        d=max(max(cos(p.y),sin(p.y)),atan(p.x))/d,
        d-=dot(p,p)/4e4;
    e-=d;
}o+=vec4(exp(p*p/e),1)*3./i;