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,1)*8.;
    p.z-=4.;
    for(j=0.;j++<4.;)
        p=abs(p-1./p)*rotate3D(t*.2,vec3(-2,0,6)),
        p.xz*=rotate2D(t*.5),
        d=max(length(p.zx),length(p.yx))/d,
    e-=d;
}o+=vec4(p,1)*3./i;