shaders 3.0.453 → 3.1.455

package/dist/core/{ParticleFlow-DfdDvN-j.js → ParticleFlow-Ck5Ta_Rj.js}

@@ -1,13 +1,13 @@
-import { At as max, Dt as floor, G as mixColorsVariants, Gt as data_exports, Ht as sqrt, Jt as texture2d, Kt as atomic, Mt as mix, Ot as fract, Tt as exp, W as mixColorsLinear, Yt as textureStorage2d, _t as abs, an as vec2u, bt as ceil, cn as vec4f, ct as textureStore, d as ZERO, dn as u32, en as struct, ft as atomicAdd, gt as select, it as std_exports, jt as min, kt as length, ln as f32, mt as atomicStore, nt as src_default, on as vec3f, ot as textureLoad, pt as atomicLoad, qt as arrayOf, rn as vec2f, un as i32, xt as clamp } from "./edges-B7Yet3LV.js";
-import { l as hash33, p as noise_exports, r as curl22z } from "./noise-B-PI7JUx.js";
-import { n as createGuardedCompute, r as createStateBuffer } from "./compute-Ax-p6Ubc.js";
-import { t as isMobileGpuViewport } from "./device-BODv_4aH.js";
-import { a as makeAgentTexelWeightFn, l as resolveAgentShape, s as orientedShapeOptions, t as AGENT_SHAPES } from "./agents-BVhq_Z2T.js";
-import { a as transformColor, o as transformColorSpace, t as colorSpaceOptions } from "./transformations-C7gn_5rF.js";
+import { At as max, Dt as floor, G as mixColorsLinear, Gt as data_exports, Ht as sqrt, Jt as texture2d, K as mixColorsVariants, Kt as atomic, Mt as mix, Ot as fract, Tt as exp, Yt as textureStorage2d, _t as abs, an as vec2u, at as std_exports, bt as ceil, cn as vec4f, dn as i32, en as struct, f as ZERO, fn as u32, gt as select, ht as atomicStore, jt as min, kt as length, lt as textureStore, mt as atomicLoad, on as vec3f, pt as atomicAdd, qt as arrayOf, rn as vec2f, rt as src_default, st as textureLoad, un as f32, xt as clamp } from "./edges-CD7gCAto.js";
+import { a as gaussianBrushSq, c as createStableFluidsPasses, l as neighbourIndex, s as buildStableFluidsKernels } from "./porters-DUtQounc.js";
+import { l as hash33, p as noise_exports, r as curl22z } from "./noise-Bd6IGS0W.js";
+import { M as createGuardedCompute, P as createStateBuffer, t as isMobileGpuViewport } from "./device-B0ZwBfOg.js";
+import { a as makeAgentTexelWeightFn, l as resolveAgentShape, s as orientedShapeOptions, t as AGENT_SHAPES } from "./agents-C_kohsYu.js";
+import { a as transformColor, o as transformColorSpace, t as colorSpaceOptions } from "./transformations-B-zNWNdd.js";
+import { i as pathStampRibbon, n as createPointerVelocityTracker, t as createIdleGate } from "./pointer-DGDGF-yO.js";
 var N = 256;
 var FLUID_COUNT = N * N;
 var JACOBI_ITERS = isMobileGpuViewport() ? 6 : 10;
-var PRESSURE_DISSIPATION = .8;
 var VEL_FORMAT = "rgba16float";
 var RES = isMobileGpuViewport() ? 512 : 1024;
 var INV_RES = 1 / RES;
@@ -97,125 +97,7 @@ var pfLayout = /* @__PURE__ */ (globalThis.__TYPEGPU_AUTONAME__ ?? ((a) => a))(s
 	params: { uniform: SimParams },
 	outTex: { storageTexture: textureStorage2d(STATE_FORMAT, "write-only") }
 }), "pfLayout");
