@mkbabb/glass-ui 4.1.0 → 4.2.0

package/dist/{uniformBridgeWGPU-CtAKl45L.js → uniformBridgeWGPU-BCse1SHU.js}

@@ -1,5 +1,5 @@
-import { a as e, n as t, o as n, r, s as i, t as a } from "./procedural-color.wgsl-m7Zgv35j.js";
-import { n as o, t as s } from "./prng-B-d9eYX8.js";
+import { a as e, n as t, o as n, r, s as i, t as a } from "./procedural-color.wgsl-C1Y5PkM6.js";
+import { n as o, t as s } from "./prng-B5tIXTtj.js";
 import { onUnmounted as c, readonly as l, ref as u, watch as d } from "vue";
 import { SpringProgress as f } from "@mkbabb/keyframes.js";
 //#region src/components/custom/goo-blob/types.ts
@@ -59,7 +59,8 @@ var p = {
 		iridSpeed: .06,
 		sssScale: .1,
 		sssPower: 2,
-		coreGlow: .06
+		coreGlow: .06,
+		fissionAmp: 0
 	},
 	interaction: {
 		pointerAttraction: .35,
@@ -120,21 +121,34 @@ var y = {
 	curious: 800,
 	sleepy: 2500,
 	excited: 600
-}, x = 6e3, S = .35, C = .85, w = .2, T = .015, E = .02, D = .01, O = .015, k = .75, A = 2e3, j = /* @__PURE__ */ "uResolution.uTime.uBaseColor.uPointer.uPointerActive.uPointerAttraction.uPointerStrength.uBodyRadius.uPulsePhase.uPulseAmp.uNoiseAmp.uNoiseFreq.uNoiseSpeed.uWarpAmp.uSmoothK.uMerge.uStage.uMaxReach.uLit.uRimColor.uLightDir.uSpecStrength.uSpecShininess.uRimPower.uRimStrength.uShadow.uShadowSoftness.uIridescence.uIridHue.uIridSpeed.uSssScale.uSssPower.uCoreGlow.uHueRange.uSatShift.uBrightnessShift.uColorNoiseFreq.uColorNoiseSpeed.uStopCount.uSatCount.uTrailCount.uVelocity.uStretch".split(".");
+}, x = 6e3, S = .35, C = .85, w = .2, T = .015, E = .02, D = .01, O = .015, k = .75, A = 2e3, j = .4, M = .18, N = 5200, P = /* @__PURE__ */ "uResolution.uTime.uBaseColor.uPointer.uPointerActive.uPointerAttraction.uPointerStrength.uBodyRadius.uPulsePhase.uPulseAmp.uNoiseAmp.uNoiseFreq.uNoiseSpeed.uWarpAmp.uSmoothK.uMerge.uStage.uMorphT.uMaxReach.uLit.uRimColor.uLightDir.uSpecStrength.uSpecShininess.uRimPower.uRimStrength.uShadow.uShadowSoftness.uIridescence.uIridHue.uIridSpeed.uSssScale.uSssPower.uCoreGlow.uHueRange.uSatShift.uBrightnessShift.uColorNoiseFreq.uColorNoiseSpeed.uStopCount.uSatCount.uTrailCount.uVelocity.uStretch".split(".");
 //#endregion
 //#region src/components/custom/goo-blob/composables/easing.ts
