Prompt:

float i,j,d=.1,e;vec3 p,q=vec3(-1);for(i=0.;i++<64.&&d>1e-4;){p=vec3((FC.xy-.5*r)/r.y,1)*e;p.z-=.5;p*=rotate3D(t,q);for(j=0.;j++<6.;)p=abs(abs(p-length(cos(q)))-.05)-.04,d=max(max(p.z,p.x),p.y),d-=dot(cos(p.xy)/d,sin(p.zy)/d)/2e4;e+=d;}o+=vec4(exp(p/d),1)*3./i;