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