Xor about 2 months ago
"Quasar"
for(float i,z,d,s;i++<1e2;o+=vec4(z,4,1,1)/d)
{vec3 p=z*normalize(FC.rgb*2.-r.xyy),a;p.z+=9.;s=length(a=mix(dot(a+=.57,p)*a,p,cos(s-=t))-sin(s)*cross(a,p));
for(d=1.;d++<9.;a+=sin(a*d-t).yzx/d);
z+=d=length(sin(a-t)+.7)*length(p.xy-p.z)/4e1;}
o=tanh(o/2e4);
