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)*6.;
    for(j=0.;j++<8.;)
        p=mod(p-.5, i/j)*rotate3D(i/j,vec3(-2)),
        p.xz*=rotate2D(t),
        p.zy*=rotate2D(t),
        d=max(length(p.xz)-.5,length(p.zy)/5.)/d,
        d-=dot(p,p)/1e4;
    e-=d;
}o+=vec4(exp(p/e),1)*3./i;