;\\n f = f + c * wave(a, r, p); c *= q;\\n f = f + c * wave(a, r, p * p); c *= q;\\n f = f + c * wave(a, r, p * p * p); c *= q;\\n f = f / 3.43;\\n return f;\\n}\\n\\n#dispatch_once init\\n@compute @workgroup_size(256)\\n#workgroup_count init COUNT_256 1 1\\nfn init(\\n @builtin(global_invocation_id) gid: vec3u,\\n @builtin(local_invocation_id) lid: vec3u,\\n @builtin(workgroup_id) wid: vec3u\\n) {\\n let random = pcg3df(vec3(gid.x, wid.x, lid.x));\\n let r = R * pow(random.x, 0.9) + .08;\\n let a = random.y * 2 * PI;\\n let x = r * sin(a);\\n let y = H * fbm(a, r) * exp(-r*r*1.5);\\n let z = r * cos(a);\\n let w = S + S * random.z; \\n let v = pow(custom.g*(.4+y*.15)/(r), 0.5);\\n let o = vec3f(custom.ox, custom.oy, custom.oz);\\n let s = Particle(\\n vec4f(x, y, z, w),\\n vec4f(v * cross(o - vec3f(x, y, z), vec3f(0,1,0)), 0)\\n );\\n particle[gid.x] = s;\\n}\\n\\n@compute @workgroup_size(1)\\n#workgroup_count clear TILES_X TILES_Y 1\\nfn clear(@builtin(global_invocation_id) wid: vec3u) {\\n atomicStore(&(tiles[wid.x][wid.y].amount), 0u);\\n}\\n\\nfn acc(r: vec3f) -> vec3f {\\n return -normalize(r) * custom.g * (max(exp(-dot(r,r)*4.), .95) + .01);\\n}\\n\\n@compute @workgroup_size(256)\\n#workgroup_count setup COUNT_256 1 1\\nfn setup(@builtin(global_invocation_id) gid: vec3u) {\\n //physics\\n var s = particle[gid.x];\\n let dt = time.delta * custom.time;\\n let m = 1.;//s.p.w*s.p.w;\\n let o = vec3f(custom.ox, custom.oy, custom.oz);\\n var a = acc(s.p.xyz - o);\\n s.p += vec4f((s.v.xyz + a * dt * 0.5) * dt, 0.);\\n a += acc(s.p.xyz - o);\\n s.v += vec4f(a * 0.5 * dt, 0.);\\n particle[gid.x] = s;\\n\\n let res = vec2f(textureDimensions(screen));\\n let sp = project(s.p, res);\\n let lb = (sp.xy - sp.z) / 16. - 1.;\\n let rt = min(res / 16., (sp.xy + sp.z) / 16.);\\n var id: u32;\\n for (var i = 0.; i < TILES_X; i += 1.) {\\n for (var j = 0.; j < TILES_Y; j += 1.) {\\n if (lb.x < i && i < rt.x && lb.y < j && j < rt.y) {\\n id = atomicAdd(&(tiles[u32(i)][u32(j)].amount), 1u);\\n tiles[u32(i)][u32(j)].ids[id % TILE_CAPACITY] = gid.x;\\n }\\n }\\n }\\n}\\n\\n@compute @workgroup_size(16, 16)\\n#workgroup_count rasterizer TILES_X TILES_Y 1\\nfn rasterizer(\\n @builtin(global_invocation_id) gid: vec3u,\\n @builtin(local_invocation_index) ldx: u32,\\n @builtin(workgroup_id) wid: vec3u\\n) {\\n // copy to wg memory\\n for (var i = 0u; i < u32(TILE_CAPACITY / 16 / 16); i++) {\\n let id = i * 16u * 16u + ldx;\\n wg_mem[id] = pack_Particle(particle[tiles[wid.x][wid.y].ids[id]]);\\n }\\n workgroupBarrier();\\n\\n let amount = max(min(atomicLoad(&tiles[wid.x][wid.y].amount), TILE_CAPACITY), 0);\\n let res = vec2f(textureDimensions(screen));\\n var col = vec3f(0.001, .0, .001);\\n //col += vec3f(1., 5., 1.)*f32(amount)/TILE_CAPACITY;\\n for (var i = 0.; i < AA; i += 1.) { for (var j = 0.; j < AA; j += 1.) {\\n let dxy = (vec2f(i, j) + .5) / AA;\\n let pix = vec2f(gid.xy) + dxy;\\n for (var id = 0u; id < amount; id++) {\\n let s = unpack_Particle(wg_mem[id]);\\n let p = project(s.p, res);\\n let L = length(s.p);\\n var size = p.z;\\n size *= exp(-L*L/12.);\\n let r = length(pix - p.xy);\\n let f = exp(-L*L/20.) * (1. - exp(-max(0, pow(.04*size/r, 1.2) - .01)));\\n col += .7*mix(vec3f(.2,.5,1.), vec3f(1.1, .5, .2), length(s.v*1.5)) * f / AA / AA;\\n }\\n }}\\n textureStore(screen, vec2u(gid.x, u32(res.y)-1-gid.y), vec4f(pow(col, vec3(2.2)), 1.));\\n}\"","uniforms":[{"name":"ox","value":0,"minRange":-1,"maxRange":1},{"name":"oy","value":0,"minRange":-1,"maxRange":1},{"name":"oz","value":0,"minRange":-1,"maxRange":1},{"name":"g","value":0.219,"minRange":0,"maxRange":1},{"name":"time","value":0.38,"minRange":0,"maxRange":2}],"textures":[{"img":"/textures/blank.png"},{"img":"/t

@compute.toys