Deltographos::Blog

float i,j,d=.01,e; vec3 p;
for(i=0.;i++<64.&&d>1e-4;){
    p=vec3((FC.xy-.5*r)*r.y,1);
    p.z-=4.;
    p*=rotate3D(t,vec3(-4,6,-4));
    for(j=0.;j++<8.;)
        p=abs(p*.08)-.015, 
        p*=rotate3D(i/j,vec3(-1)),
        p.zy*=rotate2D(t),
        d=max(cos(p.x),sin(p.z))/p.y/d,
        d-=dot(p,p)/1e4;
    e-=d;
}o+=vec4(log(p*p)/d,1)*3./i;