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)*4.;
    p.z-=4.;
    for(j=0.;j++<16.;)
        p=abs((1.-i/p))*rotate3D(t*.3,vec3(-4,6,4)),
        d=max(length(p.zx),dot(p,p))/1e4;
}
o+=vec4(exp(p),1)*3./i;