for(float i,j,g,e,s;i++<1e2;g+=e*.1){vec3 p=vec3((FC.xy-.5r)/r.yg,g-1.6);p.yz*=rotate2D(5.3);p.xz*=rotate2D(t*.5);e=p.y;for(j=9.;j-->0.;e+=min(0.,dot(cos(p.xzsPI),r/r)/s*.5))s=exp(j-fract(t/PI)3.),p.xz=rotate2D(PI/3.);o.rgb+=exp(-e1e4)/1e2(p+.5)+fwidth(p);} // https://bit.ly/2YDJqxN