-var nidx = src_default.fn([
-	i32,
-	i32,
-	i32,
-	i32
-], u32)(/* @__PURE__ */ (($) => (globalThis.__TYPEGPU_META__ ??= /* @__PURE__ */ new WeakMap()).set($.f = ((ii, ji, di, dj) => {
-	"use gpu";
-	const row = clamp(__tsover_add(ii, di), 0, __tsover_sub(N, 1));
-	const col = clamp(__tsover_add(ji, dj), 0, __tsover_sub(N, 1));
-	return u32(__tsover_add(__tsover_mul(row, N), col));
-}), {
-	v: 1,
-	name: void 0,
-	ast: {
-		"params": [
-			{
-				"type": "i",
-				"name": "ii"
-			},
-			{
-				"type": "i",
-				"name": "ji"
-			},
-			{
-				"type": "i",
-				"name": "di"
-			},
-			{
-				"type": "i",
-				"name": "dj"
-			}
-		],
-		"body": [0, [
-			[
-				13,
-				"row",
-				[
-					6,
-					[
-						7,
-						"std",
-						"clamp"
-					],
-					[
-						[
-							1,
-							"ii",
-							"+",
-							"di"
-						],
-						[5, "0"],
-						[
-							1,
-							"N",
-							"-",
-							[5, "1"]
-						]
-					]
-				]
-			],
-			[
-				13,
-				"col",
-				[
-					6,
-					[
-						7,
-						"std",
-						"clamp"
-					],
-					[
-						[
-							1,
-							"ji",
-							"+",
-							"dj"
-						],
-						[5, "0"],
-						[
-							1,
-							"N",
-							"-",
-							[5, "1"]
-						]
-					]
-				]
-			],
-			[10, [
-				6,
-				[
-					7,
-					"d",
-					"u32"
-				],
-				[[
-					1,
-					[
-						1,
-						"row",
-						"*",
-						"N"
-					],
-					"+",
-					"col"
-				]]
-			]]
-		]],
-		"externalNames": [
-			"std",
-			"N",
-			"d"
-		]
-	},
-	externals: () => ({
-		std: std_exports,
-		N,
-		d: data_exports
-	})
-}) && $.f)({}));
+var nidx = neighbourIndex(N, "clamped");
 const particleFlowForceKernel = src_default.fn([u32, u32])(/* @__PURE__ */ (($) => (globalThis.__TYPEGPU_META__ ??= /* @__PURE__ */ new WeakMap()).set($.f = ((cx, cy) => {
 	"use gpu";
 	const p = fluidLayout.$.sparams;
@@ -224,9 +106,7 @@ const particleFlowForceKernel = src_default.fn([u32, u32])(/* @__PURE__ */ (($)
 	const cellY = __tsover_add(f32(cy), .5);
 	const dx = __tsover_sub(cellX, p.posX);
 	const dy = __tsover_sub(cellY, p.posY);
-	const distSq = __tsover_add(__tsover_mul(dx, dx), __tsover_mul(dy, dy));
-	const radSq = max(__tsover_mul(p.radius, p.radius), f32(1));
-	const inf = exp(__tsover_mul(__tsover_div(distSq, radSq), -1));
+	const inf = gaussianBrushSq(__tsover_add(__tsover_mul(dx, dx), __tsover_mul(dy, dy)), max(__tsover_mul(p.radius, p.radius), f32(1)));
 	const vel = fluidLayout.$.velA[idx];
 	fluidLayout.$.velA[idx] = vec4f(__tsover_add(vel.x, __tsover_mul(p.velX, inf)), __tsover_add(vel.y, __tsover_mul(p.velY, inf)), vel.z, 0);
 }), {
@@ -400,26 +280,8 @@ const particleFlowForceKernel = src_default.fn([u32, u32])(/* @__PURE__ */ (($)
 				"inf",
 				[
 					6,
-					[
-						7,
-						"std",
-						"exp"
-					],
-					[[
-						1,
-						[
-							1,
-							"distSq",
-							"/",
-							"radSq"
-						],
-						"*",
-						[
-							4,
-							"-",
-							[5, "1"]
-						]
-					]]
+					"gaussianBrushSq",
+					["distSq", "radSq"]
 				]
 			],
 			[
@@ -515,28 +377,54 @@ const particleFlowForceKernel = src_default.fn([u32, u32])(/* @__PURE__ */ (($)
 			"fluidLayout",
 			"d",
 			"N",
-			"std"
+			"std",
+			"gaussianBrushSq"
 		]
 	},
 	externals: () => ({
 		fluidLayout,
 		d: data_exports,
 		N,
-		std: std_exports
+		std: std_exports,
+		gaussianBrushSq
 	})
 }) && $.f)({})).$name("particleFlowForce");
-const particleFlowCurlKernel = src_default.fn([u32, u32])(/* @__PURE__ */ (($) => (globalThis.__TYPEGPU_META__ ??= /* @__PURE__ */ new WeakMap()).set($.f = ((cx, cy) => {
+var solverKernels = buildStableFluidsKernels(fluidLayout, {
+	n: N,
+	namePrefix: "particleFlow",
+	boundary: "clamped",
+	dye: "none",
+	publishVelocityTexture: true
+});
+const { curl: particleFlowCurlKernel, divergence: particleFlowDivergenceKernel, jacobi: particleFlowJacobiKernel, gradSubtract: particleFlowGradSubtractKernel, advectVel: particleFlowAdvectVelKernel, copyVel: particleFlowCopyVelKernel } = solverKernels;
+const particleFlowVorticityKernel = src_default.fn([u32, u32])(/* @__PURE__ */ (($) => (globalThis.__TYPEGPU_META__ ??= /* @__PURE__ */ new WeakMap()).set($.f = ((cx, cy) => {
 	"use gpu";
+	const p = fluidLayout.$.params;
 	const idx = __tsover_add(__tsover_mul(cy, u32(N)), cx);
 	const ii = i32(cy);
 	const ji = i32(cx);
-	const vL = fluidLayout.$.velA[nidx(ii, ji, 0, -1)];
-	const vR = fluidLayout.$.velA[nidx(ii, ji, 0, 1)];
-	const vU = fluidLayout.$.velA[nidx(ii, ji, -1, 0)];
-	const vD = fluidLayout.$.velA[nidx(ii, ji, 1, 0)];
-	const curl = __tsover_mul(__tsover_sub(__tsover_sub(vR.y, vL.y), __tsover_sub(vD.x, vU.x)), .5);
+	const cL = abs(fluidLayout.$.velA[nidx(ii, ji, 0, -1)].z);
+	const cR = abs(fluidLayout.$.velA[nidx(ii, ji, 0, 1)].z);
+	const cU = abs(fluidLayout.$.velA[nidx(ii, ji, -1, 0)].z);
+	const cD = abs(fluidLayout.$.velA[nidx(ii, ji, 1, 0)].z);
+	const cc = fluidLayout.$.velA[idx].z;
+	const gx = __tsover_mul(__tsover_sub(cR, cL), .5);
+	const gy = __tsover_mul(__tsover_sub(cD, cU), .5);
+	const gl = max(sqrt(__tsover_add(__tsover_mul(gx, gx), __tsover_mul(gy, gy))), f32(1e-5));
+	const fx = __tsover_mul(__tsover_mul(__tsover_mul(__tsover_div(gy, gl), cc), p.curlStrength), p.dt);
+	const fy = __tsover_mul(__tsover_mul(__tsover_mul(__tsover_mul(__tsover_div(gx, gl), -1), cc), p.curlStrength), p.dt);
+	let ax = f32(0);
+	let ay = f32(0);
+	if (p.ambient > 0) {
+		const uvx = __tsover_div(__tsover_add(f32(cx), .5), f32(N));
+		const uvy = __tsover_div(__tsover_add(f32(cy), .5), f32(N));
+		const amb = curl22z(vec2f(__tsover_mul(uvx, p.ambientFreq), __tsover_mul(uvy, p.ambientFreq)), p.ambientTime);
+		const af = __tsover_mul(__tsover_mul(p.ambient, AMBIENT_K), p.dt);
+		ax = __tsover_mul(amb.x, af);
+		ay = __tsover_mul(amb.y, af);
+	}
 	const vel = fluidLayout.$.velA[idx];
-	fluidLayout.$.velA[idx] = vec4f(vel.x, vel.y, curl, 0);
+	fluidLayout.$.velA[idx] = vec4f(__tsover_add(__tsover_add(vel.x, fx), ax), __tsover_add(__tsover_add(vel.y, fy), ay), vel.z, 0);
 }), {
 	v: 1,
 	name: void 0,
@@ -549,6 +437,19 @@ const particleFlowCurlKernel = src_default.fn([u32, u32])(/* @__PURE__ */ (($) =
 			"name": "cy"
 		}],
 		"body": [0, [
+			[
+				13,
+				"p",
+				[
+					7,
+					[
+						7,
+						"fluidLayout",
+						"$"
+					],
+					"params"
+				]
+			],
 			[
 				13,
 				"idx",
@@ -600,339 +501,93 @@ const particleFlowCurlKernel = src_default.fn([u32, u32])(/* @__PURE__ */ (($) =
 			],
 			[
 				13,
-				"vL",
+				"cL",
 				[
-					8,
+					6,
 					[
 						7,
-						[
-							7,
-							"fluidLayout",
-							"$"
-						],
-						"velA"
+						"std",
+						"abs"
 					],
-					[
-						6,
-						"nidx",
+					[[
+						7,
 						[
-							"ii",
-							"ji",
-							[5, "0"],
+							8,
 							[
-								4,
-								"-",
-								[5, "1"]
+								7,
+								[
+									7,
+									"fluidLayout",
+									"$"
+								],
+								"velA"
+							],
+							[
+								6,
+								"nidx",
+								[
+									"ii",
+									"ji",
+									[5, "0"],
+									[
+										4,
+										"-",
+										[5, "1"]
+									]
+								]
 							]
-						]
-					]
-				]
-			],
-			[
-				13,
-				"vR",
-				[
-					8,
-					[
-						7,
-						[
-							7,
-							"fluidLayout",
-							"$"
 						],
-						"velA"
-					],
-					[
-						6,
-						"nidx",
-						[
-							"ii",
-							"ji",
-							[5, "0"],
-							[5, "1"]
-						]
-					]
+						"z"
+					]]
 				]
 			],
 			[
 				13,
-				"vU",
+				"cR",
 				[
-					8,
+					6,
 					[
 						7,
-						[
-							7,
-							"fluidLayout",
-							"$"
-						],
-						"velA"
+						"std",
+						"abs"
 					],
-					[
-						6,
-						"nidx",
+					[[
+						7,
 						[
-							"ii",
-							"ji",
+							8,
 							[
-								4,
-								"-",
-								[5, "1"]
+								7,
+								[
+									7,
+									"fluidLayout",
+									"$"
+								],
+								"velA"
 							],
-							[5, "0"]
-						]
-					]
+							[
+								6,
+								"nidx",
+								[
+									"ii",
+									"ji",
+									[5, "0"],
+									[5, "1"]
+								]
+							]
+						],
+						"z"
+					]]
 				]
 			],
 			[
 				13,
-				"vD",
+				"cU",
 				[
-					8,
+					6,
 					[
 						7,
-						[
-							7,
-							"fluidLayout",
-							"$"
-						],
-						"velA"
-					],
-					[
-						6,
-						"nidx",
-						[
-							"ii",
-							"ji",
-							[5, "1"],
-							[5, "0"]
-						]
-					]
-				]
-			],
-			[
-				13,
-				"curl",
-				[
-					1,
-					[
-						1,
-						[
-							1,
-							[
-								7,
-								"vR",
-								"y"
-							],
-							"-",
-							[
-								7,
-								"vL",
-								"y"
-							]
-						],
-						"-",
-						[
-							1,
-							[
-								7,
-								"vD",
-								"x"
-							],
-							"-",
-							[
-								7,
-								"vU",
-								"x"
-							]
-						]
-					],
-					"*",
-					[5, "0.5"]
-				]
-			],
-			[
-				13,
-				"vel",
-				[
-					8,
-					[
-						7,
-						[
-							7,
-							"fluidLayout",
-							"$"
-						],
-						"velA"
-					],
-					"idx"
-				]
-			],
-			[
-				2,
-				[
-					8,
-					[
-						7,
-						[
-							7,
-							"fluidLayout",
-							"$"
-						],
-						"velA"
-					],
-					"idx"
-				],
-				"=",
-				[
-					6,
-					[
-						7,
-						"d",
-						"vec4f"
-					],
-					[
-						[
-							7,
-							"vel",
-							"x"
-						],
-						[
-							7,
-							"vel",
-							"y"
-						],
-						"curl",
-						[5, "0"]
-					]
-				]
-			]
-		]],
-		"externalNames": [
-			"d",
-			"N",
-			"fluidLayout",
-			"nidx"
-		]
-	},
-	externals: () => ({
-		d: data_exports,
-		N,
-		fluidLayout,
-		nidx
-	})
-}) && $.f)({})).$name("particleFlowCurl");
-const particleFlowVorticityKernel = src_default.fn([u32, u32])(/* @__PURE__ */ (($) => (globalThis.__TYPEGPU_META__ ??= /* @__PURE__ */ new WeakMap()).set($.f = ((cx, cy) => {
-	"use gpu";
-	const p = fluidLayout.$.params;
-	const idx = __tsover_add(__tsover_mul(cy, u32(N)), cx);
-	const ii = i32(cy);
-	const ji = i32(cx);
-	const cL = abs(fluidLayout.$.velA[nidx(ii, ji, 0, -1)].z);
-	const cR = abs(fluidLayout.$.velA[nidx(ii, ji, 0, 1)].z);
-	const cU = abs(fluidLayout.$.velA[nidx(ii, ji, -1, 0)].z);
-	const cD = abs(fluidLayout.$.velA[nidx(ii, ji, 1, 0)].z);
-	const cc = fluidLayout.$.velA[idx].z;
-	const gx = __tsover_mul(__tsover_sub(cR, cL), .5);
-	const gy = __tsover_mul(__tsover_sub(cD, cU), .5);
-	const gl = max(sqrt(__tsover_add(__tsover_mul(gx, gx), __tsover_mul(gy, gy))), f32(1e-5));
-	const fx = __tsover_mul(__tsover_mul(__tsover_mul(__tsover_div(gy, gl), cc), p.curlStrength), p.dt);
-	const fy = __tsover_mul(__tsover_mul(__tsover_mul(__tsover_mul(__tsover_div(gx, gl), -1), cc), p.curlStrength), p.dt);
-	let ax = f32(0);
-	let ay = f32(0);
-	if (p.ambient > 0) {
-		const uvx = __tsover_div(__tsover_add(f32(cx), .5), f32(N));
-		const uvy = __tsover_div(__tsover_add(f32(cy), .5), f32(N));
-		const amb = curl22z(vec2f(__tsover_mul(uvx, p.ambientFreq), __tsover_mul(uvy, p.ambientFreq)), p.ambientTime);
-		const af = __tsover_mul(__tsover_mul(p.ambient, AMBIENT_K), p.dt);
-		ax = __tsover_mul(amb.x, af);
-		ay = __tsover_mul(amb.y, af);
-	}
-	const vel = fluidLayout.$.velA[idx];
-	fluidLayout.$.velA[idx] = vec4f(__tsover_add(__tsover_add(vel.x, fx), ax), __tsover_add(__tsover_add(vel.y, fy), ay), vel.z, 0);
-}), {
-	v: 1,
-	name: void 0,
-	ast: {
-		"params": [{
-			"type": "i",
-			"name": "cx"
-		}, {
-			"type": "i",
-			"name": "cy"
-		}],
-		"body": [0, [
-			[
-				13,
-				"p",
-				[
-					7,
-					[
-						7,
-						"fluidLayout",
-						"$"
-					],
-					"params"
-				]
-			],
-			[
-				13,
-				"idx",
-				[
-					1,
-					[
-						1,
-						"cy",
-						"*",
-						[
-							6,
-							[
-								7,
-								"d",
-								"u32"
-							],
-							["N"]
-						]
-					],
-					"+",
-					"cx"
-				]
-			],
-			[
-				13,
-				"ii",
-				[
-					6,
-					[
-						7,
-						"d",
-						"i32"
-					],
-					["cy"]
-				]
-			],
-			[
-				13,
-				"ji",
-				[
-					6,
-					[
-						7,
-						"d",
-						"i32"
-					],
-					["cx"]
-				]
-			],
-			[
-				13,
-				"cL",
-				[
-					6,
-					[
-						7,
-						"std",
-						"abs"
+						"std",
+						"abs"
 					],
 					[[
 						7,
@@ -953,12 +608,12 @@ const particleFlowVorticityKernel = src_default.fn([u32, u32])(/* @__PURE__ */ (
 								[
 									"ii",
 									"ji",
-									[5, "0"],
 									[
 										4,
 										"-",
 										[5, "1"]
-									]
+									],
+									[5, "0"]
 								]
 							]
 						],
@@ -968,7 +623,7 @@ const particleFlowVorticityKernel = src_default.fn([u32, u32])(/* @__PURE__ */ (
 			],
 			[
 				13,
-				"cR",
+				"cD",
 				[
 					6,
 					[
@@ -995,8 +650,8 @@ const particleFlowVorticityKernel = src_default.fn([u32, u32])(/* @__PURE__ */ (
 								[
 									"ii",
 									"ji",
-									[5, "0"],
-									[5, "1"]
+									[5, "1"],
+									[5, "0"]
 								]
 							]
 						],
@@ -1006,97 +661,17 @@ const particleFlowVorticityKernel = src_default.fn([u32, u32])(/* @__PURE__ */ (
 			],
 			[
 				13,
-				"cU",
+				"cc",
 				[
-					6,
+					7,
 					[
-						7,
-						"std",
-						"abs"
-					],
-					[[
-						7,
+						8,
 						[
-							8,
+							7,
 							[
 								7,
-								[
-									7,
-									"fluidLayout",
-									"$"
-								],
-								"velA"
-							],
-							[
-								6,
-								"nidx",
-								[
-									"ii",
-									"ji",
-									[
-										4,
-										"-",
-										[5, "1"]
-									],
-									[5, "0"]
-								]
-							]
-						],
-						"z"
-					]]
-				]
-			],
-			[
-				13,
-				"cD",
-				[
-					6,
-					[
-						7,
-						"std",
-						"abs"
-					],
-					[[
-						7,
-						[
-							8,
-							[
-								7,
-								[
-									7,
-									"fluidLayout",
-									"$"
-								],
-								"velA"
-							],
-							[
-								6,
-								"nidx",
-								[
-									"ii",
-									"ji",
-									[5, "1"],
-									[5, "0"]
-								]
-							]
-						],
-						"z"
-					]]
-				]
-			],
-			[
-				13,
-				"cc",
-				[
-					7,
-					[
-						8,
-						[
-							7,
-							[
-								7,
-								"fluidLayout",
-								"$"
+								"fluidLayout",
+								"$"
 							],
 							"velA"
 						],
@@ -1399,1928 +974,63 @@ const particleFlowVorticityKernel = src_default.fn([u32, u32])(/* @__PURE__ */ (
 						]
 					],
 					[
-						13,
-						"af",
-						[
-							1,
-							[
-								1,
-								[
-									7,
-									"p",
-									"ambient"
-								],
-								"*",
-								"AMBIENT_K"
-							],
-							"*",
-							[
-								7,
-								"p",
-								"dt"
-							]
-						]
-					],
-					[
-						2,
-						"ax",
-						"=",
-						[
-							1,
-							[
-								7,
-								"amb",
-								"x"
-							],
-							"*",
-							"af"
-						]
-					],
-					[
-						2,
-						"ay",
-						"=",
-						[
-							1,
-							[
-								7,
-								"amb",
-								"y"
-							],
-							"*",
-							"af"
-						]
-					]
-				]]
-			],
-			[
-				13,
-				"vel",
-				[
-					8,
-					[
-						7,
-						[
-							7,
-							"fluidLayout",
-							"$"
-						],
-						"velA"
-					],
-					"idx"
-				]
-			],
-			[
-				2,
-				[
-					8,
-					[
-						7,
-						[
-							7,
-							"fluidLayout",
-							"$"
-						],
-						"velA"
-					],
-					"idx"
-				],
-				"=",
-				[
-					6,
-					[
-						7,
-						"d",
-						"vec4f"
-					],
-					[
-						[
-							1,
-							[
-								1,
-								[
-									7,
-									"vel",
-									"x"
-								],
-								"+",
-								"fx"
-							],
-							"+",
-							"ax"
-						],
-						[
-							1,
-							[
-								1,
-								[
-									7,
-									"vel",
-									"y"
-								],
-								"+",
-								"fy"
-							],
-							"+",
-							"ay"
-						],
-						[
-							7,
-							"vel",
-							"z"
-						],
-						[5, "0"]
-					]
-				]
-			]
-		]],
-		"externalNames": [
-			"fluidLayout",
-			"d",
-			"N",
-			"std",
-			"nidx",
-			"noise",
-			"AMBIENT_K"
-		]
-	},
-	externals: () => ({
-		fluidLayout,
-		d: data_exports,
-		N,
-		std: std_exports,
-		nidx,
-		noise: noise_exports,
-		AMBIENT_K
-	})
-}) && $.f)({})).$name("particleFlowVorticity");
-const particleFlowDivergenceKernel = src_default.fn([u32, u32])(/* @__PURE__ */ (($) => (globalThis.__TYPEGPU_META__ ??= /* @__PURE__ */ new WeakMap()).set($.f = ((cx, cy) => {
-	"use gpu";
-	const idx = __tsover_add(__tsover_mul(cy, u32(N)), cx);
-	const ii = i32(cy);
-	const ji = i32(cx);
-	const jf = f32(cx);
-	const iF = f32(cy);
-	const selfVel = fluidLayout.$.velA[idx];
-	const normalL = fluidLayout.$.velA[nidx(ii, ji, 0, -1)].x;
-	const normalR = fluidLayout.$.velA[nidx(ii, ji, 0, 1)].x;
-	const normalU = fluidLayout.$.velA[nidx(ii, ji, -1, 0)].y;
-	const normalD = fluidLayout.$.velA[nidx(ii, ji, 1, 0)].y;
-	const L = select(normalL, __tsover_mul(selfVel.x, -1), jf < .5);
-	const R = select(normalR, __tsover_mul(selfVel.x, -1), jf > f32(__tsover_sub(N, 1.5)));
-	const T = select(normalU, __tsover_mul(selfVel.y, -1), iF < .5);
-	const B = select(normalD, __tsover_mul(selfVel.y, -1), iF > f32(__tsover_sub(N, 1.5)));
-	const div = __tsover_mul(__tsover_add(__tsover_sub(R, L), __tsover_sub(B, T)), .5);
-	fluidLayout.$.pressure[idx] = __tsover_mul(fluidLayout.$.pressure[idx], PRESSURE_DISSIPATION);
-	fluidLayout.$.divergence[idx] = div;
-}), {
-	v: 1,
-	name: void 0,
-	ast: {
-		"params": [{
-			"type": "i",
-			"name": "cx"
-		}, {
-			"type": "i",
-			"name": "cy"
-		}],
-		"body": [0, [
-			[
-				13,
-				"idx",
-				[
-					1,
-					[
-						1,
-						"cy",
-						"*",
-						[
-							6,
-							[
-								7,
-								"d",
-								"u32"
-							],
-							["N"]
-						]
-					],
-					"+",
-					"cx"
-				]
-			],
-			[
-				13,
-				"ii",
-				[
-					6,
-					[
-						7,
-						"d",
-						"i32"
-					],
-					["cy"]
-				]
-			],
-			[
-				13,
-				"ji",
-				[
-					6,
-					[
-						7,
-						"d",
-						"i32"
-					],
-					["cx"]
-				]
-			],
-			[
-				13,
-				"jf",
-				[
-					6,
-					[
-						7,
-						"d",
-						"f32"
-					],
-					["cx"]
-				]
-			],
-			[
-				13,
-				"iF",
-				[
-					6,
-					[
-						7,
-						"d",
-						"f32"
-					],
-					["cy"]
-				]
-			],
-			[
-				13,
-				"selfVel",
-				[
-					8,
-					[
-						7,
-						[
-							7,
-							"fluidLayout",
-							"$"
-						],
-						"velA"
-					],
-					"idx"
-				]
-			],
-			[
-				13,
-				"normalL",
-				[
-					7,
-					[
-						8,
-						[
-							7,
-							[
-								7,
-								"fluidLayout",
-								"$"
-							],
-							"velA"
-						],
-						[
-							6,
-							"nidx",
-							[
-								"ii",
-								"ji",
-								[5, "0"],
-								[
-									4,
-									"-",
-									[5, "1"]
-								]
-							]
-						]
-					],
-					"x"
-				]
-			],
-			[
-				13,
-				"normalR",
-				[
-					7,
-					[
-						8,
-						[
-							7,
-							[
-								7,
-								"fluidLayout",
-								"$"
-							],
-							"velA"
-						],
-						[
-							6,
-							"nidx",
-							[
-								"ii",
-								"ji",
-								[5, "0"],
-								[5, "1"]
-							]
-						]
-					],
-					"x"
-				]
-			],
-			[
-				13,
-				"normalU",
-				[
-					7,
-					[
-						8,
-						[
-							7,
-							[
-								7,
-								"fluidLayout",
-								"$"
-							],
-							"velA"
-						],
-						[
-							6,
-							"nidx",
-							[
-								"ii",
-								"ji",
-								[
-									4,
-									"-",
-									[5, "1"]
-								],
-								[5, "0"]
-							]
-						]
-					],
-					"y"
-				]
-			],
-			[
-				13,
-				"normalD",
-				[
-					7,
-					[
-						8,
-						[
-							7,
-							[
-								7,
-								"fluidLayout",
-								"$"
-							],
-							"velA"
-						],
-						[
-							6,
-							"nidx",
-							[
-								"ii",
-								"ji",
-								[5, "1"],
-								[5, "0"]
-							]
-						]
-					],
-					"y"
-				]
-			],
-			[
-				13,
-				"L",
-				[
-					6,
-					[
-						7,
-						"std",
-						"select"
-					],
-					[
-						"normalL",
-						[
-							1,
-							[
-								7,
-								"selfVel",
-								"x"
-							],
-							"*",
-							[
-								4,
-								"-",
-								[5, "1"]
-							]
-						],
-						[
-							1,
-							"jf",
-							"<",
-							[5, "0.5"]
-						]
-					]
-				]
-			],
-			[
-				13,
-				"R",
-				[
-					6,
-					[
-						7,
-						"std",
-						"select"
-					],
-					[
-						"normalR",
-						[
-							1,
-							[
-								7,
-								"selfVel",
-								"x"
-							],
-							"*",
-							[
-								4,
-								"-",
-								[5, "1"]
-							]
-						],
-						[
-							1,
-							"jf",
-							">",
-							[
-								6,
-								[
-									7,
-									"d",
-									"f32"
-								],
-								[[
-									1,
-									"N",
-									"-",
-									[5, "1.5"]
-								]]
-							]
-						]
-					]
-				]
-			],
-			[
-				13,
-				"T",
-				[
-					6,
-					[
-						7,
-						"std",
-						"select"
-					],
-					[
-						"normalU",
-						[
-							1,
-							[
-								7,
-								"selfVel",
-								"y"
-							],
-							"*",
-							[
-								4,
-								"-",
-								[5, "1"]
-							]
-						],
-						[
-							1,
-							"iF",
-							"<",
-							[5, "0.5"]
-						]
-					]
-				]
-			],
-			[
-				13,
-				"B",
-				[
-					6,
-					[
-						7,
-						"std",
-						"select"
-					],
-					[
-						"normalD",
-						[
-							1,
-							[
-								7,
-								"selfVel",
-								"y"
-							],
-							"*",
-							[
-								4,
-								"-",
-								[5, "1"]
-							]
-						],
-						[
-							1,
-							"iF",
-							">",
-							[
-								6,
-								[
-									7,
-									"d",
-									"f32"
-								],
-								[[
-									1,
-									"N",
-									"-",
-									[5, "1.5"]
-								]]
-							]
-						]
-					]
-				]
-			],
-			[
-				13,
-				"div",
-				[
-					1,
-					[
-						1,
-						[
-							1,
-							"R",
-							"-",
-							"L"
-						],
-						"+",
-						[
-							1,
-							"B",
-							"-",
-							"T"
-						]
-					],
-					"*",
-					[5, "0.5"]
-				]
-			],
-			[
-				2,
-				[
-					8,
-					[
-						7,
-						[
-							7,
-							"fluidLayout",
-							"$"
-						],
-						"pressure"
-					],
-					"idx"
-				],
-				"=",
-				[
-					1,
-					[
-						8,
-						[
-							7,
-							[
-								7,
-								"fluidLayout",
-								"$"
-							],
-							"pressure"
-						],
-						"idx"
-					],
-					"*",
-					"PRESSURE_DISSIPATION"
-				]
-			],
-			[
-				2,
-				[
-					8,
-					[
-						7,
-						[
-							7,
-							"fluidLayout",
-							"$"
-						],
-						"divergence"
-					],
-					"idx"
-				],
-				"=",
-				"div"
-			]
-		]],
-		"externalNames": [
-			"d",
-			"N",
-			"fluidLayout",
-			"nidx",
-			"std",
-			"PRESSURE_DISSIPATION"
-		]
-	},
-	externals: () => ({
-		d: data_exports,
-		N,
-		fluidLayout,
-		nidx,
-		std: std_exports,
-		PRESSURE_DISSIPATION
-	})
-}) && $.f)({})).$name("particleFlowDivergence");
-const particleFlowJacobiKernel = src_default.fn([u32, u32])(/* @__PURE__ */ (($) => (globalThis.__TYPEGPU_META__ ??= /* @__PURE__ */ new WeakMap()).set($.f = ((cx, cy) => {
-	"use gpu";
-	const idx = __tsover_add(__tsover_mul(cy, u32(N)), cx);
-	const ii = i32(cy);
-	const ji = i32(cx);
-	const pL = fluidLayout.$.pressure[nidx(ii, ji, 0, -1)];
-	const pR = fluidLayout.$.pressure[nidx(ii, ji, 0, 1)];
-	const pU = fluidLayout.$.pressure[nidx(ii, ji, -1, 0)];
-	const pD = fluidLayout.$.pressure[nidx(ii, ji, 1, 0)];
-	const div = fluidLayout.$.divergence[idx];
-	fluidLayout.$.pressure[idx] = __tsover_mul(__tsover_sub(__tsover_add(__tsover_add(__tsover_add(pL, pR), pU), pD), div), .25);
-}), {
-	v: 1,
-	name: void 0,
-	ast: {
-		"params": [{
-			"type": "i",
-			"name": "cx"
-		}, {
-			"type": "i",
-			"name": "cy"
-		}],
-		"body": [0, [
-			[
-				13,
-				"idx",
-				[
-					1,
-					[
-						1,
-						"cy",
-						"*",
-						[
-							6,
-							[
-								7,
-								"d",
-								"u32"
-							],
-							["N"]
-						]
-					],
-					"+",
-					"cx"
-				]
-			],
-			[
-				13,
-				"ii",
-				[
-					6,
-					[
-						7,
-						"d",
-						"i32"
-					],
-					["cy"]
-				]
-			],
-			[
-				13,
-				"ji",
-				[
-					6,
-					[
-						7,
-						"d",
-						"i32"
-					],
-					["cx"]
-				]
-			],
-			[
-				13,
-				"pL",
-				[
-					8,
-					[
-						7,
-						[
-							7,
-							"fluidLayout",
-							"$"
-						],
-						"pressure"
-					],
-					[
-						6,
-						"nidx",
-						[
-							"ii",
-							"ji",
-							[5, "0"],
-							[
-								4,
-								"-",
-								[5, "1"]
-							]
-						]
-					]
-				]
-			],
-			[
-				13,
-				"pR",
-				[
-					8,
-					[
-						7,
-						[
-							7,
-							"fluidLayout",
-							"$"
-						],
-						"pressure"
-					],
-					[
-						6,
-						"nidx",
-						[
-							"ii",
-							"ji",
-							[5, "0"],
-							[5, "1"]
-						]
-					]
-				]
-			],
-			[
-				13,
-				"pU",
-				[
-					8,
-					[
-						7,
-						[
-							7,
-							"fluidLayout",
-							"$"
-						],
-						"pressure"
-					],
-					[
-						6,
-						"nidx",
-						[
-							"ii",
-							"ji",
-							[
-								4,
-								"-",
-								[5, "1"]
-							],
-							[5, "0"]
-						]
-					]
-				]
-			],
-			[
-				13,
-				"pD",
-				[
-					8,
-					[
-						7,
-						[
-							7,
-							"fluidLayout",
-							"$"
-						],
-						"pressure"
-					],
-					[
-						6,
-						"nidx",
-						[
-							"ii",
-							"ji",
-							[5, "1"],
-							[5, "0"]
-						]
-					]
-				]
-			],
-			[
-				13,
-				"div",
-				[
-					8,
-					[
-						7,
-						[
-							7,
-							"fluidLayout",
-							"$"
-						],
-						"divergence"
-					],
-					"idx"
-				]
-			],
-			[
-				2,
-				[
-					8,
-					[
-						7,
-						[
-							7,
-							"fluidLayout",
-							"$"
-						],
-						"pressure"
-					],
-					"idx"
-				],
-				"=",
-				[
-					1,
-					[
-						1,
-						[
-							1,
-							[
-								1,
-								[
-									1,
-									"pL",
-									"+",
-									"pR"
-								],
-								"+",
-								"pU"
-							],
-							"+",
-							"pD"
-						],
-						"-",
-						"div"
-					],
-					"*",
-					[5, "0.25"]
-				]
-			]
-		]],
-		"externalNames": [
-			"d",
-			"N",
-			"fluidLayout",
-			"nidx"
-		]
-	},
-	externals: () => ({
-		d: data_exports,
-		N,
-		fluidLayout,
-		nidx
-	})
-}) && $.f)({})).$name("particleFlowJacobi");
-const particleFlowGradSubtractKernel = src_default.fn([u32, u32])(/* @__PURE__ */ (($) => (globalThis.__TYPEGPU_META__ ??= /* @__PURE__ */ new WeakMap()).set($.f = ((cx, cy) => {
-	"use gpu";
-	const idx = __tsover_add(__tsover_mul(cy, u32(N)), cx);
-	const ii = i32(cy);
-	const ji = i32(cx);
-	const pL = fluidLayout.$.pressure[nidx(ii, ji, 0, -1)];
-	const pR = fluidLayout.$.pressure[nidx(ii, ji, 0, 1)];
-	const pU = fluidLayout.$.pressure[nidx(ii, ji, -1, 0)];
-	const pD = fluidLayout.$.pressure[nidx(ii, ji, 1, 0)];
-	const vel = fluidLayout.$.velA[idx];
-	fluidLayout.$.velA[idx] = vec4f(__tsover_sub(vel.x, __tsover_mul(__tsover_sub(pR, pL), .5)), __tsover_sub(vel.y, __tsover_mul(__tsover_sub(pD, pU), .5)), vel.z, 0);
-}), {
-	v: 1,
-	name: void 0,
-	ast: {
-		"params": [{
-			"type": "i",
-			"name": "cx"
-		}, {
-			"type": "i",
-			"name": "cy"
-		}],
-		"body": [0, [
-			[
-				13,
-				"idx",
-				[
-					1,
-					[
-						1,
-						"cy",
-						"*",
-						[
-							6,
-							[
-								7,
-								"d",
-								"u32"
-							],
-							["N"]
-						]
-					],
-					"+",
-					"cx"
-				]
-			],
-			[
-				13,
-				"ii",
-				[
-					6,
-					[
-						7,
-						"d",
-						"i32"
-					],
-					["cy"]
-				]
-			],
-			[
-				13,
-				"ji",
-				[
-					6,
-					[
-						7,
-						"d",
-						"i32"
-					],
-					["cx"]
-				]
-			],
-			[
-				13,
-				"pL",
-				[
-					8,
-					[
-						7,
-						[
-							7,
-							"fluidLayout",
-							"$"
-						],
-						"pressure"
-					],
-					[
-						6,
-						"nidx",
-						[
-							"ii",
-							"ji",
-							[5, "0"],
-							[
-								4,
-								"-",
-								[5, "1"]
-							]
-						]
-					]
-				]
-			],
-			[
-				13,
-				"pR",
-				[
-					8,
-					[
-						7,
-						[
-							7,
-							"fluidLayout",
-							"$"
-						],
-						"pressure"
-					],
-					[
-						6,
-						"nidx",
-						[
-							"ii",
-							"ji",
-							[5, "0"],
-							[5, "1"]
-						]
-					]
-				]
-			],
-			[
-				13,
-				"pU",
-				[
-					8,
-					[
-						7,
-						[
-							7,
-							"fluidLayout",
-							"$"
-						],
-						"pressure"
-					],
-					[
-						6,
-						"nidx",
-						[
-							"ii",
-							"ji",
-							[
-								4,
-								"-",
-								[5, "1"]
-							],
-							[5, "0"]
-						]
-					]
-				]
-			],
-			[
-				13,
-				"pD",
-				[
-					8,
-					[
-						7,
-						[
-							7,
-							"fluidLayout",
-							"$"
-						],
-						"pressure"
-					],
-					[
-						6,
-						"nidx",
-						[
-							"ii",
-							"ji",
-							[5, "1"],
-							[5, "0"]
-						]
-					]
-				]
-			],
-			[
-				13,
-				"vel",
-				[
-					8,
-					[
-						7,
-						[
-							7,
-							"fluidLayout",
-							"$"
-						],
-						"velA"
-					],
-					"idx"
-				]
-			],
-			[
-				2,
-				[
-					8,
-					[
-						7,
-						[
-							7,
-							"fluidLayout",
-							"$"
-						],
-						"velA"
-					],
-					"idx"
-				],
-				"=",
-				[
-					6,
-					[
-						7,
-						"d",
-						"vec4f"
-					],
-					[
-						[
-							1,
-							[
-								7,
-								"vel",
-								"x"
-							],
-							"-",
-							[
-								1,
-								[
-									1,
-									"pR",
-									"-",
-									"pL"
-								],
-								"*",
-								[5, "0.5"]
-							]
-						],
-						[
-							1,
-							[
-								7,
-								"vel",
-								"y"
-							],
-							"-",
-							[
-								1,
-								[
-									1,
-									"pD",
-									"-",
-									"pU"
-								],
-								"*",
-								[5, "0.5"]
-							]
-						],
-						[
-							7,
-							"vel",
-							"z"
-						],
-						[5, "0"]
-					]
-				]
-			]
-		]],
-		"externalNames": [
-			"d",
-			"N",
-			"fluidLayout",
-			"nidx"
-		]
-	},
-	externals: () => ({
-		d: data_exports,
-		N,
-		fluidLayout,
-		nidx
-	})
-}) && $.f)({})).$name("particleFlowGradSubtract");
-const particleFlowAdvectVelKernel = src_default.fn([u32, u32])(/* @__PURE__ */ (($) => (globalThis.__TYPEGPU_META__ ??= /* @__PURE__ */ new WeakMap()).set($.f = ((cx, cy) => {
-	"use gpu";
-	const p = fluidLayout.$.params;
-	const idx = __tsover_add(__tsover_mul(cy, u32(N)), cx);
-	const vel = fluidLayout.$.velA[idx];
-	const srcX = clamp(__tsover_sub(f32(cx), __tsover_mul(vel.x, p.dt)), 0, f32(__tsover_sub(N, 1)));
-	const srcY = clamp(__tsover_sub(f32(cy), __tsover_mul(vel.y, p.dt)), 0, f32(__tsover_sub(N, 1)));
-	const x0 = clamp(floor(srcX), 0, f32(__tsover_sub(N, 2)));
-	const y0 = clamp(floor(srcY), 0, f32(__tsover_sub(N, 2)));
-	const sx = __tsover_sub(srcX, x0);
-	const sy = __tsover_sub(srcY, y0);
-	const ix0 = i32(x0);
-	const iy0 = i32(y0);
-	const v00 = fluidLayout.$.velA[u32(__tsover_add(__tsover_mul(iy0, N), ix0))];
-	const v10 = fluidLayout.$.velA[u32(__tsover_add(__tsover_add(__tsover_mul(iy0, N), ix0), 1))];
-	const v01 = fluidLayout.$.velA[u32(__tsover_add(__tsover_mul(__tsover_add(iy0, 1), N), ix0))];
-	const v11 = fluidLayout.$.velA[u32(__tsover_add(__tsover_add(__tsover_mul(__tsover_add(iy0, 1), N), ix0), 1))];
-	const decay = __tsover_add(1, __tsover_mul(p.velFade, p.dt));
-	fluidLayout.$.velB[idx] = vec4f(__tsover_div(mix(mix(v00.x, v10.x, sx), mix(v01.x, v11.x, sx), sy), decay), __tsover_div(mix(mix(v00.y, v10.y, sx), mix(v01.y, v11.y, sx), sy), decay), 0, 0);
-}), {
-	v: 1,
-	name: void 0,
-	ast: {
-		"params": [{
-			"type": "i",
-			"name": "cx"
-		}, {
-			"type": "i",
-			"name": "cy"
-		}],
-		"body": [0, [
-			[
-				13,
-				"p",
-				[
-					7,
-					[
-						7,
-						"fluidLayout",
-						"$"
-					],
-					"params"
-				]
-			],
-			[
-				13,
-				"idx",
-				[
-					1,
-					[
-						1,
-						"cy",
-						"*",
-						[
-							6,
-							[
-								7,
-								"d",
-								"u32"
-							],
-							["N"]
-						]
-					],
-					"+",
-					"cx"
-				]
-			],
-			[
-				13,
-				"vel",
-				[
-					8,
-					[
-						7,
-						[
-							7,
-							"fluidLayout",
-							"$"
-						],
-						"velA"
-					],
-					"idx"
-				]
-			],
-			[
-				13,
-				"srcX",
-				[
-					6,
-					[
-						7,
-						"std",
-						"clamp"
-					],
-					[
-						[
-							1,
-							[
-								6,
-								[
-									7,
-									"d",
-									"f32"
-								],
-								["cx"]
-							],
-							"-",
-							[
-								1,
-								[
-									7,
-									"vel",
-									"x"
-								],
-								"*",
-								[
-									7,
-									"p",
-									"dt"
-								]
-							]
-						],
-						[5, "0"],
-						[
-							6,
-							[
-								7,
-								"d",
-								"f32"
-							],
-							[[
-								1,
-								"N",
-								"-",
-								[5, "1"]
-							]]
-						]
-					]
-				]
-			],
-			[
-				13,
-				"srcY",
-				[
-					6,
-					[
-						7,
-						"std",
-						"clamp"
-					],
-					[
-						[
-							1,
-							[
-								6,
-								[
-									7,
-									"d",
-									"f32"
-								],
-								["cy"]
-							],
-							"-",
-							[
-								1,
-								[
-									7,
-									"vel",
-									"y"
-								],
-								"*",
-								[
-									7,
-									"p",
-									"dt"
-								]
-							]
-						],
-						[5, "0"],
-						[
-							6,
-							[
-								7,
-								"d",
-								"f32"
-							],
-							[[
-								1,
-								"N",
-								"-",
-								[5, "1"]
-							]]
-						]
-					]
-				]
-			],
-			[
-				13,
-				"x0",
-				[
-					6,
-					[
-						7,
-						"std",
-						"clamp"
-					],
-					[
-						[
-							6,
-							[
-								7,
-								"std",
-								"floor"
-							],
-							["srcX"]
-						],
-						[5, "0"],
-						[
-							6,
-							[
-								7,
-								"d",
-								"f32"
-							],
-							[[
-								1,
-								"N",
-								"-",
-								[5, "2"]
-							]]
-						]
-					]
-				]
-			],
-			[
-				13,
-				"y0",
-				[
-					6,
-					[
-						7,
-						"std",
-						"clamp"
-					],
-					[
-						[
-							6,
-							[
-								7,
-								"std",
-								"floor"
-							],
-							["srcY"]
-						],
-						[5, "0"],
-						[
-							6,
-							[
-								7,
-								"d",
-								"f32"
-							],
-							[[
-								1,
-								"N",
-								"-",
-								[5, "2"]
-							]]
-						]
-					]
-				]
-			],
-			[
-				13,
-				"sx",
-				[
-					1,
-					"srcX",
-					"-",
-					"x0"
-				]
-			],
-			[
-				13,
-				"sy",
-				[
-					1,
-					"srcY",
-					"-",
-					"y0"
-				]
-			],
-			[
-				13,
-				"ix0",
-				[
-					6,
-					[
-						7,
-						"d",
-						"i32"
-					],
-					["x0"]
-				]
-			],
-			[
-				13,
-				"iy0",
-				[
-					6,
-					[
-						7,
-						"d",
-						"i32"
-					],
-					["y0"]
-				]
-			],
-			[
-				13,
-				"v00",
-				[
-					8,
-					[
-						7,
-						[
-							7,
-							"fluidLayout",
-							"$"
-						],
-						"velA"
-					],
-					[
-						6,
-						[
-							7,
-							"d",
-							"u32"
-						],
-						[[
-							1,
-							[
-								1,
-								"iy0",
-								"*",
-								"N"
-							],
-							"+",
-							"ix0"
-						]]
-					]
-				]
-			],
-			[
-				13,
-				"v10",
-				[
-					8,
-					[
-						7,
-						[
-							7,
-							"fluidLayout",
-							"$"
-						],
-						"velA"
-					],
-					[
-						6,
-						[
-							7,
-							"d",
-							"u32"
-						],
-						[[
-							1,
-							[
-								1,
-								[
-									1,
-									"iy0",
-									"*",
-									"N"
-								],
-								"+",
-								"ix0"
-							],
-							"+",
-							[5, "1"]
-						]]
-					]
-				]
-			],
-			[
-				13,
-				"v01",
-				[
-					8,
-					[
-						7,
-						[
-							7,
-							"fluidLayout",
-							"$"
-						],
-						"velA"
-					],
-					[
-						6,
-						[
-							7,
-							"d",
-							"u32"
-						],
-						[[
-							1,
-							[
-								1,
-								[
-									1,
-									"iy0",
-									"+",
-									[5, "1"]
-								],
-								"*",
-								"N"
-							],
-							"+",
-							"ix0"
-						]]
-					]
-				]
-			],
-			[
-				13,
-				"v11",
-				[
-					8,
-					[
-						7,
-						[
-							7,
-							"fluidLayout",
-							"$"
-						],
-						"velA"
-					],
-					[
-						6,
-						[
-							7,
-							"d",
-							"u32"
-						],
-						[[
-							1,
-							[
-								1,
-								[
-									1,
-									[
-										1,
-										"iy0",
-										"+",
-										[5, "1"]
-									],
-									"*",
-									"N"
-								],
-								"+",
-								"ix0"
-							],
-							"+",
-							[5, "1"]
-						]]
-					]
-				]
-			],
-			[
-				13,
-				"decay",
-				[
-					1,
-					[5, "1"],
-					"+",
-					[
-						1,
-						[
-							7,
-							"p",
-							"velFade"
-						],
-						"*",
-						[
-							7,
-							"p",
-							"dt"
-						]
-					]
-				]
-			],
-			[
-				2,
-				[
-					8,
-					[
-						7,
-						[
-							7,
-							"fluidLayout",
-							"$"
-						],
-						"velB"
-					],
-					"idx"
-				],
-				"=",
-				[
-					6,
-					[
-						7,
-						"d",
-						"vec4f"
-					],
-					[
+						13,
+						"af",
 						[
 							1,
 							[
-								6,
+								1,
 								[
 									7,
-									"std",
-									"mix"
+									"p",
+									"ambient"
 								],
-								[
-									[
-										6,
-										[
-											7,
-											"std",
-											"mix"
-										],
-										[
-											[
-												7,
-												"v00",
-												"x"
-											],
-											[
-												7,
-												"v10",
-												"x"
-											],
-											"sx"
-										]
-									],
-									[
-										6,
-										[
-											7,
-											"std",
-											"mix"
-										],
-										[
-											[
-												7,
-												"v01",
-												"x"
-											],
-											[
-												7,
-												"v11",
-												"x"
-											],
-											"sx"
-										]
-									],
-									"sy"
-								]
+								"*",
+								"AMBIENT_K"
 							],
-							"/",
-							"decay"
-						],
+							"*",
+							[
+								7,
+								"p",
+								"dt"
+							]
+						]
+					],
+					[
+						2,
+						"ax",
+						"=",
 						[
 							1,
 							[
-								6,
-								[
-									7,
-									"std",
-									"mix"
-								],
-								[
-									[
-										6,
-										[
-											7,
-											"std",
-											"mix"
-										],
-										[
-											[
-												7,
-												"v00",
-												"y"
-											],
-											[
-												7,
-												"v10",
-												"y"
-											],
-											"sx"
-										]
-									],
-									[
-										6,
-										[
-											7,
-											"std",
-											"mix"
-										],
-										[
-											[
-												7,
-												"v01",
-												"y"
-											],
-											[
-												7,
-												"v11",
-												"y"
-											],
-											"sx"
-										]
-									],
-									"sy"
-								]
+								7,
+								"amb",
+								"x"
 							],
-							"/",
-							"decay"
-						],
-						[5, "0"],
-						[5, "0"]
-					]
-				]
-			]
-		]],
-		"externalNames": [
-			"fluidLayout",
-			"d",
-			"N",
-			"std"
-		]
-	},
-	externals: () => ({
-		fluidLayout,
-		d: data_exports,
-		N,
-		std: std_exports
-	})
-}) && $.f)({})).$name("particleFlowAdvectVel");
-const particleFlowCopyVelKernel = src_default.fn([u32, u32])(/* @__PURE__ */ (($) => (globalThis.__TYPEGPU_META__ ??= /* @__PURE__ */ new WeakMap()).set($.f = ((cx, cy) => {
-	"use gpu";
-	const idx = __tsover_add(__tsover_mul(cy, u32(N)), cx);
-	const v = fluidLayout.$.velB[idx];
-	fluidLayout.$.velA[idx] = vec4f(v.x, v.y, 0, 0);
-	textureStore(fluidLayout.$.velOutTex, vec2u(cx, cy), vec4f(v.x, v.y, 0, 0));
-}), {
-	v: 1,
-	name: void 0,
-	ast: {
-		"params": [{
-			"type": "i",
-			"name": "cx"
-		}, {
-			"type": "i",
-			"name": "cy"
-		}],
-		"body": [0, [
-			[
-				13,
-				"idx",
-				[
-					1,
+							"*",
+							"af"
+						]
+					],
 					[
-						1,
-						"cy",
-						"*",
+						2,
+						"ay",
+						"=",
 						[
-							6,
+							1,
 							[
 								7,
-								"d",
-								"u32"
+								"amb",
+								"y"
 							],
-							["N"]
+							"*",
+							"af"
 						]
-					],
-					"+",
-					"cx"
-				]
+					]
+				]]
 			],
 			[
 				13,
-				"v",
+				"vel",
 				[
 					8,
 					[
@@ -3330,7 +1040,7 @@ const particleFlowCopyVelKernel = src_default.fn([u32, u32])(/* @__PURE__ */ (($
 							"fluidLayout",
 							"$"
 						],
-						"velB"
+						"velA"
 					],
 					"idx"
 				]
@@ -3360,85 +1070,65 @@ const particleFlowCopyVelKernel = src_default.fn([u32, u32])(/* @__PURE__ */ (($
 					],
 					[
 						[
-							7,
-							"v",
-							"x"
-						],
-						[
-							7,
-							"v",
-							"y"
-						],
-						[5, "0"],
-						[5, "0"]
-					]
-				]
-			],
-			[
-				6,
-				[
-					7,
-					"std",
-					"textureStore"
-				],
-				[
-					[
-						7,
-						[
-							7,
-							"fluidLayout",
-							"$"
+							1,
+							[
+								1,
+								[
+									7,
+									"vel",
+									"x"
+								],
+								"+",
+								"fx"
+							],
+							"+",
+							"ax"
 						],
-						"velOutTex"
-					],
-					[
-						6,
 						[
-							7,
-							"d",
-							"vec2u"
+							1,
+							[
+								1,
+								[
+									7,
+									"vel",
+									"y"
+								],
+								"+",
+								"fy"
+							],
+							"+",
+							"ay"
 						],
-						["cx", "cy"]
-					],
-					[
-						6,
 						[
 							7,
-							"d",
-							"vec4f"
+							"vel",
+							"z"
 						],
-						[
-							[
-								7,
-								"v",
-								"x"
-							],
-							[
-								7,
-								"v",
-								"y"
-							],
-							[5, "0"],
-							[5, "0"]
-						]
+						[5, "0"]
 					]
 				]
 			]
 		]],
 		"externalNames": [
+			"fluidLayout",
 			"d",
 			"N",
-			"fluidLayout",
-			"std"
+			"std",
+			"nidx",
+			"noise",
+			"AMBIENT_K"
 		]
 	},
 	externals: () => ({
+		fluidLayout,
 		d: data_exports,
 		N,
-		fluidLayout,
-		std: std_exports
+		std: std_exports,
+		nidx,
+		noise: noise_exports,
+		AMBIENT_K
 	})
-}) && $.f)({})).$name("particleFlowCopyVel");
+}) && $.f)({})).$name("particleFlowVorticity");
 const particleFlowInitKernel = src_default.fn([u32])(/* @__PURE__ */ (($) => (globalThis.__TYPEGPU_META__ ??= /* @__PURE__ */ new WeakMap()).set($.f = ((i) => {
 	"use gpu";
 	const prm = pfLayout.$.params;
@@ -6543,9 +4233,9 @@ const componentDefinition = {
 			outTex
 		});
 		const fluidSize = [N, N];
-		const mkFluid = (kernel) => createGuardedCompute(root, /* @__PURE__ */ (($) => (globalThis.__TYPEGPU_META__ ??= /* @__PURE__ */ new WeakMap()).set($.f = ((cx, cy) => {
+		const force = createGuardedCompute(root, /* @__PURE__ */ (($) => (globalThis.__TYPEGPU_META__ ??= /* @__PURE__ */ new WeakMap()).set($.f = ((cx, cy) => {
 			"use gpu";
-			kernel(cx, cy);
+			particleFlowForceKernel(cx, cy);
 		}), {
 			v: 1,
 			name: void 0,
@@ -6559,19 +4249,19 @@ const componentDefinition = {
 				}],
 				"body": [0, [[
 					6,
-					"kernel",
+					"particleFlowForceKernel",
 					["cx", "cy"]
 				]]],
-				"externalNames": ["kernel"]
+				"externalNames": ["particleFlowForceKernel"]
 			},
-			externals: () => ({ kernel })
+			externals: () => ({ particleFlowForceKernel })
 		}) && $.f)({}), {
 			size: fluidSize,
 			bindGroup: fluidBg
 		});
-		const force = createGuardedCompute(root, /* @__PURE__ */ (($) => (globalThis.__TYPEGPU_META__ ??= /* @__PURE__ */ new WeakMap()).set($.f = ((cx, cy) => {
+		const vorticity = createGuardedCompute(root, /* @__PURE__ */ (($) => (globalThis.__TYPEGPU_META__ ??= /* @__PURE__ */ new WeakMap()).set($.f = ((cx, cy) => {
 			"use gpu";
-			particleFlowForceKernel(cx, cy);
+			particleFlowVorticityKernel(cx, cy);
 		}), {
 			v: 1,
 			name: void 0,
@@ -6585,23 +4275,20 @@ const componentDefinition = {
 				}],
 				"body": [0, [[
 					6,
-					"particleFlowForceKernel",
+					"particleFlowVorticityKernel",
 					["cx", "cy"]
 				]]],
-				"externalNames": ["particleFlowForceKernel"]
+				"externalNames": ["particleFlowVorticityKernel"]
 			},
-			externals: () => ({ particleFlowForceKernel })
+			externals: () => ({ particleFlowVorticityKernel })
 		}) && $.f)({}), {
 			size: fluidSize,
 			bindGroup: fluidBg
 		});
-		const curl = mkFluid(particleFlowCurlKernel);
-		const vorticity = mkFluid(particleFlowVorticityKernel);
-		const divergencePass = mkFluid(particleFlowDivergenceKernel);
-		const jacobi = mkFluid(particleFlowJacobiKernel);
-		const gradSubtract = mkFluid(particleFlowGradSubtractKernel);
-		const advectVel = mkFluid(particleFlowAdvectVelKernel);
-		const copyVel = mkFluid(particleFlowCopyVelKernel);
+		const solver = createStableFluidsPasses(root, solverKernels, {
+			n: N,
+			bindGroup: fluidBg
+		});
 		const trailsCpu = getCpuValue("trails");
 		const splatKernel = makeParticleFlowSplatKernel(getCpuValue("shape") || "streak");
 		const resolveKernel = makeParticleFlowResolveKernel(getCpuValue("colorSpace") ?? 2, typeof trailsCpu === "number" && trailsCpu > 0);
@@ -6689,12 +4376,10 @@ const componentDefinition = {
 			externals: () => ({ resolveKernel })
 		}) && $.f)({}), { size: [RES, RES] }).with(pfBg);
 		let lastTime = Date.now();
-		let prevMx = .5;
-		let prevMy = .5;
-		let lastActiveAt = 0;
 		let ambientTime = 0;
-		let warmup = 0;
 		let initialized = false;
+		const pointer = createPointerVelocityTracker();
+		const idle = createIdleGate({ warmupFrames: WARMUP_FRAMES + 1 });
 		const countCap = isMobileGpuViewport() ? COUNT_MAX_MOBILE : MAX_PARTICLES;
 		const num = (key, fallback) => {
 			const v = getCpuValue(key);
@@ -6718,23 +4403,12 @@ const componentDefinition = {
 				const advect = num("speed", 1);
 				const count = Math.min(Math.max(Math.round(num("count", 6e3)), 100), countCap);
 				const bodyR = Math.min(Math.max(num("size", 1.2), .3), 3) * 1.3 / SIZE_REF_RES;
-				const oldMx = prevMx;
-				const oldMy = prevMy;
-				const mx = fp.pointer?.x ?? prevMx;
-				const my = fp.pointer?.y ?? prevMy;
-				const dxp = mx - oldMx;
-				const dyp = my - oldMy;
-				prevMx = mx;
-				prevMy = my;
-				const dragDist = Math.sqrt(dxp * dxp + dyp * dyp);
-				const moving = dragDist > 6e-4 && dragDist < .25;
-				const active = moving || ambient > .001;
-				if (active) lastActiveAt = now;
+				const ptr = pointer.update(fp.pointer, dt);
+				const active = ptr.moving || ambient > .001;
+				if (active) idle.markActive(now);
 				ambientTime += dt * AMBIENT_DRIFT;
-				const idleFadeSecs = 2 + momentum * 13;
-				const settled = lastActiveAt !== 0 && (now - lastActiveAt) / 1e3 > idleFadeSecs;
-				if (initialized && warmup > WARMUP_FRAMES && !active && (settled || lastActiveAt === 0)) return null;
-				warmup++;
+				if (initialized && idle.shouldSkip(now, 2 + momentum * 13, active)) return null;
+				idle.tickFrame();
 				const velFade = VEL_FADE_MAX * Math.pow(VEL_FADE_MIN / VEL_FADE_MAX, momentum);
 				paramsU.write({
 					dt,
@@ -6761,28 +4435,32 @@ const componentDefinition = {
 					initialized = true;
 					nodes.push(() => init.dispatchThreads(MAX_PARTICLES));
 				}
-				if (moving) {
-					const velX = dxp / dt * N * forceVal * IMPULSE_K;
-					const velY = dyp / dt * N * forceVal * IMPULSE_K;
-					const stepSize = Math.max(.004, FORCE_RADIUS_GRID / N * .6);
-					const numSteps = Math.min(MAX_STEPS, Math.max(1, Math.ceil(dragDist / stepSize)));
-					for (let sIdx = 0; sIdx < numSteps; sIdx++) {
-						const t = (sIdx + .5) / numSteps;
-						const posX = (oldMx + dxp * t) * N;
-						const posY = (oldMy + dyp * t) * N;
-						nodes.push(() => splatU.write({
+				if (ptr.moving) {
+					const velX = ptr.velX * N * forceVal * IMPULSE_K;
+					const velY = ptr.velY * N * forceVal * IMPULSE_K;
+					pathStampRibbon(nodes, {
+						fromX: ptr.prevX,
+						fromY: ptr.prevY,
+						dx: ptr.dx,
+						dy: ptr.dy,
+						dragDist: ptr.dragDist,
+						stepSize: Math.max(.004, FORCE_RADIUS_GRID / N * .6),
+						maxSteps: MAX_STEPS,
+						scale: N,
+						write: (posX, posY) => splatU.write({
 							posX,
 							posY,
 							velX,
 							velY,
 							radius: FORCE_RADIUS_GRID
-						}));
-						nodes.push(force);
-					}
+						}),
+						pass: force
+					});
 				}
-				nodes.push(curl, vorticity, divergencePass);
-				for (let iter = 0; iter < JACOBI_ITERS; iter++) nodes.push(jacobi);
-				nodes.push(gradSubtract, advectVel, copyVel);
+				nodes.push(...solver.solveSteps({
+					jacobiIters: JACOBI_ITERS,
+					vorticity
+				}));
 				nodes.push(() => integrate.dispatchThreads(count), () => splat.dispatchThreads(count), resolve);
 				return nodes;
 			}