-function M(e) {
+function F(e) {
 	return e < .5 ? 2 * e * e : 1 - 2 * (1 - e) * (1 - e);
 }
-function N(e) {
+function I(e) {
 	return e * e;
 }
-function P(e) {
+function L(e) {
 	return 1 - (1 - e) * (1 - e);
 }
+function R(e) {
+	let t = e < 0 ? 0 : e > 1 ? 1 : e;
+	if (t < .3) {
+		let e = t / .3;
+		return -.06 * Math.sin(e * Math.PI);
+	}
+	if (t < .55) {
+		let e = (t - .3) / .25;
+		return (1 - (1 - e) * (1 - e) * (1 - e)) * 1.12;
+	}
+	let n = (t - .55) / .45;
+	return 1 + Math.cos(n * Math.PI) * .12 * (1 - n);
+}
 //#endregion
 //#region src/components/custom/goo-blob/composables/useBlobMood.ts
-function F(e, t, n) {
+function z(e, t, n) {
 	let r = (e, t) => e + (t - e) * n;
 	return {
 		orbitSpeedScale: r(e.orbitSpeedScale, t.orbitSpeedScale),
@@ -151,7 +165,7 @@ function F(e, t, n) {
 		iridScale: r(e.iridScale, t.iridScale)
 	};
 }
-function I() {
+function B() {
 	let e = u("idle"), t = u({ ...y.idle }), n = { ...y.idle }, r = { ...y.idle }, i = 0, a = 0, o = !1, s = 0, c = !1, d = !1, f = 0;
 	function p(s, l = {}) {
 		(l.source ?? "manual") === "manual" && (c = !0), !(s === e.value && !o) && (e.value = s, n = { ...t.value }, r = y[s], a = b[s], i = 0, o = !0);
@@ -176,8 +190,8 @@ function I() {
 	function _(e) {
 		if (s > 0 && (s = Math.max(0, s - e)), !o) return;
 		i += e;
-		let c = Math.min(i / a, 1), l = M(c);
-		t.value = F(n, r, l), c >= 1 && (o = !1);
+		let c = Math.min(i / a, 1), l = F(c);
+		t.value = z(n, r, l), c >= 1 && (o = !1);
 	}
 	return {
 		currentMood: l(e),
@@ -191,7 +205,7 @@ function I() {
 }
 //#endregion
 //#region src/components/custom/goo-blob/composables/useBlobPointer.ts
-function L(e) {
+function V(e) {
 	let t = u({
 		x: 0,
 		y: 0
@@ -291,16 +305,16 @@ function L(e) {
 }
 //#endregion
 //#region src/components/custom/goo-blob/composables/useBlobSatellites.ts
-function R(e, t, n) {
+function H(e, t, n) {
 	return t + e() * (n - t);
 }
-function z(e) {
+function U(e) {
 	return e < 0 ? 0 : e > 1 ? 1 : e;
 }
-function B(e, t, n) {
+function W(e, t, n) {
 	return e + (t - e) * n;
 }
-function V(e, t, n, r, i) {
+function G(e, t, n, r, i) {
 	let a = t / 4 * Math.PI * 2 + (e() - .5) * Math.PI, o = n * (C + e() * w), s = r * (.3 + e() * .7);
 	return {
 		phase: "orbiting",
@@ -324,27 +338,29 @@ function V(e, t, n, r, i) {
 		startX: 0,
 		startY: 0,
 		endX: 0,
-		endY: 0
+		endY: 0,
+		fissioning: !1,
+		snapFired: !1
 	};
 }
-function H(e, t, n, r) {
+function K(e, t, n, r) {
 	e.phase = t, e.phaseStart = n, e.phaseDuration = r;
 }
-function U(e, t, n, r) {
+function q(e, t, n, r) {
 	let i = (t - e.timeOrigin) / 1e3, a = e.angularSpeed * n * i + e.phaseOffset, o = r * (e.wobbleAmp1 * Math.sin(e.wobbleFreq1 * i) + e.wobbleAmp2 * Math.sin(e.wobbleFreq2 * i + 1.3));
 	return {
 		x: (e.baseRadiusX + o) * Math.cos(a) + e.pertXAmp * Math.sin(e.pertXFreq * i + e.pertXPhase),
 		y: (e.baseRadiusY + o) * Math.sin(a) + e.pertYAmp * Math.cos(e.pertYFreq * i + e.pertYPhase)
 	};
 }
-function W(e, t, n, r, i) {
+function J(e, t, n, r, i) {
 	e.timeOrigin = i, e.angularSpeed = .08 + t() * .16, e.phaseOffset = t() * Math.PI * 2;
 	let a = n * (C + t() * w), o = r * (.3 + t() * .7);
 	e.baseRadiusX = a * (1 - o), e.baseRadiusY = a * (1 + o), e.wobbleAmp1 = T + t() * E, e.wobbleFreq1 = .08 + t() * .22, e.wobbleAmp2 = D + t() * O, e.wobbleFreq2 = .18 + t() * .32, e.pertXAmp = .03 + t() * .05, e.pertXFreq = .06 + t() * .18, e.pertXPhase = t() * Math.PI * 2, e.pertYAmp = .03 + t() * .05, e.pertYFreq = .07 + t() * .17, e.pertYPhase = t() * Math.PI * 2;
 }
-function G(e, t) {
-	let n = o(s(t + "goo")), r = [], i = [], a = [], c = -Infinity, l = e.geometry.orbitRadius;
-	function u() {
+function ee(e, t) {
+	let n = o(s(t + "goo")), r = [], i = [], a = [], c = -Infinity, l = e.geometry.orbitRadius, u = -1, d = null;
+	function f() {
 		let t = e.geometry.orbitRadius;
 		if (t !== l && l > 0) {
 			let e = t / l;
@@ -352,13 +368,14 @@ function G(e, t) {
 			l = t;
 		}
 	}
-	function d(t) {
+	function p(t) {
 		let o = e.geometry.satelliteCount;
-		for (; r.length < o;) r.push(V(n, r.length, e.geometry.orbitRadius, e.geometry.eccentricity, t)), i.push({
+		for (; r.length < o;) r.push(G(n, r.length, e.geometry.orbitRadius, e.geometry.eccentricity, t)), i.push({
 			x: 0,
 			y: 0,
 			radius: e.geometry.satelliteRadius,
-			opacity: 0
+			opacity: 0,
+			fissioning: !1
 		}), a.push({
 			x: 0,
 			y: 0,
@@ -366,84 +383,96 @@ function G(e, t) {
 		});
 		r.length > o && (r.length = o, i.length = o, a.length = o);
 	}
-	function f(t, o) {
-		d(t), u();
-		let s = r.length, l = o.mergeRate, f = e.geometry.satelliteRadius, p = e.satellites.mergeDuration, m = e.satellites.emergeDuration, h = e.satellites.orbitDuration, g = e.satellites.absorbedDuration;
-		for (let u = 0; u < s; u++) {
-			let s = r[u], d = t - s.phaseStart, _ = z(s.phaseDuration > 0 ? d / s.phaseDuration : 1);
+	function m(t, o) {
+		p(t), f();
+		let s = r.length, l = o.mergeRate, m = e.geometry.satelliteRadius, h = e.satellites.mergeDuration, g = e.satellites.emergeDuration, _ = e.satellites.orbitDuration, v = e.satellites.absorbedDuration, y = e.surface.fissionAmp ?? 0, b = Math.max(0, Math.min(1, (2 - l) / 1.7)), x = Math.max(0, Math.min(1, y * (.4 + .6 * b))), S = j - M;
+		for (let f = 0; f < s; f++) {
+			let s = r[f], p = t - s.phaseStart, y = U(s.phaseDuration > 0 ? p / s.phaseDuration : 1);
 			switch (s.phase) {
 				case "orbiting":
-					if (_ >= 1) if (t - c < 3e3 * l) s.phaseStart = t, s.phaseDuration = R(n, h[0], h[1]);
+					if (y >= 1) if (t - c < 3e3 * l) s.phaseStart = t, s.phaseDuration = H(n, _[0], _[1]);
+					else if (x > 0 && u < 0 && n() < x * .75) s.phaseOffset = s.angularSpeed * o.orbitSpeedScale * ((t - s.timeOrigin) / 1e3) + s.phaseOffset, s.timeOrigin = t, s.fissioning = !0, s.snapFired = !1, u = f, K(s, "fissioning", t, N), c = t;
 					else {
-						let e = U(s, t, o.orbitSpeedScale, o.wobbleScale);
+						let e = q(s, t, o.orbitSpeedScale, o.wobbleScale);
 						s.startX = e.x, s.startY = e.y;
 						let n = Math.hypot(e.x, e.y), r = n > .01 ? .08 / n : 0;
-						s.endX = e.x * r, s.endY = e.y * r, H(s, "merging", t, p), c = t;
+						s.endX = e.x * r, s.endY = e.y * r, K(s, "merging", t, h), c = t;
 					}
 					break;
+				case "fissioning":
+					y >= 1 && (s.fissioning = !1, u === f && (u = -1), s.timeOrigin = t, K(s, "orbiting", t, H(n, _[0], _[1])));
+					break;
 				case "merging":
-					_ >= 1 && H(s, "absorbed", t, R(n, g[0], g[1]));
+					y >= 1 && K(s, "absorbed", t, H(n, v[0], v[1]));
 					break;
 				case "absorbed":
-					if (_ >= 1) {
-						W(s, n, e.geometry.orbitRadius, e.geometry.eccentricity, t);
-						let r = U(s, t + 3e3, o.orbitSpeedScale, o.wobbleScale);
+					if (y >= 1) {
+						J(s, n, e.geometry.orbitRadius, e.geometry.eccentricity, t);
+						let r = q(s, t + 3e3, o.orbitSpeedScale, o.wobbleScale);
 						s.endX = r.x, s.endY = r.y;
 						let i = Math.hypot(r.x, r.y), a = i > .01 ? .08 / i : 0;
-						s.startX = r.x * a, s.startY = r.y * a, H(s, "emerging", t, m);
+						s.startX = r.x * a, s.startY = r.y * a, K(s, "emerging", t, g);
 					}
 					break;
 				case "emerging":
-					if (_ >= 1) {
-						let e = a[u];
-						e.x = s.endX, e.y = s.endY, e.start = t, s.timeOrigin = t, s.phaseOffset = Math.atan2(s.endY, s.endX), H(s, "orbiting", t, R(n, h[0], h[1]));
+					if (y >= 1) {
+						let e = a[f];
+						e.x = s.endX, e.y = s.endY, e.start = t, s.timeOrigin = t, s.phaseOffset = Math.atan2(s.endY, s.endX), K(s, "orbiting", t, H(n, _[0], _[1]));
 					}
 					break;
 			}
-			let v, y, b, x;
+			let b, C, w, T;
 			switch (s.phase) {
 				case "orbiting": {
-					let e = U(s, t, o.orbitSpeedScale, o.wobbleScale), n = a[u], r = t - n.start;
+					let e = q(s, t, o.orbitSpeedScale, o.wobbleScale), n = a[f], r = t - n.start;
 					if (n.start > 0 && r < 2e3) {
-						let t = z(r / A), i = t * t * (3 - 2 * t);
-						v = B(n.x, e.x, i), y = B(n.y, e.y, i);
-					} else v = e.x, y = e.y;
-					b = k, x = 1;
+						let t = U(r / A), i = t * t * (3 - 2 * t);
+						b = W(n.x, e.x, i), C = W(n.y, e.y, i);
+					} else b = e.x, C = e.y;
+					w = k, T = 1;
 					break;
 				}
 				case "merging": {
-					let e = N(z((t - s.phaseStart) / s.phaseDuration));
-					v = B(s.startX, s.endX, e), y = B(s.startY, s.endY, e);
-					let n = z((t - s.phaseStart) / s.phaseDuration);
-					b = n > .85 ? B(k, 0, (n - .85) / .15) : k, x = B(1, .65, e);
+					let e = I(U((t - s.phaseStart) / s.phaseDuration));
+					b = W(s.startX, s.endX, e), C = W(s.startY, s.endY, e);
+					let n = U((t - s.phaseStart) / s.phaseDuration);
+					w = n > .85 ? W(k, 0, (n - .85) / .15) : k, T = W(1, .65, e);
 					break;
 				}
 				case "absorbed":
-					v = 0, y = 0, b = 0, x = .5;
+					b = 0, C = 0, w = 0, T = .5;
 					break;
 				case "emerging": {
-					let e = P(z((t - s.phaseStart) / s.phaseDuration));
-					v = B(s.startX, s.endX, e), y = B(s.startY, s.endY, e);
-					let n = z((t - s.phaseStart) / s.phaseDuration);
-					b = n < .1 ? B(0, k, n / .1) : k, x = B(.65, 1, e);
+					let e = L(U((t - s.phaseStart) / s.phaseDuration));
+					b = W(s.startX, s.endX, e), C = W(s.startY, s.endY, e);
+					let n = U((t - s.phaseStart) / s.phaseDuration);
+					w = n < .1 ? W(0, k, n / .1) : k, T = W(.65, 1, e);
+					break;
+				}
+				case "fissioning": {
+					let e = U((t - s.phaseStart) / s.phaseDuration), n = R(e), r = M + S * n, i = q(s, t, o.orbitSpeedScale, o.wobbleScale), a = Math.hypot(i.x, i.y), c = a > 1e-4 ? i.x / a : Math.cos(s.phaseOffset), l = a > 1e-4 ? i.y / a : Math.sin(s.phaseOffset);
+					b = c * r, C = l * r, w = k, T = W(.78, 1, U(n)), !s.snapFired && e >= .42 && (s.snapFired = !0, d?.(Math.min(.4, .18 + .22 * x)));
 					break;
 				}
 			}
-			let S = i[u];
-			S.x = v, S.y = y, S.radius = f * x, S.opacity = b;
+			let E = i[f];
+			E.x = b, E.y = C, E.radius = m * T, E.opacity = w, E.fissioning = s.fissioning;
 		}
 	}
-	function p() {
+	function h() {
 		for (let e of r) e.phaseOffset += (n() - .5) * .6, e.pertXPhase += n() * Math.PI, e.pertYPhase += n() * Math.PI;
 	}
-	function m(e) {
+	function g(e) {
 		n = o(s(e + "goo"));
 	}
-	function h() {
+	function _(e) {
+		d = e;
+	}
+	function v() {
 		for (let e of r) if (e.phase !== "orbiting") return !1;
 		return !0;
 	}
-	function g(e) {
+	function y(e) {
 		let t = Infinity;
 		for (let e of r) {
 			let n = e.phaseStart + e.phaseDuration;
@@ -453,16 +482,17 @@ function G(e, t) {
 	}
 	return {
 		sources: i,
-		tick: f,
-		nudge: p,
-		reseed: m,
-		isQuiescent: h,
-		nextEventMs: g
+		tick: m,
+		nudge: h,
+		reseed: g,
+		onPinch: _,
+		isQuiescent: v,
+		nextEventMs: y
 	};
 }
 //#endregion
 //#region src/components/custom/goo-blob/shaders/sdf-body.glsl.ts
-var K = "// --- SDF ---\n\nfloat sdCircle(vec2 p, vec2 center, float radius) {\n    return length(p - center) - radius;\n}\n\n// AX.W15 — the VALUE+GRADIENT circle: returns vec3(d, dd/dx, dd/dy). The circle\n// gradient is the UNIT radial direction (p - center)/length(p - center) (the\n// eikonal property holds EXACTLY for a circle), guarded at the degenerate centre.\n// This is the keystone of the analytic-gradient smin: carrying the gradient\n// through the scene field DELETES the per-pixel 4-tap normal (4 full sceneDist\n// evals, each running 3-octave FBM + the sat/trail loops).\nvec3 sdgCircle(vec2 p, vec2 center, float radius) {\n    vec2 d = p - center;\n    float len = length(d);\n    vec2 grad = len > 1e-6 ? d / len : vec2(0.0, 1.0);\n    return vec3(len - radius, grad);\n}\n\n// IQ 2024 quadratic-polynomial smin, NORMALIZED so k is a blend band in distance\n// units. The k *= 4.0 pre-scale makes the h*h*k*0.25 poly's effective blend width\n// equal to k (the un-normalized form's width was 0.25*k).\nfloat sminQuadratic(float a, float b, float k) {\n    k *= 4.0;\n    float h = max(k - abs(a - b), 0.0) / k;\n    return min(a, b) - h * h * k * 0.25;\n}\n\n// IQ circular smin — a quarter-circle fillet at the seam (rounder than the\n// quadratic). The k *= 1.0/(1.0-sqrt(0.5)) factor normalizes k to the same\n// blend-band units so the merge axis swaps shape without re-tuning uSmoothK.\nfloat sminCircular(float a, float b, float k) {\n    k *= 1.0 / (1.0 - sqrt(0.5));\n    float h = max(k - abs(a - b), 0.0) / k;\n    return min(a, b) - k * 0.5 * (1.0 + h - sqrt(1.0 - h * (h - 2.0)));\n}\n\n// Config-gated merge dispatch — uMerge 0 = quadratic, 1 = circular.\nfloat smin(float a, float b, float k) {\n    if (uMerge > 0.5) return sminCircular(a, b, k);\n    return sminQuadratic(a, b, k);\n}\n\n// AX.W15 — the VALUE+GRADIENT smin. a/b are vec3(dist, grad.xy); the blended\n// gradient is mix(a.yz, b.yz, w) with the interpolation weight w the distance\n// blend uses (IQ's vec2-smin insight). This propagates the analytic gradient\n// through the merge so the surface normal reads the field gradient DIRECTLY — no\n// finite-difference tap. The blend weight matches each variant.\nvec3 sminQuadraticG(vec3 a, vec3 b, float k) {\n    k *= 4.0;\n    float h = max(k - abs(a.x - b.x), 0.0) / k;\n    float m = h * h * k * 0.25;\n    // h-weight that selects which surface dominates (0 → a, 1 → b).\n    float w = (a.x < b.x) ? 0.5 * h * h : 1.0 - 0.5 * h * h;\n    float dist = min(a.x, b.x) - m;\n    vec2 grad = mix(a.yz, b.yz, w);\n    return vec3(dist, grad);\n}\n\nvec3 sminCircularG(vec3 a, vec3 b, float k) {\n    k *= 1.0 / (1.0 - sqrt(0.5));\n    float h = max(k - abs(a.x - b.x), 0.0) / k;\n    float m = k * 0.5 * (1.0 + h - sqrt(1.0 - h * (h - 2.0)));\n    float w = (a.x < b.x) ? 0.5 * h : 1.0 - 0.5 * h;\n    float dist = min(a.x, b.x) - m;\n    vec2 grad = mix(a.yz, b.yz, w);\n    return vec3(dist, grad);\n}\n\nvec3 sminG(vec3 a, vec3 b, float k) {\n    if (uMerge > 0.5) return sminCircularG(a, b, k);\n    return sminQuadraticG(a, b, k);\n}", q = "// --- Noise ---\n\nfloat hash21(vec2 p) {\n    vec3 p3 = fract(vec3(p.xyx) * 0.1031);\n    p3 += dot(p3, p3.yzx + 33.33);\n    return fract((p3.x + p3.y) * p3.z);\n}\n\n// IQ analytic-derivative gradient noise — value in .x, ANALYTIC gradient\n// (d/dx, d/dy) in .yz, from one eval. The quintic fade has a continuous 2nd\n// derivative so the gradient is smooth (a cubic smoothstep fade would crease it).\n// (Inigo Quilez, https://iquilezles.org/articles/gradientnoise/.)\nvec3 noised(vec2 p) {\n    vec2 i = floor(p);\n    vec2 f = fract(p);\n\n    // Quintic fade + its derivative.\n    vec2 u = f * f * f * (f * (f * 6.0 - 15.0) + 10.0);\n    vec2 du = 30.0 * f * f * (f * (f - 2.0) + 1.0);\n\n    // Per-corner gradient vectors from the hash (centred to [-1, 1]).\n    vec2 ga = -1.0 + 2.0 * vec2(hash21(i + vec2(0.0, 0.0)), hash21(i + vec2(7.3, 0.0)));\n    vec2 gb = -1.0 + 2.0 * vec2(hash21(i + vec2(1.0, 0.0)), hash21(i + vec2(8.3, 0.0)));\n    vec2 gc = -1.0 + 2.0 * vec2(hash21(i + vec2(0.0, 1.0)), hash21(i + vec2(7.3, 1.0)));\n    vec2 gd = -1.0 + 2.0 * vec2(hash21(i + vec2(1.0, 1.0)), hash21(i + vec2(8.3, 1.0)));\n\n    float va = dot(ga, f - vec2(0.0, 0.0));\n    float vb = dot(gb, f - vec2(1.0, 0.0));\n    float vc = dot(gc, f - vec2(0.0, 1.0));\n    float vd = dot(gd, f - vec2(1.0, 1.0));\n\n    // Value, remapped to [0,1] so it drops in for the prior valueNoise.\n    float n = va + u.x * (vb - va) + u.y * (vc - va) + u.x * u.y * (va - vb - vc + vd);\n    float value = 0.5 + 0.5 * n;\n\n    // Analytic gradient (chain rule through the bilinear blend + the fade).\n    vec2 grad =\n        ga + u.x * (gb - ga) + u.y * (gc - ga) + u.x * u.y * (ga - gb - gc + gd) +\n        du * (vec2(u.y, u.x) * (va - vb - vc + vd) + vec2(vb, vc) - va);\n    return vec3(value, 0.5 * grad);\n}\n\n// Rotated-octave FBM over the analytic gradient noise: a fixed ~0.5 rad 2x2\n// rotation between octaves breaks the axis-aligned banding. The rotation CONSTANT\n// is spliced from the shared chunk (AV.W2 — the one FBM_ROT); the loops below\n// stay blob-local (param octaves + the LIQUID lacunarity/persistence, per §3a).", J = 1.8, Y = .42, X = `float fbm(vec2 p, int octaves) {
+var te = "// --- SDF ---\n\nfloat sdCircle(vec2 p, vec2 center, float radius) {\n    return length(p - center) - radius;\n}\n\n// AX.W15 — the VALUE+GRADIENT circle: returns vec3(d, dd/dx, dd/dy). The circle\n// gradient is the UNIT radial direction (p - center)/length(p - center) (the\n// eikonal property holds EXACTLY for a circle), guarded at the degenerate centre.\n// This is the keystone of the analytic-gradient smin: carrying the gradient\n// through the scene field DELETES the per-pixel 4-tap normal (4 full sceneDist\n// evals, each running 3-octave FBM + the sat/trail loops).\nvec3 sdgCircle(vec2 p, vec2 center, float radius) {\n    vec2 d = p - center;\n    float len = length(d);\n    vec2 grad = len > 1e-6 ? d / len : vec2(0.0, 1.0);\n    return vec3(len - radius, grad);\n}\n\n// IQ 2024 quadratic-polynomial smin, NORMALIZED so k is a blend band in distance\n// units. The k *= 4.0 pre-scale makes the h*h*k*0.25 poly's effective blend width\n// equal to k (the un-normalized form's width was 0.25*k).\nfloat sminQuadratic(float a, float b, float k) {\n    k *= 4.0;\n    float h = max(k - abs(a - b), 0.0) / k;\n    return min(a, b) - h * h * k * 0.25;\n}\n\n// IQ circular smin — a quarter-circle fillet at the seam (rounder than the\n// quadratic). The k *= 1.0/(1.0-sqrt(0.5)) factor normalizes k to the same\n// blend-band units so the merge axis swaps shape without re-tuning uSmoothK.\nfloat sminCircular(float a, float b, float k) {\n    k *= 1.0 / (1.0 - sqrt(0.5));\n    float h = max(k - abs(a - b), 0.0) / k;\n    return min(a, b) - k * 0.5 * (1.0 + h - sqrt(1.0 - h * (h - 2.0)));\n}\n\n// Config-gated merge dispatch — uMerge 0 = quadratic, 1 = circular.\nfloat smin(float a, float b, float k) {\n    if (uMerge > 0.5) return sminCircular(a, b, k);\n    return sminQuadratic(a, b, k);\n}\n\n// AX.W15 — the VALUE+GRADIENT smin. a/b are vec3(dist, grad.xy); the blended\n// gradient is mix(a.yz, b.yz, w) with the interpolation weight w the distance\n// blend uses (IQ's vec2-smin insight). This propagates the analytic gradient\n// through the merge so the surface normal reads the field gradient DIRECTLY — no\n// finite-difference tap. The blend weight matches each variant.\nvec3 sminQuadraticG(vec3 a, vec3 b, float k) {\n    k *= 4.0;\n    float h = max(k - abs(a.x - b.x), 0.0) / k;\n    float m = h * h * k * 0.25;\n    // h-weight that selects which surface dominates (0 → a, 1 → b).\n    float w = (a.x < b.x) ? 0.5 * h * h : 1.0 - 0.5 * h * h;\n    float dist = min(a.x, b.x) - m;\n    vec2 grad = mix(a.yz, b.yz, w);\n    return vec3(dist, grad);\n}\n\nvec3 sminCircularG(vec3 a, vec3 b, float k) {\n    k *= 1.0 / (1.0 - sqrt(0.5));\n    float h = max(k - abs(a.x - b.x), 0.0) / k;\n    float m = k * 0.5 * (1.0 + h - sqrt(1.0 - h * (h - 2.0)));\n    float w = (a.x < b.x) ? 0.5 * h : 1.0 - 0.5 * h;\n    float dist = min(a.x, b.x) - m;\n    vec2 grad = mix(a.yz, b.yz, w);\n    return vec3(dist, grad);\n}\n\nvec3 sminG(vec3 a, vec3 b, float k) {\n    if (uMerge > 0.5) return sminCircularG(a, b, k);\n    return sminQuadraticG(a, b, k);\n}", ne = "// --- Noise ---\n\nfloat hash21(vec2 p) {\n    vec3 p3 = fract(vec3(p.xyx) * 0.1031);\n    p3 += dot(p3, p3.yzx + 33.33);\n    return fract((p3.x + p3.y) * p3.z);\n}\n\n// IQ analytic-derivative gradient noise — value in .x, ANALYTIC gradient\n// (d/dx, d/dy) in .yz, from one eval. The quintic fade has a continuous 2nd\n// derivative so the gradient is smooth (a cubic smoothstep fade would crease it).\n// (Inigo Quilez, https://iquilezles.org/articles/gradientnoise/.)\nvec3 noised(vec2 p) {\n    vec2 i = floor(p);\n    vec2 f = fract(p);\n\n    // Quintic fade + its derivative.\n    vec2 u = f * f * f * (f * (f * 6.0 - 15.0) + 10.0);\n    vec2 du = 30.0 * f * f * (f * (f - 2.0) + 1.0);\n\n    // Per-corner gradient vectors from the hash (centred to [-1, 1]).\n    vec2 ga = -1.0 + 2.0 * vec2(hash21(i + vec2(0.0, 0.0)), hash21(i + vec2(7.3, 0.0)));\n    vec2 gb = -1.0 + 2.0 * vec2(hash21(i + vec2(1.0, 0.0)), hash21(i + vec2(8.3, 0.0)));\n    vec2 gc = -1.0 + 2.0 * vec2(hash21(i + vec2(0.0, 1.0)), hash21(i + vec2(7.3, 1.0)));\n    vec2 gd = -1.0 + 2.0 * vec2(hash21(i + vec2(1.0, 1.0)), hash21(i + vec2(8.3, 1.0)));\n\n    float va = dot(ga, f - vec2(0.0, 0.0));\n    float vb = dot(gb, f - vec2(1.0, 0.0));\n    float vc = dot(gc, f - vec2(0.0, 1.0));\n    float vd = dot(gd, f - vec2(1.0, 1.0));\n\n    // Value, remapped to [0,1] so it drops in for the prior valueNoise.\n    float n = va + u.x * (vb - va) + u.y * (vc - va) + u.x * u.y * (va - vb - vc + vd);\n    float value = 0.5 + 0.5 * n;\n\n    // Analytic gradient (chain rule through the bilinear blend + the fade).\n    vec2 grad =\n        ga + u.x * (gb - ga) + u.y * (gc - ga) + u.x * u.y * (ga - gb - gc + gd) +\n        du * (vec2(u.y, u.x) * (va - vb - vc + vd) + vec2(vb, vc) - va);\n    return vec3(value, 0.5 * grad);\n}\n\n// Rotated-octave FBM over the analytic gradient noise: a fixed ~0.5 rad 2x2\n// rotation between octaves breaks the axis-aligned banding. The rotation CONSTANT\n// is spliced from the shared chunk (AV.W2 — the one FBM_ROT); the loops below\n// stay blob-local (param octaves + the LIQUID lacunarity/persistence, per §3a).", Y = 1.8, X = .42, re = `float fbm(vec2 p, int octaves) {
     float value = 0.0;
     float amp = 0.5;
     float freq = 1.0;
@@ -470,8 +500,8 @@ var K = "// --- SDF ---\n\nfloat sdCircle(vec2 p, vec2 center, float radius) {\n
         if (i >= octaves) break;
         value += amp * noised(p * freq).x;
         p = FBM_ROT * p;
-        freq *= ${J.toFixed(1)};
-        amp *= ${Y.toFixed(2)};
+        freq *= ${Y.toFixed(1)};
+        amp *= ${X.toFixed(2)};
     }
     return value;
 }
@@ -498,8 +528,8 @@ vec3 fbmG(vec2 p, int octaves) {
         // d(value)/dp = amp · freq · (R^i)^T · noised.yz  (the rotated chain rule).
         grad += amp * freq * (transpose(rotAccum) * n.yz);
         pp = FBM_ROT * pp;
-        freq *= ${J.toFixed(1)};
-        amp *= ${Y.toFixed(2)};
+        freq *= ${Y.toFixed(1)};
+        amp *= ${X.toFixed(2)};
         rotAccum = FBM_ROT * rotAccum;
     }
     return vec3(value, grad);
@@ -514,9 +544,9 @@ vec3 fbmWarpedG(vec2 p, int octaves, float warpAmp) {
     vec2 q = vec2(fbm(p + vec2(0.0, 0.0), octaves), fbm(p + vec2(5.2, 1.3), octaves));
     return fbmG(p + warpAmp * (q - 0.5), octaves);
 }
-`, ee = "// Edit #4's gamut step — hue-preserving inward chroma clamp. If the linear result\n// is out of [0,1], bisect chroma toward 0 (L + h fixed) until in gamut. Preserves\n// hue exactly (only C shrinks). 16 steps resolve C to < 2^-16 of its span.\nbool inGamut(vec3 lin) {\n    return all(greaterThanEqual(lin, vec3(0.0))) && all(lessThanEqual(lin, vec3(1.0)));\n}\nvec3 gamutClampOklch(vec3 lch) {\n    vec3 lin = oklabToLinearSrgb(oklchToOklab(lch));\n    if (inGamut(lin)) return lch;\n    float lo = 0.0;\n    float hi = lch.y;\n    for (int i = 0; i < 16; i++) {\n        float mid = 0.5 * (lo + hi);\n        vec3 test = vec3(lch.x, mid, lch.z);\n        if (inGamut(oklabToLinearSrgb(oklchToOklab(test)))) {\n            lo = mid;\n        } else {\n            hi = mid;\n        }\n    }\n    return vec3(lch.x, lo, lch.z);\n}", te = "#version 300 es\nprecision highp float;\n\nin vec2 vUv;\nout vec4 fragColor;\n\nuniform vec2 uResolution;\nuniform float uTime;\nuniform vec3 uBaseColor;\n\n// Multi-stop palette (W11.b) — 2-4 in-family OKLCh stops (uploaded GAMMA sRGB, like\n// uBaseColor). uStopCount <= 1 falls back to uBaseColor (the single-color default,\n// zero regression). MAX_BLOB_STOPS is a compile-time #define.\n#define MAX_BLOB_STOPS 4\nuniform int uStopCount;\nuniform vec3 uPalette[MAX_BLOB_STOPS];\n\n// Pointer\nuniform vec2 uPointer;\nuniform float uPointerActive;\nuniform float uPointerAttraction;\nuniform float uPointerStrength;\n\n// Main body\nuniform float uBodyRadius;\n\n// Pulsation\nuniform float uPulsePhase;\nuniform float uPulseAmp;\n\n// Surface noise\nuniform float uNoiseAmp;\nuniform float uNoiseFreq;\nuniform float uNoiseSpeed;\n\n// Gooey\nuniform float uSmoothK;\nuniform float uMerge;   // 0 = quadratic smin, 1 = circular smin (rounder menisci)\n\n// BC.W-GOOBLOB-PLAIN — the STAGE-1 gate. 1.0 = the plain shadowless lightless\n// fill-only floor (variant=\"blob\"); 0.0 = the full lit pipeline (variant=\"meatball\" /\n// STAGE 2). The WebGL2 fallback transcribes the SAME stripped branch the WGSL primary\n// gates (research/viz/goo-blob.md §5 — the GLSL twin).\nuniform float uStage;\n\n// AX.W16 (arm 5) — the PRE-FBM bounding-discard radius (UV space). main() early-outs\n// to a transparent write for any fragment OUTSIDE this radius BEFORE the two 3-octave\n// FBM calls + the OKLCh round-trip. The renderer uploads it PADDED by every\n// outward-expanding term (body + eccentric orbit + sat radius + smin band + FBM edge\n// amplitude + click-pulse + a lean/squash safety pad) so it NEVER clips the wet\n// meniscus (IQ: inflate the bound to match any outward op).\nuniform float uMaxReach;\n\n// Membrane — domain-warp strength on the FBM edge displacement (0 = plain fbm).\nuniform float uWarpAmp;\n\n// Lit glass surface (W9.b) — gated behind uLit so the flat fill stays default.\nuniform float uLit;             // 0 = flat fill (default), 1 = lit droplet\nuniform vec3 uRimColor;         // Fresnel rim tint (GAMMA sRGB, via the /color leaf)\nuniform vec3 uLightDir;         // normalized light direction (x, y, z)\nuniform float uSpecStrength;    // Blinn-Phong specular weight\nuniform float uSpecShininess;   // specular exponent (16-64, tight glint)\nuniform float uRimPower;        // Fresnel/Schlick exponent (~2-3)\nuniform float uRimStrength;     // Fresnel rim weight\n\n// BC.W-GOOBLOB-MEATBALL — the procedural 2D SDF soft-shadow march (T2, the GLSL twin of\n// the WGSL res.z/res.w lanes). uShadow gates the IQ rmshadows soft contact shadow that\n// FOLLOWS the irregular silhouette; uShadowSoftness is the penumbra hardness (4-48).\nuniform float uShadow;          // 0 = no shadow, 1 = the soft contact shadow ON\nuniform float uShadowSoftness;  // penumbra hardness (inverse light-source size, 4-48)\n\n// Iridescence + fake subsurface (W11.a) — translucent-gel read.\nuniform float uIridescence;     // warm-pearl rim sheen weight (0 = off)\nuniform float uIridHue;         // base hue (radians) the cosine palette biases toward (warm)\nuniform float uIridSpeed;       // animated-thickness scroll for the iridescent t\nuniform float uSssScale;        // fast-SSS back-light weight\nuniform float uSssPower;        // fast-SSS exponent\nuniform float uCoreGlow;        // thickness-driven inner-luminosity lift\n\n// Color perturbation (perceptually uniform — OKLCh L/C/h)\nuniform float uHueRange;        // degrees of hue swing (converted to radians)\nuniform float uSatShift;        // OKLCh chroma swing\nuniform float uBrightnessShift; // OKLCh lightness bias\nuniform float uColorNoiseFreq;\nuniform float uColorNoiseSpeed;\n\n// Satellites (max 4)\n#define MAX_SATS 4\nuniform int uSatCount;\nuniform vec2 uSatPos[MAX_SATS];\nuniform float uSatRadius[MAX_SATS];\nuniform float uSatOpacity[MAX_SATS];\n\n// Pointer trail (W10) — a decaying-radius pseudopod toward the cursor. uTrail is a\n// COMPILE-TIME-SIZED array (GLSL ES 3.00 forbids a uniform-sized array) with a\n// DYNAMIC break on uTrailCount, mirroring the satellite loop. Later samples paint\n// smaller metaballs (r *= 1 - i/N) so the trail tapers into a tail.\n#define TRAIL_N 15\nuniform int uTrailCount;\nuniform vec2 uTrailPos[TRAIL_N];\nuniform float uTrailRadius[TRAIL_N];\n\n// Velocity-driven volume-preserving squash-and-stretch (W10).\nuniform vec2 uVelocity;   // smoothed pointer velocity (motion direction + speed)\nuniform float uStretch;   // squash-stretch magnitude (0 = off)\n\nconst float PI = 3.141592653589793;\n", Z = "\n", ne = te + Z + q + Z + `${e}
+`, ie = "// Edit #4's gamut step — hue-preserving inward chroma clamp. If the linear result\n// is out of [0,1], bisect chroma toward 0 (L + h fixed) until in gamut. Preserves\n// hue exactly (only C shrinks). 16 steps resolve C to < 2^-16 of its span.\nbool inGamut(vec3 lin) {\n    return all(greaterThanEqual(lin, vec3(0.0))) && all(lessThanEqual(lin, vec3(1.0)));\n}\nvec3 gamutClampOklch(vec3 lch) {\n    vec3 lin = oklabToLinearSrgb(oklchToOklab(lch));\n    if (inGamut(lin)) return lch;\n    float lo = 0.0;\n    float hi = lch.y;\n    for (int i = 0; i < 16; i++) {\n        float mid = 0.5 * (lo + hi);\n        vec3 test = vec3(lch.x, mid, lch.z);\n        if (inGamut(oklabToLinearSrgb(oklchToOklab(test)))) {\n            lo = mid;\n        } else {\n            hi = mid;\n        }\n    }\n    return vec3(lch.x, lo, lch.z);\n}", ae = "#version 300 es\nprecision highp float;\n\nin vec2 vUv;\nout vec4 fragColor;\n\nuniform vec2 uResolution;\nuniform float uTime;\nuniform vec3 uBaseColor;\n\n// Multi-stop palette (W11.b) — 2-4 in-family OKLCh stops (uploaded GAMMA sRGB, like\n// uBaseColor). uStopCount <= 1 falls back to uBaseColor (the single-color default,\n// zero regression). MAX_BLOB_STOPS is a compile-time #define.\n#define MAX_BLOB_STOPS 4\nuniform int uStopCount;\nuniform vec3 uPalette[MAX_BLOB_STOPS];\n\n// Pointer\nuniform vec2 uPointer;\nuniform float uPointerActive;\nuniform float uPointerAttraction;\nuniform float uPointerStrength;\n\n// Main body\nuniform float uBodyRadius;\n\n// Pulsation\nuniform float uPulsePhase;\nuniform float uPulseAmp;\n\n// Surface noise\nuniform float uNoiseAmp;\nuniform float uNoiseFreq;\nuniform float uNoiseSpeed;\n\n// Gooey\nuniform float uSmoothK;\nuniform float uMerge;   // 0 = quadratic smin, 1 = circular smin (rounder menisci)\n\n// BC.W-GOOBLOB-PLAIN — the STAGE-1 gate. 1.0 = the plain shadowless lightless\n// fill-only floor (variant=\"blob\"); 0.0 = the full lit pipeline (variant=\"meatball\" /\n// STAGE 2). The WebGL2 fallback transcribes the SAME stripped branch the WGSL primary\n// gates (research/viz/goo-blob.md §5 — the GLSL twin).\n// BD.W-GOO-CAROUSEL-DECK — uStage is now DERIVED from uMorphT (uStage = uMorphT<=0 ?\n// 1 : 0) at upload; the SHADING morph (below) reads uMorphT directly.\nuniform float uStage;\n\n// BD.W-GOO-CAROUSEL-DECK — the blob to meatball SHADING-MORPH scalar. 0.0 = the flat blob\n// fill (byte-identical to the STAGE-1 floor), 1.0 = the fully-dressed lit meatball; the\n// in-between LERPS the surface shading (NOT the geometry — the smin field is shared). The\n// renderer resolves it from config.morphT (variant to endpoint back-compat: blob to 0,\n// meatball to 1). A consumer ANIMATING morphT 0 to 1 gets the live morph.\nuniform float uMorphT;\n\n// AX.W16 (arm 5) — the PRE-FBM bounding-discard radius (UV space). main() early-outs\n// to a transparent write for any fragment OUTSIDE this radius BEFORE the two 3-octave\n// FBM calls + the OKLCh round-trip. The renderer uploads it PADDED by every\n// outward-expanding term (body + eccentric orbit + sat radius + smin band + FBM edge\n// amplitude + click-pulse + a lean/squash safety pad) so it NEVER clips the wet\n// meniscus (IQ: inflate the bound to match any outward op).\nuniform float uMaxReach;\n\n// Membrane — domain-warp strength on the FBM edge displacement (0 = plain fbm).\nuniform float uWarpAmp;\n\n// Lit glass surface (W9.b) — gated behind uLit so the flat fill stays default.\nuniform float uLit;             // 0 = flat fill (default), 1 = lit droplet\nuniform vec3 uRimColor;         // Fresnel rim tint (GAMMA sRGB, via the /color leaf)\nuniform vec3 uLightDir;         // normalized light direction (x, y, z)\nuniform float uSpecStrength;    // Blinn-Phong specular weight\nuniform float uSpecShininess;   // specular exponent (16-64, tight glint)\nuniform float uRimPower;        // Fresnel/Schlick exponent (~2-3)\nuniform float uRimStrength;     // Fresnel rim weight\n\n// BC.W-GOOBLOB-MEATBALL — the procedural 2D SDF soft-shadow march (T2, the GLSL twin of\n// the WGSL res.z/res.w lanes). uShadow gates the IQ rmshadows soft contact shadow that\n// FOLLOWS the irregular silhouette; uShadowSoftness is the penumbra hardness (4-48).\nuniform float uShadow;          // 0 = no shadow, 1 = the soft contact shadow ON\nuniform float uShadowSoftness;  // penumbra hardness (inverse light-source size, 4-48)\n\n// Iridescence + fake subsurface (W11.a) — translucent-gel read.\nuniform float uIridescence;     // warm-pearl rim sheen weight (0 = off)\nuniform float uIridHue;         // base hue (radians) the cosine palette biases toward (warm)\nuniform float uIridSpeed;       // animated-thickness scroll for the iridescent t\nuniform float uSssScale;        // fast-SSS back-light weight\nuniform float uSssPower;        // fast-SSS exponent\nuniform float uCoreGlow;        // thickness-driven inner-luminosity lift\n\n// Color perturbation (perceptually uniform — OKLCh L/C/h)\nuniform float uHueRange;        // degrees of hue swing (converted to radians)\nuniform float uSatShift;        // OKLCh chroma swing\nuniform float uBrightnessShift; // OKLCh lightness bias\nuniform float uColorNoiseFreq;\nuniform float uColorNoiseSpeed;\n\n// Satellites (max 4)\n#define MAX_SATS 4\nuniform int uSatCount;\nuniform vec2 uSatPos[MAX_SATS];\nuniform float uSatRadius[MAX_SATS];\nuniform float uSatOpacity[MAX_SATS];\n\n// Pointer trail (W10) — a decaying-radius pseudopod toward the cursor. uTrail is a\n// COMPILE-TIME-SIZED array (GLSL ES 3.00 forbids a uniform-sized array) with a\n// DYNAMIC break on uTrailCount, mirroring the satellite loop. Later samples paint\n// smaller metaballs (r *= 1 - i/N) so the trail tapers into a tail.\n#define TRAIL_N 15\nuniform int uTrailCount;\nuniform vec2 uTrailPos[TRAIL_N];\nuniform float uTrailRadius[TRAIL_N];\n\n// Velocity-driven volume-preserving squash-and-stretch (W10).\nuniform vec2 uVelocity;   // smoothed pointer velocity (motion direction + speed)\nuniform float uStretch;   // squash-stretch magnitude (0 = off)\n\nconst float PI = 3.141592653589793;\n", Z = "\n", oe = ae + Z + ne + Z + `${e}
 
-` + X + Z + K + Z + `
+` + re + Z + te + Z + `
 // --- Color (Ottosson OKLab/OKLCh — value.js EXACT constants, transposed) ---
 //
 // The OETF (srgbToLinear/linearToSrgb) + the four Ottosson matrices + their space
@@ -528,10 +558,10 @@ vec3 fbmWarpedG(vec2 p, int octaves, float warpAmp) {
 ${n}
 ${i}
 
-` + ee + Z + "\n// AX.W15 — de-synced breath. The body throb was a single sin(uPulsePhase) — a\n// mechanical pulse that re-syncs every cycle. Three DETUNED sines at IRRATIONAL\n// frequency ratios (1, 0.13/0.09 ≈ 1.444, 0.31) never re-phase, so the membrane\n// breathes like a living creature (an asymmetric slower-exhale calm band) rather\n// than a metronome. Normalized so the composite stays in ~[-1, 1].\nfloat breath(float phase) {\n    return (sin(phase)\n          + 0.5 * sin(phase * 1.4444 + 1.7)\n          + 0.28 * sin(phase * 0.31 + 4.1)) / 1.78;\n}\n\n// AX.W15 — the composite SDF field WITH its ANALYTIC GRADIENT, returned as\n// vec3(dist, ∂d/∂x, ∂d/∂y). The domain-warped body membrane smin-merged with the\n// satellites + trail; the gradient propagates through sminG's mix(a.yz, b.yz, h)\n// so the surface normal reads the field gradient DIRECTLY (the 4-tap finite\n// difference is DELETED — its 4 full evals per lit pixel are gone).\n//\n// The ANISOTROPIC SQUASH (W10 volume-preserving stretch) is a linear map\n// bodyUv = D·R·uv (R rotates world→(motion-axis, perp); D = diag(1/sa, sa)). The\n// body field's gradient is computed in bodyUv space, so it MUST be transformed back\n// to uv space by M^T = R^T·D (the inverse-transpose handling the SOTA-deepening [2]\n// warns the 4-tap did implicitly by sampling in screen space).\nvec3 sceneDistG(vec2 uv) {\n    // De-synced pulsing body radius.\n    float bodyR = uBodyRadius + breath(uPulsePhase) * uPulseAmp;\n\n    // Velocity-driven VOLUME-PRESERVING squash-and-stretch (W10), SATURATED (AX.W46\n    // D5). sa is the stretch along motion, EXACTLY 1/sa perpendicular (area-preserving).\n    // The basis (ax, perp) is captured so the gradient can be transformed back below.\n    //\n    // The W10 form sa = 1 + speed*uStretch rode an UNBOUNDED critically-damped spring\n    // velocity (O(5-8)/s on a fast flick) -> sa ~ 1.78-2.25x elongation, a violent\n    // taffy-pull. A tanh saturation sa = 1 + tanh(speed*k)*uStretch caps the elongation\n    // at 1 + uStretch (~1.5x at the default stretch 0.5) no matter how fast the flick —\n    // the fastest flick reads as a lively bounce, never a rubber-band snap (the standard\n    // volume-preserving squash-and-stretch restraint). k = 1.6 keeps the tanh in its\n    // near-linear regime for a slow drag (the gentle squash still reads) and saturates\n    // only the fast flick.\n    vec2 bodyUv = uv;\n    vec2 ax = vec2(1.0, 0.0);\n    vec2 perp = vec2(0.0, 1.0);\n    float sa = 1.0;\n    float speed = length(uVelocity);\n    bool squashed = uStretch > 0.0 && speed > 1e-4;\n    if (squashed) {\n        ax = uVelocity / speed;\n        perp = vec2(-ax.y, ax.x);\n        sa = 1.0 + tanh(speed * 1.6) * uStretch;\n        bodyUv = vec2(dot(uv, ax) / sa, dot(uv, perp) * sa);\n    }\n\n    // Domain-warped FBM displacement WITH its analytic gradient (the warp-Jacobian\n    // approximation per [2]). disp = (noiseVal - 0.5) · uNoiseAmp; its gradient in\n    // bodyUv space is uNoiseAmp · uNoiseFreq · fbmGrad (the uNoiseFreq chain factor\n    // from the bodyUv·uNoiseFreq argument).\n    vec3 fbmv = fbmWarpedG(bodyUv * uNoiseFreq + uTime * uNoiseSpeed, 3, uWarpAmp);\n    float bodyDisplacement = (fbmv.x - 0.5) * uNoiseAmp;\n    vec2 dispGradBody = uNoiseAmp * uNoiseFreq * fbmv.yz;\n\n    // Body circle WITH gradient (in bodyUv space). The displaced radius subtracts\n    // the displacement, so the body gradient is circleGrad - dispGrad (sign per\n    // [2]: d = circleDist - displacement ⇒ grad = circleGrad - dispGrad).\n    vec3 bodyG = sdgCircle(bodyUv, vec2(0.0), bodyR + bodyDisplacement);\n    vec2 bodyGradBody = bodyG.yz - dispGradBody;\n\n    // Transform the body gradient back to uv space: g_uv = R^T · D · g_body, i.e.\n    // scale by D = diag(1/sa, sa) in the (ax, perp) frame then reconstruct in world.\n    vec2 bodyGradUv;\n    if (squashed) {\n        vec2 scaled = vec2(bodyGradBody.x / sa, bodyGradBody.y * sa);\n        bodyGradUv = scaled.x * ax + scaled.y * perp;\n    } else {\n        bodyGradUv = bodyGradBody;\n    }\n    vec3 d = vec3(bodyG.x, bodyGradUv);\n\n    // Satellites — smin-merged into the body (uMerge selects quadratic/circular).\n    for (int i = 0; i < MAX_SATS; i++) {\n        if (i >= uSatCount) break;\n        vec3 satG = sdgCircle(uv, uSatPos[i], uSatRadius[i]);\n        satG.x += (1.0 - uSatOpacity[i]) * 0.3;\n        d = sminG(d, satG, uSmoothK);\n    }\n\n    // Pointer trail — a decaying-radius pseudopod reaching toward the cursor,\n    // smin-merged so it stretches an elastic limb and snaps back.\n    for (int i = 0; i < TRAIL_N; i++) {\n        if (i >= uTrailCount) break;\n        vec3 trailG = sdgCircle(uv, uTrailPos[i], uTrailRadius[i]);\n        d = sminG(d, trailG, uSmoothK);\n    }\n    return d;\n}\n\n// AX.W15 — the analytic surface normal. The field gradient grad2d arrives\n// DIRECTLY from sceneDistG (no 4-tap). The smin field is sub-unit (the CD family\n// |grad| ≤ 1) so the gradient is NOT unit-length — normalize() before the dome\n// lift (the unit contract is on the FINAL lifted normal N, not the raw field\n// gradient). The Z dome: 'interior' is the normalized depth inward from the rim\n// (-d/bodyR); z = sqrt(1 - (1 - interior)^2) is the unit half-sphere profile (flat\n// centre, steep rim). N tilts outward along the field gradient near the rim and\n// points at the viewer deep inside.\nvec3 surfaceNormalFromGrad(vec2 grad2d, float d, float bodyR) {\n    vec2 g = normalize(grad2d + vec2(1e-6));\n    float interior = clamp(-d / max(bodyR, 1e-4), 0.0, 1.0);\n    float z = sqrt(max(0.0, 1.0 - (1.0 - interior) * (1.0 - interior)));\n    return normalize(vec3(g * (1.0 - z), z) + vec3(0.0, 0.0, 1e-6));\n}\n\n// BC.W-GOOBLOB-MEATBALL — the 2D SDF soft-shadow march (IQ rmshadows improved-penumbra,\n// research/viz/goo-blob.md §2.3, §5). The GLSL twin of the WGSL softShadow2D — a procedural\n// soft contact shadow FOLLOWING the irregular metaball silhouette (NOT a hard disc/box\n// shadow). Marches from ro along rd (the in-plane projection of uLightDir),\n// accumulating the closest miss via the Aaltonen penumbra (y = h*h/(2*ph) + d =\n// sqrt(h*h - y*y)) that kills the banding the naive res = min(res, h/(w*t)) shows. CHEAP:\n// 24 steps, re-uses the SAME sceneDistG().x field (NO separate field re-eval, NO FBO/\n// multi-pass — the procedural field needs none). w is the penumbra hardness (inverse\n// light-source size -> uShadowSoftness). Derivative-FREE. Returns 1.0 (lit) -> 0.0 (occluded).\nfloat softShadow2D(vec2 ro, vec2 rd, float mint, float maxt, float w) {\n    float res = 1.0;\n    float ph = 1e20;\n    float t = mint;\n    for (int i = 0; i < 24; i++) {\n        if (t >= maxt) break;\n        float h = sceneDistG(ro + rd * t).x;\n        if (h < 0.001) return 0.0;\n        float y = h * h / (2.0 * ph);\n        float dd = sqrt(max(h * h - y * y, 0.0));\n        res = min(res, dd / (w * max(0.001, t - y)));\n        ph = h;\n        t += h;\n    }\n    return clamp(res, 0.0, 1.0);\n}\n\n// W11.b — sample the multi-stop palette at t in [0,1], interpolating adjacent\n// stops in OKLab with a MIDPOINT CHROMA-BUMP (a linear OKLab mix of a vivid pair\n// dips chroma through grey; the bump lifts it back). Falls back to uBaseColor when\n// uStopCount <= 1. Returns an OKLCh stop [L, C, h(rad)].\nvec3 samplePaletteOklch(float t) {\n    if (uStopCount <= 1) return oklabToOklch(srgbToOklab(uBaseColor));\n    float n = float(uStopCount);\n    float ft = clamp(t, 0.0, 1.0) * (n - 1.0);\n    int i0 = int(floor(ft));\n    int i1 = min(i0 + 1, uStopCount - 1);\n    float f = ft - float(i0);\n    // Fetch adjacent stops (uPalette is gamma sRGB) → OKLab.\n    vec3 labA = srgbToOklab(uPalette[i0]);\n    vec3 labB = srgbToOklab(uPalette[i1]);\n    vec3 lab = mix(labA, labB, f);\n    vec3 lch = oklabToOklch(lab);\n    // Midpoint chroma-bump: a bell (peaks at f=0.5) lifts C off the grey midpoint.\n    lch.y += 0.03 * sin(PI * f);\n    return lch;\n}\n\nvoid main() {\n    vec2 uv = vUv - 0.5;\n\n    // AX.W16 (arm 5) — PRE-FBM bounding early-out. The oversized canvas (1.6x the\n    // wrapper) runs the full fragment ALU — two 3-octave FBM evals + the OKLCh\n    // round-trip + the lit/iridescence/SSS block — on a large transparent border\n    // otherwise (~60% of the canvas is outside the droplet reach). Any fragment beyond\n    // uMaxReach (the PADDED worst-case painted reach, uploaded so it never clips the\n    // wet meniscus) writes transparent and returns BEFORE any of that work. This is a\n    // transparent WRITE, NOT a GLSL discard: the pass is a pure premultiplied blend\n    // with no depth, and discard disables the tiled-renderer fast path on mobile GPUs.\n    // No fwidth is computed before the return, so there is no\n    // derivative-in-non-uniform-control-flow hazard (the kept pixels compute their own\n    // fwidth normally below). The pad covers the pointer-lean uv shift, so testing the\n    // pre-deform uv is safe.\n    if (dot(uv, uv) > uMaxReach * uMaxReach) {\n        fragColor = vec4(0.0);\n        return;\n    }\n\n    // Pointer deformation — honor the SIGN of uPointerAttraction (W10): a positive\n    // attraction leans the body IN toward the cursor, a negative shies it AWAY. The\n    // signed influence flows straight into the UV shift (no hardcoded repulsion).\n    //\n    // AY.W-BLOB-CONFIG D2 — the SIGN was INVERTED in the rendered result. The\n    // sample-shift sign here was uv -= dir*influence; the live readback (an attraction\n    // sweep over a HELD pointer) measured the body shifting MORE toward the cursor as\n    // attraction went MORE NEGATIVE (a=-1 -> +0.081, a=0 -> +0.063, a=+1 -> +0.042) —\n    // exactly inverted: a shy-away -1 LUNGED at the cursor, a lean-in +1 pulled back.\n    // pointerDir = uPointer - uv is the toward-pointer direction in SAMPLE space; to\n    // MOVE the rendered body TOWARD the pointer the sample must shift ALONG pointerDir\n    // (sample the field from where the pointer is), so the sign is uv += dir*influence.\n    // A positive attraction now shifts the body toward the cursor (leans in), a negative\n    // shifts it away (shies) — the sign the comment + the demo slider always promised,\n    // now true in the render.\n    //\n    // AX.W46 D5 — the falloff radius NARROWS 0.65 → 0.5 (the calm-lean reconciliation).\n    // The W15 REDRESS widened it to 0.65 to clear a synthetic floor; the live π-lane\n    // read the result as a LUNGE. The drama lives in the STRENGTH (pointerStrength,\n    // dropped to 0.18 in types.ts), NOT the falloff — the falloff's job is purely to keep\n    // the lean COHERENT across the creature so the body, satellites and trail all tilt\n    // toward the cursor as ONE. At 0.5 the body (sitting ≈0.2 uv from the cursor on a\n    // hover-flick) is still well inside the cutoff, so the whole creature tilts as one —\n    // a gentle, coherent \"it notices you\" lean rather than the body lunging a body-width.\n    if (uPointerActive > 0.5) {\n        vec2 pointerDir = uPointer - uv;\n        float pointerDist = length(pointerDir);\n        float influence = smoothstep(0.5, 0.0, pointerDist) * uPointerAttraction * uPointerStrength;\n        uv += normalize(pointerDir + 1e-6) * influence;\n    }\n\n    // The de-synced pulsing body radius — also reused below for the inner-glow scale.\n    float bodyR = uBodyRadius + breath(uPulsePhase) * uPulseAmp;\n\n    // Composite domain-warped membrane field WITH its analytic gradient (AX.W15).\n    vec3 scene = sceneDistG(uv);\n    float d = scene.x;\n    vec2 fieldGrad = scene.yz;\n\n    // Edit #1 — fwidth-based anti-aliased edge: derive the smoothstep half-width\n    // from the SDF screen-space gradient so the edge stays ~1px regardless of zoom.\n    float aa = max(fwidth(d), 1e-6);\n    float alpha = 1.0 - smoothstep(-aa, aa, d);\n\n    if (alpha < 0.001) {\n        fragColor = vec4(0.0);\n        return;\n    }\n\n    // Edit #4 — perceptually-uniform OKLCh variation. Gamma base → linear → OKLab\n    // → OKLCh; perturb L/C/h; OKLCh → OKLab → linear; hue-preserving gamut clamp.\n    // The base is the multi-stop palette sample (W11.b) — distributed across the\n    // body/satellites by the color noise field — or uBaseColor when single-stop.\n    float colorNoise = fbm(uv * uColorNoiseFreq + uTime * uColorNoiseSpeed, 3);\n    vec3 oklch = samplePaletteOklch(colorNoise);\n\n    oklch.z += (colorNoise - 0.5) * uHueRange * (PI / 180.0);     // hue swing, radians\n    oklch.y = max(oklch.y + (colorNoise - 0.5) * uSatShift, 0.0); // chroma swing\n    oklch.x = clamp(oklch.x + uBrightnessShift, 0.0, 1.0);        // lightness bias\n\n    // ── BC.W-GOOBLOB-PLAIN — the STAGE-1 stripped floor (variant=blob). ──\n    // The minimal verifiable floor: SDF + smin (the meatball field, already in\n    // d/alpha above) + fwidth-AA (the crisp alpha) + warm-cream fill. NO surface\n    // normal, NO Fresnel, NO lit glint, NO iridescence, NO fake-SSS, NO shadow —\n    // deliberately FLAT. It returns BEFORE the STAGE-2 dressing, the GLSL twin of\n    // the WGSL primary's uStage > 0.5 branch. The teaching contrast: the \"it\n    // renders, it meatballs\" floor STAGE 2 layers the lit/shadow onto via the SAME\n    // uStage/uLit uniforms.\n    if (uStage > 0.5) {\n        oklch = gamutClampOklch(oklch);\n        vec3 lin1 = oklabToLinearSrgb(oklchToOklab(oklch));\n        vec3 rgb1 = clamp(linearToSrgb(lin1), 0.0, 1.0);\n        float ign1 = fract(52.9829189 * fract(dot(gl_FragCoord.xy, vec2(0.06711056, 0.00583715))));\n        rgb1 = clamp(rgb1 + (ign1 - 0.5) / 255.0, 0.0, 1.0);\n        fragColor = vec4(rgb1 * alpha, alpha);\n        return;\n    }\n\n    // Surface normal — from the ANALYTIC field gradient (AX.W15: the 4-tap is\n    // deleted). Computed ONCE here and reused by the iridescence (W11.a), the\n    // fake-SSS (W11.a), and the lit glass block (W9.b).\n    vec3 N = surfaceNormalFromGrad(fieldGrad, d, bodyR);\n    vec3 V = vec3(0.0, 0.0, 1.0);\n    float thickness = clamp(-d / max(bodyR, 1e-4), 0.0, 1.0); // 0 at rim, ~1 deep in\n    float fres = pow(1.0 - max(dot(N, V), 0.0), 2.5);          // rim-weighted angle\n\n    // ── W11.a iridescence — warm-biased IQ cosine palette driving OKLCh HUE ──\n    //\n    // A thin-film-like sheen on the RIM (fres-weighted), NOT the cold-blue default:\n    // an IQ cosine palette (a + b*cos(2pi(c*t + d))) maps the Fresnel/edge angle +\n    // the FBM color field + an animated thickness onto a HUE OFFSET biased to the\n    // warm arc (uIridHue), with the chroma clamped to a warm-pearl band so the rim\n    // reads cream/gold (congruent with the warm-cream-glass identity). Mixed into\n    // OKLCh BEFORE the gamut clamp. Default LOW.\n    if (uIridescence > 0.0) {\n        float t = fres + 0.3 * colorNoise + uTime * uIridSpeed;\n        // Cosine palette in HUE space, centred on the warm uIridHue (~0.18 turn of\n        // warm sweep so the sheen shimmers gold→peach, never swinging to cold blue).\n        float iridHue = uIridHue + 0.18 * PI * cos(2.0 * PI * t);\n        float w = fres * uIridescence;          // rim-weighted blend\n        oklch.z = mix(oklch.z, iridHue, w);\n        // Warm-pearl chroma cap: lift chroma toward a pale-pearl band, clamped.\n        oklch.y = min(oklch.y + 0.04 * w, oklch.y + 0.08);\n        oklch.x = min(oklch.x + 0.05 * w, 1.0); // a touch brighter at the sheen\n    }\n\n    // ── W11.a fake subsurface translucency — thickness inner-glow + back-light ──\n    //\n    // A bright translucent CORE fading to a light-leaking EDGE (the -d ramp) plus\n    // the fast-SSS back-light (light wrapping through the thin rim). Both lift OKLCh\n    // L and warm the hue, consuming the W9 normal. In OKLCh before the gamut clamp.\n    if (uCoreGlow > 0.0 || uSssScale > 0.0) {\n        // Beer-Lambert inner luminosity (AX.W15): a SATURATING 1 - exp(-k·thickness)\n        // curve — a flat thick core with a fast warm rim falloff — NOT the linear\n        // coreGlow·thickness ramp (which over-brightens deep interiors). k ≈ 3 reads\n        // as glass depth.\n        oklch.x = min(oklch.x + uCoreGlow * (1.0 - exp(-3.0 * thickness)), 1.0);\n        // Fast-SSS back-light: light wrapping through the thin (low-thickness) rim.\n        vec3 L = normalize(uLightDir + vec3(1e-6));\n        float back = pow(clamp(dot(V, -(L + N * thickness)), 0.0, 1.0), uSssPower);\n        float sss = back * uSssScale * (1.0 - thickness);\n        oklch.x = min(oklch.x + sss, 1.0);\n        // Warm the THIN leaking rim only (scale by 1 - thickness): the SSS hue-warm\n        // shift rides the rim, not the core (AX.W15 — only the thin rim warms).\n        oklch.z += sss * 0.1 * (1.0 - thickness);\n    }\n\n    // ── BC.W-GOOBLOB-MEATBALL — the soft contact shadow (T2, the GLSL twin). A procedural\n    //    soft shadow FOLLOWING the irregular silhouette: march from this fragment's uv\n    //    toward the in-plane projection of uLightDir; where the metaball field occludes the\n    //    light, DARKEN the OKLCh L — a soft grounded contact band under the dome\n    //    (PROCEDURAL, silhouette-following, unlike the static CSS gel-dome drop-shadow).\n    //    Rim-weighted (1 - thickness) so the lit interior stays bright. Gated uShadow > 0.5;\n    //    variant=blob (STAGE 1) ships shadowless via the uStage early-return above. ──\n    if (uShadow > 0.5) {\n        vec2 L2 = normalize(uLightDir.xy + vec2(1e-6));\n        float sh = softShadow2D(uv, L2, max(aa, 0.004), bodyR * 1.5, max(uShadowSoftness, 4.0));\n        // The contact band sits on the LOWER rim (the side AWAY from the light), gated by\n        // the away-facing rim weight + a thin-rim Gaussian band (peak thickness ≈ 0.18) so\n        // the deep-lit interior + the light-facing edge stay BRIGHT (the warm-cream body\n        // mean + the proof:blob-warm-default floor hold). The 0.20 cap is a soft grounding.\n        float away = clamp(-dot(normalize(fieldGrad + vec2(1e-6)), L2), 0.0, 1.0);\n        float band = exp(-pow((thickness - 0.18) / 0.16, 2.0));\n        float shadowDarken = (1.0 - sh) * away * band * 0.20;\n        oklch.x = max(oklch.x - shadowDarken, 0.0);\n    }\n\n    oklch = gamutClampOklch(oklch);\n\n    vec3 lin = oklabToLinearSrgb(oklchToOklab(oklch));\n\n    // ── W9.b lit glass surface — Blinn-Phong glint + Fresnel rim, in LINEAR ──\n    //\n    // The lit terms enter lin in LINEAR space BEFORE the OETF (linearToSrgb)\n    // and BEFORE the * alpha premultiply — a post-OETF apply double-gammas the\n    // highlight and fringes the premultiplied edge (the named A5/A2 trap). The\n    // prior flat inner-glow lift is SUBSUMED by the Fresnel rim (the curvature\n    // read now comes from the surface normal, not a flat smoothstep). Gated behind\n    // uLit so the flat fill stays the default (zero regression).\n    if (uLit > 0.5) {\n        // N, V, thickness are already computed in main() (reused by iridescence/SSS).\n        vec3 L = normalize(uLightDir + vec3(1e-6)); // light direction\n        vec3 H = normalize(L + V);                  // Blinn-Phong half-vector\n\n        // Warm-cream specular glint — a near-white OKLCh TINT (L~0.97, C~0.03,\n        // hue~85°) routed through the SAME spliced OKLCh matrices, NOT hardcoded\n        // sRGB white (pure white reads cheap-CG). Linearized for the linear add.\n        vec3 warmCream = oklabToLinearSrgb(oklchToOklab(vec3(0.97, 0.03, radians(85.0))));\n\n        // ENERGY-CONSERVING Blinn-Phong (AX.W15 [9]): the (shininess+2)/8 factor\n        // DECOUPLES shininess from strength (without it, raising shininess dims the\n        // glint, so the two knobs fight). SPECULAR ANTIALIASING: the FBM membrane\n        // makes the normal vary fast, so a tight glint (16-64) STROBES on small\n        // dock-grid blobs. Widen the effective lobe where the normal varies\n        // (Toksvig-style fwidth clamp) — drop shininess, NOT raise it. nVar reads\n        // the screen-space normal derivative; a high-variance pixel softens the\n        // exponent so the glint stays stable.\n        float nVar = length(fwidth(N));\n        float shininess = uSpecShininess / (1.0 + 24.0 * nVar);\n        float energyNorm = (shininess + 2.0) / 8.0;\n        float spec = pow(max(dot(N, H), 0.0), shininess) * uSpecStrength * energyNorm;\n\n        // Fresnel/Schlick rim — fed uRimColor (the --foreground warm rim via the\n        // /color leaf). FOREGROUND-AWARE MIN-CONTRAST GUARD (AX.W15 [10]):\n        // a var(--primary) blob in dark mode resolves a rim near the BODY color, so\n        // the rim washes out (no contrast). When the rim's luminance sits too close\n        // to the body's, the guard CHROMA-REDUCES and L-LIFTS the rim stop (a\n        // perceptual move, NOT a re-tint) so the curved rim always reads. Computed in\n        // OKLCh so the lift is perceptual.\n        vec3 rimOkl = oklabToOklch(srgbToOklab(uRimColor));\n        float bodyL = oklch.x; // the post-perturb body lightness\n        float dL = abs(rimOkl.x - bodyL);\n        // If the rim is within 0.22 L of the body, push it AWAY (lift if the body is\n        // dark, the dark-mode primary case; the chroma-reduce keeps it neutral-warm).\n        float lack = clamp((0.22 - dL) / 0.22, 0.0, 1.0);\n        rimOkl.x = clamp(mix(rimOkl.x, (bodyL < 0.5 ? 0.92 : 0.18), lack), 0.0, 1.0);\n        rimOkl.y *= mix(1.0, 0.5, lack); // chroma-reduce toward a pale pearl\n        vec3 rimLin = oklabToLinearSrgb(oklchToOklab(rimOkl));\n\n        float rimFres = pow(1.0 - max(dot(N, V), 0.0), uRimPower);\n        float rim = rimFres * uRimStrength * (1.0 - 0.6 * thickness);\n\n        // Combine the two warm highlights and add in LINEAR. max(spec, rim*scale)\n        // keeps the glint from stacking on the rim into a blown hotspot.\n        //\n        // AX.W46 D4 — CLAMP the linear highlight below unity before the add (the\n        // belt-and-braces guard the re-derived specStrength makes rarely-binding but\n        // never lets a worst-case normal blow to a hard white spot). The energy-\n        // conserving Blinn-Phong CAN still spike on a grazing normal where the Toksvig\n        // widen has not kicked in; capping the highlight magnitude at 0.85 (a contained\n        // warm gleam, sub-unity) means the OETF clamp below can never crush a\n        // pure-white pixel over the dome. The cap is on the highlight ALONE (not the\n        // composited lin), so the body color underneath is untouched — it caps only\n        // the additive catch-light's energy.\n        vec3 highlight = max(warmCream * spec, rimLin * rim);\n        highlight = min(highlight, vec3(0.85));\n        lin += highlight;\n    }\n\n    vec3 rgb = clamp(linearToSrgb(lin), 0.0, 1.0); // MANDATORY OETF — closes the seam\n\n    // IGN dither (AX.W15 [2][9]) — the low-chroma warm-cream dome BANDS on 8-bit\n    // panels (visible Mach steps in the smooth L roll). Interleaved-gradient noise at\n    // 1/255, applied AFTER linearToSrgb and BEFORE the *alpha premultiply (aurora\n    // ships the same splice). A triangular ±0.5 LSB dither decorrelates the steps.\n    float ign = fract(52.9829189 * fract(dot(gl_FragCoord.xy, vec2(0.06711056, 0.00583715))));\n    rgb = clamp(rgb + (ign - 0.5) / 255.0, 0.0, 1.0);\n\n    // Edit #8 — premultiply AFTER the OETF: straight-alpha gamma → premultiplied.\n    fragColor = vec4(rgb * alpha, alpha);\n}";
+` + ie + Z + "\n// AX.W15 — de-synced breath. The body throb was a single sin(uPulsePhase) — a\n// mechanical pulse that re-syncs every cycle. Three DETUNED sines at IRRATIONAL\n// frequency ratios (1, 0.13/0.09 ≈ 1.444, 0.31) never re-phase, so the membrane\n// breathes like a living creature (an asymmetric slower-exhale calm band) rather\n// than a metronome. Normalized so the composite stays in ~[-1, 1].\nfloat breath(float phase) {\n    return (sin(phase)\n          + 0.5 * sin(phase * 1.4444 + 1.7)\n          + 0.28 * sin(phase * 0.31 + 4.1)) / 1.78;\n}\n\n// AX.W15 — the composite SDF field WITH its ANALYTIC GRADIENT, returned as\n// vec3(dist, ∂d/∂x, ∂d/∂y). The domain-warped body membrane smin-merged with the\n// satellites + trail; the gradient propagates through sminG's mix(a.yz, b.yz, h)\n// so the surface normal reads the field gradient DIRECTLY (the 4-tap finite\n// difference is DELETED — its 4 full evals per lit pixel are gone).\n//\n// The ANISOTROPIC SQUASH (W10 volume-preserving stretch) is a linear map\n// bodyUv = D·R·uv (R rotates world→(motion-axis, perp); D = diag(1/sa, sa)). The\n// body field's gradient is computed in bodyUv space, so it MUST be transformed back\n// to uv space by M^T = R^T·D (the inverse-transpose handling the SOTA-deepening [2]\n// warns the 4-tap did implicitly by sampling in screen space).\nvec3 sceneDistG(vec2 uv) {\n    // De-synced pulsing body radius.\n    float bodyR = uBodyRadius + breath(uPulsePhase) * uPulseAmp;\n\n    // Velocity-driven VOLUME-PRESERVING squash-and-stretch (W10), SATURATED (AX.W46\n    // D5). sa is the stretch along motion, EXACTLY 1/sa perpendicular (area-preserving).\n    // The basis (ax, perp) is captured so the gradient can be transformed back below.\n    //\n    // The W10 form sa = 1 + speed*uStretch rode an UNBOUNDED critically-damped spring\n    // velocity (O(5-8)/s on a fast flick) -> sa ~ 1.78-2.25x elongation, a violent\n    // taffy-pull. A tanh saturation sa = 1 + tanh(speed*k)*uStretch caps the elongation\n    // at 1 + uStretch (~1.5x at the default stretch 0.5) no matter how fast the flick —\n    // the fastest flick reads as a lively bounce, never a rubber-band snap (the standard\n    // volume-preserving squash-and-stretch restraint). k = 1.6 keeps the tanh in its\n    // near-linear regime for a slow drag (the gentle squash still reads) and saturates\n    // only the fast flick.\n    vec2 bodyUv = uv;\n    vec2 ax = vec2(1.0, 0.0);\n    vec2 perp = vec2(0.0, 1.0);\n    float sa = 1.0;\n    float speed = length(uVelocity);\n    bool squashed = uStretch > 0.0 && speed > 1e-4;\n    if (squashed) {\n        ax = uVelocity / speed;\n        perp = vec2(-ax.y, ax.x);\n        sa = 1.0 + tanh(speed * 1.6) * uStretch;\n        bodyUv = vec2(dot(uv, ax) / sa, dot(uv, perp) * sa);\n    }\n\n    // Domain-warped FBM displacement WITH its analytic gradient (the warp-Jacobian\n    // approximation per [2]). disp = (noiseVal - 0.5) · uNoiseAmp; its gradient in\n    // bodyUv space is uNoiseAmp · uNoiseFreq · fbmGrad (the uNoiseFreq chain factor\n    // from the bodyUv·uNoiseFreq argument).\n    vec3 fbmv = fbmWarpedG(bodyUv * uNoiseFreq + uTime * uNoiseSpeed, 3, uWarpAmp);\n    float bodyDisplacement = (fbmv.x - 0.5) * uNoiseAmp;\n    vec2 dispGradBody = uNoiseAmp * uNoiseFreq * fbmv.yz;\n\n    // Body circle WITH gradient (in bodyUv space). The displaced radius subtracts\n    // the displacement, so the body gradient is circleGrad - dispGrad (sign per\n    // [2]: d = circleDist - displacement ⇒ grad = circleGrad - dispGrad).\n    vec3 bodyG = sdgCircle(bodyUv, vec2(0.0), bodyR + bodyDisplacement);\n    vec2 bodyGradBody = bodyG.yz - dispGradBody;\n\n    // Transform the body gradient back to uv space: g_uv = R^T · D · g_body, i.e.\n    // scale by D = diag(1/sa, sa) in the (ax, perp) frame then reconstruct in world.\n    vec2 bodyGradUv;\n    if (squashed) {\n        vec2 scaled = vec2(bodyGradBody.x / sa, bodyGradBody.y * sa);\n        bodyGradUv = scaled.x * ax + scaled.y * perp;\n    } else {\n        bodyGradUv = bodyGradBody;\n    }\n    vec3 d = vec3(bodyG.x, bodyGradUv);\n\n    // Satellites — smin-merged into the body (uMerge selects quadratic/circular).\n    for (int i = 0; i < MAX_SATS; i++) {\n        if (i >= uSatCount) break;\n        vec3 satG = sdgCircle(uv, uSatPos[i], uSatRadius[i]);\n        satG.x += (1.0 - uSatOpacity[i]) * 0.3;\n        d = sminG(d, satG, uSmoothK);\n    }\n\n    // Pointer trail — a decaying-radius pseudopod reaching toward the cursor,\n    // smin-merged so it stretches an elastic limb and snaps back.\n    for (int i = 0; i < TRAIL_N; i++) {\n        if (i >= uTrailCount) break;\n        vec3 trailG = sdgCircle(uv, uTrailPos[i], uTrailRadius[i]);\n        d = sminG(d, trailG, uSmoothK);\n    }\n    return d;\n}\n\n// AX.W15 — the analytic surface normal. The field gradient grad2d arrives\n// DIRECTLY from sceneDistG (no 4-tap). The smin field is sub-unit (the CD family\n// |grad| ≤ 1) so the gradient is NOT unit-length — normalize() before the dome\n// lift (the unit contract is on the FINAL lifted normal N, not the raw field\n// gradient). The Z dome: 'interior' is the normalized depth inward from the rim\n// (-d/bodyR); z = sqrt(1 - (1 - interior)^2) is the unit half-sphere profile (flat\n// centre, steep rim). N tilts outward along the field gradient near the rim and\n// points at the viewer deep inside.\nvec3 surfaceNormalFromGrad(vec2 grad2d, float d, float bodyR) {\n    vec2 g = normalize(grad2d + vec2(1e-6));\n    float interior = clamp(-d / max(bodyR, 1e-4), 0.0, 1.0);\n    float z = sqrt(max(0.0, 1.0 - (1.0 - interior) * (1.0 - interior)));\n    return normalize(vec3(g * (1.0 - z), z) + vec3(0.0, 0.0, 1e-6));\n}\n\n// BC.W-GOOBLOB-MEATBALL — the 2D SDF soft-shadow march (IQ rmshadows improved-penumbra,\n// research/viz/goo-blob.md §2.3, §5). The GLSL twin of the WGSL softShadow2D — a procedural\n// soft contact shadow FOLLOWING the irregular metaball silhouette (NOT a hard disc/box\n// shadow). Marches from ro along rd (the in-plane projection of uLightDir),\n// accumulating the closest miss via the Aaltonen penumbra (y = h*h/(2*ph) + d =\n// sqrt(h*h - y*y)) that kills the banding the naive res = min(res, h/(w*t)) shows. CHEAP:\n// 24 steps, re-uses the SAME sceneDistG().x field (NO separate field re-eval, NO FBO/\n// multi-pass — the procedural field needs none). w is the penumbra hardness (inverse\n// light-source size -> uShadowSoftness). Derivative-FREE. Returns 1.0 (lit) -> 0.0 (occluded).\nfloat softShadow2D(vec2 ro, vec2 rd, float mint, float maxt, float w) {\n    float res = 1.0;\n    float ph = 1e20;\n    float t = mint;\n    for (int i = 0; i < 24; i++) {\n        if (t >= maxt) break;\n        float h = sceneDistG(ro + rd * t).x;\n        if (h < 0.001) return 0.0;\n        float y = h * h / (2.0 * ph);\n        float dd = sqrt(max(h * h - y * y, 0.0));\n        res = min(res, dd / (w * max(0.001, t - y)));\n        ph = h;\n        t += h;\n    }\n    return clamp(res, 0.0, 1.0);\n}\n\n// W11.b — sample the multi-stop palette at t in [0,1], interpolating adjacent\n// stops in OKLab with a MIDPOINT CHROMA-BUMP (a linear OKLab mix of a vivid pair\n// dips chroma through grey; the bump lifts it back). Falls back to uBaseColor when\n// uStopCount <= 1. Returns an OKLCh stop [L, C, h(rad)].\nvec3 samplePaletteOklch(float t) {\n    if (uStopCount <= 1) return oklabToOklch(srgbToOklab(uBaseColor));\n    float n = float(uStopCount);\n    float ft = clamp(t, 0.0, 1.0) * (n - 1.0);\n    int i0 = int(floor(ft));\n    int i1 = min(i0 + 1, uStopCount - 1);\n    float f = ft - float(i0);\n    // Fetch adjacent stops (uPalette is gamma sRGB) → OKLab.\n    vec3 labA = srgbToOklab(uPalette[i0]);\n    vec3 labB = srgbToOklab(uPalette[i1]);\n    vec3 lab = mix(labA, labB, f);\n    vec3 lch = oklabToOklch(lab);\n    // Midpoint chroma-bump: a bell (peaks at f=0.5) lifts C off the grey midpoint.\n    lch.y += 0.03 * sin(PI * f);\n    return lch;\n}\n\nvoid main() {\n    vec2 uv = vUv - 0.5;\n\n    // AX.W16 (arm 5) — PRE-FBM bounding early-out. The oversized canvas (1.6x the\n    // wrapper) runs the full fragment ALU — two 3-octave FBM evals + the OKLCh\n    // round-trip + the lit/iridescence/SSS block — on a large transparent border\n    // otherwise (~60% of the canvas is outside the droplet reach). Any fragment beyond\n    // uMaxReach (the PADDED worst-case painted reach, uploaded so it never clips the\n    // wet meniscus) writes transparent and returns BEFORE any of that work. This is a\n    // transparent WRITE, NOT a GLSL discard: the pass is a pure premultiplied blend\n    // with no depth, and discard disables the tiled-renderer fast path on mobile GPUs.\n    // No fwidth is computed before the return, so there is no\n    // derivative-in-non-uniform-control-flow hazard (the kept pixels compute their own\n    // fwidth normally below). The pad covers the pointer-lean uv shift, so testing the\n    // pre-deform uv is safe.\n    if (dot(uv, uv) > uMaxReach * uMaxReach) {\n        fragColor = vec4(0.0);\n        return;\n    }\n\n    // Pointer deformation — honor the SIGN of uPointerAttraction (W10): a positive\n    // attraction leans the body IN toward the cursor, a negative shies it AWAY. The\n    // signed influence flows straight into the UV shift (no hardcoded repulsion).\n    //\n    // AY.W-BLOB-CONFIG D2 — the SIGN was INVERTED in the rendered result. The\n    // sample-shift sign here was uv -= dir*influence; the live readback (an attraction\n    // sweep over a HELD pointer) measured the body shifting MORE toward the cursor as\n    // attraction went MORE NEGATIVE (a=-1 -> +0.081, a=0 -> +0.063, a=+1 -> +0.042) —\n    // exactly inverted: a shy-away -1 LUNGED at the cursor, a lean-in +1 pulled back.\n    // pointerDir = uPointer - uv is the toward-pointer direction in SAMPLE space; to\n    // MOVE the rendered body TOWARD the pointer the sample must shift ALONG pointerDir\n    // (sample the field from where the pointer is), so the sign is uv += dir*influence.\n    // A positive attraction now shifts the body toward the cursor (leans in), a negative\n    // shifts it away (shies) — the sign the comment + the demo slider always promised,\n    // now true in the render.\n    //\n    // AX.W46 D5 — the falloff radius NARROWS 0.65 → 0.5 (the calm-lean reconciliation).\n    // The W15 REDRESS widened it to 0.65 to clear a synthetic floor; the live π-lane\n    // read the result as a LUNGE. The drama lives in the STRENGTH (pointerStrength,\n    // dropped to 0.18 in types.ts), NOT the falloff — the falloff's job is purely to keep\n    // the lean COHERENT across the creature so the body, satellites and trail all tilt\n    // toward the cursor as ONE. At 0.5 the body (sitting ≈0.2 uv from the cursor on a\n    // hover-flick) is still well inside the cutoff, so the whole creature tilts as one —\n    // a gentle, coherent \"it notices you\" lean rather than the body lunging a body-width.\n    if (uPointerActive > 0.5) {\n        vec2 pointerDir = uPointer - uv;\n        float pointerDist = length(pointerDir);\n        float influence = smoothstep(0.5, 0.0, pointerDist) * uPointerAttraction * uPointerStrength;\n        uv += normalize(pointerDir + 1e-6) * influence;\n    }\n\n    // The de-synced pulsing body radius — also reused below for the inner-glow scale.\n    float bodyR = uBodyRadius + breath(uPulsePhase) * uPulseAmp;\n\n    // Composite domain-warped membrane field WITH its analytic gradient (AX.W15).\n    vec3 scene = sceneDistG(uv);\n    float d = scene.x;\n    vec2 fieldGrad = scene.yz;\n\n    // Edit #1 — fwidth-based anti-aliased edge: derive the smoothstep half-width\n    // from the SDF screen-space gradient so the edge stays ~1px regardless of zoom.\n    float aa = max(fwidth(d), 1e-6);\n    float alpha = 1.0 - smoothstep(-aa, aa, d);\n\n    if (alpha < 0.001) {\n        fragColor = vec4(0.0);\n        return;\n    }\n\n    // Edit #4 — perceptually-uniform OKLCh variation. Gamma base → linear → OKLab\n    // → OKLCh; perturb L/C/h; OKLCh → OKLab → linear; hue-preserving gamut clamp.\n    // The base is the multi-stop palette sample (W11.b) — distributed across the\n    // body/satellites by the color noise field — or uBaseColor when single-stop.\n    float colorNoise = fbm(uv * uColorNoiseFreq + uTime * uColorNoiseSpeed, 3);\n    vec3 oklch = samplePaletteOklch(colorNoise);\n\n    oklch.z += (colorNoise - 0.5) * uHueRange * (PI / 180.0);     // hue swing, radians\n    oklch.y = max(oklch.y + (colorNoise - 0.5) * uSatShift, 0.0); // chroma swing\n    oklch.x = clamp(oklch.x + uBrightnessShift, 0.0, 1.0);        // lightness bias\n\n    // ── BD.W-GOO-CAROUSEL-DECK — the blob to meatball SHADING MORPH (uMorphT). ──\n    // The user #1 ask: the goo should MORPH BLOB and MEATBALL from one to another\n    // — a CONTINUOUS in-between, not a hard variant cut. The KEY FACT: the body\n    // GEOMETRY is IDENTICAL across blob/meatball (the smin SDF field, the satellites,\n    // the fwidth-AA — all computed ABOVE, before this gate). ONLY the SURFACE shading\n    // (flat fill vs lit/shadow/iridescence/SSS) differs. So the morph is a uMorphT\n    // scalar LERPING the SHADING — NOT a geometry rebuild, NOT a path interpolation.\n    //\n    //   uMorphT = 0  -> the flat warm-cream blob (the STAGE-1 floor, byte-identical to\n    //                   the old uStage > 0.5 early-return);\n    //   uMorphT = 1  -> the lit meatball (the full STAGE-2 dressing, byte-identical to\n    //                   the old else path);\n    //   0 < uMorphT < 1 -> the CONTINUOUS in-between (the morph the user asked for).\n    //\n    // uMorphT resolves from variant for back-compat (blob -> 0, meatball -> 1); a\n    // consumer ANIMATING config.morphT 0..1 gets the live morph. The smin field is\n    // byte-UNTOUCHED (the geometry is shared — that is the whole point).\n    //\n    // The flat color is computed ALWAYS (cheap); the EXPENSIVE dressing (the soft-shadow\n    // secondary march) is gated uMorphT > 0.0 so a pure blob pays ZERO shadow-march\n    // cost (the SLOW-fix discipline — a flat blob is as cheap as the old early-return).\n    float morphT = clamp(uMorphT, 0.0, 1.0);\n\n    // STAGE-1 flat fill — the warm-cream blob color, off the PRE-dressing oklch (the\n    // dressing mutates oklch below, so the flat read must snapshot it here).\n    vec3 flatOkl = gamutClampOklch(oklch);\n    vec3 flatLin = oklabToLinearSrgb(oklchToOklab(flatOkl));\n    vec3 flatRgb = clamp(linearToSrgb(flatLin), 0.0, 1.0);\n\n    // If we are fully flat (a pure blob, morphT == 0), emit the flat fill + dither and\n    // return — the byte-identical STAGE-1 floor, zero dressing cost.\n    if (morphT <= 0.0) {\n        float ign1 = fract(52.9829189 * fract(dot(gl_FragCoord.xy, vec2(0.06711056, 0.00583715))));\n        flatRgb = clamp(flatRgb + (ign1 - 0.5) / 255.0, 0.0, 1.0);\n        fragColor = vec4(flatRgb * alpha, alpha);\n        return;\n    }\n\n    // Surface normal — from the ANALYTIC field gradient (AX.W15: the 4-tap is\n    // deleted). Computed ONCE here and reused by the iridescence (W11.a), the\n    // fake-SSS (W11.a), and the lit glass block (W9.b).\n    vec3 N = surfaceNormalFromGrad(fieldGrad, d, bodyR);\n    vec3 V = vec3(0.0, 0.0, 1.0);\n    float thickness = clamp(-d / max(bodyR, 1e-4), 0.0, 1.0); // 0 at rim, ~1 deep in\n    float fres = pow(1.0 - max(dot(N, V), 0.0), 2.5);          // rim-weighted angle\n\n    // ── W11.a iridescence — warm-biased IQ cosine palette driving OKLCh HUE ──\n    //\n    // A thin-film-like sheen on the RIM (fres-weighted), NOT the cold-blue default:\n    // an IQ cosine palette (a + b*cos(2pi(c*t + d))) maps the Fresnel/edge angle +\n    // the FBM color field + an animated thickness onto a HUE OFFSET biased to the\n    // warm arc (uIridHue), with the chroma clamped to a warm-pearl band so the rim\n    // reads cream/gold (congruent with the warm-cream-glass identity). Mixed into\n    // OKLCh BEFORE the gamut clamp. Default LOW.\n    if (uIridescence > 0.0) {\n        float t = fres + 0.3 * colorNoise + uTime * uIridSpeed;\n        // Cosine palette in HUE space, centred on the warm uIridHue (~0.18 turn of\n        // warm sweep so the sheen shimmers gold→peach, never swinging to cold blue).\n        float iridHue = uIridHue + 0.18 * PI * cos(2.0 * PI * t);\n        float w = fres * uIridescence;          // rim-weighted blend\n        oklch.z = mix(oklch.z, iridHue, w);\n        // Warm-pearl chroma cap: lift chroma toward a pale-pearl band, clamped.\n        oklch.y = min(oklch.y + 0.04 * w, oklch.y + 0.08);\n        oklch.x = min(oklch.x + 0.05 * w, 1.0); // a touch brighter at the sheen\n    }\n\n    // ── W11.a fake subsurface translucency — thickness inner-glow + back-light ──\n    //\n    // A bright translucent CORE fading to a light-leaking EDGE (the -d ramp) plus\n    // the fast-SSS back-light (light wrapping through the thin rim). Both lift OKLCh\n    // L and warm the hue, consuming the W9 normal. In OKLCh before the gamut clamp.\n    if (uCoreGlow > 0.0 || uSssScale > 0.0) {\n        // Beer-Lambert inner luminosity (AX.W15): a SATURATING 1 - exp(-k·thickness)\n        // curve — a flat thick core with a fast warm rim falloff — NOT the linear\n        // coreGlow·thickness ramp (which over-brightens deep interiors). k ≈ 3 reads\n        // as glass depth.\n        oklch.x = min(oklch.x + uCoreGlow * (1.0 - exp(-3.0 * thickness)), 1.0);\n        // Fast-SSS back-light: light wrapping through the thin (low-thickness) rim.\n        vec3 L = normalize(uLightDir + vec3(1e-6));\n        float back = pow(clamp(dot(V, -(L + N * thickness)), 0.0, 1.0), uSssPower);\n        float sss = back * uSssScale * (1.0 - thickness);\n        oklch.x = min(oklch.x + sss, 1.0);\n        // Warm the THIN leaking rim only (scale by 1 - thickness): the SSS hue-warm\n        // shift rides the rim, not the core (AX.W15 — only the thin rim warms).\n        oklch.z += sss * 0.1 * (1.0 - thickness);\n    }\n\n    // ── BC.W-GOOBLOB-MEATBALL — the soft contact shadow (T2, the GLSL twin). A procedural\n    //    soft shadow FOLLOWING the irregular silhouette: march from this fragment's uv\n    //    toward the in-plane projection of uLightDir; where the metaball field occludes the\n    //    light, DARKEN the OKLCh L — a soft grounded contact band under the dome\n    //    (PROCEDURAL, silhouette-following, unlike the static CSS gel-dome drop-shadow).\n    //    Rim-weighted (1 - thickness) so the lit interior stays bright. Gated uShadow > 0.5;\n    //    variant=blob (STAGE 1) ships shadowless via the uStage early-return above. ──\n    if (uShadow > 0.5) {\n        vec2 L2 = normalize(uLightDir.xy + vec2(1e-6));\n        float sh = softShadow2D(uv, L2, max(aa, 0.004), bodyR * 1.5, max(uShadowSoftness, 4.0));\n        // The contact band sits on the LOWER rim (the side AWAY from the light), gated by\n        // the away-facing rim weight + a thin-rim Gaussian band (peak thickness ≈ 0.18) so\n        // the deep-lit interior + the light-facing edge stay BRIGHT (the warm-cream body\n        // mean + the proof:blob-warm-default floor hold). The 0.20 cap is a soft grounding.\n        float away = clamp(-dot(normalize(fieldGrad + vec2(1e-6)), L2), 0.0, 1.0);\n        float band = exp(-pow((thickness - 0.18) / 0.16, 2.0));\n        float shadowDarken = (1.0 - sh) * away * band * 0.20;\n        oklch.x = max(oklch.x - shadowDarken, 0.0);\n    }\n\n    oklch = gamutClampOklch(oklch);\n\n    vec3 lin = oklabToLinearSrgb(oklchToOklab(oklch));\n\n    // ── W9.b lit glass surface — Blinn-Phong glint + Fresnel rim, in LINEAR ──\n    //\n    // The lit terms enter lin in LINEAR space BEFORE the OETF (linearToSrgb)\n    // and BEFORE the * alpha premultiply — a post-OETF apply double-gammas the\n    // highlight and fringes the premultiplied edge (the named A5/A2 trap). The\n    // prior flat inner-glow lift is SUBSUMED by the Fresnel rim (the curvature\n    // read now comes from the surface normal, not a flat smoothstep). Gated behind\n    // uLit so the flat fill stays the default (zero regression).\n    if (uLit > 0.5) {\n        // N, V, thickness are already computed in main() (reused by iridescence/SSS).\n        vec3 L = normalize(uLightDir + vec3(1e-6)); // light direction\n        vec3 H = normalize(L + V);                  // Blinn-Phong half-vector\n\n        // Warm-cream specular glint — a near-white OKLCh TINT (L~0.97, C~0.03,\n        // hue~85°) routed through the SAME spliced OKLCh matrices, NOT hardcoded\n        // sRGB white (pure white reads cheap-CG). Linearized for the linear add.\n        vec3 warmCream = oklabToLinearSrgb(oklchToOklab(vec3(0.97, 0.03, radians(85.0))));\n\n        // ENERGY-CONSERVING Blinn-Phong (AX.W15 [9]): the (shininess+2)/8 factor\n        // DECOUPLES shininess from strength (without it, raising shininess dims the\n        // glint, so the two knobs fight). SPECULAR ANTIALIASING: the FBM membrane\n        // makes the normal vary fast, so a tight glint (16-64) STROBES on small\n        // dock-grid blobs. Widen the effective lobe where the normal varies\n        // (Toksvig-style fwidth clamp) — drop shininess, NOT raise it. nVar reads\n        // the screen-space normal derivative; a high-variance pixel softens the\n        // exponent so the glint stays stable.\n        float nVar = length(fwidth(N));\n        float shininess = uSpecShininess / (1.0 + 24.0 * nVar);\n        float energyNorm = (shininess + 2.0) / 8.0;\n        float spec = pow(max(dot(N, H), 0.0), shininess) * uSpecStrength * energyNorm;\n\n        // Fresnel/Schlick rim — fed uRimColor (the --foreground warm rim via the\n        // /color leaf). FOREGROUND-AWARE MIN-CONTRAST GUARD (AX.W15 [10]):\n        // a var(--primary) blob in dark mode resolves a rim near the BODY color, so\n        // the rim washes out (no contrast). When the rim's luminance sits too close\n        // to the body's, the guard CHROMA-REDUCES and L-LIFTS the rim stop (a\n        // perceptual move, NOT a re-tint) so the curved rim always reads. Computed in\n        // OKLCh so the lift is perceptual.\n        vec3 rimOkl = oklabToOklch(srgbToOklab(uRimColor));\n        float bodyL = oklch.x; // the post-perturb body lightness\n        float dL = abs(rimOkl.x - bodyL);\n        // If the rim is within 0.22 L of the body, push it AWAY (lift if the body is\n        // dark, the dark-mode primary case; the chroma-reduce keeps it neutral-warm).\n        float lack = clamp((0.22 - dL) / 0.22, 0.0, 1.0);\n        rimOkl.x = clamp(mix(rimOkl.x, (bodyL < 0.5 ? 0.92 : 0.18), lack), 0.0, 1.0);\n        rimOkl.y *= mix(1.0, 0.5, lack); // chroma-reduce toward a pale pearl\n        vec3 rimLin = oklabToLinearSrgb(oklchToOklab(rimOkl));\n\n        float rimFres = pow(1.0 - max(dot(N, V), 0.0), uRimPower);\n        float rim = rimFres * uRimStrength * (1.0 - 0.6 * thickness);\n\n        // Combine the two warm highlights and add in LINEAR. max(spec, rim*scale)\n        // keeps the glint from stacking on the rim into a blown hotspot.\n        //\n        // AX.W46 D4 — CLAMP the linear highlight below unity before the add (the\n        // belt-and-braces guard the re-derived specStrength makes rarely-binding but\n        // never lets a worst-case normal blow to a hard white spot). The energy-\n        // conserving Blinn-Phong CAN still spike on a grazing normal where the Toksvig\n        // widen has not kicked in; capping the highlight magnitude at 0.85 (a contained\n        // warm gleam, sub-unity) means the OETF clamp below can never crush a\n        // pure-white pixel over the dome. The cap is on the highlight ALONE (not the\n        // composited lin), so the body color underneath is untouched — it caps only\n        // the additive catch-light's energy.\n        vec3 highlight = max(warmCream * spec, rimLin * rim);\n        highlight = min(highlight, vec3(0.85));\n        lin += highlight;\n    }\n\n    vec3 rgb = clamp(linearToSrgb(lin), 0.0, 1.0); // MANDATORY OETF — closes the seam\n\n    // ── BD.W-GOO-CAROUSEL-DECK — the blob to meatball SHADING MORPH compose. ──\n    // rgb is the fully-dressed meatball surface; flatRgb (snapshotted before the\n    // dressing) is the flat blob fill. The morph is a CONTINUOUS lerp between them on\n    // morphT — at morphT==1 it is byte-identical to the dressed meatball, and the\n    // morphT==0 case already returned the flat fill above. The geometry (alpha/silhouette)\n    // is shared, so only the SURFACE shading interpolates (the morph from one to\n    // another). The dither below applies to the composited result.\n    rgb = mix(flatRgb, rgb, morphT);\n\n    // IGN dither (AX.W15 [2][9]) — the low-chroma warm-cream dome BANDS on 8-bit\n    // panels (visible Mach steps in the smooth L roll). Interleaved-gradient noise at\n    // 1/255, applied AFTER linearToSrgb and BEFORE the *alpha premultiply (aurora\n    // ships the same splice). A triangular ±0.5 LSB dither decorrelates the steps.\n    float ign = fract(52.9829189 * fract(dot(gl_FragCoord.xy, vec2(0.06711056, 0.00583715))));\n    rgb = clamp(rgb + (ign - 0.5) / 255.0, 0.0, 1.0);\n\n    // Edit #8 — premultiply AFTER the OETF: straight-alpha gamma → premultiplied.\n    fragColor = vec4(rgb * alpha, alpha);\n}";
 //#endregion
 //#region src/components/custom/goo-blob/composables/uploadBlobUniforms.ts
-function re(e, t, n, r, i, a, o, s, c) {
+function se(e, t, n, r, i, a, o, s, c) {
 	let { U: l, satPosLocs: u, satRadLocs: d, satOpLocs: f, trailPosLocs: p, trailRadLocs: m, paletteLocs: g } = r, { params: _, rgb: v, simTimeMs: y, timeSec: b, resolveColor: x, rimColor: S, paletteStops: C } = c, { geometry: w, membrane: T, color: E, surface: D, interaction: O } = a;
 	e.useProgram(t), e.bindVertexArray(n), e.uniform2f(l.uResolution, i.width, i.height), e.uniform1f(l.uTime, y / 1e3), e.uniform3f(l.uBaseColor, v[0], v[1], v[2]);
 	let k = C, A = Math.min(k.length, 4);
@@ -543,38 +573,40 @@ function re(e, t, n, r, i, a, o, s, c) {
 			e.uniform3f(n, t[0], t[1], t[2]);
 		} else e.uniform3f(n, v[0], v[1], v[2]);
 	}
-	let j = o.pointer.value, M = O.pointerAttraction + _.pointerAttraction;
-	e.uniform2f(l.uPointer, j.x * .5 * h, j.y * .5 * h), e.uniform1f(l.uPointerActive, +!!o.active.value), e.uniform1f(l.uPointerAttraction, M), e.uniform1f(l.uPointerStrength, O.pointerStrength * h);
-	let N = Math.max(0, Math.min(1, M)), P = o.velocity.value;
-	e.uniform2f(l.uVelocity, P.x * .5 * h, P.y * .5 * h), e.uniform1f(l.uStretch, O.stretch);
-	let F = o.trailSources(w.satelliteRadius * .7 * N);
-	e.uniform1i(l.uTrailCount, F.count);
+	let M = o.pointer.value, N = O.pointerAttraction + _.pointerAttraction;
+	e.uniform2f(l.uPointer, M.x * .5 * h, M.y * .5 * h), e.uniform1f(l.uPointerActive, +!!o.active.value), e.uniform1f(l.uPointerAttraction, N), e.uniform1f(l.uPointerStrength, O.pointerStrength * h);
+	let P = Math.max(0, Math.min(1, N)), F = o.velocity.value;
+	e.uniform2f(l.uVelocity, F.x * .5 * h, F.y * .5 * h), e.uniform1f(l.uStretch, O.stretch);
+	let I = o.trailSources(w.satelliteRadius * .7 * P);
+	e.uniform1i(l.uTrailCount, I.count);
 	for (let t = 0; t < 15; t++) {
-		let n = p[t] ?? null, r = m[t] ?? null, i = F.sources[t];
+		let n = p[t] ?? null, r = m[t] ?? null, i = I.sources[t];
 		i ? (e.uniform2f(n, i.x * .5 * h, i.y * .5 * h), e.uniform1f(r, i.radius * h)) : (e.uniform2f(n, 0, 0), e.uniform1f(r, 0));
 	}
 	e.uniform1f(l.uBodyRadius, w.bodyRadius * h), e.uniform1f(l.uPulsePhase, b * T.pulseFreq * _.pulseFreq * Math.PI * 2);
-	let I = o.pulse.value * O.clickImpulse;
-	e.uniform1f(l.uPulseAmp, (T.pulseAmp * _.pulseAmp / .015 + I) * h), e.uniform1f(l.uNoiseAmp, T.noiseAmp * _.noiseAmp / .025 * h), e.uniform1f(l.uNoiseFreq, T.noiseFreq), e.uniform1f(l.uNoiseSpeed, T.noiseSpeed), e.uniform1f(l.uWarpAmp, T.warpAmp);
-	let L = w.orbitRadius * 1.2 * (1 + w.eccentricity), R = T.smoothK * _.smoothK, z = Math.max(0, L - w.satelliteRadius - w.bodyRadius), B = R > 0 ? Math.min(1.25, 1 + z / R) : 1;
-	e.uniform1f(l.uSmoothK, R * B * h), e.uniform1f(l.uMerge, +(T.merge === "circular")), e.uniform1f(l.uStage, +(a.variant === "blob"));
-	let V = L + w.satelliteRadius, H = (Math.max(w.bodyRadius, V) + R * B + T.noiseAmp * _.noiseAmp / .025 + Math.abs(o.pulse.value) * O.clickImpulse) * h + .1;
-	e.uniform1f(l.uMaxReach, H), e.uniform1f(l.uHueRange, E.hueRange + _.hueRange), e.uniform1f(l.uSatShift, E.satShift + _.satShift), e.uniform1f(l.uBrightnessShift, E.brightnessShift + _.brightnessShift), e.uniform1f(l.uColorNoiseFreq, E.colorNoiseFreq), e.uniform1f(l.uColorNoiseSpeed, E.colorNoiseSpeed);
-	let U = a.variant !== "blob";
-	e.uniform1f(l.uLit, U && D.lit ? 1 : 0);
-	let W = x(S);
-	e.uniform3f(l.uRimColor, W[0], W[1], W[2]), e.uniform3f(l.uLightDir, D.lightDir[0], D.lightDir[1], D.lightDir[2]), e.uniform1f(l.uSpecStrength, D.specStrength), e.uniform1f(l.uSpecShininess, D.specShininess), e.uniform1f(l.uRimPower, D.rimPower), e.uniform1f(l.uRimStrength, D.rimStrength), e.uniform1f(l.uShadow, U && D.shadow ? 1 : 0), e.uniform1f(l.uShadowSoftness, D.shadowSoftness), e.uniform1f(l.uIridescence, D.iridescence * _.iridScale), e.uniform1f(l.uIridHue, D.iridHue * (Math.PI / 180)), e.uniform1f(l.uIridSpeed, D.iridSpeed), e.uniform1f(l.uSssScale, D.sssScale * _.iridScale), e.uniform1f(l.uSssPower, D.sssPower), e.uniform1f(l.uCoreGlow, D.coreGlow);
-	let G = s.sources;
-	e.uniform1i(l.uSatCount, G.length);
+	let L = o.pulse.value * O.clickImpulse;
+	e.uniform1f(l.uPulseAmp, (T.pulseAmp * _.pulseAmp / .015 + L) * h), e.uniform1f(l.uNoiseAmp, T.noiseAmp * _.noiseAmp / .025 * h), e.uniform1f(l.uNoiseFreq, T.noiseFreq), e.uniform1f(l.uNoiseSpeed, T.noiseSpeed), e.uniform1f(l.uWarpAmp, T.warpAmp);
+	let R = w.orbitRadius * 1.2 * (1 + w.eccentricity), z = T.smoothK * _.smoothK, B = Math.max(0, R - w.satelliteRadius - w.bodyRadius), V = z > 0 ? Math.min(1.25, 1 + B / z) : 1;
+	e.uniform1f(l.uSmoothK, z * V * h), e.uniform1f(l.uMerge, +(T.merge === "circular"));
+	let H = typeof a.morphT == "number" ? Math.max(0, Math.min(1, a.morphT)) : a.variant === "blob" ? 0 : 1;
+	e.uniform1f(l.uMorphT, H), e.uniform1f(l.uStage, +(H <= 0));
+	let U = (a.surface.fissionAmp ?? 0) > 0 ? j * 1.12 : 0, W = Math.max(R, U) + w.satelliteRadius, G = (Math.max(w.bodyRadius, W) + z * V + T.noiseAmp * _.noiseAmp / .025 + Math.abs(o.pulse.value) * O.clickImpulse) * h + .1;
+	e.uniform1f(l.uMaxReach, G), e.uniform1f(l.uHueRange, E.hueRange + _.hueRange), e.uniform1f(l.uSatShift, E.satShift + _.satShift), e.uniform1f(l.uBrightnessShift, E.brightnessShift + _.brightnessShift), e.uniform1f(l.uColorNoiseFreq, E.colorNoiseFreq), e.uniform1f(l.uColorNoiseSpeed, E.colorNoiseSpeed);
+	let K = H > 0;
+	e.uniform1f(l.uLit, K && D.lit ? 1 : 0);
+	let q = x(S);
+	e.uniform3f(l.uRimColor, q[0], q[1], q[2]), e.uniform3f(l.uLightDir, D.lightDir[0], D.lightDir[1], D.lightDir[2]), e.uniform1f(l.uSpecStrength, D.specStrength), e.uniform1f(l.uSpecShininess, D.specShininess), e.uniform1f(l.uRimPower, D.rimPower), e.uniform1f(l.uRimStrength, D.rimStrength), e.uniform1f(l.uShadow, K && D.shadow ? 1 : 0), e.uniform1f(l.uShadowSoftness, D.shadowSoftness), e.uniform1f(l.uIridescence, D.iridescence * _.iridScale), e.uniform1f(l.uIridHue, D.iridHue * (Math.PI / 180)), e.uniform1f(l.uIridSpeed, D.iridSpeed), e.uniform1f(l.uSssScale, D.sssScale * _.iridScale), e.uniform1f(l.uSssPower, D.sssPower), e.uniform1f(l.uCoreGlow, D.coreGlow);
+	let J = s.sources;
+	e.uniform1i(l.uSatCount, J.length);
 	for (let t = 0; t < 4; t++) {
-		let n = u[t] ?? null, r = d[t] ?? null, i = f[t] ?? null, a = G[t];
+		let n = u[t] ?? null, r = d[t] ?? null, i = f[t] ?? null, a = J[t];
 		a ? (e.uniform2f(n, a.x * h, a.y * h), e.uniform1f(r, a.radius * h), e.uniform1f(i, a.opacity)) : (e.uniform2f(n, 0, 0), e.uniform1f(r, 0), e.uniform1f(i, 0));
 	}
 	e.clearColor(0, 0, 0, 0), e.clear(e.COLOR_BUFFER_BIT), e.drawArrays(e.TRIANGLES, 0, 6), e.bindVertexArray(null);
 }
 //#endregion
 //#region src/components/custom/goo-blob/shaders/metaball.wgsl.ts
-var ie = `
+var ce = `
 const MAX_SATS: i32 = 4;
 const TRAIL_N: i32 = 15;
 const MAX_BLOB_STOPS: i32 = 4;
@@ -628,7 +660,10 @@ struct Uniforms {
 
 @group(0) @binding(0) var<uniform> u: Uniforms;
 
-// Noise + FBM (BC.W-CARVE6 → metaball-noise.wgsl) — spliced AFTER ${a}.
+// Noise + FBM (BC.W-CARVE6 -> metaball-noise.wgsl) — spliced AFTER the FBM_ROT_WGSL chunk
+// below. NB: do NOT name the chunk with a template-interpolation form in a WGSL comment —
+// it injects the whole multi-line chunk into this // comment line and breaks the parse
+// (BD.W-VIZ-BROKEN-FIX D3c).
 ${a}
 
 const FBM_LACUNARITY: f32 = 1.8;
@@ -754,8 +789,12 @@ fn sminG(a: vec3<f32>, b: vec3<f32>, k: f32) -> vec3<f32> {
 ${t}
 ${r}
 
-// OKLCh gamut clamp + palette sample + breath (BC.W-CARVE6 → metaball-palette.wgsl)
-// — spliced AFTER ${r} (its matrix source).
+// OKLCh gamut clamp + palette sample + breath (BC.W-CARVE6 -> metaball-palette.wgsl)
+// — spliced AFTER the OKLCH_MATRICES_WGSL chunk above (its matrix source). NB: do NOT
+// write the chunk's identifier with a template-interpolation form in a WGSL comment — it
+// would inject the whole multi-line matrices chunk INTO this // comment line and break the
+// parse at the first injected newline (BD.W-VIZ-BROKEN-FIX D3c — the live WGSL syntax error
+// that silently fell every metaball/goo-dot to the WebGL2 net, masking the WGSL primary).
 
 // ── OKLCh gamut clamp (hue-preserving inward chroma clamp — blob-local) ──
 fn inGamut(lin: vec3<f32>) -> bool {
@@ -933,6 +972,10 @@ fn fs_main(in: VSOut) -> @location(0) vec4<f32> {
   let uPointerStrength = u.s7.z;
   // BC.W-GOOBLOB-PLAIN — the STAGE-1 gate (1.0 = the plain shadowless fill-only floor).
   let uStage = u.s7.w;
+  // BD.W-GOO-CAROUSEL-DECK — the blob to meatball SHADING-MORPH scalar (the base.w pad
+  // lane; the typed-struct SoT writes it). 0 = flat blob, 1 = lit meatball, in-between
+  // lerps the surface shading (the smin geometry is shared).
+  let uMorphT = clamp(u.base.w, 0.0, 1.0);
   let uPointer = u.ptr.xy;
   let uHueRange = u.s5.z;
   let uSatShift = u.s5.w;
@@ -957,10 +1000,19 @@ fn fs_main(in: VSOut) -> @location(0) vec4<f32> {
 
   var uv = in.uv - 0.5;
 
-  // PRE-FBM bounding early-out — transparent WRITE (not a discard).
-  if (dot(uv, uv) > uMaxReach * uMaxReach) {
-    return vec4<f32>(0.0);
-  }
+  // PRE-FBM bounding test (BD.W-VIZ-BROKEN-FIX D3c — the uniformity fix). The GLSL fallback
+  // EARLY-RETURNS transparent here (a perf skip of the FBM border); WGSL CANNOT — an early
+  // return predicated on the per-fragment uv makes the fwidth(d)/fwidth(Nh) calls below
+  // reachable only through non-uniform control flow, which the WGSL uniformity analyzer
+  // REJECTS (fwidth must only be called from uniform control flow) so the metaball WGSL
+  // module never compiled and every GooBlob + goo-dot fell silently to the WebGL2 net (the
+  // "broken TOTALLY" the user reports, the un-exercised Safari-primary surface). So the
+  // bound is folded into the FINAL alpha (a flag, NOT a control-flow return): the derivative
+  // builtins stay in UNIFORM control flow, and an outside fragment writes transparent at the
+  // end — gestalt-identical to the GLSL early-out, just without the WGSL-illegal pre-
+  // derivative return. sceneDistG is derivative-FREE so running it on the border costs only
+  // ALU (no correctness hazard); the real GPU swallows it.
+  let inBounds = dot(uv, uv) <= uMaxReach * uMaxReach;
 
   // Pointer deformation — sign per AY.W-BLOB-CONFIG D2.
   if (uPointerActive > 0.5) {
@@ -978,9 +1030,12 @@ fn fs_main(in: VSOut) -> @location(0) vec4<f32> {
   let fieldGrad = scene.yz;
 
   // ── fwidth SITE #1 (metaball.frag.ts line 266) — AA-edge half-width from the SDF
-  //    screen-space gradient. The WGSL fragment-stage fwidth() transcription. ──
+  //    screen-space gradient. The WGSL fragment-stage fwidth() transcription. In UNIFORM
+  //    control flow (no early return precedes it — the bounding test above is a flag). ──
   let aa = max(fwidth(d), 1e-6);
-  let alpha = 1.0 - smoothstep(-aa, aa, d);
+  // The bounding test (D3c) folds in here: an out-of-bounds fragment is forced transparent
+  // (the GLSL early-out's transparent WRITE, expressed as an alpha gate, not a return).
+  let alpha = (1.0 - smoothstep(-aa, aa, d)) * select(0.0, 1.0, inBounds);
 
   // ── BC.W-GOOBLOB-MEATBALL — the WGSL uniformity STRUCTURAL fix (the meatball-armer).
   //    The Toksvig screen-space derivative (fwidth SITE #2) is hoisted HERE, into the SAME
@@ -1027,13 +1082,16 @@ fn fs_main(in: VSOut) -> @location(0) vec4<f32> {
   // (the BB residual that kept STAGE 2 falling to the WebGL2 net is CLOSED). The
   // teaching contrast: this is the "it renders, it meatballs, it works on Safari"
   // floor STAGE 2 layers the lit/shadow onto via the SAME uniforms.
-  if (uStage > 0.5) {
-    let okl1 = gamutClampOklch(oklch);
-    var lin1 = oklabToLinearSrgb(oklchToOklab(okl1));
-    var rgb1 = clamp(linearToSrgb(lin1), vec3<f32>(0.0), vec3<f32>(1.0));
+  // BD.W-GOO-CAROUSEL-DECK — the flat blob fill is snapshotted from the PRE-dressing
+  // oklch (the dressing mutates oklch below). At morphT <= 0 (a pure blob) emit the
+  // flat fill + dither and return — the byte-identical STAGE-1 floor, zero dressing cost.
+  let flatOkl = gamutClampOklch(oklch);
+  let flatLin = oklabToLinearSrgb(oklchToOklab(flatOkl));
+  let flatRgbBase = clamp(linearToSrgb(flatLin), vec3<f32>(0.0), vec3<f32>(1.0));
+  if (uMorphT <= 0.0) {
     let fc1 = in.pos.xy;
     let ign1 = fract(52.9829189 * fract(dot(fc1, vec2<f32>(0.06711056, 0.00583715))));
-    rgb1 = clamp(rgb1 + (ign1 - 0.5) / 255.0, vec3<f32>(0.0), vec3<f32>(1.0));
+    let rgb1 = clamp(flatRgbBase + (ign1 - 0.5) / 255.0, vec3<f32>(0.0), vec3<f32>(1.0));
     return vec4<f32>(rgb1 * alpha, alpha);
   }
 
@@ -1137,6 +1195,12 @@ fn fs_main(in: VSOut) -> @location(0) vec4<f32> {
 
   var rgb = clamp(linearToSrgb(lin), vec3<f32>(0.0), vec3<f32>(1.0));
 
+  // BD.W-GOO-CAROUSEL-DECK — the blob↔meatball SHADING-MORPH compose: lerp the flat blob
+  // fill -> the dressed meatball surface on uMorphT (the geometry/alpha is shared; only
+  // the surface shading interpolates). morphT == 1 is byte-identical to the dressed
+  // meatball; the morphT == 0 case already returned the flat fill above.
+  rgb = mix(flatRgbBase, rgb, uMorphT);
+
   // IGN dither — 1-LSB triangular dither in DISPLAY space, AFTER linearToSrgb.
   let fragCoord = in.pos.xy;
   let ign = fract(52.9829189 * fract(dot(fragCoord, vec2<f32>(0.06711056, 0.00583715))));
@@ -1164,7 +1228,7 @@ fn fs_main(in: VSOut) -> @location(0) vec4<f32> {
 	satPos: 72,
 	trailPos: 88
 };
-function ae() {
+function le() {
 	let e = new ArrayBuffer(Q);
 	return {
 		buffer: e,
@@ -1172,33 +1236,33 @@ function ae() {
 		i32: new Int32Array(e)
 	};
 }
-function oe(e, t, n, r, i, a) {
-	let { f32: o, i32: s } = e, { params: c, rgb: l, simTimeMs: u, timeSec: d, resolveColor: f, rimColor: p, paletteStops: m } = a, { geometry: g, membrane: _, color: v, surface: y, interaction: b } = n, x = r.pointer.value, S = b.pointerAttraction + c.pointerAttraction, C = Math.max(0, Math.min(1, S)), w = r.velocity.value, T = r.pulse.value * b.clickImpulse, E = (_.pulseAmp * c.pulseAmp / .015 + T) * h, D = d * _.pulseFreq * c.pulseFreq * Math.PI * 2, O = _.noiseAmp * c.noiseAmp / .025 * h, k = g.orbitRadius * 1.2 * (1 + g.eccentricity), A = _.smoothK * c.smoothK, j = Math.max(0, k - g.satelliteRadius - g.bodyRadius), M = A > 0 ? Math.min(1.25, 1 + j / A) : 1, N = A * M * h, P = k + g.satelliteRadius, F = (Math.max(g.bodyRadius, P) + A * M + _.noiseAmp * c.noiseAmp / .025 + Math.abs(r.pulse.value) * b.clickImpulse) * h + .1;
+function ue(e, t, n, r, i, a) {
+	let { f32: o, i32: s } = e, { params: c, rgb: l, simTimeMs: u, timeSec: d, resolveColor: f, rimColor: p, paletteStops: m } = a, { geometry: g, membrane: _, color: v, surface: y, interaction: b } = n, x = r.pointer.value, S = b.pointerAttraction + c.pointerAttraction, C = Math.max(0, Math.min(1, S)), w = r.velocity.value, T = r.pulse.value * b.clickImpulse, E = (_.pulseAmp * c.pulseAmp / .015 + T) * h, D = d * _.pulseFreq * c.pulseFreq * Math.PI * 2, O = _.noiseAmp * c.noiseAmp / .025 * h, k = g.orbitRadius * 1.2 * (1 + g.eccentricity), A = _.smoothK * c.smoothK, M = Math.max(0, k - g.satelliteRadius - g.bodyRadius), N = A > 0 ? Math.min(1.25, 1 + M / A) : 1, P = A * N * h, F = (n.surface.fissionAmp ?? 0) > 0 ? j * 1.12 : 0, I = Math.max(k, F) + g.satelliteRadius, L = (Math.max(g.bodyRadius, I) + A * N + _.noiseAmp * c.noiseAmp / .025 + Math.abs(r.pulse.value) * b.clickImpulse) * h + .1;
 	o[$.s0 + 0] = u / 1e3, o[$.s0 + 1] = g.bodyRadius * h, o[$.s0 + 2] = D, o[$.s0 + 3] = E, o[$.s1 + 0] = O, o[$.s1 + 1] = _.noiseFreq, o[$.s1 + 2] = _.noiseSpeed, o[$.s1 + 3] = _.warpAmp;
-	let I = n.variant !== "blob";
-	o[$.s2 + 0] = N, o[$.s2 + 1] = +(_.merge === "circular"), o[$.s2 + 2] = F, o[$.s2 + 3] = I && y.lit ? 1 : 0, o[$.s3 + 0] = y.specStrength, o[$.s3 + 1] = y.specShininess, o[$.s3 + 2] = y.rimPower, o[$.s3 + 3] = y.rimStrength, o[$.s4 + 0] = y.iridescence * c.iridScale, o[$.s4 + 1] = y.iridHue * (Math.PI / 180), o[$.s4 + 2] = y.iridSpeed, o[$.s4 + 3] = y.sssScale * c.iridScale, o[$.s5 + 0] = y.sssPower, o[$.s5 + 1] = y.coreGlow, o[$.s5 + 2] = v.hueRange + c.hueRange, o[$.s5 + 3] = v.satShift + c.satShift, o[$.s6 + 0] = v.brightnessShift + c.brightnessShift, o[$.s6 + 1] = v.colorNoiseFreq, o[$.s6 + 2] = v.colorNoiseSpeed, o[$.s6 + 3] = b.stretch, o[$.s7 + 0] = +!!r.active.value, o[$.s7 + 1] = S, o[$.s7 + 2] = b.pointerStrength * h, o[$.s7 + 3] = +(n.variant === "blob"), o[$.ptr + 0] = x.x * .5 * h, o[$.ptr + 1] = x.y * .5 * h, o[$.ptr + 2] = w.x * .5 * h, o[$.ptr + 3] = w.y * .5 * h, o[$.base + 0] = l[0], o[$.base + 1] = l[1], o[$.base + 2] = l[2], o[$.base + 3] = 0;
-	let L = f(p);
-	o[$.rim + 0] = L[0], o[$.rim + 1] = L[1], o[$.rim + 2] = L[2], o[$.rim + 3] = 0, o[$.light + 0] = y.lightDir[0], o[$.light + 1] = y.lightDir[1], o[$.light + 2] = y.lightDir[2], o[$.light + 3] = 0, o[$.res + 0] = t.width, o[$.res + 1] = t.height, o[$.res + 2] = I && y.shadow ? 1 : 0, o[$.res + 3] = y.shadowSoftness;
-	let R = m, z = Math.min(R.length, 4);
+	let R = typeof n.morphT == "number" ? Math.max(0, Math.min(1, n.morphT)) : n.variant === "blob" ? 0 : 1, z = R > 0;
+	o[$.s2 + 0] = P, o[$.s2 + 1] = +(_.merge === "circular"), o[$.s2 + 2] = L, o[$.s2 + 3] = z && y.lit ? 1 : 0, o[$.s3 + 0] = y.specStrength, o[$.s3 + 1] = y.specShininess, o[$.s3 + 2] = y.rimPower, o[$.s3 + 3] = y.rimStrength, o[$.s4 + 0] = y.iridescence * c.iridScale, o[$.s4 + 1] = y.iridHue * (Math.PI / 180), o[$.s4 + 2] = y.iridSpeed, o[$.s4 + 3] = y.sssScale * c.iridScale, o[$.s5 + 0] = y.sssPower, o[$.s5 + 1] = y.coreGlow, o[$.s5 + 2] = v.hueRange + c.hueRange, o[$.s5 + 3] = v.satShift + c.satShift, o[$.s6 + 0] = v.brightnessShift + c.brightnessShift, o[$.s6 + 1] = v.colorNoiseFreq, o[$.s6 + 2] = v.colorNoiseSpeed, o[$.s6 + 3] = b.stretch, o[$.s7 + 0] = +!!r.active.value, o[$.s7 + 1] = S, o[$.s7 + 2] = b.pointerStrength * h, o[$.s7 + 3] = +(R <= 0), o[$.ptr + 0] = x.x * .5 * h, o[$.ptr + 1] = x.y * .5 * h, o[$.ptr + 2] = w.x * .5 * h, o[$.ptr + 3] = w.y * .5 * h, o[$.base + 0] = l[0], o[$.base + 1] = l[1], o[$.base + 2] = l[2], o[$.base + 3] = R;
+	let B = f(p);
+	o[$.rim + 0] = B[0], o[$.rim + 1] = B[1], o[$.rim + 2] = B[2], o[$.rim + 3] = 0, o[$.light + 0] = y.lightDir[0], o[$.light + 1] = y.lightDir[1], o[$.light + 2] = y.lightDir[2], o[$.light + 3] = 0, o[$.res + 0] = t.width, o[$.res + 1] = t.height, o[$.res + 2] = z && y.shadow ? 1 : 0, o[$.res + 3] = y.shadowSoftness;
+	let V = m, H = Math.min(V.length, 4);
 	for (let e = 0; e < 4; e++) {
-		let t = $.palette + e * 4, n = R[e];
+		let t = $.palette + e * 4, n = V[e];
 		if (n) {
 			let e = f(n);
 			o[t + 0] = e[0], o[t + 1] = e[1], o[t + 2] = e[2];
 		} else o[t + 0] = l[0], o[t + 1] = l[1], o[t + 2] = l[2];
 		o[t + 3] = 0;
 	}
-	let B = i.sources;
+	let U = i.sources;
 	for (let e = 0; e < 4; e++) {
-		let t = $.satPos + e * 4, n = B[e];
+		let t = $.satPos + e * 4, n = U[e];
 		n ? (o[t + 0] = n.x * h, o[t + 1] = n.y * h, o[t + 2] = n.radius * h, o[t + 3] = n.opacity) : (o[t + 0] = 0, o[t + 1] = 0, o[t + 2] = 0, o[t + 3] = 0);
 	}
-	let V = r.trailSources(g.satelliteRadius * .7 * C);
+	let W = r.trailSources(g.satelliteRadius * .7 * C);
 	for (let e = 0; e < 15; e++) {
-		let t = $.trailPos + e * 4, n = V.sources[e];
+		let t = $.trailPos + e * 4, n = W.sources[e];
 		n ? (o[t + 0] = n.x * .5 * h, o[t + 1] = n.y * .5 * h, o[t + 2] = n.radius * h) : (o[t + 0] = 0, o[t + 1] = 0, o[t + 2] = 0), o[t + 3] = 0;
 	}
-	return s[$.ints + 0] = z, s[$.ints + 1] = B.length, s[$.ints + 2] = V.count, s[$.ints + 3] = 0, e;
+	return s[$.ints + 0] = H, s[$.ints + 1] = U.length, s[$.ints + 2] = W.count, s[$.ints + 3] = 0, e;
 }
 //#endregion
-export { re as a, L as c, p as d, m as f, ie as i, I as l, ae as n, ne as o, oe as r, G as s, Q as t, j as u };
+export { se as a, V as c, p as d, m as f, ce as i, B as l, le as n, oe as o, ue as r, ee as s, Q as t, P as u };