orage[idx*4+1]));\\n let z = float(atomicLoad(&atomic_storage[idx*4+2]));\\n return float3(x,y,z) / 255.;\\n}\\n\\nstruct Particle {\\n position: float4,\\n velocity: float4,\\n}\\n\\nfn LoadParticle(id: int2) -> Particle {\\n var p = Particle();\\n p.position = passLoad(0, id, 0); \\n p.velocity = passLoad(1, id, 0);\\n return p;\\n}\\n\\nfn SaveParticle(id: int2, p: Particle) {\\n passStore(0, id, p.position); \\n passStore(1, id, p.velocity); \\n}\\n\\n@stage(compute) @workgroup_size(16, 16)\\nfn SimulatePrey(@builtin(global_invocation_id) id: uint3) {\\n var p = LoadParticle(int2(id.xy));\\n\\n if(int(id.x) > PREY || int(id.y) > PREY) { \\n return;\\n }\\n\\n if(time.frame == 0u) {\\n p.position = hash42(float2(id.xy));\\n }\\n\\n var separation = float4();\\n var mean = Particle();\\n var count = 0;\\n\\n // other prey\\n for(var i = 0; i < PREY; i += 1) {\\n for(var j = 0; j < PREY; j += 1) {\\n let q = LoadParticle(int2(i,j));\\n let d = p.position - q.position;\\n if (length(d) < 1e-3) {\\n continue;\\n }\\n if (length(d) < custom.DistSeparation * 30.) {\\n separation += d / pow(length(d), 2.);\\n }\\n if (length(d) < custom.DistCohesion * 30.) {\\n mean.position += q.position;\\n mean.velocity += q.velocity;\\n count += 1;\\n }\\n }\\n }\\n\\n // predators\\n for(var i = PREY; i < PREY + PREDATORS; i += 1) {\\n for(var j = 0; j < PREDATORS; j += 1) {\\n let q = LoadParticle(int2(i,j));\\n let d = p.position - q.position;\\n if (length(d) < 1e-3) {\\n continue;\\n }\\n if (length(d) < custom.DistPrey * 30.)\\n {\\n separation += d / pow(length(d), 1.);\\n }\\n }\\n }\\n\\n if (count > 0) {\\n mean.position /= float(count);\\n mean.velocity /= float(count);\\n }\\n\\n var cohesion = mean.position - p.position;\\n var alignment = mean.velocity;\\n\\n let v = p.velocity;\\n p.velocity += clamp_length(normz(separation) - v, custom.MaxForce * .1) * 1.5;\\n p.velocity += clamp_length(normz(cohesion) - v, custom.MaxForce * .1);\\n p.velocity += clamp_length(normz(alignment) - v, custom.MaxForce * .1);\\n\\n if (length(p.position) > 20.) {\\n p.velocity -= pow(10., -10. * custom.Homing) * normalize(p.position);\\n }\\n\\n p.velocity = clamp_length(p.velocity, custom.MaxSpeed);\\n\\n p.position += custom.Timescale * p.velocity;\\n\\n SaveParticle(int2(id.xy), p);\\n}\\n\\n@stage(compute) @workgroup_size(16, 16)\\nfn SimulatePredators(@builtin(global_invocation_id) id: uint3) {\\n var p = LoadParticle(int2(id.xy));\\n\\n if(int(id.x) <= PREY || int(id.x) > PREY + PREDATORS || int(id.y) > PREDATORS) { \\n return;\\n }\\n\\n if(time.frame == 0u) {\\n p.position = hash42(float2(id.xy));\\n }\\n\\n var separation = float4();\\n\\n // other predators\\n for(var i = PREY; i < PREY + PREDATORS; i += 1) {\\n for(var j = 0; j < PREDATORS; j += 1) {\\n let q = LoadParticle(int2(i,j));\\n let d = p.position - q.position;\\n if (length(d) < 1e-3) {\\n continue;\\n }\\n if (length(d) < custom.DistSeparationPred * 30.) {\\n separation += d / pow(length(d), 1.);\\n }\\n }\\n }\\n\\n // prey\\n for(var i = 0; i < PREY; i += 1) {\\n for(var j = 0; j < PREY; j += 1) {\\n let q = LoadParticle(int2(i,j));\\n let d = p.position - q.position;\\n if (length(d) < 1e-3) {\\n continue;\\n }\\n if (length(d) < custom.DistPredation * 30.) {\\n separation -= d / pow(length(d), 1.);\\n }\\n }\\n }\\n\\n p.velocity = clamp_length(p.velocity +

@compute.toys