https://twigl.app?ol=true&ss=-OKVZVCZDJ_JS7Qk5lDh

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)*4.;
    p.z-=4.;
    p.xy*=rotate2D(11.);
    for(j=0.;j++<8.;)
        p=abs(abs(p)-.405)-.4,
        p.xy*=rotate2D(PI/12.),
        d=max(max(p.z,p.y),p.x-cos(t)),
        d-=dot(p,p)/1e4;
    e-=d;
}o+=vec4(exp(p/d/i),1)*5./i;