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)*7.;
    p.z-=4.;
    for(j=0.;j++<5.;)
        p=abs(mod(p-7., 12.)-.05)*rotate3D(i/j-t,vec3(-4)),
        d=max(dot(p.zy,p.zx),length(p.xy))-5.,
        d-=dot(p,p)/1e4;
    e-=d;
}o+=vec4(exp(-p/i-.5),1)*8./i;