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,e)/1.2;
    p*=rotate3D(t,vec3(-8,6,4));
    for(j=0.;j++<8.;)
        p=abs(p/d)*rotate3D(i/j*2.,vec3(-1)),
        d=max(sin(length(-p)),cos(length(p/d)))/d,
        d-=dot(p,p)/1e4;
    e-=log(-d);
}o+=vec4(exp(p/d),1)*20./i;