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-.8);
    p*=rotate3D(t,vec3(0,4,2));
    for(j=0.;j++<4.;)
        p=abs(p*p)*rotate3D(i/j,vec3(-2,4,2)),
        d=max(sin(length(p.xy)),(length(p.yz)))/d,
        d-=dot(p/i,-p/i);
    e-=sqrt(d);
}o+=vec4(1./p,1)*3./i;