2022-01-19 22:46:57 +00:00
|
|
|
@stage(vertex)
|
|
|
|
fn vert_main(@location(0) a_particlePos : vec2<f32>, @location(1) a_particleVel : vec2<f32>, @location(2) a_pos : vec2<f32>) -> @builtin(position) vec4<f32> {
|
2021-05-18 09:24:18 +00:00
|
|
|
var angle : f32 = -(atan2(a_particleVel.x, a_particleVel.y));
|
|
|
|
var pos : vec2<f32> = vec2<f32>(((a_pos.x * cos(angle)) - (a_pos.y * sin(angle))), ((a_pos.x * sin(angle)) + (a_pos.y * cos(angle))));
|
|
|
|
return vec4<f32>((pos + a_particlePos), 0.0, 1.0);
|
|
|
|
}
|
|
|
|
|
2022-01-19 22:46:57 +00:00
|
|
|
@stage(fragment)
|
|
|
|
fn frag_main() -> @location(0) vec4<f32> {
|
2021-05-18 09:24:18 +00:00
|
|
|
return vec4<f32>(1.0, 1.0, 1.0, 1.0);
|
|
|
|
}
|
|
|
|
|
|
|
|
struct Particle {
|
|
|
|
pos : vec2<f32>;
|
|
|
|
vel : vec2<f32>;
|
2022-01-19 18:11:17 +00:00
|
|
|
}
|
2021-05-18 09:24:18 +00:00
|
|
|
|
|
|
|
struct SimParams {
|
|
|
|
deltaT : f32;
|
|
|
|
rule1Distance : f32;
|
|
|
|
rule2Distance : f32;
|
|
|
|
rule3Distance : f32;
|
|
|
|
rule1Scale : f32;
|
|
|
|
rule2Scale : f32;
|
|
|
|
rule3Scale : f32;
|
2022-01-19 18:11:17 +00:00
|
|
|
}
|
2021-05-18 09:24:18 +00:00
|
|
|
|
|
|
|
struct Particles {
|
|
|
|
particles : array<Particle, 5>;
|
2022-01-19 18:11:17 +00:00
|
|
|
}
|
2021-05-18 09:24:18 +00:00
|
|
|
|
2022-01-19 22:46:57 +00:00
|
|
|
@binding(0) @group(0) var<uniform> params : SimParams;
|
2021-05-18 09:24:18 +00:00
|
|
|
|
2022-01-19 22:46:57 +00:00
|
|
|
@binding(1) @group(0) var<storage, read_write> particlesA : Particles;
|
2021-05-18 09:24:18 +00:00
|
|
|
|
2022-01-19 22:46:57 +00:00
|
|
|
@binding(2) @group(0) var<storage, read_write> particlesB : Particles;
|
2021-05-18 09:24:18 +00:00
|
|
|
|
2022-01-19 22:46:57 +00:00
|
|
|
@stage(compute) @workgroup_size(1)
|
|
|
|
fn comp_main(@builtin(global_invocation_id) gl_GlobalInvocationID : vec3<u32>) {
|
2021-05-18 09:24:18 +00:00
|
|
|
var index : u32 = gl_GlobalInvocationID.x;
|
|
|
|
if ((index >= 5u)) {
|
|
|
|
return;
|
|
|
|
}
|
|
|
|
var vPos : vec2<f32> = particlesA.particles[index].pos;
|
|
|
|
var vVel : vec2<f32> = particlesA.particles[index].vel;
|
|
|
|
var cMass : vec2<f32> = vec2<f32>(0.0, 0.0);
|
|
|
|
var cVel : vec2<f32> = vec2<f32>(0.0, 0.0);
|
|
|
|
var colVel : vec2<f32> = vec2<f32>(0.0, 0.0);
|
|
|
|
var cMassCount : i32 = 0;
|
|
|
|
var cVelCount : i32 = 0;
|
|
|
|
var pos : vec2<f32>;
|
|
|
|
var vel : vec2<f32>;
|
2021-07-08 21:23:33 +00:00
|
|
|
for(var i : u32 = 0u; (i < 5u); i = (i + 1u)) {
|
|
|
|
if ((i == index)) {
|
|
|
|
continue;
|
|
|
|
}
|
|
|
|
pos = particlesA.particles[i].pos.xy;
|
|
|
|
vel = particlesA.particles[i].vel.xy;
|
|
|
|
if ((distance(pos, vPos) < params.rule1Distance)) {
|
|
|
|
cMass = (cMass + pos);
|
|
|
|
cMassCount = (cMassCount + 1);
|
|
|
|
}
|
|
|
|
if ((distance(pos, vPos) < params.rule2Distance)) {
|
|
|
|
colVel = (colVel - (pos - vPos));
|
|
|
|
}
|
|
|
|
if ((distance(pos, vPos) < params.rule3Distance)) {
|
|
|
|
cVel = (cVel + vel);
|
|
|
|
cVelCount = (cVelCount + 1);
|
2021-05-18 09:24:18 +00:00
|
|
|
}
|
|
|
|
}
|
|
|
|
if ((cMassCount > 0)) {
|
|
|
|
cMass = ((cMass / vec2<f32>(f32(cMassCount), f32(cMassCount))) - vPos);
|
|
|
|
}
|
|
|
|
if ((cVelCount > 0)) {
|
|
|
|
cVel = (cVel / vec2<f32>(f32(cVelCount), f32(cVelCount)));
|
|
|
|
}
|
|
|
|
vVel = (((vVel + (cMass * params.rule1Scale)) + (colVel * params.rule2Scale)) + (cVel * params.rule3Scale));
|
|
|
|
vVel = (normalize(vVel) * clamp(length(vVel), 0.0, 0.100000001));
|
|
|
|
vPos = (vPos + (vVel * params.deltaT));
|
|
|
|
if ((vPos.x < -1.0)) {
|
|
|
|
vPos.x = 1.0;
|
|
|
|
}
|
|
|
|
if ((vPos.x > 1.0)) {
|
|
|
|
vPos.x = -1.0;
|
|
|
|
}
|
|
|
|
if ((vPos.y < -1.0)) {
|
|
|
|
vPos.y = 1.0;
|
|
|
|
}
|
|
|
|
if ((vPos.y > 1.0)) {
|
|
|
|
vPos.y = -1.0;
|
|
|
|
}
|
|
|
|
particlesB.particles[index].pos = vPos;
|
|
|
|
particlesB.particles[index].vel = vVel;
|
|
|
|
}
|