@mkbabb/glass-ui 3.13.0 → 4.0.1

package/dist/aurora.js

@@ -1,13 +1,15 @@
-import { t as e } from "./_plugin-vue_export-helper-ok8pgSl5.js";
-import { a as t, c as n, i as r, l as i, n as a, o, r as s, s as c, t as l, u } from "./compile-B7NGOBcP.js";
-import { n as d, r as f, t as p } from "./presets-BKVyQfbt.js";
-import { a as m, c as h, i as g, l as _, o as v, t as y } from "./color-CUYeFbgd.js";
-import { t as ee } from "./useIntersectionPause-Cc-H9ees.js";
-import { t as te } from "./useScrollProgress-Dlvx-L3v.js";
-import { computed as b, createElementBlock as x, createElementVNode as S, defineComponent as ne, normalizeClass as re, normalizeStyle as C, onBeforeUnmount as w, onMounted as T, openBlock as ie, readonly as E, ref as D, unref as O, watch as k } from "vue";
-import { interpolateHue as ae, rawOklabToOklch as oe, srgbToOKLab as se } from "@mkbabb/value.js";
+import { t as e } from "./_plugin-vue_export-helper-C1je1s0u.js";
+import { a as t, i as n, n as r, r as i, t as a } from "./procedural-color.glsl-CIFsJhwI.js";
+import { a as o, c as s, i as c, l, n as u, r as d, s as f, t as p } from "./procedural-color.wgsl-CToRHugY.js";
+import { n as m, t as h } from "./compile-CQEFRTmx.js";
+import { i as g, n as _, r as v, t as y } from "./presets-BW2b-fr1.js";
+import { a as ee, c as b, i as te, l as ne, o as x, t as S } from "./color-DweYl7pE.js";
+import { t as re } from "./useIntersectionPause-Bj4sPkMh.js";
+import { t as ie } from "./useScrollProgress-D7bYU2RK.js";
+import { computed as C, createElementBlock as ae, createElementVNode as oe, defineComponent as se, normalizeClass as ce, normalizeStyle as le, onBeforeUnmount as w, onMounted as T, openBlock as ue, readonly as de, ref as E, unref as D, watch as O } from "vue";
+import { interpolateHue as fe, rawOklabToOklch as pe, srgbToOKLab as me } from "@mkbabb/value.js";
 //#region src/components/custom/aurora/constants/shaders/aurora.vert.ts
-var ce = "#version 300 es\nin vec2 aPos;\nout vec2 vUv;\nvoid main() {\n  vUv = aPos * 0.5 + 0.5;\n  gl_Position = vec4(aPos, 0.0, 1.0);\n}", le = "// ── Palette LUT (W5 — OKLCh/OKLab interpolation) ──────────────────────────────\n// The endpoints stay CPU-baked to LINEAR sRGB (the Aras precompute pattern); only\n// the BETWEEN-endpoint interpolation runs in perceptual space. AX.W11 — the ramp\n// itself (the smoothstep ease + the OKLab-rectangular-vs-OKLCh-hue-arc dispatch on\n// uHuePath) is the shared samplePaletteRamp from procedural-color.glsl, the SAME\n// source the WebGPU twin splices; samplePalette here only selects the bracketing\n// stop pair + the raw inter-stop t and hands them to the shared ramp.\nvec3 samplePalette(float id) {\n  id = clamp(id, 0.0, 1.0);\n  float scaled = id * float(uStopCount - 1);\n  int i0 = int(floor(scaled));\n  int i1 = min(i0 + 1, uStopCount - 1);\n  float t = fract(scaled);\n  return samplePaletteRamp(uPalette[i0], uPalette[i1], t, uHuePath);\n}\n\n// ── Nuclei field ──────────────────────────────────────────────────────────\nvoid nucleiField(vec2 p, float t, out float paletteId, out float valueMod) {\n  float accumBias  = 0.0;\n  float accumValue = 0.0;\n  float accumW     = 0.0;\n  for (int i = 0; i < MAX_NUCLEI; i++) {\n    if (i >= uNucleiCount) break;\n    vec2 posI = uNucleiPos[i]\n              + uNucleiDriftRadius[i] * vec2(\n                  cos(t * uNucleiDrift * K_NUCLEI + uNucleiDriftPhase[i]),\n                  sin(t * uNucleiDrift * K_NUCLEI + uNucleiDriftPhase[i] * 1.13)\n                );\n    vec2 diff = p - posI;\n    // Anisotropic Gaussian: rotate diff into the nucleus's local frame\n    // (major axis along uNucleiAngle), then scale the major-axis component by\n    // 1/elongation so the squared distance describes an ellipse. Defaults\n    // 1.0/0.0 reduce to the isotropic dot(diff, diff).\n    float ca = cos(uNucleiAngle[i]);\n    float sa = sin(uNucleiAngle[i]);\n    vec2  local = vec2( ca * diff.x + sa * diff.y,\n                       -sa * diff.x + ca * diff.y);\n    float along  = local.x / max(uNucleiElong[i], 0.01);\n    float across = local.y;\n    float d2 = along * along + across * across;\n    float r = max(uNucleiRadius[i], 0.01);\n    float w = exp(-uSoftmaxBeta * d2 / (r * r));\n    accumBias  += w * uNucleiPaletteBias[i];\n    accumValue += w * uNucleiValueBias[i];\n    accumW     += w;\n  }\n  paletteId = accumBias  / max(accumW, 1e-4);\n  valueMod  = accumValue / max(accumW, 1e-4);\n\n  // Palette drift — slow global paletteId breathe between adjacent stops. The\n  // rate rides K_PAL (perceptible ~30–60s hue cycle at the default coefficient);\n  // the amplitude lifts the authored coefficient into a 0.03..0.06 paletteId band\n  // so the palette visibly travels rather than dithering within one stop.\n  float palAmp = clamp(uPaletteDrift * 6.0, 0.0, 0.06);\n  paletteId += palAmp * sin(t * uPaletteDrift * K_PAL);\n  paletteId = clamp(paletteId, 0.0, 1.0);\n}\n", ue = "// ── Flow field ────────────────────────────────────────────────────────────\n// Structure-tensor prototype (body defined post-sampleBase in the mediums block).\n// Returns vec3(dir.x, dir.y, coherence) — the minor-eigenvector tangent + A.\nvec3 structureTensorField(vec2 p, float t, vec2 fallbackDir);\n\nvec2 flowField(vec2 p, float t) {\n  vec2 dir = vec2(1.0, 0.0);\n  if (uFlowPattern == 1) {\n    // radial\n    dir = normalize(p - uFlowFocal + 1e-4);\n  } else if (uFlowPattern == 2) {\n    // swirl — tangent to radial\n    vec2 rad = normalize(p - uFlowFocal + 1e-4);\n    dir = vec2(-rad.y, rad.x);\n  } else if (uFlowPattern == 3) {\n    // diagonal\n    float a = radians(uFlowAngle);\n    dir = vec2(cos(a), sin(a));\n  } else if (uFlowPattern == 4) {\n    // multi — curl-noise driven\n    float n = fbm(p * 2.0 + t * 0.02);\n    float a = n * 6.2831;\n    dir = vec2(cos(a), sin(a));\n  } else if (uFlowPattern == 5) {\n    // tensor / ETF — the color field's OWN edge-tangent flow (AW.W4.1). The\n    // diagonal angle seeds the smooth fallback the coherence-weighted blend\n    // relaxes toward in low-structure zones. Returns directly (no curl/cursor\n    // post-mix — the tensor field already carries the structure orientation;\n    // the cursor warp still drives the underlying color field in domainWarp).\n    float a = radians(uFlowAngle);\n    vec2 fallback = vec2(cos(a), sin(a));\n    return structureTensorField(p, t, fallback).xy;\n  }\n  // curl — perturb by noise. Radial gets much less curl so rays stay clean.\n  if (uFlowCurl > 0.0) {\n    float n = fbm(p * 3.0) - 0.5;\n    float curlAmt = uFlowCurl;\n    if (uFlowPattern == 1) curlAmt *= 0.25;  // radial: preserve ray clarity\n    if (uFlowPattern == 2) curlAmt *= 0.55;  // swirl: allow moderate curl\n    float a = atan(dir.y, dir.x) + n * 3.14159 * curlAmt;\n    dir = vec2(cos(a), sin(a));\n  }\n  // cursor influence — swirl around the cursor position\n  if (uCursorStrength > 0.001) {\n    vec2 toCur = uCursor - p;\n    float d = length(toCur);\n    float r = max(uCursorRadius, 0.01);\n    float w = exp(-(d * d) / (r * r * 0.5));\n    // swirl tangent\n    vec2 tangent = vec2(-toCur.y, toCur.x) / max(d, 1e-4);\n    float a0 = atan(dir.y, dir.x);\n    float a1 = atan(tangent.y, tangent.x);\n    // blend angle\n    float da = a1 - a0;\n    da = atan(sin(da), cos(da)); // wrap to [-pi, pi]\n    float a = a0 + da * w * uCursorStrength;\n    dir = vec2(cos(a), sin(a));\n  }\n  // AW.W8.1 — velocity-reactive swirl-BURST: a fast flick (uCursorBurst, decaying\n  // over ~1s) injects a transient flow bias along the pointer velocity, distinct from\n  // the steady attraction above (the stateless swirl is SUPERSEDED by this form). The\n  // burst is gated by the master tempo scalar on the CPU side (it decays to 0 under\n  // PRM/pause), so this term cannot leak motion under reduced-motion.\n  if (uCursorBurst > 0.001) {\n    vec2 toCur = uCursor - p;\n    float d = length(toCur);\n    float r = max(uCursorRadius, 0.01);\n    float w = exp(-(d * d) / (r * r * 0.7));\n    // Bias the flow toward the velocity direction, scaled by the decaying burst.\n    vec2 velDir = uCursorVelocity / max(length(uCursorVelocity), 1e-4);\n    float a0 = atan(dir.y, dir.x);\n    float a1 = atan(velDir.y, velDir.x);\n    float da = atan(sin(a1 - a0), cos(a1 - a0));\n    float a = a0 + da * w * uCursorBurst * 0.8;\n    dir = vec2(cos(a), sin(a));\n  }\n  return dir;\n}\n", de = "// Khronos PBR-Neutral tonemap (hue + saturation preserving)\nvec3 aces(vec3 color) {\n  const float startCompression = 0.8 - 0.04;\n  const float desaturation = 0.15;\n  float x = min(color.r, min(color.g, color.b));\n  float offset = x < 0.08 ? x - 6.25 * x * x : 0.04;\n  color -= offset;\n  float peak = max(color.r, max(color.g, color.b));\n  if (peak < startCompression) return clamp(color, 0.0, 1.0);\n  float d = 1.0 - startCompression;\n  float newPeak = 1.0 - d * d / (peak + d - startCompression);\n  color *= newPeak / peak;\n  float g = 1.0 - 1.0 / (desaturation * (peak - newPeak) + 1.0);\n  return clamp(mix(color, newPeak * vec3(1.0), g), 0.0, 1.0);\n}\n", fe = "// ── Curved swept brushstroke primitive ────────────────────────────────────\n// A stroke is a curved spine from A to B with quadratic bulge K (sideways offset\n// at midpoint). Width varies along length via a shape profile. Edge is ragged\n// (bristle-modulated), ends are rounded blobs, inside has streak modulation.\n//\n// For pixel p, we invert an approximate projection onto the curved spine:\n//   1) project p onto straight AB axis -> along0, cross0\n//   2) at parameter along0, spine offset = 4·K·along0·(1-along0) perpendicular\n//   3) cross = cross0 - spineOffset, refine along one iteration\n// Coverage blends paint atop col based on per-stroke color sampled from the\n// underlying base at spine midpoint — so overlapping strokes from different\n// palette regions meet at hard, ragged demarcations.\n\nstruct StrokeHit {\n  float coverage;  // 0..1\n  vec3  color;\n  float alongT;    // 0..1 along spine, for internal modulation\n  float crossN;    // cross distance / current half-width, signed\n  float edgeN;     // distance to edge in half-widths (0 = at edge, 1 = at spine)\n};\n\nStrokeHit noHit() { return StrokeHit(0.0, vec3(0.0), 0.0, 0.0, 0.0); }\n\nvec2 rotateDir(vec2 dir, float angle) {\n  float ca = cos(angle), sa = sin(angle);\n  return vec2(dir.x * ca - dir.y * sa, dir.x * sa + dir.y * ca);\n}\n\nvec2 safeDir(vec2 dir) {\n  float len = length(dir);\n  return len > 1e-4 ? dir / len : vec2(1.0, 0.0);\n}\n\n// Shape profile along the stroke. type:\n//   0 tapered    — thin-fat-thin (classic brush)\n//   1 load-drag  — fat start, tapers to a point\n//   2 dab        — ellipse-like, round center\n//   3 even       — near-constant, slight end taper\n//   4 comma      — van-Gogh atomic dab: ASYMMETRIC taper, fat near one end and\n//                  drawn to a fine point at the other (the divisionist crescent\n//                  mark). The width peak sits off-centre (~0.32) so the two ends\n//                  taper at DIFFERENT rates — a loaded round head + a dragged tail.\nfloat strokeShape(float t, int type) {\n  t = clamp(t, 0.0, 1.0);\n  if (type == 1) {\n    // loaded at t=0, tapering\n    return pow(1.0 - t, 0.55) * smoothstep(0.0, 0.08, t);\n  } else if (type == 2) {\n    // dab / blob\n    float d = t - 0.5;\n    return exp(-d * d * 12.0);\n  } else if (type == 3) {\n    // mostly even, slight end softening\n    return smoothstep(0.0, 0.08, t) * smoothstep(1.0, 0.92, t) * 0.95 + 0.05;\n  } else if (type == 4) {\n    // comma / crescent — width peaks at ~0.32 then drags to a fine tail. The two\n    // sides of the peak fall at different exponents (a fat round head, a long thin\n    // tail), so the mark reads as a directional van-Gogh comma, not a symmetric oval.\n    float peak = 0.32;\n    float head = smoothstep(0.0, peak, t);                            // rounded loaded head\n    float tail = pow(clamp((1.0 - t) / (1.0 - peak), 0.0, 1.0), 1.6); // dragged thin tail\n    return clamp(min(head, tail) * 1.04, 0.0, 1.0);\n  }\n  // 0 tapered\n  return smoothstep(0.0, 0.22, t) * smoothstep(1.0, 0.78, t);\n}\n\n// A single curved stroke.\n//   a, b           — endpoints\n//   halfW          — base half-width (world units)\n//   bulge          — signed perpendicular midpoint offset, in world units\n//   shapeType      — 0..3\n//   bristleFreq    — spatial frequency of edge raggedness\n//   bristleAmp     — 0..0.5 fraction of halfW chewed away at edge extrema\n//   streakSeed     — uniqueness seed for internal streaks\n//   colAtMid       — pre-sampled color at midpoint\nStrokeHit curvedStroke(vec2 p, vec2 a, vec2 b, float halfW,\n                       float bulge, int shapeType,\n                       float bristleFreq, float bristleAmp,\n                       float streakSeed, vec3 colAtMid) {\n  vec2 ab = b - a;\n  float L = length(ab);\n  if (L < 1e-5) return noHit();\n  vec2 tang = ab / L;\n  vec2 norm = vec2(-tang.y, tang.x);\n\n  vec2 rel = p - a;\n  float along0 = dot(rel, tang) / L;      // straight projection 0..1 along AB\n  float cross0 = dot(rel, norm);          // signed cross\n\n  // Spine sideways offset due to bulge (parabolic)\n  float bend = 4.0 * bulge * along0 * (1.0 - along0);\n  float cross1 = cross0 - bend;\n\n  // One refinement step: as the spine bends, the closest-point along shifts.\n  // Spine tangent differs from straight tangent; first-order correction:\n  float dBend = 4.0 * bulge * (1.0 - 2.0 * along0);  // d(bend)/d(along)\n  float along1 = along0 + (cross1 * dBend) / (L * (1.0 + dBend * dBend / (L * L)));\n\n  // Bristle-ragged edge: half-width gets chewed by 1D noise along the spine\n  float edgeNoise = fbm(vec2(along1 * bristleFreq, streakSeed * 7.3)) - 0.5;\n  float edgeNoise2 = fbm(vec2(along1 * bristleFreq * 2.3, streakSeed * 3.1 + 17.0)) - 0.5;\n  float edgeMod = 1.0 - bristleAmp * (0.6 + 0.4 * edgeNoise) * (0.5 + edgeNoise2);\n\n  // Width profile along the stroke\n  float shape = strokeShape(along1, shapeType);\n  float halfWNow = halfW * shape * edgeMod;\n\n  // Inside-segment coverage\n  float cov = 0.0;\n  if (along1 >= 0.0 && along1 <= 1.0 && halfWNow > 1e-6) {\n    float cn = abs(cross1) / halfWNow;\n    // fwidth-derived AA: half-width from the field's screen-space gradient so the\n    // edge holds ~1px regardless of zoom/DPR (mirrors metaball.frag.ts:252-255).\n    float aaC = max(fwidth(cn), 1e-6);\n    cov = 1.0 - smoothstep(1.0 - aaC, 1.0 + aaC, cn);\n  }\n\n  // End-cap blobs (rounded tips, not perpendicular cuts)\n  float capA = 0.0, capB = 0.0;\n  {\n    vec2 capCenterA = a + norm * (bulge * 0.0); // at a\n    vec2 capCenterB = a + tang * L + norm * bend * 0.0 + norm * 0.0; // at b; bend is 0 at endpoints\n    capCenterB = b;\n    float endShapeA = strokeShape(0.02, shapeType);\n    float endShapeB = strokeShape(0.98, shapeType);\n    float rA = halfW * endShapeA * (0.9 + 0.2 * fbm(vec2(streakSeed * 13.0)));\n    float rB = halfW * endShapeB * (0.9 + 0.2 * fbm(vec2(streakSeed * 19.0 + 5.0)));\n    float dA = length(p - capCenterA);\n    float dB = length(p - capCenterB);\n    // fwidth-derived AA on the clean Euclidean cap distances (mirrors metaball.frag.ts:252-255).\n    float aaA = max(fwidth(dA), 1e-6);\n    float aaB = max(fwidth(dB), 1e-6);\n    capA = 1.0 - smoothstep(rA - aaA, rA + aaA, dA);\n    capB = 1.0 - smoothstep(rB - aaB, rB + aaB, dB);\n    // Only apply caps if we're BEYOND the segment; otherwise the segment wins.\n    capA *= (along1 < 0.05) ? 1.0 : 0.0;\n    capB *= (along1 > 0.95) ? 1.0 : 0.0;\n  }\n\n  float coverage = max(cov, max(capA, capB));\n  if (coverage < 0.002) return noHit();\n\n  // Edge distance (for impasto later): 0 at edge, 1 at spine\n  float edgeDist = halfWNow > 1e-6 ? clamp(1.0 - abs(cross1) / halfWNow, 0.0, 1.0) : 0.0;\n\n  return StrokeHit(\n    coverage,\n    colAtMid,\n    clamp(along1, 0.0, 1.0),\n    halfWNow > 1e-6 ? cross1 / halfWNow : 0.0,\n    edgeDist\n  );\n}\n\n// The painterly STROKE mediums (those that deposit via paintStrokeLayers/paintOver):\n// oil (3), van-Gogh (5), oil-pastel (6). For these the stroke OVER-composite runs in\n// OKLab so overlapping complementaries mix through a chromatic path, not the linear-\n// RGB grey mud; the smooth/atmospheric pole never reaches paintOver so it pays nothing.\nbool isPainterlyStroke() {\n  return uMedium == 3 || uMedium == 5 || uMedium == 6;\n}\n\n// OKLab OVER-composite: lerp L,a,b of the over-color toward the under-color (the\n// Ottosson matrices spliced from procedural-color.glsl.ts — no new color math). Two\n// complementary-hued strokes overlapping transition through OKLab, where the midpoint\n// holds chroma, instead of the linear mix() that desaturates toward grey (slice 8 F2).\nvec3 paintOverOklab(vec3 under, vec3 over, float alpha) {\n  vec3 labU = linOklab(under);\n  vec3 labO = linOklab(over);\n  return max(oklabToLinearSrgb(mix(labU, labO, alpha)), vec3(0.0));\n}\n\n// Composite stroke over col with internal streaking + ACCUMULATE paint HEIGHT.\n//   streakAmp    — 0..0.2 how much internal streaks darken/lighten\n//   impastoAmp   — 0..1 paint-thickness contribution (the height the relight reads)\n//   hardness     — 0..1 how crisp the compositing transition is (1 = razor, 0 = creamy)\n//\n// AW.W4.2 — the fixed-RGB edge rim (the phantom upper-left light) is RETIRED. Each\n// stroke now deposits PAINT HEIGHT into the height accumulator (coverage x per-stroke\n// thickness, perturbed by the bristle/streak fbm for ridges/grooves); mediumOil\n// derives a normal from the accumulated height gradient (dFdx/dFdy) and relights it\n// once with the movable uLightDir source (relightImpasto), in linear before aces().\n//\n// AX.W13 — the internal streak now perturbs HUE + CHROMA in OKLCh (not value-only),\n// so a single stroke carries a small hue gradient (broken color at the ATOM level —\n// adjacent strokes shimmer like real impasto, not flat stamped swatches; slice 8 F5);\n// and the OVER-composite runs in OKLab on the painterly stroke mediums so overlapping\n// complementaries mix as pigment, not the linear-mix grey (slice 8 F2).\n//   impastoFloor — the per-stroke height-crown floor: oil's 0.4+0.6*edgeN reads\n//                  floor=0.4 (the crown tapers to 40% at the stroke edge); the van-Gogh\n//                  profile passes floor=1.0 so each atomic dab catches a FULL-height\n//                  glint (scope item 1d). The crown magnitude is floor + (1-floor)*edgeN.\nvoid paintOver(inout vec3 col, inout float height, StrokeHit s,\n               float streakFreq, float streakAmp,\n               float impastoAmp, float hardness, float streakSeed,\n               float impastoFloor) {\n  if (s.coverage < 0.002) return;\n  float strokeOpacity = clamp(uStrokeAmount, 0.0, 1.0);\n  if (strokeOpacity <= 0.001) return;\n  vec3 c = s.color;\n\n  // Internal streaks — fbm along spine, modulated by cross position. Gives\n  // the loaded-brush look: some bristles carry more pigment than others.\n  float streakA = fbm(vec2(s.alongT * streakFreq, s.crossN * 2.7 + streakSeed));\n  float streakB = fbm(vec2(s.alongT * streakFreq * 0.6 + streakSeed * 3.7, s.crossN * 4.1));\n  float streak = 0.6 * (streakA - 0.5) + 0.4 * (streakB - 0.5);\n\n  // ── WITHIN-stroke OKLCh broken color (AX.W13, slice 8 F5). The streak fBm perturbs\n  // the stroke's HUE + CHROMA along alongT/crossN, seeded per-stroke, so a single\n  // stroke reads a small hue gradient (the atom-level broken-color shimmer). Value is\n  // carried by the same streak so the loaded-brush light/dark survives. Gated to the\n  // painterly stroke mediums + uBrokenColor (the smooth pole keeps the cheap value-only\n  // streak). A second decorrelated fBm drives the hue axis so hue ≠ a scaled value.\n  float hueStreak = fbm(vec2(s.alongT * streakFreq * 0.45 + streakSeed * 5.3,\n                             s.crossN * 1.9 + streakSeed * 1.7)) - 0.5;\n  if (isPainterlyStroke() && uBrokenColor > 0.001) {\n    vec3 lch = oklabToOklch(linOklab(c));\n    // ±~10° within-stroke hue swing + ±~9% chroma, scaled by uBrokenColor. Smaller\n    // than the per-stroke cell jitter (brokenColorJitter) — this is the intra-stroke\n    // gradient layered on top of the per-cell pigment patch.\n    lch.z += hueStreak * (20.0 * PI / 180.0) * uBrokenColor;\n    lch.y = max(lch.y * (1.0 + (streak) * 0.18 * uBrokenColor), 0.0);\n    c = max(oklabToLinearSrgb(oklchToOklab(lch)), vec3(0.0));\n  }\n  // Value streak (the loaded-brush light/dark across the bristles) — kept for every\n  // medium that reaches paintOver (the smooth pole never does). The along-spine value\n  // modulation is the loaded-brush light/dark that gives each stroke its OWN directional\n  // luminance band — the structure-tensor signal a low-contrast color field (a uniform-hue\n  // oil field) has no other source for (the §4.2 anisotropy lift on the flat-palette case).\n  c *= 1.0 + streak * streakAmp * 3.2;\n\n  // Cross-stroke ridge/valley shading — a catch-light along the spine centre, darkening to\n  // the stroke edges. This is the strongest DIRECTIONAL luminance cue: it carves a bright\n  // ridge down the stroke's length and dark grooves between adjacent strokes, so a field of\n  // parallel strokes reads as parallel light/dark bands (the high-A van-Gogh signature) even\n  // where the underlying colour is flat. s.edgeN is 0 at the edge, 1 at the spine.\n  float ridgeShade = smoothstep(0.0, 0.55, s.edgeN);\n  c *= 1.0 + (ridgeShade - 0.5) * 0.30;\n\n  float softLimit = mix(0.35, 0.98, hardness);\n  float alpha = smoothstep(0.0, 1.0 - softLimit, s.coverage) * strokeOpacity;\n  // OKLab OVER-composite on the painterly stroke mediums (overlap mixes as pigment,\n  // not grey); the smooth pole never reaches paintOver, so it keeps the cheap path.\n  col = isPainterlyStroke() ? paintOverOklab(col, c, alpha) : mix(col, c, alpha);\n\n  // ── Accumulate paint HEIGHT (AW.W4.2). Thickness peaks at the stroke spine and\n  // falls toward the impastoFloor at the edge (edgeN: 0=edge, 1=spine), perturbed by\n  // the bristle streak so ridges/grooves form. The relight reads the gradient of this\n  // field. The van-Gogh profile's floor=1.0 gives every atomic dab a full-height crown.\n  float ridge = 0.7 + 0.3 * (streak + 0.5);\n  float crown = impastoFloor + (1.0 - impastoFloor) * s.edgeN;\n  height += s.coverage * impastoAmp * crown * ridge * alpha;\n}\n\n// AW.W4.2 — relight the accumulated paint height. Derive a normal from the height\n// gradient (dFdx/dFdy — derivatives are in-pattern; fwidth is already used for AA),\n// then apply diffuse + Blinn specular from the movable uLightDir/uLightColor source.\n// LINEAR-space lighting BEFORE aces() (the tonemap/OETF stay locked); thin strokes\n// inherit a low shininess, thick impasto raises it (high-height ridges catch a sharp\n// glint). canvasBase is the tooth height so the canvas weave also catches the light.\nvec3 relightImpasto(vec3 col, float height, float canvasBase) {\n  float h = height + canvasBase;\n  // Height-gradient normal. heightScale trades relief depth for surface flatness.\n  float heightScale = 6.0;\n  vec3 N = normalize(vec3(-dFdx(h) * 40.0, -dFdy(h) * 40.0, 1.0 / heightScale));\n  vec3 L = normalize(uLightDir);\n  vec3 V = vec3(0.0, 0.0, 1.0);            // view straight on\n  vec3 H = normalize(L + V);\n  float diff = max(dot(N, L), 0.0);\n  // Thick impasto → high shininess (sharp glint); thin paint → broad soft sheen.\n  float shininess = mix(8.0, 64.0, clamp(height * 2.0, 0.0, 1.0));\n  float spec = pow(max(dot(N, H), 0.0), shininess);\n  // Modulate the lit terms by how much paint sits here (no relief on bare ground).\n  float relief = clamp(h * 1.5, 0.0, 1.0);\n  vec3 lit = col * (1.0 + uImpasto * relief * (diff - 0.5) * 0.5 * uLightColor);\n  // The specular glint is gated to ACCUMULATED-thick paint (height²) — a per-fine-stroke\n  // sharp glint on every thin fill dab read as bright isotropic SPECKLE that buried the\n  // directional stroke bands (the §4.2 anisotropy / §4.3 slope sink on the low-contrast\n  // oil field). Only where real impasto has built up does the raked glint fire, so the\n  // highlight tracks the macro ridges, not the fine speckle.\n  float specGate = clamp(height, 0.0, 1.0);\n  lit += uImpasto * relief * spec * specGate * 0.5 * uLightColor;\n  return lit;\n}\n\n// Best-of-9-neighbor placement: sample 3x3 surrounding cells and take the\n// stroke that covers this pixel most. Breaks the grid by per-cell jitter and\n// sparse density via noise thresholding.\n//   energyGrade — the SBR energy-grade magnitude (0 = none, the oil/oil-pastel\n//                 baseline; 1 = the full van-Gogh Starry-Night cascade — long\n//                 confident strokes in the lights/coherent zones, short dabs in\n//                 the darks/flat zones). The van-Gogh profile carries 1.0 (the\n//                 energy grade is a PROFILE field, not a buried uMedium==5 branch).\nStrokeHit bestOil(vec2 p, float cellSize, float lenMul, float halfWMul,\n                  float jitterAmt, float density, int shapeType,\n                  float bristleAmp, vec2 flow, float t, float seed,\n                  float energyGrade) {\n  vec2 cell = floor(p / cellSize);\n\n  StrokeHit best = noHit();\n  // 3x3 neighborhood\n  for (int dy = -1; dy <= 1; dy++) {\n    for (int dx = -1; dx <= 1; dx++) {\n      vec2 cc = cell + vec2(float(dx), float(dy));\n      vec2 hh = hash22(cc + seed) - 0.5;\n      // Density gate — noise threshold; sparse placement.\n      float gate = hash21(cc * 1.7 + seed * 0.3);\n      if (gate > density) continue;\n\n      vec2 center = (cc + 0.5 + hh * jitterAmt) * cellSize;\n\n      // Per-stroke direction source — the AW.W4.1 strokeOrient switch:\n      //   uStrokeOrient==0 (flow)   — the layer-provided hand-authored flowField.\n      //   uStrokeOrient==1 (tensor) — the structure-tensor minor eigenvector at the\n      //                               cell center (the color field's edge-tangent),\n      //                               so strokes hug the color zones (van Gogh).\n      // coh carries the per-cell coherence A for the W4.3 energy grading (long\n      // confident strokes where structure is coherent, stubby dabs in flat zones).\n      vec2 f;\n      float coh = 1.0;\n      if (uStrokeOrient == 1) {\n        vec3 tf = structureTensorField(center, t, safeDir(flow));\n        f = safeDir(tf.xy);\n        coh = tf.z;\n      } else {\n        f = safeDir(flow);\n      }\n      // Per-stroke deterministic local perturbation so alternate stroke layers stay\n      // live. Tensor strokes get FAR less angular jitter (the field already orients them\n      // — over-jittering the tensor strokes scatters their directions and the rendered\n      // field reads isotropic, the §4.2 anisotropy floor). The flow-oriented (non-tensor)\n      // strokes keep the full jitter (they have no per-cell orientation source).\n      float jScale = (uStrokeOrient == 1) ? 0.15 : 1.0;\n      float angJ = (hash21(cc + seed + 11.0) - 0.5) * 0.9 * jScale;  // +/- 0.45/0.07 rad\n      float curlScale = (uStrokeOrient == 1) ? 0.30 : 1.0;\n      float localCurl = (fbm(center * (2.6 + seed * 0.11) + seed * 1.9) - 0.5) * 0.55 * uFlowCurl * curlScale;\n      vec2 dir = rotateDir(f, angJ + localCurl);\n\n      // AX.W13 — the SBR energy grade is a PROFILE field (energyGrade), not a buried\n      // uMedium==5 branch. The van-Gogh profile passes energyGrade=1.0 (the full\n      // Starry-Night cascade); oil/oil-pastel pass 0.0 (uniform length). The grade\n      // modulates stroke length by local luminance (bright → long confident strokes,\n      // dark → short dabs, the Kolmogorov/Batchelor congruence) AND by coherence\n      // (coherent zones → long, flat → stubby), lerped in by energyGrade so a profile\n      // can dial the cascade from off (0) to full (1).\n      float energy = 1.0;\n      if (energyGrade > 0.001) {\n        float luma = clamp(dot(sampleBase(center, t), W_LUMA), 0.0, 1.0);\n        float graded = mix(0.5, 1.5, luma) * mix(0.6, 1.0, coh);\n        energy = mix(1.0, graded, energyGrade);\n      }\n\n      float lenV = cellSize * lenMul * (0.65 + 0.55 * hash21(cc + seed + 23.0)) * energy;\n      float halfW = cellSize * halfWMul * (0.70 + 0.55 * hash21(cc + seed + 41.0));\n      float bulge = (hash21(cc + seed + 53.0) - 0.5) * lenV * 0.35;\n\n      vec2 a = center - dir * (lenV * 0.5);\n      vec2 b = center + dir * (lenV * 0.5);\n\n      vec3 colMid = brokenColorJitter(\n        sampleBase(center, t),\n        hash21(cc + seed + 89.0),\n        hash21(cc * 2.3 + seed + 97.0),\n        1.0\n      );\n\n      StrokeHit h = curvedStroke(p, a, b, halfW, bulge, shapeType,\n                                 7.0, bristleAmp,\n                                 hash21(cc + seed + 67.0), colMid);\n      if (h.coverage > best.coverage) best = h;\n    }\n  }\n  return best;\n}\n", pe = "// ── Van-Gogh — DEDICATED atomic-dab body (W-AUR-VANGOGH-REBUILD) ─────────────\n// REBUILT FROM FIRST PRINCIPLES. The prior path routed van-Gogh through the shared\n// oil paintStrokeMedium cascade — four DENSE layers of LONG-THIN tensor-oriented\n// strokes (lenMul up to 5.2, widMul 0.22) hugging a continuous coherent field. Packed\n// that tight, thin strokes that follow continuous iso-bands MERGE into smooth flowing\n// filaments — the \"marbled flow-bands\" read (a spirograph/fingerprint, not impasto),\n// AND the per-cell 8-tap structure tensor × 4 layers × 9 cells = ~324 sampleBase/pixel\n// made it run at ~4fps. Van Gogh's sky is NOT continuous flow — it is DISCRETE, SHORT,\n// LOADED comma dabs, separated by visible ground, queued into swirl ROWS.\n//\n// The rebuild is a purpose-built two-layer DAB pass, disjoint from the oil cascade:\n//   • orientation from the analytic SWIRL flowField (van-Gogh skies ARE swirl rows) +\n//     per-cell angular jitter — NOT the expensive per-cell tensor (the lag source, and\n//     the marble source: a jittered discrete dab can NOT merge into a continuous\n//     filament the way a tensor-locked thin stroke does).\n//   • SHORT comma dabs (len/wid ≈ 2.6 : 1 — a loaded crescent mark, not a 10:1\n//     filament), so each reads as a SEPARABLE atom.\n//   • SPARSE placement with visible inter-dab ground (density ≈ 0.42), over a darkened\n//     swirl underpainting (the Starry-Night visible ground the bright dabs read over).\n//   • a single sampleBase per cell + OKLCh broken-color jitter (pigment shimmer);\n//     full-height impasto crown per dab (each catches its own raked glint).\n// This is ~2 layers × 9 cells × ~1 sampleBase ≈ 20/pixel (a ~16× cost cut, ~60fps), and\n// jittered short dabs CANNOT marble.\n\n// One atomic comma dab — a CLEAN analytic crescent SDF (NOT the oil curvedStroke).\n// The oil curvedStroke's spine-projection + rounded end-caps + fwidth-AA + bristle\n// edge each leave a thin sub-pixel skirt, and a field of those skirts over open ground\n// abuts into a faint network of SEAM LINES. A van-Gogh dab needs none of that\n// machinery — it is a short oriented capsule (rounded line segment) whose half-width\n// tapers along its length (the loaded-head → dragged-tail comma). The capsule SDF is\n// derivative-free (smoothstep AA against the analytic distance), so no fwidth skirt, no\n// best-of-9 seam. Returns a StrokeHit (coverage/color/alongT/edgeN) for paintOver.\n//   cellSize — dab cell pitch (UV units); lenMul/widMul — dab aspect off the cell;\n//   density  — placement gate (sparse → visible gaps); seed — layer decorrelation.\nStrokeHit vangoghDab(vec2 p, float cellSize, float lenMul, float widMul,\n                     float density, float jitterAmt, float angJitter,\n                     vec2 flow, float t, float seed) {\n  vec2 cell = floor(p / cellSize);\n  StrokeHit best = noHit();\n  for (int dy = -1; dy <= 1; dy++) {\n    for (int dx = -1; dx <= 1; dx++) {\n      vec2 cc = cell + vec2(float(dx), float(dy));\n      // Sparse density gate — leaves visible canvas ground between dabs (atomicity).\n      if (hash21(cc * 1.7 + seed * 0.3) > density) continue;\n      vec2 jit = (hash22(cc + seed) - 0.5) * jitterAmt;\n      vec2 center = (cc + 0.5 + jit) * cellSize;\n\n      // Orientation: the analytic swirl flow at the cell center + per-cell angular\n      // jitter (so adjacent dabs FAN — the divisionist scatter that keeps the dabs from\n      // locking into perfectly-concentric rows) + a flow-noise curl wobble.\n      float ang = (hash21(cc + seed + 11.0) - 0.5) * angJitter;\n      float curl = (fbm(center * 3.1 + seed * 1.9) - 0.5) * (0.18 + 0.3 * uFlowCurl);\n      vec2 dir = rotateDir(safeDir(flow), ang + curl);\n\n      // SHORT comma — len ≈ the cell pitch (a contained dab, not a tiling filament).\n      float len = cellSize * lenMul * (0.82 + 0.30 * hash21(cc + seed + 23.0));\n      float halfW = max(cellSize * widMul * (0.86 + 0.28 * hash21(cc + seed + 41.0)), 0.010);\n\n      // Local dab frame: along = projection onto the dab axis (centred 0), perp = cross.\n      vec2 rel = p - center;\n      float along = dot(rel, dir);\n      float perp  = dot(rel, vec2(-dir.y, dir.x));\n      // The comma TAPER: half-width peaks ~1/3 from the loaded head and drags to a fine\n      // tail (the asymmetric crescent). u in [0,1] along the dab from head to tail.\n      float halfLen = len * 0.5;\n      float u = clamp((along + halfLen) / max(len, 1e-5), 0.0, 1.0);\n      float head = smoothstep(0.0, 0.30, u);                 // rounded loaded head\n      float tail = pow(clamp((1.0 - u) / 0.70, 0.0, 1.0), 1.3); // dragged thin tail\n      float prof = clamp(min(head, tail) * 1.05, 0.0, 1.0);\n      float wNow = halfW * prof;\n      // A gentle comma CURVE — bow the spine sideways so the mark reads as a crescent.\n      float bowAmt = (hash21(cc + seed + 53.0) - 0.5) * halfW * 1.4;\n      float spine = bowAmt * sin(u * PI);                    // 0 at ends, max mid\n      float d = abs(perp - spine);\n      // Derivative-free AA: a fixed fraction of the dab width (no fwidth skirt). Coverage\n      // is 1 in the dab body, smoothly 0 at the edge; outside the [head,tail] span wNow→0\n      // so the dab self-terminates with no end-cap blob.\n      float aa = max(halfW * 0.45, 0.002);\n      float cov = (wNow > 1e-5) ? (1.0 - smoothstep(wNow - aa, wNow + aa, d)) : 0.0;\n      if (cov < 0.004) continue;\n\n      // One sampleBase per cell + OKLCh broken-color jitter (pigment patch shimmer). The\n      // dab pigment is chroma-boosted (saturate3 1.22) AND value-lifted (×1.14) at the\n      // SOURCE — the loaded brush mark deposits BRIGHTER + more saturated pigment than the\n      // base field, so the loaded-mark pixels lift BOTH the §4.1 colourfulness AND the field\n      // mean-luma. The higher mean raises the gap-fraction threshold (meanLum·0.6) so the\n      // moderate-depth inter-dab ground falls below it (the gap-fraction atomicity read is\n      // DECOUPLED from the ground floor depth — the ground can stay only moderately dark, so\n      // the §4.3 slope is NOT over-steepened, while the bright dabs manufacture the gap). The\n      // dark ground is untouched, so this lifts C + gap WITHOUT a band-conflicting deep ground.\n      vec3 colMid = brokenColorJitter(\n        saturate3(sampleBase(center, t), 1.22),\n        hash21(cc + seed + 89.0),\n        hash21(cc * 2.3 + seed + 97.0),\n        1.0\n      );\n\n      // edgeN: 1 at the spine, 0 at the edge (the impasto crown reads it).\n      float edgeN = (wNow > 1e-5) ? clamp(1.0 - d / wNow, 0.0, 1.0) : 0.0;\n      if (cov > best.coverage) {\n        best = StrokeHit(cov, colMid, u, (wNow > 1e-5 ? (perp - spine) / wNow : 0.0), edgeN);\n      }\n    }\n  }\n  return best;\n}\n\nvec3 mediumVangogh(vec3 col, vec2 p, float t) {\n  vec3 result = col;\n  float height = 0.0;\n\n  // The swirl flow basis (van-Gogh sky = swirl rows). Computed ONCE per pixel.\n  vec2 flow = flowField(p, t);\n\n  // Dab cell pitch in UV [0,1]. CHUNKY, SEPARABLE marks — the big layer is ~1/14 of\n  // the frame so it reads as a distinct loaded brush dab (not a fine speckle, not a\n  // screen-spanning filament). uStrokeScale (authored ~110) modulates the pitch\n  // INVERSELY (higher = finer dabs) around the ~0.07 chunky default; clamped so the\n  // dab never collapses to noise or balloons past one screen. uStrokeAnisotropy\n  // stretches the comma aspect (longer, more directional marks).\n  float pitch = clamp(7.0 / max(uStrokeScale, 20.0), 0.044, 0.12);\n  // The dab LENGTH ≈ the cell pitch (lenAniso ~1.2-1.7) so a dab sits WITHIN its cell\n  // and leaves visible ground around it — a longer dab (lenMul > 2) spans neighbour\n  // cells and tiles the field into continuous coverage (no gaps, the gap-fraction floor\n  // miss + the marble). Narrow width keeps it a directional crescent, not a wide shard.\n  // lenAniso/widAniso — the comma aspect. The W-AUR-VANGOGH-CONFIG reconcile lengthened\n  // the comma (1.25-1.75 → 1.5-2.0) and THINNED it (0.40-0.30 → 0.32-0.24): the rebuilt\n  // dabs measured §4.2 structure-tensor anisotropy A=0.699, BELOW the reference band\n  // floor 0.732 (the dabs read slightly too short-fat → locally isotropic, dragging mean\n  // A down). A LONGER, THINNER comma hugs the swirl-flow tangent more strongly so the\n  // rendered field reads coherent the way the tensor metric measures it (A lift), while\n  // the thinner width LEAVES MORE inter-dab ground (the gap-fraction atomicity read is\n  // preserved, not closed). The length stays ≈ the cell pitch so each dab remains a\n  // CONTAINED separable mark, not a tiling filament.\n  float lenAniso = mix(2.1, 2.7, uStrokeAnisotropy);   // comma length / cell pitch (a contained dab)\n  float widAniso = mix(0.23, 0.16, uStrokeAnisotropy); // comma width / cell pitch (a thin loaded crescent)\n\n  // Layer 1 — the big swirl-row dabs: the loaded gestural marks, SPARSE so the dark\n  // underpainting shows BETWEEN them (the atomicity read — visible inter-dab ground).\n  // hardness 0.74 — the loaded brush deposits with a SOFT shoulder, so adjacent dabs\n  // feather instead of razor-clipping into hard polygon shards. jitterAmt 0.85 fully\n  // scatters the cell centres so the placement GRID never reads as a periodic seam.\n  // angJitter EASED 0.55 → 0.38 (the reconcile A lift): a tighter angular spread keeps\n  // adjacent dabs more flow-aligned (coherent) while the per-cell curl wobble still\n  // breaks the perfectly-concentric-rows read (the dabs FAN gently, not lock to a grid).\n  // angJitter eased 0.38 → 0.26 + hardness 0.74 → 0.84 (the §4.2 A lift): a tighter angular\n  // spread keeps adjacent dabs strongly flow-aligned (the structure-tensor coherence) and a\n  // sharper edge gives each dab a STRONGER directional gradient (a soft dab's faint edge\n  // reads weakly anisotropic). The per-cell curl wobble still breaks the perfectly-\n  // concentric-rows read, so the dabs FAN gently rather than locking to a grid.\n  StrokeHit dBig = vangoghDab(p, pitch, lenAniso, widAniso,\n                              0.30, 0.50, 0.08, flow, t, 1.3);\n  // streakFreq raised 14 → 34 (the §4.3 slope fix): a high-frequency streak runs ALONG the\n  // dab spine, injecting fine-scale luma energy that flattens the too-steep β WITHOUT\n  // adding new isotropic marks (the along-spine streak stays directional, so the §4.2\n  // anisotropy the long-thin dab carries is preserved — the band-conflict resolution).\n  paintOver(result, height, dBig, 48.0, 0.34,\n            uImpasto * 0.75 * uStrokeAmount, 0.84, 1.3, 1.0 /* full crown */);\n\n  // Layer 2 — two-thirds-pitch accent dabs (the divisionist sparkle between the big\n  // marks — denser placement but still contained dabs, so dark canvas survives between).\n  // angJitter eased 0.5 → 0.34 (same A-lift reasoning — the accent dabs queue along the\n  // flow rather than scattering isotropically); hardness 0.78 → 0.86 sharpens the gradient.\n  StrokeHit dSml = vangoghDab(p + vec2(-5.1, 8.4), pitch * 0.62, lenAniso * 0.95,\n                              widAniso * 0.95, 0.36, 0.60, 0.20, flow, t, 4.1);\n  paintOver(result, height, dSml, 42.0, 0.22,\n            uImpasto * 0.55 * uStrokeAmount, 0.86, 4.1, 1.0);\n\n  // Layer 3 — the FINE divisionist flecks (W-AUR-VANGOGH-CONFIG reconcile). The two-layer\n  // field read as a few BIG cobalt dabs over a smooth empty teal ground (the side-by-side\n  // against starry-night-crop.png: the painting is DENSE small varied flecks, not sparse\n  // large marks). The flecks are LONG·THIN directional commas hugging the swirl tangent\n  // (the A-preserving fleck — short fat flecks read isotropic and drag §4.2 down). Sparse\n  // density 0.30 keeps the inter-fleck ground open (the gap-fraction atomicity read).\n  // The flecks are FINE (≈0.42× pitch) directional commas, tightly flow-aligned (angJitter\n  // 0.22) and SPARSE (0.30) — the divisionist sparkle between the loaded marks that keeps\n  // the inter-dab ground open (atomicity) while adding the small-mark mid-band energy. Thin\n  // + flow-aligned so the §4.2 anisotropy the long dabs carry survives the added small marks.\n  StrokeHit dFleck = vangoghDab(p + vec2(8.7, -3.2), pitch * 0.42, lenAniso * 1.0,\n                                widAniso * 0.68, 0.30, 0.85, 0.22, flow, t, 7.7);\n  paintOver(result, height, dFleck, 30.0, 0.16,\n            uImpasto * 0.40 * uStrokeAmount, 0.80, 7.7, 1.0);\n\n  // ── Flow-aligned GROUND swirl-GROOVING (W-AUR-VANGOGH-CONFIG reconcile) ──\n  // The bare ground between the dabs was a SMOOTH low-frequency teal gradient — the §4.3\n  // power-spectrum slope ran too STEEP (β≈−2.0: too much low-frequency nuclei-field energy\n  // vs the high band). Van Gogh's actual sky has the swirl BRUSHWORK everywhere, so groove\n  // the OPEN GROUND with a flow-aligned ripple. It is SUBTRACTIVE-biased (it only darkens,\n  // never brightens the ground), so it DEEPENS the inter-dab valleys (MORE gap-fraction —\n  // atomicity preserved/strengthened) while injecting the across-flow mid-band luma energy\n  // that flattens β. Gated to the open ground (low height) via the (1-height) gate so it\n  // never lifts/erodes the loaded dab bodies (the dab read is untouched). Across the swirl\n  // FLOW axis so the grooves run WITH the sky's brushwork direction.\n  // The ground swirl-combing — a FIXED-global-angle directional grain on the open ground.\n  // A flow-FOLLOWING comb varies its direction with the swirl across the readback window,\n  // so its blurred structure tensor averages toward isotropic (no §4.2 lift). A comb at one\n  // GLOBAL angle (van-Gogh's dominant diagonal stroke grain) is COHERENT over the whole\n  // window ⇒ every open-ground pixel reads the SAME direction (the §4.2 mean-anisotropy lift\n  // the swirl comb could not reach). Two resolved bands (≈22 + ≈44 across the global axis ⇒\n  // ~21px + ~11px at 464px, both in the §4.3 fit window) also feed the β flatten. Gated to\n  // the open ground (low height) so the loaded dabs are untouched; SUBTRACTIVE-biased so the\n  // inter-dab valleys deepen (gap-fraction preserved, not filled).\n  vec2 gdir = normalize(vec2(0.82, 0.57));   // the dominant global stroke grain axis\n  float gacross = dot(p, vec2(-gdir.y, gdir.x));\n  float groundGate = 1.0 - smoothstep(0.05, 0.28, height);\n  float groove = sin(gacross * 22.0) + 0.5 * sin(gacross * 44.0 + 1.3);\n  result *= 1.0 - 0.16 * groundGate * (0.5 - 0.5 * groove);\n\n  // Canvas tooth (fine linen weave — base relief the raking light catches).\n  float tooth1 = vnoise(p * 220.0);\n  float tooth2 = vnoise(p * 220.0 * vec2(0.6, 2.4) + 37.0);\n  float tooth  = (0.6 * tooth1 + 0.4 * tooth2) - 0.5;\n  result *= 1.0 + tooth * 0.06 * uCanvasGrain;\n\n  // Relight the accumulated impasto height (raking light catches each dab's crown —\n  // the visible impasto ridge that makes a dab read as a RAISED loaded mark).\n  float canvasBase = tooth * 0.06 * 0.5;\n  result = relightImpasto(result, height, canvasBase);\n\n  // The visible Starry-Night GROUND — the open canvas between dabs (low accumulated\n  // height) darkens toward the deep underpainting, so the bright impasto dabs read as\n  // SEPARABLE atomic marks over a darker field (the dab↔ground contrast the user's\n  // \"reads as strokes\" bar needs, AND the painterly-statistics gap-fraction floor — the\n  // dark inter-dab ground must register as a measurable fraction of the field). 0.30\n  // floor over a WIDER smoothstep (0..0.40): deep enough that the open gaps read as the\n  // dark Starry-Night sky between marks (clearing the 4% gap floor) while the gradual\n  // ramp keeps the gaps reading as flowing underpainting, not the hard black ink LINES a\n  // narrow near-black step (0.26 / 0.0..0.32) produced.\n  // The smoothstep FLOOR starts at 0.10 (not 0): a dab's faint AA coverage TAIL deposits\n  // a sliver of height, and where many tails align along the swirl flow those slivers\n  // accumulate into a faint continuous height ridge that read as thin curved LINES in the\n  // ground. Holding the ground full-dark until a real dab CORE (height > 0.10) lifts out\n  // keeps the open field clean — only the loaded dab bodies emerge, the tails stay ground.\n  result *= mix(0.115, 1.0, smoothstep(0.085, 0.31, height));\n\n  // The van-Gogh chroma — a light OKLCh punch. EASED to 0.92 (off 1.04) in the\n  // W-AUR-VANGOGH-CONFIG reconcile: the rebuilt discrete-dab register's per-dab OKLCh\n  // broken-color jitter (the ±16° hue + ±chroma swing per cell) lifts the field's\n  // colourfulness on its own, so the prior 1.04 punch on TOP overshot the §4.1\n  // reference band (rendered C=100.5 vs the painting's measured 70.67 ± 15/+25 →\n  // ceiling 95.67 — neon, not Starry Night's rich-but-not-garish blue/gold). The 0.92\n  // demote pulls the rendered colourfulness back IN-band while the broken-color jitter\n  // keeps the per-dab pigment variation (the discrete-stroke read is untouched — this\n  // is a chroma scale, not a dab-geometry change). Band-anchored on the real painting\n  // (proof:aurora-arresting), NOT a widened band.\n  result = saturate3(result, 1.0);\n  return result;\n}\n\n", me = "// ── Medium overlays ───────────────────────────────────────────────────────\n\n// A quick re-computation of base color for edge-mask sampling\nvec3 sampleBase(vec2 p, float t) {\n  vec2 pw = domainWarp(p, t);\n  float id; float vm;\n  nucleiField(pw, t, id, vm);\n  vec3 c = samplePalette(id);\n  c *= 1.0 + uValueVariance * vm;\n  return c;\n}\n\n// ── Structure-tensor / edge-tangent-flow (AW.W4.1) ─────────────────────────\n// The keystone of the painterly engine: derive stroke orientation from the COLOR\n// FIELD's OWN structure, not a hand-authored flow pattern. A 3x3 Sobel over\n// luma(sampleBase) yields the gradient (Gx,Gy); the 2x2 structure tensor\n// J = [[Gx², GxGy],[GxGy, Gy²]] eigen-decomposes closed-form into the MINOR\n// eigenvector (the edge-TANGENT — least color change, the stroke direction; the\n// MAJOR eigenvector is the gradient/edge-NORMAL and would make strokes cross the\n// bands) and the coherence A = (λ1-λ2)/(λ1+λ2). Low-coherence (flat) zones relax\n// toward the smooth fallbackDir by (1-A) so tensor noise never reads as jitter.\n// Single-pass small-tap WebGL2 form; the Gaussian-smoothed multi-tap + LIC smear\n// is the AW.W7 WebGPU multi-pass fold. Returns vec3(tangent.x, tangent.y, A).\nvec3 structureTensorField(vec2 p, float t, vec2 fallbackDir) {\n  float e = 0.0035; // small-tap neighborhood (edge-mask scale)\n  // 3x3 luma samples.\n  float l00 = dot(sampleBase(p + vec2(-e, -e), t), W_LUMA);\n  float l10 = dot(sampleBase(p + vec2( 0.0, -e), t), W_LUMA);\n  float l20 = dot(sampleBase(p + vec2( e, -e), t), W_LUMA);\n  float l01 = dot(sampleBase(p + vec2(-e, 0.0), t), W_LUMA);\n  float l21 = dot(sampleBase(p + vec2( e, 0.0), t), W_LUMA);\n  float l02 = dot(sampleBase(p + vec2(-e, e), t), W_LUMA);\n  float l12 = dot(sampleBase(p + vec2( 0.0, e), t), W_LUMA);\n  float l22 = dot(sampleBase(p + vec2( e, e), t), W_LUMA);\n  // Sobel derivatives.\n  float Gx = (l20 + 2.0 * l21 + l22) - (l00 + 2.0 * l01 + l02);\n  float Gy = (l02 + 2.0 * l12 + l22) - (l00 + 2.0 * l10 + l20);\n  // Structure tensor components.\n  float Jxx = Gx * Gx;\n  float Jyy = Gy * Gy;\n  float Jxy = Gx * Gy;\n  // Closed-form 2x2 eigenvalues: λ = 0.5(tr ± sqrt((Jxx-Jyy)² + 4Jxy²)).\n  float tr = Jxx + Jyy;\n  float disc = sqrt(max((Jxx - Jyy) * (Jxx - Jyy) + 4.0 * Jxy * Jxy, 0.0));\n  float lambdaMaj = 0.5 * (tr + disc);\n  float lambdaMin = 0.5 * (tr - disc);\n  // The MAJOR eigenvector points along the gradient (the edge NORMAL); the MINOR\n  // is its perpendicular — the edge TANGENT (the stroke direction). The principal\n  // angle is θ = 0.5·atan2(2·Jxy, Jxx-Jyy); this form is robust at Jxy≈0 (where the\n  // (λ-Jyy, Jxy) eigenvector formula collapses to a zero vector). The tangent is at\n  // θ + 90°. Guard the fully-isotropic (no-gradient) degenerate case → fallback dir.\n  vec2 tangent;\n  if (disc < 1e-7) {\n    tangent = fallbackDir;\n  } else {\n    float theta = 0.5 * atan(2.0 * Jxy, Jxx - Jyy); // major (gradient) angle\n    tangent = vec2(-sin(theta), cos(theta));          // perpendicular = tangent\n  }\n  // Coherence A in [0,1] — high where one eigenvalue dominates (a strong edge),\n  // ~0 in flat zones. Relax toward the smooth fallback by (1 - A).\n  float A = (lambdaMaj - lambdaMin) / (lambdaMaj + lambdaMin + 1e-6);\n  A = clamp(A, 0.0, 1.0);\n  // Coherence-weighted blend: full tensor in coherent zones, smooth in flat ones.\n  // Resolve the 180° sign ambiguity toward the fallback so the blend is continuous.\n  if (dot(tangent, fallbackDir) < 0.0) tangent = -tangent;\n  // The blend weight is a low-power curve (A^0.35) so even MODERATE-coherence zones queue\n  // strongly along the tangent (a flowing van-Gogh band reads coherent the way the\n  // structure-tensor metric measures it; a linear A under-orients the mid-band and the\n  // rendered field reads too isotropic — the §4.2 anisotropy lift). The flat-zone (A≈0)\n  // relax to the fallback survives (0^0.35=0), so tensor noise never reads as jitter.\n  float blendW = pow(A, 0.28);\n  vec2 dir = normalize(mix(fallbackDir, tangent, blendW) + 1e-6);\n  return vec3(dir, A);\n}\n\nvec3 mediumWatercolor(vec3 col, vec2 p, float t) {\n  // Wet-edge cauliflowers via luma-gradient magnitude\n  float eps = 0.004;\n  vec3 cx1 = sampleBase(p + vec2(eps, 0.0), t);\n  vec3 cx2 = sampleBase(p - vec2(eps, 0.0), t);\n  vec3 cy1 = sampleBase(p + vec2(0.0, eps), t);\n  vec3 cy2 = sampleBase(p - vec2(0.0, eps), t);\n  float gx = dot(cx1 - cx2, W_LUMA);\n  float gy = dot(cy1 - cy2, W_LUMA);\n  float edge = sqrt(gx * gx + gy * gy) / (2.0 * eps);\n  float mask = smoothstep(0.0, 2.5, edge);\n  col *= mix(1.0, 0.78, mask * uWetEdge);\n\n  // Granulation — pigment settles in paper tooth\n  float paper = 0.5 * vnoise(p * 160.0) + 0.5 * vnoise(p * 360.0);\n  float pigLoad = 1.0 - dot(col, W_LUMA);\n  col *= 1.0 - uGranulation * pigLoad * (paper - 0.5);\n\n  // Wash banding — faint horizontal wet gradient\n  float band = fbm(vec2(p.x * 1.5, p.y * 0.4));\n  col *= 1.0 + 0.04 * (band - 0.5);\n  return col;\n}\n\nvec3 mediumPastel(vec3 col, vec2 p, float t) {\n  vec2 flow = flowField(p, t);\n  vec2 perp = vec2(-flow.y, flow.x);\n  float along  = dot(p, flow) * uStrokeScale;\n  float across = dot(p, perp) * uStrokeScale;\n  across += 0.03 * (vnoise(p * 260.0) - 0.5);\n\n  float aniso = mix(1.0, 0.18, uStrokeAnisotropy);\n  float stroke = fbm(vec2(along * aniso, across));\n  col *= mix(1.0, 0.82 + 0.32 * stroke, uStrokeAmount);\n\n  // Pastel tooth — tiny high-frequency grain\n  float tooth = vnoise(p * 800.0);\n  col *= 1.0 - 0.08 * uStrokeAmount * (tooth - 0.5);\n  return col;\n}\n", he = "// ── Medium id constants (mirror uniformBridge.ts MEDIUM_ID — the uMedium ladder) ──\n// The StrokeProfile selector (AX.W12 + AX.W13) dispatches on the medium id: oil,\n// van-Gogh, and oil-pastel each AUTHOR a first-class profile in profileFor.\n#define MEDIUM_OIL 3\n#define MEDIUM_CRAYON 4\n#define MEDIUM_VANGOGH 5\n#define MEDIUM_OILPASTEL 6\n\n// ── Crayon — DRY wax pigment on paper tooth (first-class medium, uMedium==4) ──\n// Crayon is not strokes. It's dry pigment crumbs dragged across paper tooth — the\n// DRY tooth-multiply model: anisotropic tooth noise multiplied into the base, a hard\n// SCUMBLE that lets the lower color show through the broken upper layer, and stable\n// OKLCh broken-color pigment patches. NO sheen, NO burnish film — that waxy gloss is\n// the OIL-PASTEL deposition's signature (mediumOilPastel), distinct from dry crayon\n// (AX.W13 slice 8 F1 — the split). Crayon shares the SUBSTRATE (the structure-tensor\n// orientation, the tooth noise, OKLCh brokenColorJitter), not the dispatch body.\nvec3 mediumCrayon(vec3 col, vec2 p, float t) {\n  // Orient the deposition along the structure-tensor edge-tangent (the bridge forces\n  // tensor for crayon, so read the tensor field directly).\n  vec2 flow = structureTensorField(p, t, flowField(p, t)).xy;\n  float ang = atan(flow.y, flow.x);\n  float ca = cos(-ang), sa = sin(-ang);\n  vec2 pr = vec2(p.x * ca - p.y * sa, p.x * sa + p.y * ca);\n\n  float aniso = mix(0.45, 0.95, uStrokeAnisotropy);\n  float scale = max(uStrokeScale * 1.6, 180.0);\n\n  // ── Paper TOOTH height field — the relief the dry pigment rides. Anisotropic\n  // (squished along flow, coarse cross-flow). Peaks = ridges, valleys = paper pits.\n  float t1 = vnoise(vec2(pr.x * scale * aniso, pr.y * scale));\n  float t2 = vnoise(vec2(pr.x * scale * aniso * 0.4, pr.y * scale * 0.4) + 11.0);\n  float t3 = vnoise(vec2(pr.x * scale * aniso * 2.1, pr.y * scale * 2.1) + 23.0);\n  float tooth = 0.55 * t1 + 0.30 * t2 + 0.15 * t3;\n\n  // ── DRY TOOTH-MULTIPLY — the crayon multiplies the tooth relief into the base\n  // color (dark in the pits, bright on the ridges). The pressure (uStrokeAmount) sets\n  // how hard the multiply bites. This is the DRY model — no pigment build-up, no sheen.\n  float pressure = clamp(uStrokeAmount, 0.0, 1.0);\n  float bite = mix(0.18, 0.46, pressure);        // multiply strength\n  vec3 result = col * (1.0 - bite * (0.5 - tooth));\n\n  // ── HARD SCUMBLE — a broken upper layer that lets the lower color show through the\n  // gaps. A coarse mask gates a HARD (high-contrast) dragged pass so the dry crayon\n  // reads as broken pigment crumbs, not a smooth wash. Drag follows the tensor flow.\n  float scumbleMask = smoothstep(0.42, 0.78, vnoise(vec2(pr.x * scale * 0.55, pr.y * scale * 0.22) + 7.0));\n  float scumbleCoverage = scumbleMask * pressure * 0.7;     // hard broken upper layer\n  vec3 paper = col * mix(0.80, 0.94, pressure);             // the paper showing through\n  result = mix(paper, result, scumbleCoverage);\n\n  // ── Broken-color pigment: stable wax/pigment patches in OKLCh (W5), not flicker.\n  vec2 pigmentCell = floor(pr * max(scale * 0.18, 32.0));\n  float pigmentMask = smoothstep(0.28, 0.82, vnoise(pr * scale * 0.21 + 19.0));\n  result = brokenColorJitter(\n    result,\n    hash21(pigmentCell + 17.0),\n    hash21(pigmentCell * 2.1 + 31.0),\n    0.45 + 0.55 * pigmentMask\n  );\n\n  // Crayon is saturation-amplified (the wax chroma punch — but DRY, no sheen).\n  result = saturate3(result, 1.08);\n\n  return result;\n}\n\n// ── StrokeProfile — the oil-stroke parameter vector (AX.W12, slice 8 F6) ──────\n// The SBR \"stroke = parameter vector\" canon (facet 11): the per-mode if-ladder knobs\n// become struct FIELDS — logic-as-DATA, not an imperative branch. profileFor(medium,\n// mode) populates the profile for a (medium, mode) pair; paintStrokeLayers(profile) runs\n// the four-layer bestOil/paintOver cascade off it. A new medium (the W13 van-Gogh /\n// oil-pastel profiles) AUTHORS a profile entry — it never edits a monolith.\nstruct StrokeProfile {\n  int   shapeType;   // 0 tapered, 1 load-drag, 2 dab, 3 even, 4 comma/crescent\n  float bristleAmp;  // edge raggedness 0..0.5\n  float streakFreq;  // internal-streak spatial frequency\n  float streakAmp;   // internal-streak amplitude\n  float impastoAmp;  // paint-thickness contribution (the relight reads the height)\n  float hardness;    // edge-compositing crispness 0..1\n  float toothScale;  // canvas-tooth (linen weave) spatial frequency\n  float toothAmp;    // canvas-tooth amplitude\n  float pigmentSat;  // pigment saturation boost (OKLCh chroma scale)\n  float densityBig;  // layer-1 (big gestural) placement density gate\n  float densityMed;  // layer-2 (medium body) placement density gate\n  float densitySml;  // layer-3 (small dabs) placement density gate\n  float energyGrade; // AX.W13 — SBR energy-grade magnitude (0 off; 1 the van-Gogh\n                     // Starry-Night length cascade). The energy grade is a PROFILE\n                     // field, NOT a buried uMedium==5 bestOil branch (slice 8 F0).\n  float impastoFloor;// AX.W13 — the per-stroke height-crown floor (0.4 oil falloff;\n                     // 1.0 van-Gogh FULL-height crown so each dab catches its glint).\n  float densityFill; // AX.W13 — layer-4 fill-dab density. Oil/oil-pastel fill the\n                     // bald spots (0.95, full coverage); van-Gogh sets it LOW so the\n                     // sparse atomic dabs keep their visible inter-stroke ground gaps\n                     // (the atomicity read — a 0.95 fill would close every gap to a smear).\n  float groundFloor; // AX.W13 — the bare-ground darken floor. The gaps between strokes\n                     // (low accumulated height) multiply toward this floor, so van-Gogh's\n                     // sparse atomic dabs read over a visibly darker underpainting (the\n                     // Starry-Night visible ground). Oil/oil-pastel keep 1.0 (no darken).\n};\n\n// The (medium, mode) -> StrokeProfile selector. The if-ladder's knobs are the profile's\n// fields; a new medium adds a CASE here (never an edit to a dispatch body). mode is the\n// uStrokeMode oil-stroke sub-mode (0 oil, 1 palette-knife, 3 brushwork). crayon\n// (uMedium==4) is a PEER tooth-multiply medium that never reaches the stroke cascade.\n// The painterly STROKE mediums are oil (MEDIUM_OIL), van-Gogh (MEDIUM_VANGOGH) and\n// oil-pastel (MEDIUM_OILPASTEL) — each AUTHORS its own profile entry below.\nStrokeProfile profileFor(int medium, int mode) {\n  // The oil baseline (mode 0) — balanced modern-abstract/palette-knife hybrid.\n  StrokeProfile prof = StrokeProfile(\n    0,      // shapeType — tapered\n    0.25,   // bristleAmp\n    9.0,    // streakFreq\n    0.09,   // streakAmp\n    0.9,    // impastoAmp\n    0.80,   // hardness\n    340.0,  // toothScale — W-AUR-VANGOGH-CONFIG reconcile: FINER 240 → 340. Oil's §4.3 β sat\n            //               right AT the −1.45 band ceiling (±0.02 frame variance straddling\n            //               it, surfaced once the vangogh bands went green and the spec reached\n            //               the oil β assert). MEASURED on the live GPU: ADDING fine high-\n            //               frequency tooth energy steepens oil's power-spectrum slope (more\n            //               negative — coarsening it the other way pushed β SHALLOWER to\n            //               −1.33). A finer, stronger tooth moves β to a safe distance inside\n            //               the −5/3 band [−1.85, −1.45]. Oil-only (oil-pastel sets 280).\n    0.14,   // toothAmp — raised 0.09 → 0.14 (same oil-β steepen: more fine-grain high band).\n    1.03,   // pigmentSat\n    0.65,   // densityBig\n    0.78,   // densityMed\n    0.90,   // densitySml\n    0.0,    // energyGrade — oil is uniform-length (no Starry-Night cascade)\n    0.4,    // impastoFloor — oil's 0.4+0.6·edgeN crown falloff\n    0.92,   // densityFill — oil fills DENSELY with the directional low-relief fill so the\n            //               flat-colour regions between the macro strokes are TEXTURED with\n            //               flow-aligned fine strokes (the §4.2 anisotropy needs every pixel\n            //               to carry directional gradient; a flat untextured region reads as\n            //               an isotropic/degenerate tensor and drags mean A down). The fill\n            //               is now directional + low-streak, so dense coverage textures\n            //               WITHOUT the prior round-dab speckle.\n    1.0     // groundFloor — oil keeps the full base (no ground darken; fills everywhere)\n  );\n\n  // W-AUR-VANGOGH-REBUILD — the van-Gogh medium NO LONGER routes through this oil\n  // StrokeProfile cascade. It is a DEDICATED atomic-dab body (mediumVangogh / vangoghDab\n  // below) — the prior profile-driven cascade (four dense layers of long-thin\n  // tensor-oriented strokes) merged into the \"marbled flow-bands\" read and ran at ~4fps.\n  // profileFor now serves ONLY the oil + oil-pastel stroke mediums; the MEDIUM_VANGOGH\n  // case is removed (no caller). MEDIUM_VANGOGH stays #define'd (the main() dispatch +\n  // isPainterlyStroke still key off uMedium==5).\n\n  // ── Oil-pastel — broad smeared directional strokes deposited via the brush engine\n  // (slice 8 F1). Creamy soft edges (low hardness), heavy build-up, and a chroma punch.\n  // Distinct from the dry-crayon tooth-multiply (mediumCrayon): oil-pastel LAYS strokes.\n  if (medium == MEDIUM_OILPASTEL) {\n    prof.shapeType   = 0;     // tapered directional smear — a round dab reads locally\n                              // isotropic (steep β / low A); the creamy oil-pastel mark\n                              // is a broad DIRECTIONAL smear, which carries the §4.2\n                              // coherence + adds the directional mid-band the §4.3 slope\n                              // needs (the round-dab profile rolled β off to −2.8)\n    prof.bristleAmp  = 0.10;  // smooth creamy edges (waxy pastel, not bristle)\n    prof.streakFreq  = 18.0;  // finer internal smear — adds the mid/high-band the smooth\n                              // round-dab profile lacked (the β-too-steep fix); the\n                              // directional along-spine streak flattens the §4.3 slope\n                              // without crossing the strokes\n    prof.streakAmp   = 0.20;\n    prof.impastoAmp  = 0.7;   // build-up, but flatter than oil impasto\n    prof.hardness    = 0.42;  // CREAMY — soft compositing, strokes blend on overlap\n    prof.toothScale  = 280.0;\n    prof.toothAmp    = 0.12;  // the pastel tooth reads stronger\n    prof.pigmentSat  = 0.80;  // the waxy chroma — eased off 1.16 so the rendered C clears\n                              // the §4.1 band ceiling with margin while the K-M subtractive\n                              // overlap path (paintOverOklab) keeps the chroma OFF the grey floor\n    prof.hardness    = 0.58;  // creamy but the strokes still REGISTER — a too-soft blend\n                              // (0.42) smeared the marks into the smooth colour zones and\n                              // the field rolled off too steep (β≈−2.9, a near-gradient);\n                              // a crisper deposit keeps the mid/high-band the −5/3 slope needs\n    prof.densityBig  = 0.62;  // broad directional coverage — but NOT over-dense: too many\n                              // overlapping creamy strokes point in slightly-varied\n                              // directions and the local §4.2 coherence DROPS (over-dense\n                              // measured A=0.57 < the 0.62-density A=0.63). The mid-density\n                              // flow-aligned deposit is the coherence sweet spot.\n    prof.densityMed  = 0.78;\n    prof.densitySml  = 0.90;\n    prof.energyGrade = 0.80;  // a strong Starry-Night length cascade — varies stroke length\n                              // by luma/coherence so the dab-size energy spreads across\n                              // scales (the −5/3 cascade the uniform-length 0.0 grade lacked,\n                              // the §4.3 β-too-steep fix). Not the full van-Gogh 1.0 (oil-\n                              // pastel keeps its creamy broad character).\n    prof.impastoFloor = 0.55; // a softer crown than van-Gogh's full-height glint\n    return prof;\n  }\n\n  // ── Oil sub-modes (mode dispatch on the oil medium).\n  if (mode == 1) {           // palette knife — razor edges, heavy impasto\n    prof.shapeType  = 3;     // flat, even\n    prof.bristleAmp = 0.12;\n    prof.streakFreq = 4.0;  prof.streakAmp = 0.05;\n    prof.impastoAmp = 1.6;\n    prof.hardness   = 0.95;\n    prof.toothAmp   = 0.04;\n    prof.densityBig = 0.80; prof.densityMed = 0.88; prof.densitySml = 0.70;\n  } else if (mode == 3) {    // thick brushwork — heavy bristle brush\n    prof.shapeType  = 0;     // tapered\n    prof.bristleAmp = 0.32;\n    prof.streakFreq = 14.0; prof.streakAmp = 0.14;\n    prof.impastoAmp = 1.2;\n    prof.hardness   = 0.85;\n    prof.toothAmp   = 0.07;\n  }\n  return prof;\n}\n\n// The single parameterized four-layer stroke cascade (AX.W12). The four hand-unrolled\n// bestOil/paintOver invocations collapse into ONE body driven by the profile + the\n// uniform-derived per-layer scale/anisotropy multipliers (which are mode-INVARIANT — the\n// per-layer offsets, seeds, and len/wid muls are fixed across modes, so they stay here as\n// the cascade's structure; only the profile's knobs differentiate the medium). The\n// mode passes through for the two mode-special-cased shapes (knife layer-3 dabs +\n// layer-4 fill shape), which are placement details of the cascade, not profile knobs.\nvoid paintStrokeLayers(inout vec3 col, inout float height, StrokeProfile prof,\n                       int mode, vec2 p, float t) {\n  // Scales & multipliers from uniforms (mode-invariant cascade structure).\n  float baseScale = max(uStrokeScale * 0.006, 0.008);\n  // Three primary layers: big gestural, medium body, small dabs.\n  // W-AUR-STUDIO D5 — the −5/3 radii respacing was MEASURED + REVERTED (recorded, not\n  // silently struck). The φ-adjacent geometric candidate (2.4 / 1.45 / 0.87, a fixed\n  // ≈1.66× step) moved oil-pastel β −2.534 → −2.413 — TOWARD the −5/3 band but NOT into\n  // it (floor −1.85), so the keep-iff-into-band gate failed. The single-pass WebGL2 path\n  // is the true limiter (HC-aurora §5): the oil-pastel β + the oil/oil-pastel anisotropy\n  // residual is the named AY.W-AUR-T5 multi-pass anisotropic-Kuwahara candidate. The\n  // hand-set 2.4 / 1.1 / 0.45 stays — it is van-Gogh's landed-band spacing.\n  float sBig = baseScale * 2.4;\n  float sMed = baseScale * 1.1;\n  float sSml = baseScale * 0.45;\n\n  // Strokes are elongated HARD (long·thin) so the rendered field reads strongly\n  // directional — the §4.2 structure-tensor anisotropy lift. A short fat dab is locally\n  // isotropic and drags the measured A down; a long thin stroke that hugs the tensor\n  // tangent is what gives Starry Night its A≈0.83. The width muls are pulled tighter to\n  // keep the long·thin aspect.\n  float lenMulBig = mix(3.0, 5.2, uStrokeAnisotropy);\n  float widMulBig = mix(0.42, 0.22, uStrokeAnisotropy);\n  float lenMulMed = mix(2.8, 4.8, uStrokeAnisotropy);\n  float widMulMed = mix(0.40, 0.21, uStrokeAnisotropy);\n  float lenMulSml = mix(2.2, 3.6, uStrokeAnisotropy);\n  float widMulSml = mix(0.36, 0.23, uStrokeAnisotropy);\n\n  float jitterAmt = 0.75;   // large jitter — no grid\n  vec2 flow = flowField(p, t);\n\n  // Layer 1 — big gestural strokes (sparse, shaping)\n  StrokeHit hBig = bestOil(p, sBig, lenMulBig, widMulBig, jitterAmt * 0.55,\n                           prof.densityBig, prof.shapeType, prof.bristleAmp, flow, t, 1.3,\n                           prof.energyGrade);\n  paintOver(col, height, hBig, prof.streakFreq * 0.7, prof.streakAmp,\n            uImpasto * prof.impastoAmp * uStrokeAmount, prof.hardness, 1.3, prof.impastoFloor);\n\n  // Layer 2 — medium body strokes\n  StrokeHit hMed = bestOil(p + vec2(11.3, 3.7), sMed, lenMulMed, widMulMed,\n                           jitterAmt, prof.densityMed, prof.shapeType, prof.bristleAmp, flow, t, 2.7,\n                           prof.energyGrade);\n  paintOver(col, height, hMed, prof.streakFreq, prof.streakAmp,\n            uImpasto * prof.impastoAmp * uStrokeAmount, prof.hardness, 2.7, prof.impastoFloor);\n\n  // Layer 3 — small strokes (more frequent, smaller). Even the knife layer-3 stays the\n  // medium's DIRECTIONAL shape, not a round dab — the round-dab sparkle was a fine-scale\n  // isotropic speckle that dragged the §4.2 anisotropy down + rolled the §4.3 slope shallow.\n  int smlShape = prof.shapeType;\n  StrokeHit hSml = bestOil(p + vec2(-5.1, 8.4), sSml, lenMulSml, widMulSml,\n                           jitterAmt * 0.9, prof.densitySml, smlShape,\n                           prof.bristleAmp * 0.85, flow, t, 4.1, prof.energyGrade);\n  paintOver(col, height, hSml, prof.streakFreq * 1.4, prof.streakAmp * 0.6,\n            uImpasto * prof.impastoAmp * 0.45 * uStrokeAmount, prof.hardness, 4.1, prof.impastoFloor);\n\n  // Layer 4 — fill strokes (dense, fine) — cover bald spots. The density is\n  // PROFILE-driven (prof.densityFill): oil/oil-pastel fill to full coverage; van-Gogh\n  // keeps it sparse so the bare-ground gaps between its atomic dabs survive. The fill is\n  // ELONGATED + flow-oriented (NOT a round dab) — a dense round-dab fill is locally\n  // isotropic and was the dominant §4.2 anisotropy sink; an oriented fine stroke keeps the\n  // fill directional so the whole field reads coherent.\n  float sFill = baseScale * 0.22;\n  float lenMulFill = mix(2.0, 3.2, uStrokeAnisotropy);\n  float widMulFill = mix(0.34, 0.24, uStrokeAnisotropy);\n  int fillShape = (mode == 1) ? 3 : prof.shapeType; // knife=even; others inherit the medium's directional shape\n  StrokeHit hFill = bestOil(p + vec2(3.9, -6.2), sFill, lenMulFill, widMulFill,\n                            jitterAmt * 0.8, prof.densityFill, fillShape,\n                            prof.bristleAmp * 0.6, flow, t, 8.9, prof.energyGrade);\n  // The fill is a low-relief, low-streak base coat — it covers bald spots WITHOUT adding\n  // fine bright speckle (a high-impasto/high-streak fill caught a per-dab glint that read\n  // as isotropic noise over the directional macro strokes — the §4.2/§4.3 sink). Its\n  // positional jitter is tamed too so the fill strokes stay flow-aligned, not scattered.\n  paintOver(col, height, hFill, prof.streakFreq * 1.2, prof.streakAmp * 0.35,\n            uImpasto * prof.impastoAmp * 0.18 * uStrokeAmount, prof.hardness * 0.9, 8.9, prof.impastoFloor);\n\n  // Optional crosshatch layer\n  if (uStrokeLayers == 2) {\n    vec2 flow2 = vec2(-flow.y, flow.x);\n    StrokeHit hX = bestOil(p + vec2(7.3, -2.1), sMed, lenMulMed * 0.9, widMulMed,\n                           jitterAmt, prof.densityMed * 0.7, prof.shapeType, prof.bristleAmp, flow2, t, 6.5,\n                           prof.energyGrade);\n    paintOver(col, height, hX, prof.streakFreq, prof.streakAmp * 0.85,\n              uImpasto * prof.impastoAmp * 0.55 * uStrokeAmount, prof.hardness, 6.5, prof.impastoFloor);\n  }\n}\n\n// The shared stroke-medium SUBSTRATE finish (AX.W12 + AX.W13): paint the four-layer\n// cascade off the profile, then canvas-tooth + relight + saturation. The three painterly\n// STROKE mediums (oil / van-Gogh / oil-pastel) each AUTHOR their own profile and call\n// this — the substrate is shared, the PROFILE (not a dispatch-body fork) differentiates\n// the medium (slice 8 F6). mode passes through for the cascade's mode-special-cased shapes.\nvec3 paintStrokeMedium(vec3 col, vec2 p, float t, StrokeProfile prof, int mode) {\n  vec3 result = col;\n  // AW.W4.2 — accumulated paint HEIGHT across the stroke layers. The relight reads\n  // its gradient for the normal; the canvas tooth seeds the base relief.\n  float height = 0.0;\n  paintStrokeLayers(result, height, prof, mode, p, t);\n\n  // Canvas tooth — linen weave\n  float tooth1 = vnoise(p * prof.toothScale);\n  float tooth2 = vnoise(p * prof.toothScale * vec2(0.6, 2.4) + 37.0);\n  float tooth  = (0.6 * tooth1 + 0.4 * tooth2) - 0.5;\n  result *= 1.0 + tooth * prof.toothAmp * uCanvasGrain;\n\n  // AW.W4.2 — relight the accumulated paint height with the movable uLightDir\n  // source (diffuse + Blinn specular, in LINEAR before aces()). The canvas tooth\n  // is the base relief term so the weave also catches the raking light.\n  float canvasBase = tooth * prof.toothAmp * 0.5;\n  result = relightImpasto(result, height, canvasBase);\n\n  // AX.W13 — the visible ground. The bare-ground gaps (low accumulated stroke height)\n  // multiply toward prof.groundFloor, so van-Gogh's sparse atomic dabs read as separable\n  // marks over a darker underpainting (oil/oil-pastel keep groundFloor=1.0 → no-op). The\n  // smoothstep keeps the impasto dabs (high height) at full value; only the open ground\n  // between them darkens.\n  result *= mix(prof.groundFloor, 1.0, smoothstep(0.0, 0.38, height));\n\n  // Pigment saturation boost\n  result = saturate3(result, prof.pigmentSat);\n\n  return result;\n}\n\nvec3 mediumOil(vec3 col, vec2 p, float t) {\n  // The oil medium is a thin body: fetch the oil profile (mode-dispatched), paint via\n  // the shared stroke substrate.  uStrokeMode: 0 oil (gestural), 1 knife, 3 brushwork.\n  int mode = uStrokeMode;\n  StrokeProfile prof = profileFor(MEDIUM_OIL, mode);\n  return paintStrokeMedium(col, p, t, prof, mode);\n}\n\n" + pe + "// AX.W13 — the oil-pastel medium (uMedium==6) is a DISTINCT stroke-deposition body,\n// split out of the dry-crayon tooth-multiply (mediumCrayon) (slice 8 F1). It DEPOSITS\n// broad smeared directional strokes via the same brush engine with a creamy soft\n// hardness, heavy pigment build-up where strokes overlap, and a chroma punch — the\n// stroke-deposition model that gives oil-pastel its depth and less-uniform read.\n// Shares the SUBSTRATE (the stroke cascade + tooth + relight) with oil/van-Gogh, NOT\n// the dispatch body of dry crayon. The bridge forces tensor orientation.\nvec3 mediumOilPastel(vec3 col, vec2 p, float t) {\n  StrokeProfile prof = profileFor(MEDIUM_OILPASTEL, 0);\n  return paintStrokeMedium(col, p, t, prof, 0);\n}", A = "\n", ge = `#version 300 es
+var he = "#version 300 es\nin vec2 aPos;\nout vec2 vUv;\nvoid main() {\n  vUv = aPos * 0.5 + 0.5;\n  gl_Position = vec4(aPos, 0.0, 1.0);\n}", ge = "\nconst float CURL_EPS = 0.012;\n\n// Scalar fbm potential psi — the host's noise basis, wrapped so the curl operator is\n// basis-agnostic. A host carrying its own fbm defines this as: return fbm(p);\nfloat potentialFBM(vec2 p);\n\n// ∇×ψ in 2D — the divergence-free flow direction. Central finite differences of the\n// scalar potential give the partials; the cross-pairing (∂y, −∂x) is the 2D curl.\nvec2 curlFBM(vec2 p) {\n  float dx = potentialFBM(p + vec2(CURL_EPS, 0.0)) - potentialFBM(p - vec2(CURL_EPS, 0.0));\n  float dy = potentialFBM(p + vec2(0.0, CURL_EPS)) - potentialFBM(p - vec2(0.0, CURL_EPS));\n  // ψ-gradient is (dx, dy)/(2·eps); the curl rotates it 90°: (∂ψ/∂y, −∂ψ/∂x).\n  vec2 g = vec2(dx, dy) / (2.0 * CURL_EPS);\n  return vec2(g.y, -g.x);\n}\n", _e = "// ── Palette LUT (W5 — OKLCh/OKLab interpolation) ──────────────────────────────\n// The endpoints stay CPU-baked to LINEAR sRGB (the Aras precompute pattern); only\n// the BETWEEN-endpoint interpolation runs in perceptual space. AX.W11 — the ramp\n// itself (the smoothstep ease + the OKLab-rectangular-vs-OKLCh-hue-arc dispatch on\n// uHuePath) is the shared samplePaletteRamp from procedural-color.glsl, the SAME\n// source the WebGPU twin splices; samplePalette here only selects the bracketing\n// stop pair + the raw inter-stop t and hands them to the shared ramp.\nvec3 samplePalette(float id) {\n  id = clamp(id, 0.0, 1.0);\n  float scaled = id * float(uStopCount - 1);\n  int i0 = int(floor(scaled));\n  int i1 = min(i0 + 1, uStopCount - 1);\n  float t = fract(scaled);\n  return samplePaletteRamp(uPalette[i0], uPalette[i1], t, uHuePath);\n}\n\n// ── Nuclei field ──────────────────────────────────────────────────────────\nvoid nucleiField(vec2 p, float t, out float paletteId, out float valueMod) {\n  float accumBias  = 0.0;\n  float accumValue = 0.0;\n  float accumW     = 0.0;\n  for (int i = 0; i < MAX_NUCLEI; i++) {\n    if (i >= uNucleiCount) break;\n    vec2 posI = uNucleiPos[i]\n              + uNucleiDriftRadius[i] * vec2(\n                  cos(t * uNucleiDrift * K_NUCLEI + uNucleiDriftPhase[i]),\n                  sin(t * uNucleiDrift * K_NUCLEI + uNucleiDriftPhase[i] * 1.13)\n                );\n    vec2 diff = p - posI;\n    // Anisotropic Gaussian: rotate diff into the nucleus's local frame\n    // (major axis along uNucleiAngle), then scale the major-axis component by\n    // 1/elongation so the squared distance describes an ellipse. Defaults\n    // 1.0/0.0 reduce to the isotropic dot(diff, diff).\n    float ca = cos(uNucleiAngle[i]);\n    float sa = sin(uNucleiAngle[i]);\n    vec2  local = vec2( ca * diff.x + sa * diff.y,\n                       -sa * diff.x + ca * diff.y);\n    float along  = local.x / max(uNucleiElong[i], 0.01);\n    float across = local.y;\n    float d2 = along * along + across * across;\n    float r = max(uNucleiRadius[i], 0.01);\n    float w = exp(-uSoftmaxBeta * d2 / (r * r));\n    accumBias  += w * uNucleiPaletteBias[i];\n    accumValue += w * uNucleiValueBias[i];\n    accumW     += w;\n  }\n  paletteId = accumBias  / max(accumW, 1e-4);\n  valueMod  = accumValue / max(accumW, 1e-4);\n\n  // Palette drift — slow global paletteId breathe between adjacent stops. The\n  // rate rides K_PAL (perceptible ~30–60s hue cycle at the default coefficient);\n  // the amplitude lifts the authored coefficient into a 0.03..0.06 paletteId band\n  // so the palette visibly travels rather than dithering within one stop.\n  float palAmp = clamp(uPaletteDrift * 6.0, 0.0, 0.06);\n  paletteId += palAmp * sin(t * uPaletteDrift * K_PAL);\n  paletteId = clamp(paletteId, 0.0, 1.0);\n}\n", ve = "// ── Flow field ────────────────────────────────────────────────────────────\n// Structure-tensor prototype (body defined post-sampleBase in the mediums block).\n// Returns vec3(dir.x, dir.y, coherence) — the minor-eigenvector tangent + A.\nvec3 structureTensorField(vec2 p, float t, vec2 fallbackDir);\n\nvec2 flowField(vec2 p, float t) {\n  vec2 dir = vec2(1.0, 0.0);\n  if (uFlowPattern == 1) {\n    // radial\n    dir = normalize(p - uFlowFocal + 1e-4);\n  } else if (uFlowPattern == 2) {\n    // swirl — tangent to radial\n    vec2 rad = normalize(p - uFlowFocal + 1e-4);\n    dir = vec2(-rad.y, rad.x);\n  } else if (uFlowPattern == 3) {\n    // diagonal\n    float a = radians(uFlowAngle);\n    dir = vec2(cos(a), sin(a));\n  } else if (uFlowPattern == 4) {\n    // multi — curl-noise driven\n    float n = fbm(p * 2.0 + t * 0.02);\n    float a = n * 6.2831;\n    dir = vec2(cos(a), sin(a));\n  } else if (uFlowPattern == 5) {\n    // tensor / ETF — the color field's OWN edge-tangent flow (AW.W4.1). The\n    // diagonal angle seeds the smooth fallback the coherence-weighted blend\n    // relaxes toward in low-structure zones. Returns directly (no curl/cursor\n    // post-mix — the tensor field already carries the structure orientation;\n    // the cursor warp still drives the underlying color field in domainWarp).\n    float a = radians(uFlowAngle);\n    vec2 fallback = vec2(cos(a), sin(a));\n    return structureTensorField(p, t, fallback).xy;\n  }\n  // curl — perturb by noise. Radial gets much less curl so rays stay clean.\n  if (uFlowCurl > 0.0) {\n    float n = fbm(p * 3.0) - 0.5;\n    float curlAmt = uFlowCurl;\n    if (uFlowPattern == 1) curlAmt *= 0.25;  // radial: preserve ray clarity\n    if (uFlowPattern == 2) curlAmt *= 0.55;  // swirl: allow moderate curl\n    float a = atan(dir.y, dir.x) + n * 3.14159 * curlAmt;\n    dir = vec2(cos(a), sin(a));\n  }\n  // cursor influence — swirl around the cursor position\n  if (uCursorStrength > 0.001) {\n    vec2 toCur = uCursor - p;\n    float d = length(toCur);\n    float r = max(uCursorRadius, 0.01);\n    float w = exp(-(d * d) / (r * r * 0.5));\n    // swirl tangent\n    vec2 tangent = vec2(-toCur.y, toCur.x) / max(d, 1e-4);\n    float a0 = atan(dir.y, dir.x);\n    float a1 = atan(tangent.y, tangent.x);\n    // blend angle\n    float da = a1 - a0;\n    da = atan(sin(da), cos(da)); // wrap to [-pi, pi]\n    float a = a0 + da * w * uCursorStrength;\n    dir = vec2(cos(a), sin(a));\n  }\n  // AW.W8.1 — velocity-reactive swirl-BURST: a fast flick (uCursorBurst, decaying\n  // over ~1s) injects a transient flow bias along the pointer velocity, distinct from\n  // the steady attraction above (the stateless swirl is SUPERSEDED by this form). The\n  // burst is gated by the master tempo scalar on the CPU side (it decays to 0 under\n  // PRM/pause), so this term cannot leak motion under reduced-motion.\n  if (uCursorBurst > 0.001) {\n    vec2 toCur = uCursor - p;\n    float d = length(toCur);\n    float r = max(uCursorRadius, 0.01);\n    float w = exp(-(d * d) / (r * r * 0.7));\n    // Bias the flow toward the velocity direction, scaled by the decaying burst.\n    vec2 velDir = uCursorVelocity / max(length(uCursorVelocity), 1e-4);\n    float a0 = atan(dir.y, dir.x);\n    float a1 = atan(velDir.y, velDir.x);\n    float da = atan(sin(a1 - a0), cos(a1 - a0));\n    float a = a0 + da * w * uCursorBurst * 0.8;\n    dir = vec2(cos(a), sin(a));\n  }\n  return dir;\n}\n", ye = "// Khronos PBR-Neutral tonemap (hue + saturation preserving)\nvec3 aces(vec3 color) {\n  const float startCompression = 0.8 - 0.04;\n  const float desaturation = 0.15;\n  float x = min(color.r, min(color.g, color.b));\n  float offset = x < 0.08 ? x - 6.25 * x * x : 0.04;\n  color -= offset;\n  float peak = max(color.r, max(color.g, color.b));\n  if (peak < startCompression) return clamp(color, 0.0, 1.0);\n  float d = 1.0 - startCompression;\n  float newPeak = 1.0 - d * d / (peak + d - startCompression);\n  color *= newPeak / peak;\n  float g = 1.0 - 1.0 / (desaturation * (peak - newPeak) + 1.0);\n  return clamp(mix(color, newPeak * vec3(1.0), g), 0.0, 1.0);\n}\n", be = "// ── Curved swept brushstroke primitive ────────────────────────────────────\n// A stroke is a curved spine from A to B with quadratic bulge K (sideways offset\n// at midpoint). Width varies along length via a shape profile. Edge is ragged\n// (bristle-modulated), ends are rounded blobs, inside has streak modulation.\n//\n// For pixel p, we invert an approximate projection onto the curved spine:\n//   1) project p onto straight AB axis -> along0, cross0\n//   2) at parameter along0, spine offset = 4·K·along0·(1-along0) perpendicular\n//   3) cross = cross0 - spineOffset, refine along one iteration\n// Coverage blends paint atop col based on per-stroke color sampled from the\n// underlying base at spine midpoint — so overlapping strokes from different\n// palette regions meet at hard, ragged demarcations.\n\nstruct StrokeHit {\n  float coverage;  // 0..1\n  vec3  color;\n  float alongT;    // 0..1 along spine, for internal modulation\n  float crossN;    // cross distance / current half-width, signed\n  float edgeN;     // distance to edge in half-widths (0 = at edge, 1 = at spine)\n};\n\nStrokeHit noHit() { return StrokeHit(0.0, vec3(0.0), 0.0, 0.0, 0.0); }\n\nvec2 rotateDir(vec2 dir, float angle) {\n  float ca = cos(angle), sa = sin(angle);\n  return vec2(dir.x * ca - dir.y * sa, dir.x * sa + dir.y * ca);\n}\n\nvec2 safeDir(vec2 dir) {\n  float len = length(dir);\n  return len > 1e-4 ? dir / len : vec2(1.0, 0.0);\n}\n\n// Shape profile along the stroke. type:\n//   0 tapered    — thin-fat-thin (classic brush)\n//   1 load-drag  — fat start, tapers to a point\n//   2 dab        — ellipse-like, round center\n//   3 even       — near-constant, slight end taper\n//   4 comma      — van-Gogh atomic dab: ASYMMETRIC taper, fat near one end and\n//                  drawn to a fine point at the other (the divisionist crescent\n//                  mark). The width peak sits off-centre (~0.32) so the two ends\n//                  taper at DIFFERENT rates — a loaded round head + a dragged tail.\nfloat strokeShape(float t, int type) {\n  t = clamp(t, 0.0, 1.0);\n  if (type == 1) {\n    // loaded at t=0, tapering\n    return pow(1.0 - t, 0.55) * smoothstep(0.0, 0.08, t);\n  } else if (type == 2) {\n    // dab / blob\n    float d = t - 0.5;\n    return exp(-d * d * 12.0);\n  } else if (type == 3) {\n    // mostly even, slight end softening\n    return smoothstep(0.0, 0.08, t) * smoothstep(1.0, 0.92, t) * 0.95 + 0.05;\n  } else if (type == 4) {\n    // comma / crescent — width peaks at ~0.32 then drags to a fine tail. The two\n    // sides of the peak fall at different exponents (a fat round head, a long thin\n    // tail), so the mark reads as a directional van-Gogh comma, not a symmetric oval.\n    float peak = 0.32;\n    float head = smoothstep(0.0, peak, t);                            // rounded loaded head\n    float tail = pow(clamp((1.0 - t) / (1.0 - peak), 0.0, 1.0), 1.6); // dragged thin tail\n    return clamp(min(head, tail) * 1.04, 0.0, 1.0);\n  }\n  // 0 tapered\n  return smoothstep(0.0, 0.22, t) * smoothstep(1.0, 0.78, t);\n}\n\n// A single curved stroke.\n//   a, b           — endpoints\n//   halfW          — base half-width (world units)\n//   bulge          — signed perpendicular midpoint offset, in world units\n//   shapeType      — 0..3\n//   bristleFreq    — spatial frequency of edge raggedness\n//   bristleAmp     — 0..0.5 fraction of halfW chewed away at edge extrema\n//   streakSeed     — uniqueness seed for internal streaks\n//   colAtMid       — pre-sampled color at midpoint\nStrokeHit curvedStroke(vec2 p, vec2 a, vec2 b, float halfW,\n                       float bulge, int shapeType,\n                       float bristleFreq, float bristleAmp,\n                       float streakSeed, vec3 colAtMid) {\n  vec2 ab = b - a;\n  float L = length(ab);\n  if (L < 1e-5) return noHit();\n  vec2 tang = ab / L;\n  vec2 norm = vec2(-tang.y, tang.x);\n\n  vec2 rel = p - a;\n  float along0 = dot(rel, tang) / L;      // straight projection 0..1 along AB\n  float cross0 = dot(rel, norm);          // signed cross\n\n  // Spine sideways offset due to bulge (parabolic)\n  float bend = 4.0 * bulge * along0 * (1.0 - along0);\n  float cross1 = cross0 - bend;\n\n  // One refinement step: as the spine bends, the closest-point along shifts.\n  // Spine tangent differs from straight tangent; first-order correction:\n  float dBend = 4.0 * bulge * (1.0 - 2.0 * along0);  // d(bend)/d(along)\n  float along1 = along0 + (cross1 * dBend) / (L * (1.0 + dBend * dBend / (L * L)));\n\n  // Bristle-ragged edge: half-width gets chewed by 1D noise along the spine\n  float edgeNoise = fbm(vec2(along1 * bristleFreq, streakSeed * 7.3)) - 0.5;\n  float edgeNoise2 = fbm(vec2(along1 * bristleFreq * 2.3, streakSeed * 3.1 + 17.0)) - 0.5;\n  float edgeMod = 1.0 - bristleAmp * (0.6 + 0.4 * edgeNoise) * (0.5 + edgeNoise2);\n\n  // Width profile along the stroke\n  float shape = strokeShape(along1, shapeType);\n  float halfWNow = halfW * shape * edgeMod;\n\n  // Inside-segment coverage\n  float cov = 0.0;\n  if (along1 >= 0.0 && along1 <= 1.0 && halfWNow > 1e-6) {\n    float cn = abs(cross1) / halfWNow;\n    // fwidth-derived AA: half-width from the field's screen-space gradient so the\n    // edge holds ~1px regardless of zoom/DPR (mirrors metaball.frag.ts:252-255).\n    float aaC = max(fwidth(cn), 1e-6);\n    cov = 1.0 - smoothstep(1.0 - aaC, 1.0 + aaC, cn);\n  }\n\n  // End-cap blobs (rounded tips, not perpendicular cuts)\n  float capA = 0.0, capB = 0.0;\n  {\n    vec2 capCenterA = a + norm * (bulge * 0.0); // at a\n    vec2 capCenterB = a + tang * L + norm * bend * 0.0 + norm * 0.0; // at b; bend is 0 at endpoints\n    capCenterB = b;\n    float endShapeA = strokeShape(0.02, shapeType);\n    float endShapeB = strokeShape(0.98, shapeType);\n    float rA = halfW * endShapeA * (0.9 + 0.2 * fbm(vec2(streakSeed * 13.0)));\n    float rB = halfW * endShapeB * (0.9 + 0.2 * fbm(vec2(streakSeed * 19.0 + 5.0)));\n    float dA = length(p - capCenterA);\n    float dB = length(p - capCenterB);\n    // fwidth-derived AA on the clean Euclidean cap distances (mirrors metaball.frag.ts:252-255).\n    float aaA = max(fwidth(dA), 1e-6);\n    float aaB = max(fwidth(dB), 1e-6);\n    capA = 1.0 - smoothstep(rA - aaA, rA + aaA, dA);\n    capB = 1.0 - smoothstep(rB - aaB, rB + aaB, dB);\n    // Only apply caps if we're BEYOND the segment; otherwise the segment wins.\n    capA *= (along1 < 0.05) ? 1.0 : 0.0;\n    capB *= (along1 > 0.95) ? 1.0 : 0.0;\n  }\n\n  float coverage = max(cov, max(capA, capB));\n  if (coverage < 0.002) return noHit();\n\n  // Edge distance (for impasto later): 0 at edge, 1 at spine\n  float edgeDist = halfWNow > 1e-6 ? clamp(1.0 - abs(cross1) / halfWNow, 0.0, 1.0) : 0.0;\n\n  return StrokeHit(\n    coverage,\n    colAtMid,\n    clamp(along1, 0.0, 1.0),\n    halfWNow > 1e-6 ? cross1 / halfWNow : 0.0,\n    edgeDist\n  );\n}\n\n// The painterly STROKE mediums (those that deposit via paintStrokeLayers/paintOver):\n// oil (3), van-Gogh (5), oil-pastel (6). For these the stroke OVER-composite runs in\n// OKLab so overlapping complementaries mix through a chromatic path, not the linear-\n// RGB grey mud; the smooth/atmospheric pole never reaches paintOver so it pays nothing.\nbool isPainterlyStroke() {\n  return uMedium == 3 || uMedium == 5 || uMedium == 6;\n}\n\n// OKLab OVER-composite: lerp L,a,b of the over-color toward the under-color (the\n// Ottosson matrices spliced from procedural-color.glsl.ts — no new color math). Two\n// complementary-hued strokes overlapping transition through OKLab, where the midpoint\n// holds chroma, instead of the linear mix() that desaturates toward grey (slice 8 F2).\nvec3 paintOverOklab(vec3 under, vec3 over, float alpha) {\n  vec3 labU = linOklab(under);\n  vec3 labO = linOklab(over);\n  return max(oklabToLinearSrgb(mix(labU, labO, alpha)), vec3(0.0));\n}\n\n// Composite stroke over col with internal streaking + ACCUMULATE paint HEIGHT.\n//   streakAmp    — 0..0.2 how much internal streaks darken/lighten\n//   impastoAmp   — 0..1 paint-thickness contribution (the height the relight reads)\n//   hardness     — 0..1 how crisp the compositing transition is (1 = razor, 0 = creamy)\n//\n// AW.W4.2 — the fixed-RGB edge rim (the phantom upper-left light) is RETIRED. Each\n// stroke now deposits PAINT HEIGHT into the height accumulator (coverage x per-stroke\n// thickness, perturbed by the bristle/streak fbm for ridges/grooves); mediumOil\n// derives a normal from the accumulated height gradient (dFdx/dFdy) and relights it\n// once with the movable uLightDir source (relightImpasto), in linear before aces().\n//\n// AX.W13 — the internal streak now perturbs HUE + CHROMA in OKLCh (not value-only),\n// so a single stroke carries a small hue gradient (broken color at the ATOM level —\n// adjacent strokes shimmer like real impasto, not flat stamped swatches; slice 8 F5);\n// and the OVER-composite runs in OKLab on the painterly stroke mediums so overlapping\n// complementaries mix as pigment, not the linear-mix grey (slice 8 F2).\n//   impastoFloor — the per-stroke height-crown floor: oil's 0.4+0.6*edgeN reads\n//                  floor=0.4 (the crown tapers to 40% at the stroke edge); the van-Gogh\n//                  profile passes floor=1.0 so each atomic dab catches a FULL-height\n//                  glint (scope item 1d). The crown magnitude is floor + (1-floor)*edgeN.\nvoid paintOver(inout vec3 col, inout float height, StrokeHit s,\n               float streakFreq, float streakAmp,\n               float impastoAmp, float hardness, float streakSeed,\n               float impastoFloor) {\n  if (s.coverage < 0.002) return;\n  float strokeOpacity = clamp(uStrokeAmount, 0.0, 1.0);\n  if (strokeOpacity <= 0.001) return;\n  vec3 c = s.color;\n\n  // Internal streaks — fbm along spine, modulated by cross position. Gives\n  // the loaded-brush look: some bristles carry more pigment than others.\n  float streakA = fbm(vec2(s.alongT * streakFreq, s.crossN * 2.7 + streakSeed));\n  float streakB = fbm(vec2(s.alongT * streakFreq * 0.6 + streakSeed * 3.7, s.crossN * 4.1));\n  float streak = 0.6 * (streakA - 0.5) + 0.4 * (streakB - 0.5);\n\n  // ── WITHIN-stroke OKLCh broken color (AX.W13, slice 8 F5). The streak fBm perturbs\n  // the stroke's HUE + CHROMA along alongT/crossN, seeded per-stroke, so a single\n  // stroke reads a small hue gradient (the atom-level broken-color shimmer). Value is\n  // carried by the same streak so the loaded-brush light/dark survives. Gated to the\n  // painterly stroke mediums + uBrokenColor (the smooth pole keeps the cheap value-only\n  // streak). A second decorrelated fBm drives the hue axis so hue ≠ a scaled value.\n  float hueStreak = fbm(vec2(s.alongT * streakFreq * 0.45 + streakSeed * 5.3,\n                             s.crossN * 1.9 + streakSeed * 1.7)) - 0.5;\n  if (isPainterlyStroke() && uBrokenColor > 0.001) {\n    vec3 lch = oklabToOklch(linOklab(c));\n    // ±~10° within-stroke hue swing + ±~9% chroma, scaled by uBrokenColor. Smaller\n    // than the per-stroke cell jitter (brokenColorJitter) — this is the intra-stroke\n    // gradient layered on top of the per-cell pigment patch.\n    lch.z += hueStreak * (20.0 * PI / 180.0) * uBrokenColor;\n    lch.y = max(lch.y * (1.0 + (streak) * 0.18 * uBrokenColor), 0.0);\n    c = max(oklabToLinearSrgb(oklchToOklab(lch)), vec3(0.0));\n  }\n  // Value streak (the loaded-brush light/dark across the bristles) — kept for every\n  // medium that reaches paintOver (the smooth pole never does). The along-spine value\n  // modulation is the loaded-brush light/dark that gives each stroke its OWN directional\n  // luminance band — the structure-tensor signal a low-contrast color field (a uniform-hue\n  // oil field) has no other source for (the §4.2 anisotropy lift on the flat-palette case).\n  c *= 1.0 + streak * streakAmp * 3.2;\n\n  // Cross-stroke ridge/valley shading — a catch-light along the spine centre, darkening to\n  // the stroke edges. This is the strongest DIRECTIONAL luminance cue: it carves a bright\n  // ridge down the stroke's length and dark grooves between adjacent strokes, so a field of\n  // parallel strokes reads as parallel light/dark bands (the high-A van-Gogh signature) even\n  // where the underlying colour is flat. s.edgeN is 0 at the edge, 1 at the spine.\n  float ridgeShade = smoothstep(0.0, 0.55, s.edgeN);\n  c *= 1.0 + (ridgeShade - 0.5) * 0.30;\n\n  float softLimit = mix(0.35, 0.98, hardness);\n  float alpha = smoothstep(0.0, 1.0 - softLimit, s.coverage) * strokeOpacity;\n  // OKLab OVER-composite on the painterly stroke mediums (overlap mixes as pigment,\n  // not grey); the smooth pole never reaches paintOver, so it keeps the cheap path.\n  col = isPainterlyStroke() ? paintOverOklab(col, c, alpha) : mix(col, c, alpha);\n\n  // ── Accumulate paint HEIGHT (AW.W4.2). Thickness peaks at the stroke spine and\n  // falls toward the impastoFloor at the edge (edgeN: 0=edge, 1=spine), perturbed by\n  // the bristle streak so ridges/grooves form. The relight reads the gradient of this\n  // field. The van-Gogh profile's floor=1.0 gives every atomic dab a full-height crown.\n  float ridge = 0.7 + 0.3 * (streak + 0.5);\n  float crown = impastoFloor + (1.0 - impastoFloor) * s.edgeN;\n  height += s.coverage * impastoAmp * crown * ridge * alpha;\n}\n\n// AW.W4.2 — relight the accumulated paint height. Derive a normal from the height\n// gradient (dFdx/dFdy — derivatives are in-pattern; fwidth is already used for AA),\n// then apply diffuse + Blinn specular from the movable uLightDir/uLightColor source.\n// LINEAR-space lighting BEFORE aces() (the tonemap/OETF stay locked); thin strokes\n// inherit a low shininess, thick impasto raises it (high-height ridges catch a sharp\n// glint). canvasBase is the tooth height so the canvas weave also catches the light.\nvec3 relightImpasto(vec3 col, float height, float canvasBase) {\n  float h = height + canvasBase;\n  // Height-gradient normal. heightScale trades relief depth for surface flatness.\n  float heightScale = 6.0;\n  vec3 N = normalize(vec3(-dFdx(h) * 40.0, -dFdy(h) * 40.0, 1.0 / heightScale));\n  vec3 L = normalize(uLightDir);\n  vec3 V = vec3(0.0, 0.0, 1.0);            // view straight on\n  vec3 H = normalize(L + V);\n  float diff = max(dot(N, L), 0.0);\n  // Thick impasto → high shininess (sharp glint); thin paint → broad soft sheen.\n  float shininess = mix(8.0, 64.0, clamp(height * 2.0, 0.0, 1.0));\n  float spec = pow(max(dot(N, H), 0.0), shininess);\n  // Modulate the lit terms by how much paint sits here (no relief on bare ground).\n  float relief = clamp(h * 1.5, 0.0, 1.0);\n  vec3 lit = col * (1.0 + uImpasto * relief * (diff - 0.5) * 0.5 * uLightColor);\n  // The specular glint is gated to ACCUMULATED-thick paint (height²) — a per-fine-stroke\n  // sharp glint on every thin fill dab read as bright isotropic SPECKLE that buried the\n  // directional stroke bands (the §4.2 anisotropy / §4.3 slope sink on the low-contrast\n  // oil field). Only where real impasto has built up does the raked glint fire, so the\n  // highlight tracks the macro ridges, not the fine speckle.\n  float specGate = clamp(height, 0.0, 1.0);\n  lit += uImpasto * relief * spec * specGate * 0.5 * uLightColor;\n  return lit;\n}\n\n// Best-of-9-neighbor placement: sample 3x3 surrounding cells and take the\n// stroke that covers this pixel most. Breaks the grid by per-cell jitter and\n// sparse density via noise thresholding.\n//   energyGrade — the SBR energy-grade magnitude (0 = none, the oil/oil-pastel\n//                 baseline; 1 = the full van-Gogh Starry-Night cascade — long\n//                 confident strokes in the lights/coherent zones, short dabs in\n//                 the darks/flat zones). The van-Gogh profile carries 1.0 (the\n//                 energy grade is a PROFILE field, not a buried uMedium==5 branch).\nStrokeHit bestOil(vec2 p, float cellSize, float lenMul, float halfWMul,\n                  float jitterAmt, float density, int shapeType,\n                  float bristleAmp, vec2 flow, float t, float seed,\n                  float energyGrade) {\n  vec2 cell = floor(p / cellSize);\n\n  StrokeHit best = noHit();\n  // 3x3 neighborhood\n  for (int dy = -1; dy <= 1; dy++) {\n    for (int dx = -1; dx <= 1; dx++) {\n      vec2 cc = cell + vec2(float(dx), float(dy));\n      vec2 hh = hash22(cc + seed) - 0.5;\n      // Density gate — noise threshold; sparse placement.\n      float gate = hash21(cc * 1.7 + seed * 0.3);\n      if (gate > density) continue;\n\n      vec2 center = (cc + 0.5 + hh * jitterAmt) * cellSize;\n\n      // Per-stroke direction source — the AW.W4.1 strokeOrient switch:\n      //   uStrokeOrient==0 (flow)   — the layer-provided hand-authored flowField.\n      //   uStrokeOrient==1 (tensor) — the structure-tensor minor eigenvector at the\n      //                               cell center (the color field's edge-tangent),\n      //                               so strokes hug the color zones (van Gogh).\n      // coh carries the per-cell coherence A for the W4.3 energy grading (long\n      // confident strokes where structure is coherent, stubby dabs in flat zones).\n      vec2 f;\n      float coh = 1.0;\n      if (uStrokeOrient == 1) {\n        vec3 tf = structureTensorField(center, t, safeDir(flow));\n        f = safeDir(tf.xy);\n        coh = tf.z;\n      } else {\n        f = safeDir(flow);\n      }\n      // Per-stroke deterministic local perturbation so alternate stroke layers stay\n      // live. Tensor strokes get FAR less angular jitter (the field already orients them\n      // — over-jittering the tensor strokes scatters their directions and the rendered\n      // field reads isotropic, the §4.2 anisotropy floor). The flow-oriented (non-tensor)\n      // strokes keep the full jitter (they have no per-cell orientation source).\n      float jScale = (uStrokeOrient == 1) ? 0.15 : 1.0;\n      float angJ = (hash21(cc + seed + 11.0) - 0.5) * 0.9 * jScale;  // +/- 0.45/0.07 rad\n      float curlScale = (uStrokeOrient == 1) ? 0.30 : 1.0;\n      float localCurl = (fbm(center * (2.6 + seed * 0.11) + seed * 1.9) - 0.5) * 0.55 * uFlowCurl * curlScale;\n      vec2 dir = rotateDir(f, angJ + localCurl);\n\n      // AX.W13 — the SBR energy grade is a PROFILE field (energyGrade), not a buried\n      // uMedium==5 branch. The van-Gogh profile passes energyGrade=1.0 (the full\n      // Starry-Night cascade); oil/oil-pastel pass 0.0 (uniform length). The grade\n      // modulates stroke length by local luminance (bright → long confident strokes,\n      // dark → short dabs, the Kolmogorov/Batchelor congruence) AND by coherence\n      // (coherent zones → long, flat → stubby), lerped in by energyGrade so a profile\n      // can dial the cascade from off (0) to full (1).\n      float energy = 1.0;\n      if (energyGrade > 0.001) {\n        float luma = clamp(dot(sampleBase(center, t), W_LUMA), 0.0, 1.0);\n        float graded = mix(0.5, 1.5, luma) * mix(0.6, 1.0, coh);\n        energy = mix(1.0, graded, energyGrade);\n      }\n\n      float lenV = cellSize * lenMul * (0.65 + 0.55 * hash21(cc + seed + 23.0)) * energy;\n      float halfW = cellSize * halfWMul * (0.70 + 0.55 * hash21(cc + seed + 41.0));\n      float bulge = (hash21(cc + seed + 53.0) - 0.5) * lenV * 0.35;\n\n      vec2 a = center - dir * (lenV * 0.5);\n      vec2 b = center + dir * (lenV * 0.5);\n\n      vec3 colMid = brokenColorJitter(\n        sampleBase(center, t),\n        hash21(cc + seed + 89.0),\n        hash21(cc * 2.3 + seed + 97.0),\n        1.0\n      );\n\n      StrokeHit h = curvedStroke(p, a, b, halfW, bulge, shapeType,\n                                 7.0, bristleAmp,\n                                 hash21(cc + seed + 67.0), colMid);\n      if (h.coverage > best.coverage) best = h;\n    }\n  }\n  return best;\n}\n", xe = "// ── Van-Gogh — DEDICATED atomic-dab body (W-AUR-VANGOGH-REBUILD) ─────────────\n// REBUILT FROM FIRST PRINCIPLES. The prior path routed van-Gogh through the shared\n// oil paintStrokeMedium cascade — four DENSE layers of LONG-THIN tensor-oriented\n// strokes (lenMul up to 5.2, widMul 0.22) hugging a continuous coherent field. Packed\n// that tight, thin strokes that follow continuous iso-bands MERGE into smooth flowing\n// filaments — the \"marbled flow-bands\" read (a spirograph/fingerprint, not impasto),\n// AND the per-cell 8-tap structure tensor × 4 layers × 9 cells = ~324 sampleBase/pixel\n// made it run at ~4fps. Van Gogh's sky is NOT continuous flow — it is DISCRETE, SHORT,\n// LOADED comma dabs, separated by visible ground, queued into swirl ROWS.\n//\n// The rebuild is a purpose-built two-layer DAB pass, disjoint from the oil cascade:\n//   • orientation from the analytic SWIRL flowField (van-Gogh skies ARE swirl rows) +\n//     per-cell angular jitter — NOT the expensive per-cell tensor (the lag source, and\n//     the marble source: a jittered discrete dab can NOT merge into a continuous\n//     filament the way a tensor-locked thin stroke does).\n//   • SHORT comma dabs (len/wid ≈ 2.6 : 1 — a loaded crescent mark, not a 10:1\n//     filament), so each reads as a SEPARABLE atom.\n//   • SPARSE placement with visible inter-dab ground (density ≈ 0.42), over a darkened\n//     swirl underpainting (the Starry-Night visible ground the bright dabs read over).\n//   • a single sampleBase per cell + OKLCh broken-color jitter (pigment shimmer);\n//     full-height impasto crown per dab (each catches its own raked glint).\n// This is ~2 layers × 9 cells × ~1 sampleBase ≈ 20/pixel (a ~16× cost cut, ~60fps), and\n// jittered short dabs CANNOT marble.\n\n// One atomic comma dab — a CLEAN analytic crescent SDF (NOT the oil curvedStroke).\n// The oil curvedStroke's spine-projection + rounded end-caps + fwidth-AA + bristle\n// edge each leave a thin sub-pixel skirt, and a field of those skirts over open ground\n// abuts into a faint network of SEAM LINES. A van-Gogh dab needs none of that\n// machinery — it is a short oriented capsule (rounded line segment) whose half-width\n// tapers along its length (the loaded-head → dragged-tail comma). The capsule SDF is\n// derivative-free (smoothstep AA against the analytic distance), so no fwidth skirt, no\n// best-of-9 seam. Returns a StrokeHit (coverage/color/alongT/edgeN) for paintOver.\n//   cellSize — dab cell pitch (UV units); lenMul/widMul — dab aspect off the cell;\n//   density  — placement gate (sparse → visible gaps); seed — layer decorrelation.\nStrokeHit vangoghDab(vec2 p, float cellSize, float lenMul, float widMul,\n                     float density, float jitterAmt, float angJitter,\n                     vec2 flow, float t, float seed) {\n  vec2 cell = floor(p / cellSize);\n  StrokeHit best = noHit();\n  for (int dy = -1; dy <= 1; dy++) {\n    for (int dx = -1; dx <= 1; dx++) {\n      vec2 cc = cell + vec2(float(dx), float(dy));\n      // Sparse density gate — leaves visible canvas ground between dabs (atomicity).\n      if (hash21(cc * 1.7 + seed * 0.3) > density) continue;\n      vec2 jit = (hash22(cc + seed) - 0.5) * jitterAmt;\n      vec2 center = (cc + 0.5 + jit) * cellSize;\n\n      // Orientation: the analytic swirl flow at the cell center + per-cell angular\n      // jitter (so adjacent dabs FAN — the divisionist scatter that keeps the dabs from\n      // locking into perfectly-concentric rows) + a flow-noise curl wobble.\n      float ang = (hash21(cc + seed + 11.0) - 0.5) * angJitter;\n      float curl = (fbm(center * 3.1 + seed * 1.9) - 0.5) * (0.18 + 0.3 * uFlowCurl);\n      vec2 dir = rotateDir(safeDir(flow), ang + curl);\n\n      // SHORT comma — len ≈ the cell pitch (a contained dab, not a tiling filament).\n      float len = cellSize * lenMul * (0.82 + 0.30 * hash21(cc + seed + 23.0));\n      float halfW = max(cellSize * widMul * (0.86 + 0.28 * hash21(cc + seed + 41.0)), 0.010);\n\n      // Local dab frame: along = projection onto the dab axis (centred 0), perp = cross.\n      vec2 rel = p - center;\n      float along = dot(rel, dir);\n      float perp  = dot(rel, vec2(-dir.y, dir.x));\n      // The comma TAPER: half-width peaks ~1/3 from the loaded head and drags to a fine\n      // tail (the asymmetric crescent). u in [0,1] along the dab from head to tail.\n      float halfLen = len * 0.5;\n      float u = clamp((along + halfLen) / max(len, 1e-5), 0.0, 1.0);\n      float head = smoothstep(0.0, 0.30, u);                 // rounded loaded head\n      float tail = pow(clamp((1.0 - u) / 0.70, 0.0, 1.0), 1.3); // dragged thin tail\n      float prof = clamp(min(head, tail) * 1.05, 0.0, 1.0);\n      float wNow = halfW * prof;\n      // A gentle comma CURVE — bow the spine sideways so the mark reads as a crescent.\n      float bowAmt = (hash21(cc + seed + 53.0) - 0.5) * halfW * 1.4;\n      float spine = bowAmt * sin(u * PI);                    // 0 at ends, max mid\n      float d = abs(perp - spine);\n      // Derivative-free AA: a fixed fraction of the dab width (no fwidth skirt). Coverage\n      // is 1 in the dab body, smoothly 0 at the edge; outside the [head,tail] span wNow→0\n      // so the dab self-terminates with no end-cap blob.\n      float aa = max(halfW * 0.45, 0.002);\n      float cov = (wNow > 1e-5) ? (1.0 - smoothstep(wNow - aa, wNow + aa, d)) : 0.0;\n      if (cov < 0.004) continue;\n\n      // One sampleBase per cell + OKLCh broken-color jitter (pigment patch shimmer). The\n      // dab pigment is chroma-boosted (saturate3 1.22) AND value-lifted (×1.14) at the\n      // SOURCE — the loaded brush mark deposits BRIGHTER + more saturated pigment than the\n      // base field, so the loaded-mark pixels lift BOTH the §4.1 colourfulness AND the field\n      // mean-luma. The higher mean raises the gap-fraction threshold (meanLum·0.6) so the\n      // moderate-depth inter-dab ground falls below it (the gap-fraction atomicity read is\n      // DECOUPLED from the ground floor depth — the ground can stay only moderately dark, so\n      // the §4.3 slope is NOT over-steepened, while the bright dabs manufacture the gap). The\n      // dark ground is untouched, so this lifts C + gap WITHOUT a band-conflicting deep ground.\n      vec3 colMid = brokenColorJitter(\n        saturate3(sampleBase(center, t), 1.22),\n        hash21(cc + seed + 89.0),\n        hash21(cc * 2.3 + seed + 97.0),\n        1.0\n      );\n\n      // edgeN: 1 at the spine, 0 at the edge (the impasto crown reads it).\n      float edgeN = (wNow > 1e-5) ? clamp(1.0 - d / wNow, 0.0, 1.0) : 0.0;\n      if (cov > best.coverage) {\n        best = StrokeHit(cov, colMid, u, (wNow > 1e-5 ? (perp - spine) / wNow : 0.0), edgeN);\n      }\n    }\n  }\n  return best;\n}\n\nvec3 mediumVangogh(vec3 col, vec2 p, float t) {\n  vec3 result = col;\n  float height = 0.0;\n\n  // The swirl flow basis (van-Gogh sky = swirl rows). Computed ONCE per pixel.\n  vec2 flow = flowField(p, t);\n\n  // Dab cell pitch in UV [0,1]. CHUNKY, SEPARABLE marks — the big layer is ~1/14 of\n  // the frame so it reads as a distinct loaded brush dab (not a fine speckle, not a\n  // screen-spanning filament). uStrokeScale (authored ~110) modulates the pitch\n  // INVERSELY (higher = finer dabs) around the ~0.07 chunky default; clamped so the\n  // dab never collapses to noise or balloons past one screen. uStrokeAnisotropy\n  // stretches the comma aspect (longer, more directional marks).\n  float pitch = clamp(7.0 / max(uStrokeScale, 20.0), 0.044, 0.12);\n  // The dab LENGTH ≈ the cell pitch (lenAniso ~1.2-1.7) so a dab sits WITHIN its cell\n  // and leaves visible ground around it — a longer dab (lenMul > 2) spans neighbour\n  // cells and tiles the field into continuous coverage (no gaps, the gap-fraction floor\n  // miss + the marble). Narrow width keeps it a directional crescent, not a wide shard.\n  // lenAniso/widAniso — the comma aspect. The W-AUR-VANGOGH-CONFIG reconcile lengthened\n  // the comma (1.25-1.75 → 1.5-2.0) and THINNED it (0.40-0.30 → 0.32-0.24): the rebuilt\n  // dabs measured §4.2 structure-tensor anisotropy A=0.699, BELOW the reference band\n  // floor 0.732 (the dabs read slightly too short-fat → locally isotropic, dragging mean\n  // A down). A LONGER, THINNER comma hugs the swirl-flow tangent more strongly so the\n  // rendered field reads coherent the way the tensor metric measures it (A lift), while\n  // the thinner width LEAVES MORE inter-dab ground (the gap-fraction atomicity read is\n  // preserved, not closed). The length stays ≈ the cell pitch so each dab remains a\n  // CONTAINED separable mark, not a tiling filament.\n  float lenAniso = mix(2.1, 2.7, uStrokeAnisotropy);   // comma length / cell pitch (a contained dab)\n  float widAniso = mix(0.23, 0.16, uStrokeAnisotropy); // comma width / cell pitch (a thin loaded crescent)\n\n  // Layer 1 — the big swirl-row dabs: the loaded gestural marks, SPARSE so the dark\n  // underpainting shows BETWEEN them (the atomicity read — visible inter-dab ground).\n  // hardness 0.74 — the loaded brush deposits with a SOFT shoulder, so adjacent dabs\n  // feather instead of razor-clipping into hard polygon shards. jitterAmt 0.85 fully\n  // scatters the cell centres so the placement GRID never reads as a periodic seam.\n  // angJitter EASED 0.55 → 0.38 (the reconcile A lift): a tighter angular spread keeps\n  // adjacent dabs more flow-aligned (coherent) while the per-cell curl wobble still\n  // breaks the perfectly-concentric-rows read (the dabs FAN gently, not lock to a grid).\n  // angJitter eased 0.38 → 0.26 + hardness 0.74 → 0.84 (the §4.2 A lift): a tighter angular\n  // spread keeps adjacent dabs strongly flow-aligned (the structure-tensor coherence) and a\n  // sharper edge gives each dab a STRONGER directional gradient (a soft dab's faint edge\n  // reads weakly anisotropic). The per-cell curl wobble still breaks the perfectly-\n  // concentric-rows read, so the dabs FAN gently rather than locking to a grid.\n  StrokeHit dBig = vangoghDab(p, pitch, lenAniso, widAniso,\n                              0.30, 0.50, 0.08, flow, t, 1.3);\n  // streakFreq raised 14 → 34 (the §4.3 slope fix): a high-frequency streak runs ALONG the\n  // dab spine, injecting fine-scale luma energy that flattens the too-steep β WITHOUT\n  // adding new isotropic marks (the along-spine streak stays directional, so the §4.2\n  // anisotropy the long-thin dab carries is preserved — the band-conflict resolution).\n  paintOver(result, height, dBig, 48.0, 0.34,\n            uImpasto * 0.75 * uStrokeAmount, 0.84, 1.3, 1.0 /* full crown */);\n\n  // Layer 2 — two-thirds-pitch accent dabs (the divisionist sparkle between the big\n  // marks — denser placement but still contained dabs, so dark canvas survives between).\n  // angJitter eased 0.5 → 0.34 (same A-lift reasoning — the accent dabs queue along the\n  // flow rather than scattering isotropically); hardness 0.78 → 0.86 sharpens the gradient.\n  StrokeHit dSml = vangoghDab(p + vec2(-5.1, 8.4), pitch * 0.62, lenAniso * 0.95,\n                              widAniso * 0.95, 0.36, 0.60, 0.20, flow, t, 4.1);\n  paintOver(result, height, dSml, 42.0, 0.22,\n            uImpasto * 0.55 * uStrokeAmount, 0.86, 4.1, 1.0);\n\n  // Layer 3 — the FINE divisionist flecks (W-AUR-VANGOGH-CONFIG reconcile). The two-layer\n  // field read as a few BIG cobalt dabs over a smooth empty teal ground (the side-by-side\n  // against starry-night-crop.png: the painting is DENSE small varied flecks, not sparse\n  // large marks). The flecks are LONG·THIN directional commas hugging the swirl tangent\n  // (the A-preserving fleck — short fat flecks read isotropic and drag §4.2 down). Sparse\n  // density 0.30 keeps the inter-fleck ground open (the gap-fraction atomicity read).\n  // The flecks are FINE (≈0.42× pitch) directional commas, tightly flow-aligned (angJitter\n  // 0.22) and SPARSE (0.30) — the divisionist sparkle between the loaded marks that keeps\n  // the inter-dab ground open (atomicity) while adding the small-mark mid-band energy. Thin\n  // + flow-aligned so the §4.2 anisotropy the long dabs carry survives the added small marks.\n  StrokeHit dFleck = vangoghDab(p + vec2(8.7, -3.2), pitch * 0.42, lenAniso * 1.0,\n                                widAniso * 0.68, 0.30, 0.85, 0.22, flow, t, 7.7);\n  paintOver(result, height, dFleck, 30.0, 0.16,\n            uImpasto * 0.40 * uStrokeAmount, 0.80, 7.7, 1.0);\n\n  // ── Flow-aligned GROUND swirl-GROOVING (W-AUR-VANGOGH-CONFIG reconcile) ──\n  // The bare ground between the dabs was a SMOOTH low-frequency teal gradient — the §4.3\n  // power-spectrum slope ran too STEEP (β≈−2.0: too much low-frequency nuclei-field energy\n  // vs the high band). Van Gogh's actual sky has the swirl BRUSHWORK everywhere, so groove\n  // the OPEN GROUND with a flow-aligned ripple. It is SUBTRACTIVE-biased (it only darkens,\n  // never brightens the ground), so it DEEPENS the inter-dab valleys (MORE gap-fraction —\n  // atomicity preserved/strengthened) while injecting the across-flow mid-band luma energy\n  // that flattens β. Gated to the open ground (low height) via the (1-height) gate so it\n  // never lifts/erodes the loaded dab bodies (the dab read is untouched). Across the swirl\n  // FLOW axis so the grooves run WITH the sky's brushwork direction.\n  // The ground swirl-combing — a FIXED-global-angle directional grain on the open ground.\n  // A flow-FOLLOWING comb varies its direction with the swirl across the readback window,\n  // so its blurred structure tensor averages toward isotropic (no §4.2 lift). A comb at one\n  // GLOBAL angle (van-Gogh's dominant diagonal stroke grain) is COHERENT over the whole\n  // window ⇒ every open-ground pixel reads the SAME direction (the §4.2 mean-anisotropy lift\n  // the swirl comb could not reach). Two resolved bands (≈22 + ≈44 across the global axis ⇒\n  // ~21px + ~11px at 464px, both in the §4.3 fit window) also feed the β flatten. Gated to\n  // the open ground (low height) so the loaded dabs are untouched; SUBTRACTIVE-biased so the\n  // inter-dab valleys deepen (gap-fraction preserved, not filled).\n  vec2 gdir = normalize(vec2(0.82, 0.57));   // the dominant global stroke grain axis\n  float gacross = dot(p, vec2(-gdir.y, gdir.x));\n  float groundGate = 1.0 - smoothstep(0.05, 0.28, height);\n  float groove = sin(gacross * 22.0) + 0.5 * sin(gacross * 44.0 + 1.3);\n  result *= 1.0 - 0.16 * groundGate * (0.5 - 0.5 * groove);\n\n  // Canvas tooth (fine linen weave — base relief the raking light catches).\n  float tooth1 = vnoise(p * 220.0);\n  float tooth2 = vnoise(p * 220.0 * vec2(0.6, 2.4) + 37.0);\n  float tooth  = (0.6 * tooth1 + 0.4 * tooth2) - 0.5;\n  result *= 1.0 + tooth * 0.06 * uCanvasGrain;\n\n  // Relight the accumulated impasto height (raking light catches each dab's crown —\n  // the visible impasto ridge that makes a dab read as a RAISED loaded mark).\n  float canvasBase = tooth * 0.06 * 0.5;\n  result = relightImpasto(result, height, canvasBase);\n\n  // The visible Starry-Night GROUND — the open canvas between dabs (low accumulated\n  // height) darkens toward the deep underpainting, so the bright impasto dabs read as\n  // SEPARABLE atomic marks over a darker field (the dab↔ground contrast the user's\n  // \"reads as strokes\" bar needs, AND the painterly-statistics gap-fraction floor — the\n  // dark inter-dab ground must register as a measurable fraction of the field). 0.30\n  // floor over a WIDER smoothstep (0..0.40): deep enough that the open gaps read as the\n  // dark Starry-Night sky between marks (clearing the 4% gap floor) while the gradual\n  // ramp keeps the gaps reading as flowing underpainting, not the hard black ink LINES a\n  // narrow near-black step (0.26 / 0.0..0.32) produced.\n  // The smoothstep FLOOR starts at 0.10 (not 0): a dab's faint AA coverage TAIL deposits\n  // a sliver of height, and where many tails align along the swirl flow those slivers\n  // accumulate into a faint continuous height ridge that read as thin curved LINES in the\n  // ground. Holding the ground full-dark until a real dab CORE (height > 0.10) lifts out\n  // keeps the open field clean — only the loaded dab bodies emerge, the tails stay ground.\n  result *= mix(0.115, 1.0, smoothstep(0.085, 0.31, height));\n\n  // The van-Gogh chroma — a light OKLCh punch. EASED to 0.92 (off 1.04) in the\n  // W-AUR-VANGOGH-CONFIG reconcile: the rebuilt discrete-dab register's per-dab OKLCh\n  // broken-color jitter (the ±16° hue + ±chroma swing per cell) lifts the field's\n  // colourfulness on its own, so the prior 1.04 punch on TOP overshot the §4.1\n  // reference band (rendered C=100.5 vs the painting's measured 70.67 ± 15/+25 →\n  // ceiling 95.67 — neon, not Starry Night's rich-but-not-garish blue/gold). The 0.92\n  // demote pulls the rendered colourfulness back IN-band while the broken-color jitter\n  // keeps the per-dab pigment variation (the discrete-stroke read is untouched — this\n  // is a chroma scale, not a dab-geometry change). Band-anchored on the real painting\n  // (proof:aurora-arresting), NOT a widened band.\n  result = saturate3(result, 1.0);\n  return result;\n}\n\n", Se = "StrokeProfile profileFor(int medium, int mode) {\n  // The oil baseline (mode 0) — balanced modern-abstract/palette-knife hybrid.\n  StrokeProfile prof = StrokeProfile(\n    0,      // shapeType — tapered\n    0.25,   // bristleAmp\n    9.0,    // streakFreq\n    0.09,   // streakAmp\n    0.9,    // impastoAmp\n    0.80,   // hardness\n    340.0,  // toothScale — W-AUR-VANGOGH-CONFIG reconcile: FINER 240 → 340. Oil's §4.3 β sat\n            //               right AT the −1.45 band ceiling (±0.02 frame variance straddling\n            //               it, surfaced once the vangogh bands went green and the spec reached\n            //               the oil β assert). MEASURED on the live GPU: ADDING fine high-\n            //               frequency tooth energy steepens oil's power-spectrum slope (more\n            //               negative — coarsening it the other way pushed β SHALLOWER to\n            //               −1.33). A finer, stronger tooth moves β to a safe distance inside\n            //               the −5/3 band [−1.85, −1.45]. Oil-only (oil-pastel sets 280).\n    0.14,   // toothAmp — raised 0.09 → 0.14 (same oil-β steepen: more fine-grain high band).\n    1.03,   // pigmentSat\n    0.65,   // densityBig\n    0.78,   // densityMed\n    0.90,   // densitySml\n    0.0,    // energyGrade — oil is uniform-length (no Starry-Night cascade)\n    0.4,    // impastoFloor — oil's 0.4+0.6·edgeN crown falloff\n    0.92,   // densityFill — oil fills DENSELY with the directional low-relief fill so the\n            //               flat-colour regions between the macro strokes are TEXTURED with\n            //               flow-aligned fine strokes (the §4.2 anisotropy needs every pixel\n            //               to carry directional gradient; a flat untextured region reads as\n            //               an isotropic/degenerate tensor and drags mean A down). The fill\n            //               is now directional + low-streak, so dense coverage textures\n            //               WITHOUT the prior round-dab speckle.\n    1.0     // groundFloor — oil keeps the full base (no ground darken; fills everywhere)\n  );\n\n  // W-AUR-VANGOGH-REBUILD — the van-Gogh medium NO LONGER routes through this oil\n  // StrokeProfile cascade. It is a DEDICATED atomic-dab body (mediumVangogh / vangoghDab\n  // below) — the prior profile-driven cascade (four dense layers of long-thin\n  // tensor-oriented strokes) merged into the \"marbled flow-bands\" read and ran at ~4fps.\n  // profileFor now serves ONLY the oil + oil-pastel stroke mediums; the MEDIUM_VANGOGH\n  // case is removed (no caller). MEDIUM_VANGOGH stays #define'd (the main() dispatch +\n  // isPainterlyStroke still key off uMedium==5).\n\n  // ── Oil-pastel — broad smeared directional strokes deposited via the brush engine\n  // (slice 8 F1). Creamy soft edges (low hardness), heavy build-up, and a chroma punch.\n  // Distinct from the dry-crayon tooth-multiply (mediumCrayon): oil-pastel LAYS strokes.\n  if (medium == MEDIUM_OILPASTEL) {\n    prof.shapeType   = 0;     // tapered directional smear — a round dab reads locally\n                              // isotropic (steep β / low A); the creamy oil-pastel mark\n                              // is a broad DIRECTIONAL smear, which carries the §4.2\n                              // coherence + adds the directional mid-band the §4.3 slope\n                              // needs (the round-dab profile rolled β off to −2.8)\n    prof.bristleAmp  = 0.10;  // smooth creamy edges (waxy pastel, not bristle)\n    prof.streakFreq  = 18.0;  // finer internal smear — adds the mid/high-band the smooth\n                              // round-dab profile lacked (the β-too-steep fix); the\n                              // directional along-spine streak flattens the §4.3 slope\n                              // without crossing the strokes\n    prof.streakAmp   = 0.20;\n    prof.impastoAmp  = 0.7;   // build-up, but flatter than oil impasto\n    prof.hardness    = 0.42;  // CREAMY — soft compositing, strokes blend on overlap\n    prof.toothScale  = 280.0;\n    prof.toothAmp    = 0.12;  // the pastel tooth reads stronger\n    prof.pigmentSat  = 0.80;  // the waxy chroma — eased off 1.16 so the rendered C clears\n                              // the §4.1 band ceiling with margin while the K-M subtractive\n                              // overlap path (paintOverOklab) keeps the chroma OFF the grey floor\n    prof.hardness    = 0.58;  // creamy but the strokes still REGISTER — a too-soft blend\n                              // (0.42) smeared the marks into the smooth colour zones and\n                              // the field rolled off too steep (β≈−2.9, a near-gradient);\n                              // a crisper deposit keeps the mid/high-band the −5/3 slope needs\n    prof.densityBig  = 0.62;  // broad directional coverage — but NOT over-dense: too many\n                              // overlapping creamy strokes point in slightly-varied\n                              // directions and the local §4.2 coherence DROPS (over-dense\n                              // measured A=0.57 < the 0.62-density A=0.63). The mid-density\n                              // flow-aligned deposit is the coherence sweet spot.\n    prof.densityMed  = 0.78;\n    prof.densitySml  = 0.90;\n    prof.energyGrade = 0.80;  // a strong Starry-Night length cascade — varies stroke length\n                              // by luma/coherence so the dab-size energy spreads across\n                              // scales (the −5/3 cascade the uniform-length 0.0 grade lacked,\n                              // the §4.3 β-too-steep fix). Not the full van-Gogh 1.0 (oil-\n                              // pastel keeps its creamy broad character).\n    prof.impastoFloor = 0.55; // a softer crown than van-Gogh's full-height glint\n    return prof;\n  }\n\n  // ── Oil sub-modes (mode dispatch on the oil medium).\n  if (mode == 1) {           // palette knife — razor edges, heavy impasto\n    prof.shapeType  = 3;     // flat, even\n    prof.bristleAmp = 0.12;\n    prof.streakFreq = 4.0;  prof.streakAmp = 0.05;\n    prof.impastoAmp = 1.6;\n    prof.hardness   = 0.95;\n    prof.toothAmp   = 0.04;\n    prof.densityBig = 0.80; prof.densityMed = 0.88; prof.densitySml = 0.70;\n  } else if (mode == 3) {    // thick brushwork — heavy bristle brush\n    prof.shapeType  = 0;     // tapered\n    prof.bristleAmp = 0.32;\n    prof.streakFreq = 14.0; prof.streakAmp = 0.14;\n    prof.impastoAmp = 1.2;\n    prof.hardness   = 0.85;\n    prof.toothAmp   = 0.07;\n  }\n  return prof;\n}\n", Ce = "// ── Medium overlays ───────────────────────────────────────────────────────\n\n// A quick re-computation of base color for edge-mask sampling\nvec3 sampleBase(vec2 p, float t) {\n  vec2 pw = domainWarp(p, t);\n  float id; float vm;\n  nucleiField(pw, t, id, vm);\n  vec3 c = samplePalette(id);\n  c *= 1.0 + uValueVariance * vm;\n  return c;\n}\n\n// ── Structure-tensor / edge-tangent-flow (AW.W4.1) ─────────────────────────\n// The keystone of the painterly engine: derive stroke orientation from the COLOR\n// FIELD's OWN structure, not a hand-authored flow pattern. A 3x3 Sobel over\n// luma(sampleBase) yields the gradient (Gx,Gy); the 2x2 structure tensor\n// J = [[Gx², GxGy],[GxGy, Gy²]] eigen-decomposes closed-form into the MINOR\n// eigenvector (the edge-TANGENT — least color change, the stroke direction; the\n// MAJOR eigenvector is the gradient/edge-NORMAL and would make strokes cross the\n// bands) and the coherence A = (λ1-λ2)/(λ1+λ2). Low-coherence (flat) zones relax\n// toward the smooth fallbackDir by (1-A) so tensor noise never reads as jitter.\n// Single-pass small-tap WebGL2 form; the Gaussian-smoothed multi-tap + LIC smear\n// is the AW.W7 WebGPU multi-pass fold. Returns vec3(tangent.x, tangent.y, A).\nvec3 structureTensorField(vec2 p, float t, vec2 fallbackDir) {\n  float e = 0.0035; // small-tap neighborhood (edge-mask scale)\n  // 3x3 luma samples.\n  float l00 = dot(sampleBase(p + vec2(-e, -e), t), W_LUMA);\n  float l10 = dot(sampleBase(p + vec2( 0.0, -e), t), W_LUMA);\n  float l20 = dot(sampleBase(p + vec2( e, -e), t), W_LUMA);\n  float l01 = dot(sampleBase(p + vec2(-e, 0.0), t), W_LUMA);\n  float l21 = dot(sampleBase(p + vec2( e, 0.0), t), W_LUMA);\n  float l02 = dot(sampleBase(p + vec2(-e, e), t), W_LUMA);\n  float l12 = dot(sampleBase(p + vec2( 0.0, e), t), W_LUMA);\n  float l22 = dot(sampleBase(p + vec2( e, e), t), W_LUMA);\n  // Sobel derivatives.\n  float Gx = (l20 + 2.0 * l21 + l22) - (l00 + 2.0 * l01 + l02);\n  float Gy = (l02 + 2.0 * l12 + l22) - (l00 + 2.0 * l10 + l20);\n  // Structure tensor components.\n  float Jxx = Gx * Gx;\n  float Jyy = Gy * Gy;\n  float Jxy = Gx * Gy;\n  // Closed-form 2x2 eigenvalues: λ = 0.5(tr ± sqrt((Jxx-Jyy)² + 4Jxy²)).\n  float tr = Jxx + Jyy;\n  float disc = sqrt(max((Jxx - Jyy) * (Jxx - Jyy) + 4.0 * Jxy * Jxy, 0.0));\n  float lambdaMaj = 0.5 * (tr + disc);\n  float lambdaMin = 0.5 * (tr - disc);\n  // The MAJOR eigenvector points along the gradient (the edge NORMAL); the MINOR\n  // is its perpendicular — the edge TANGENT (the stroke direction). The principal\n  // angle is θ = 0.5·atan2(2·Jxy, Jxx-Jyy); this form is robust at Jxy≈0 (where the\n  // (λ-Jyy, Jxy) eigenvector formula collapses to a zero vector). The tangent is at\n  // θ + 90°. Guard the fully-isotropic (no-gradient) degenerate case → fallback dir.\n  vec2 tangent;\n  if (disc < 1e-7) {\n    tangent = fallbackDir;\n  } else {\n    float theta = 0.5 * atan(2.0 * Jxy, Jxx - Jyy); // major (gradient) angle\n    tangent = vec2(-sin(theta), cos(theta));          // perpendicular = tangent\n  }\n  // Coherence A in [0,1] — high where one eigenvalue dominates (a strong edge),\n  // ~0 in flat zones. Relax toward the smooth fallback by (1 - A).\n  float A = (lambdaMaj - lambdaMin) / (lambdaMaj + lambdaMin + 1e-6);\n  A = clamp(A, 0.0, 1.0);\n  // Coherence-weighted blend: full tensor in coherent zones, smooth in flat ones.\n  // Resolve the 180° sign ambiguity toward the fallback so the blend is continuous.\n  if (dot(tangent, fallbackDir) < 0.0) tangent = -tangent;\n  // The blend weight is a low-power curve (A^0.35) so even MODERATE-coherence zones queue\n  // strongly along the tangent (a flowing van-Gogh band reads coherent the way the\n  // structure-tensor metric measures it; a linear A under-orients the mid-band and the\n  // rendered field reads too isotropic — the §4.2 anisotropy lift). The flat-zone (A≈0)\n  // relax to the fallback survives (0^0.35=0), so tensor noise never reads as jitter.\n  float blendW = pow(A, 0.28);\n  vec2 dir = normalize(mix(fallbackDir, tangent, blendW) + 1e-6);\n  return vec3(dir, A);\n}\n\nvec3 mediumWatercolor(vec3 col, vec2 p, float t) {\n  // Wet-edge cauliflowers via luma-gradient magnitude\n  float eps = 0.004;\n  vec3 cx1 = sampleBase(p + vec2(eps, 0.0), t);\n  vec3 cx2 = sampleBase(p - vec2(eps, 0.0), t);\n  vec3 cy1 = sampleBase(p + vec2(0.0, eps), t);\n  vec3 cy2 = sampleBase(p - vec2(0.0, eps), t);\n  float gx = dot(cx1 - cx2, W_LUMA);\n  float gy = dot(cy1 - cy2, W_LUMA);\n  float edge = sqrt(gx * gx + gy * gy) / (2.0 * eps);\n  float mask = smoothstep(0.0, 2.5, edge);\n  col *= mix(1.0, 0.78, mask * uWetEdge);\n\n  // Granulation — pigment settles in paper tooth\n  float paper = 0.5 * vnoise(p * 160.0) + 0.5 * vnoise(p * 360.0);\n  float pigLoad = 1.0 - dot(col, W_LUMA);\n  col *= 1.0 - uGranulation * pigLoad * (paper - 0.5);\n\n  // Wash banding — faint horizontal wet gradient\n  float band = fbm(vec2(p.x * 1.5, p.y * 0.4));\n  col *= 1.0 + 0.04 * (band - 0.5);\n  return col;\n}\n\nvec3 mediumPastel(vec3 col, vec2 p, float t) {\n  vec2 flow = flowField(p, t);\n  vec2 perp = vec2(-flow.y, flow.x);\n  float along  = dot(p, flow) * uStrokeScale;\n  float across = dot(p, perp) * uStrokeScale;\n  across += 0.03 * (vnoise(p * 260.0) - 0.5);\n\n  float aniso = mix(1.0, 0.18, uStrokeAnisotropy);\n  float stroke = fbm(vec2(along * aniso, across));\n  col *= mix(1.0, 0.82 + 0.32 * stroke, uStrokeAmount);\n\n  // Pastel tooth — tiny high-frequency grain\n  float tooth = vnoise(p * 800.0);\n  col *= 1.0 - 0.08 * uStrokeAmount * (tooth - 0.5);\n  return col;\n}\n", we = "// ── Medium id constants (mirror uniformBridge.ts MEDIUM_ID — the uMedium ladder) ──\n// The StrokeProfile selector (AX.W12 + AX.W13) dispatches on the medium id: oil,\n// van-Gogh, and oil-pastel each AUTHOR a first-class profile in profileFor.\n#define MEDIUM_OIL 3\n#define MEDIUM_CRAYON 4\n#define MEDIUM_VANGOGH 5\n#define MEDIUM_OILPASTEL 6\n\n// ── Crayon — DRY wax pigment on paper tooth (first-class medium, uMedium==4) ──\n// Crayon is not strokes. It's dry pigment crumbs dragged across paper tooth — the\n// DRY tooth-multiply model: anisotropic tooth noise multiplied into the base, a hard\n// SCUMBLE that lets the lower color show through the broken upper layer, and stable\n// OKLCh broken-color pigment patches. NO sheen, NO burnish film — that waxy gloss is\n// the OIL-PASTEL deposition's signature (mediumOilPastel), distinct from dry crayon\n// (AX.W13 slice 8 F1 — the split). Crayon shares the SUBSTRATE (the structure-tensor\n// orientation, the tooth noise, OKLCh brokenColorJitter), not the dispatch body.\nvec3 mediumCrayon(vec3 col, vec2 p, float t) {\n  // Orient the deposition along the structure-tensor edge-tangent (the bridge forces\n  // tensor for crayon, so read the tensor field directly).\n  vec2 flow = structureTensorField(p, t, flowField(p, t)).xy;\n  float ang = atan(flow.y, flow.x);\n  float ca = cos(-ang), sa = sin(-ang);\n  vec2 pr = vec2(p.x * ca - p.y * sa, p.x * sa + p.y * ca);\n\n  float aniso = mix(0.45, 0.95, uStrokeAnisotropy);\n  float scale = max(uStrokeScale * 1.6, 180.0);\n\n  // ── Paper TOOTH height field — the relief the dry pigment rides. Anisotropic\n  // (squished along flow, coarse cross-flow). Peaks = ridges, valleys = paper pits.\n  float t1 = vnoise(vec2(pr.x * scale * aniso, pr.y * scale));\n  float t2 = vnoise(vec2(pr.x * scale * aniso * 0.4, pr.y * scale * 0.4) + 11.0);\n  float t3 = vnoise(vec2(pr.x * scale * aniso * 2.1, pr.y * scale * 2.1) + 23.0);\n  float tooth = 0.55 * t1 + 0.30 * t2 + 0.15 * t3;\n\n  // ── DRY TOOTH-MULTIPLY — the crayon multiplies the tooth relief into the base\n  // color (dark in the pits, bright on the ridges). The pressure (uStrokeAmount) sets\n  // how hard the multiply bites. This is the DRY model — no pigment build-up, no sheen.\n  float pressure = clamp(uStrokeAmount, 0.0, 1.0);\n  float bite = mix(0.18, 0.46, pressure);        // multiply strength\n  vec3 result = col * (1.0 - bite * (0.5 - tooth));\n\n  // ── HARD SCUMBLE — a broken upper layer that lets the lower color show through the\n  // gaps. A coarse mask gates a HARD (high-contrast) dragged pass so the dry crayon\n  // reads as broken pigment crumbs, not a smooth wash. Drag follows the tensor flow.\n  float scumbleMask = smoothstep(0.42, 0.78, vnoise(vec2(pr.x * scale * 0.55, pr.y * scale * 0.22) + 7.0));\n  float scumbleCoverage = scumbleMask * pressure * 0.7;     // hard broken upper layer\n  vec3 paper = col * mix(0.80, 0.94, pressure);             // the paper showing through\n  result = mix(paper, result, scumbleCoverage);\n\n  // ── Broken-color pigment: stable wax/pigment patches in OKLCh (W5), not flicker.\n  vec2 pigmentCell = floor(pr * max(scale * 0.18, 32.0));\n  float pigmentMask = smoothstep(0.28, 0.82, vnoise(pr * scale * 0.21 + 19.0));\n  result = brokenColorJitter(\n    result,\n    hash21(pigmentCell + 17.0),\n    hash21(pigmentCell * 2.1 + 31.0),\n    0.45 + 0.55 * pigmentMask\n  );\n\n  // Crayon is saturation-amplified (the wax chroma punch — but DRY, no sheen).\n  result = saturate3(result, 1.08);\n\n  return result;\n}\n\n// ── StrokeProfile — the oil-stroke parameter vector (AX.W12, slice 8 F6) ──────\n// The SBR \"stroke = parameter vector\" canon (facet 11): the per-mode if-ladder knobs\n// become struct FIELDS — logic-as-DATA, not an imperative branch. profileFor(medium,\n// mode) populates the profile for a (medium, mode) pair; paintStrokeLayers(profile) runs\n// the four-layer bestOil/paintOver cascade off it. A new medium (the W13 van-Gogh /\n// oil-pastel profiles) AUTHORS a profile entry — it never edits a monolith.\nstruct StrokeProfile {\n  int   shapeType;   // 0 tapered, 1 load-drag, 2 dab, 3 even, 4 comma/crescent\n  float bristleAmp;  // edge raggedness 0..0.5\n  float streakFreq;  // internal-streak spatial frequency\n  float streakAmp;   // internal-streak amplitude\n  float impastoAmp;  // paint-thickness contribution (the relight reads the height)\n  float hardness;    // edge-compositing crispness 0..1\n  float toothScale;  // canvas-tooth (linen weave) spatial frequency\n  float toothAmp;    // canvas-tooth amplitude\n  float pigmentSat;  // pigment saturation boost (OKLCh chroma scale)\n  float densityBig;  // layer-1 (big gestural) placement density gate\n  float densityMed;  // layer-2 (medium body) placement density gate\n  float densitySml;  // layer-3 (small dabs) placement density gate\n  float energyGrade; // AX.W13 — SBR energy-grade magnitude (0 off; 1 the van-Gogh\n                     // Starry-Night length cascade). The energy grade is a PROFILE\n                     // field, NOT a buried uMedium==5 bestOil branch (slice 8 F0).\n  float impastoFloor;// AX.W13 — the per-stroke height-crown floor (0.4 oil falloff;\n                     // 1.0 van-Gogh FULL-height crown so each dab catches its glint).\n  float densityFill; // AX.W13 — layer-4 fill-dab density. Oil/oil-pastel fill the\n                     // bald spots (0.95, full coverage); van-Gogh sets it LOW so the\n                     // sparse atomic dabs keep their visible inter-stroke ground gaps\n                     // (the atomicity read — a 0.95 fill would close every gap to a smear).\n  float groundFloor; // AX.W13 — the bare-ground darken floor. The gaps between strokes\n                     // (low accumulated height) multiply toward this floor, so van-Gogh's\n                     // sparse atomic dabs read over a visibly darker underpainting (the\n                     // Starry-Night visible ground). Oil/oil-pastel keep 1.0 (no darken).\n};\n\n// The (medium, mode) -> StrokeProfile selector. The if-ladder's knobs are the profile's\n// fields; a new medium adds a CASE here (never an edit to a dispatch body). mode is the\n// uStrokeMode oil-stroke sub-mode (0 oil, 1 palette-knife, 3 brushwork). crayon\n// (uMedium==4) is a PEER tooth-multiply medium that never reaches the stroke cascade.\n// The painterly STROKE mediums are oil (MEDIUM_OIL), van-Gogh (MEDIUM_VANGOGH) and\n// oil-pastel (MEDIUM_OILPASTEL) — each AUTHORS its own profile entry below.\n" + Se + "\n// The single parameterized four-layer stroke cascade (AX.W12). The four hand-unrolled\n// bestOil/paintOver invocations collapse into ONE body driven by the profile + the\n// uniform-derived per-layer scale/anisotropy multipliers (which are mode-INVARIANT — the\n// per-layer offsets, seeds, and len/wid muls are fixed across modes, so they stay here as\n// the cascade's structure; only the profile's knobs differentiate the medium). The\n// mode passes through for the two mode-special-cased shapes (knife layer-3 dabs +\n// layer-4 fill shape), which are placement details of the cascade, not profile knobs.\nvoid paintStrokeLayers(inout vec3 col, inout float height, StrokeProfile prof,\n                       int mode, vec2 p, float t) {\n  // Scales & multipliers from uniforms (mode-invariant cascade structure).\n  float baseScale = max(uStrokeScale * 0.006, 0.008);\n  // Three primary layers: big gestural, medium body, small dabs.\n  // W-AUR-STUDIO D5 — the −5/3 radii respacing was MEASURED + REVERTED (recorded, not\n  // silently struck). The φ-adjacent geometric candidate (2.4 / 1.45 / 0.87, a fixed\n  // ≈1.66× step) moved oil-pastel β −2.534 → −2.413 — TOWARD the −5/3 band but NOT into\n  // it (floor −1.85), so the keep-iff-into-band gate failed. The single-pass WebGL2 path\n  // is the true limiter (HC-aurora §5): the oil-pastel β + the oil/oil-pastel anisotropy\n  // residual stays the oil/oil-pastel stroke-cascade ceiling. BB.W-AUR-KUWAHARA DECIDED the\n  // booking: the SOFT anisotropic-Kuwahara finish ships as the OPT-IN mediumKuwahara() body\n  // (uMedium==7) the consumer selects — NOT a re-tune of these oil/oil-pastel knobs. The\n  // hand-set 2.4 / 1.1 / 0.45 stays — it is van-Gogh's landed-band spacing.\n  float sBig = baseScale * 2.4;\n  float sMed = baseScale * 1.1;\n  float sSml = baseScale * 0.45;\n\n  // Strokes are elongated HARD (long·thin) so the rendered field reads strongly\n  // directional — the §4.2 structure-tensor anisotropy lift. A short fat dab is locally\n  // isotropic and drags the measured A down; a long thin stroke that hugs the tensor\n  // tangent is what gives Starry Night its A≈0.83. The width muls are pulled tighter to\n  // keep the long·thin aspect.\n  float lenMulBig = mix(3.0, 5.2, uStrokeAnisotropy);\n  float widMulBig = mix(0.42, 0.22, uStrokeAnisotropy);\n  float lenMulMed = mix(2.8, 4.8, uStrokeAnisotropy);\n  float widMulMed = mix(0.40, 0.21, uStrokeAnisotropy);\n  float lenMulSml = mix(2.2, 3.6, uStrokeAnisotropy);\n  float widMulSml = mix(0.36, 0.23, uStrokeAnisotropy);\n\n  float jitterAmt = 0.75;   // large jitter — no grid\n  vec2 flow = flowField(p, t);\n\n  // Layer 1 — big gestural strokes (sparse, shaping)\n  StrokeHit hBig = bestOil(p, sBig, lenMulBig, widMulBig, jitterAmt * 0.55,\n                           prof.densityBig, prof.shapeType, prof.bristleAmp, flow, t, 1.3,\n                           prof.energyGrade);\n  paintOver(col, height, hBig, prof.streakFreq * 0.7, prof.streakAmp,\n            uImpasto * prof.impastoAmp * uStrokeAmount, prof.hardness, 1.3, prof.impastoFloor);\n\n  // Layer 2 — medium body strokes\n  StrokeHit hMed = bestOil(p + vec2(11.3, 3.7), sMed, lenMulMed, widMulMed,\n                           jitterAmt, prof.densityMed, prof.shapeType, prof.bristleAmp, flow, t, 2.7,\n                           prof.energyGrade);\n  paintOver(col, height, hMed, prof.streakFreq, prof.streakAmp,\n            uImpasto * prof.impastoAmp * uStrokeAmount, prof.hardness, 2.7, prof.impastoFloor);\n\n  // Layer 3 — small strokes (more frequent, smaller). Even the knife layer-3 stays the\n  // medium's DIRECTIONAL shape, not a round dab — the round-dab sparkle was a fine-scale\n  // isotropic speckle that dragged the §4.2 anisotropy down + rolled the §4.3 slope shallow.\n  int smlShape = prof.shapeType;\n  StrokeHit hSml = bestOil(p + vec2(-5.1, 8.4), sSml, lenMulSml, widMulSml,\n                           jitterAmt * 0.9, prof.densitySml, smlShape,\n                           prof.bristleAmp * 0.85, flow, t, 4.1, prof.energyGrade);\n  paintOver(col, height, hSml, prof.streakFreq * 1.4, prof.streakAmp * 0.6,\n            uImpasto * prof.impastoAmp * 0.45 * uStrokeAmount, prof.hardness, 4.1, prof.impastoFloor);\n\n  // Layer 4 — fill strokes (dense, fine) — cover bald spots. The density is\n  // PROFILE-driven (prof.densityFill): oil/oil-pastel fill to full coverage; van-Gogh\n  // keeps it sparse so the bare-ground gaps between its atomic dabs survive. The fill is\n  // ELONGATED + flow-oriented (NOT a round dab) — a dense round-dab fill is locally\n  // isotropic and was the dominant §4.2 anisotropy sink; an oriented fine stroke keeps the\n  // fill directional so the whole field reads coherent.\n  float sFill = baseScale * 0.22;\n  float lenMulFill = mix(2.0, 3.2, uStrokeAnisotropy);\n  float widMulFill = mix(0.34, 0.24, uStrokeAnisotropy);\n  int fillShape = (mode == 1) ? 3 : prof.shapeType; // knife=even; others inherit the medium's directional shape\n  StrokeHit hFill = bestOil(p + vec2(3.9, -6.2), sFill, lenMulFill, widMulFill,\n                            jitterAmt * 0.8, prof.densityFill, fillShape,\n                            prof.bristleAmp * 0.6, flow, t, 8.9, prof.energyGrade);\n  // The fill is a low-relief, low-streak base coat — it covers bald spots WITHOUT adding\n  // fine bright speckle (a high-impasto/high-streak fill caught a per-dab glint that read\n  // as isotropic noise over the directional macro strokes — the §4.2/§4.3 sink). Its\n  // positional jitter is tamed too so the fill strokes stay flow-aligned, not scattered.\n  paintOver(col, height, hFill, prof.streakFreq * 1.2, prof.streakAmp * 0.35,\n            uImpasto * prof.impastoAmp * 0.18 * uStrokeAmount, prof.hardness * 0.9, 8.9, prof.impastoFloor);\n\n  // Optional crosshatch layer\n  if (uStrokeLayers == 2) {\n    vec2 flow2 = vec2(-flow.y, flow.x);\n    StrokeHit hX = bestOil(p + vec2(7.3, -2.1), sMed, lenMulMed * 0.9, widMulMed,\n                           jitterAmt, prof.densityMed * 0.7, prof.shapeType, prof.bristleAmp, flow2, t, 6.5,\n                           prof.energyGrade);\n    paintOver(col, height, hX, prof.streakFreq, prof.streakAmp * 0.85,\n              uImpasto * prof.impastoAmp * 0.55 * uStrokeAmount, prof.hardness, 6.5, prof.impastoFloor);\n  }\n}\n\n// The shared stroke-medium SUBSTRATE finish (AX.W12 + AX.W13): paint the four-layer\n// cascade off the profile, then canvas-tooth + relight + saturation. The three painterly\n// STROKE mediums (oil / van-Gogh / oil-pastel) each AUTHOR their own profile and call\n// this — the substrate is shared, the PROFILE (not a dispatch-body fork) differentiates\n// the medium (slice 8 F6). mode passes through for the cascade's mode-special-cased shapes.\nvec3 paintStrokeMedium(vec3 col, vec2 p, float t, StrokeProfile prof, int mode) {\n  vec3 result = col;\n  // AW.W4.2 — accumulated paint HEIGHT across the stroke layers. The relight reads\n  // its gradient for the normal; the canvas tooth seeds the base relief.\n  float height = 0.0;\n  paintStrokeLayers(result, height, prof, mode, p, t);\n\n  // Canvas tooth — linen weave\n  float tooth1 = vnoise(p * prof.toothScale);\n  float tooth2 = vnoise(p * prof.toothScale * vec2(0.6, 2.4) + 37.0);\n  float tooth  = (0.6 * tooth1 + 0.4 * tooth2) - 0.5;\n  result *= 1.0 + tooth * prof.toothAmp * uCanvasGrain;\n\n  // AW.W4.2 — relight the accumulated paint height with the movable uLightDir\n  // source (diffuse + Blinn specular, in LINEAR before aces()). The canvas tooth\n  // is the base relief term so the weave also catches the raking light.\n  float canvasBase = tooth * prof.toothAmp * 0.5;\n  result = relightImpasto(result, height, canvasBase);\n\n  // AX.W13 — the visible ground. The bare-ground gaps (low accumulated stroke height)\n  // multiply toward prof.groundFloor, so van-Gogh's sparse atomic dabs read as separable\n  // marks over a darker underpainting (oil/oil-pastel keep groundFloor=1.0 → no-op). The\n  // smoothstep keeps the impasto dabs (high height) at full value; only the open ground\n  // between them darkens.\n  result *= mix(prof.groundFloor, 1.0, smoothstep(0.0, 0.38, height));\n\n  // Pigment saturation boost\n  result = saturate3(result, prof.pigmentSat);\n\n  return result;\n}\n\nvec3 mediumOil(vec3 col, vec2 p, float t) {\n  // The oil medium is a thin body: fetch the oil profile (mode-dispatched), paint via\n  // the shared stroke substrate.  uStrokeMode: 0 oil (gestural), 1 knife, 3 brushwork.\n  int mode = uStrokeMode;\n  StrokeProfile prof = profileFor(MEDIUM_OIL, mode);\n  return paintStrokeMedium(col, p, t, prof, mode);\n}\n\n#define MEDIUM_KUWAHARA 7\n\n// ── Kuwahara — the anisotropic generalized Kuwahara painterly finish (uMedium==7,\n// BB.W-AUR-KUWAHARA) ───────────────────────────────────────────────────────────────\n// The SOTA edge-preserving painterly smoothing: the generalized/anisotropic Kuwahara\n// of Kyprianidis 2010, the SOFT polynomial-weighted variant (NOT the pre-2010 hard\n// argmin, which BANDS the flat field into an 8-spoke pinwheel — the §4.2 anti-regression).\n// Aurora is a PROCEDURAL field (no input texture), so the operator runs over the\n// sampleBase() color field directly: an elliptical kernel oriented along the structure\n// tensor (elongated along the edge-tangent, narrowed across by the anisotropy), divided\n// into 8 overlapping sectors via smooth weights; each sector accumulates a weighted mean\n// + variance, and the output is the variance-weighted blend of the sector means\n// (1/(1+var^q) — the low-variance sectors dominate, the SOFT criterion). The result is\n// flat painterly oil-paint patches with crisp colour-zone boundaries — the §4.3 −5/3\n// cascade carried by the multi-scale patch structure, the §4.2 anisotropy carried by the\n// ellipse hugging the tensor. Single-pass small-radius WebGL2 form: 4 rings × 8 angular\n// taps = 32 procedural neighbourhood samples (cost-bounded; the offscreen-park / PRM\n// freeze the substrate owns is untouched — a parked frame attaches zero of these taps).\nvec3 mediumKuwahara(vec3 col, vec2 p, float t) {\n  // Orient the kernel along the structure-tensor edge-tangent + read the coherence A.\n  vec3 stf = structureTensorField(p, t, flowField(p, t));\n  vec2 tangent = stf.xy;\n  float A = stf.z;\n  // The ellipse axes: major along the tangent, minor across. Anisotropy squeezes the\n  // minor axis (coherent zones → long-thin patch along the colour band; flat zones →\n  // near-circular so the operator degrades gracefully to the isotropic Kuwahara). The\n  // base radius rides uStrokeScale so the painterly patch size tracks the studio knob.\n  float radius = clamp(0.010 + uStrokeScale * 0.00018, 0.006, 0.024);\n  float aniso = mix(1.0, 0.34, clamp(A, 0.0, 1.0)); // minor/major axis ratio\n  vec2 ax = tangent;                  // major-axis unit vector (edge-tangent)\n  vec2 ay = vec2(-tangent.y, tangent.x); // minor-axis unit vector (edge-normal)\n\n  // 8 sectors, soft polynomial weights. Each sector i centers on angle i·45°; a sample's\n  // sector weight is a smooth gaussian-of-angular-distance so sectors OVERLAP (no hard\n  // argmin → no pinwheel banding). The center sample seeds all sectors equally.\n  const int SECTORS = 8;\n  vec3  m[8];   // per-sector weighted colour sum\n  vec3  s[8];   // per-sector weighted colour² sum (for variance)\n  float w[8];   // per-sector weight sum\n  for (int i = 0; i < SECTORS; i++) { m[i] = vec3(0.0); s[i] = vec3(0.0); w[i] = 0.0; }\n\n  // 4 rings × 8 angular taps over the elliptical disc. Each tap: sample the procedural\n  // field, then deposit into the sectors by the soft angular weight × the radial falloff.\n  const int RINGS = 4;\n  const int TAPS  = 8;\n  for (int r = 1; r <= RINGS; r++) {\n    float rr = float(r) / float(RINGS);          // 0.25 .. 1.0\n    float radial = 1.0 - 0.6 * rr;               // inner taps weigh more (gaussian-ish)\n    for (int a = 0; a < TAPS; a++) {\n      float ang = 6.2831853 * float(a) / float(TAPS);\n      // Elliptical offset: major axis full radius, minor axis squeezed by aniso.\n      vec2 dirE = ax * cos(ang) + ay * (sin(ang) * aniso);\n      vec2 off  = dirE * (radius * rr);\n      vec3 c = sampleBase(p + off, t);\n      // The sample's angular bin (its sector) — soft, so it spills into neighbours.\n      for (int k = 0; k < SECTORS; k++) {\n        float sectorAng = 6.2831853 * float(k) / float(SECTORS);\n        // Smallest signed angular distance to the sector center.\n        float d = ang - sectorAng;\n        d = atan(sin(d), cos(d));                 // wrap to [-π, π]\n        float sw = exp(-d * d * 2.2) * radial;    // gaussian angular weight × radial\n        m[k] += c * sw;\n        s[k] += c * c * sw;\n        w[k] += sw;\n      }\n    }\n  }\n\n  // Per-sector mean + scalar luma-variance, then the SOFT variance-weighted blend.\n  vec3  accum = vec3(0.0);\n  float accumW = 0.0;\n  for (int k = 0; k < SECTORS; k++) {\n    float wk = max(w[k], 1e-4);\n    vec3  mean = m[k] / wk;\n    vec3  var3 = max(s[k] / wk - mean * mean, vec3(0.0));\n    float variance = dot(var3, W_LUMA);           // perceptual scalar variance\n    // SOFT criterion: low-variance sectors (flat, inside a colour zone) dominate; a\n    // high-variance sector (straddling an edge) is suppressed. q controls the sharpness;\n    // a finite q (not argmin) keeps the blend continuous → no pinwheel spoke.\n    float sectorW = 1.0 / (1.0 + pow(variance * 320.0, 4.0));\n    accum  += mean * sectorW;\n    accumW += sectorW;\n  }\n  vec3 result = accum / max(accumW, 1e-4);\n\n  // The center sample anchors the result in near-flat zones (where every sector reads the\n  // same colour) so the operator is a no-op there — the smoothing only BITES at edges.\n  vec3 center = sampleBase(p, t);\n  result = mix(center, result, clamp(0.35 + 0.65 * A, 0.0, 1.0));\n\n  // A faint canvas tooth + pigment-saturation lift so the flattened patches read as\n  // oil-paint, not a posterized blur (the painterly identity the smoothing alone lacks).\n  float tooth = vnoise(p * 300.0) - 0.5;\n  result *= 1.0 + tooth * 0.06 * uCanvasGrain;\n  result = saturate3(result, 1.04);\n\n  // Blend toward the incoming col by uStrokeAmount so the studio amount knob controls\n  // the finish strength (uStrokeAmount=0 → the raw field passes through; default → full).\n  return mix(col, result, clamp(uStrokeAmount, 0.0, 1.0));\n}\n\n" + xe + "// AX.W13 — the oil-pastel medium (uMedium==6) is a DISTINCT stroke-deposition body,\n// split out of the dry-crayon tooth-multiply (mediumCrayon) (slice 8 F1). It DEPOSITS\n// broad smeared directional strokes via the same brush engine with a creamy soft\n// hardness, heavy pigment build-up where strokes overlap, and a chroma punch — the\n// stroke-deposition model that gives oil-pastel its depth and less-uniform read.\n// Shares the SUBSTRATE (the stroke cascade + tooth + relight) with oil/van-Gogh, NOT\n// the dispatch body of dry crayon. The bridge forces tensor orientation.\nvec3 mediumOilPastel(vec3 col, vec2 p, float t) {\n  StrokeProfile prof = profileFor(MEDIUM_OILPASTEL, 0);\n  return paintStrokeMedium(col, p, t, prof, 0);\n}", k = "\n", Te = `#version 300 es
 precision highp float;
 
 in vec2 vUv;
@@ -42,7 +44,7 @@ uniform float uValueVariance;
 uniform float uWarpAmount;
 uniform float uWarpScale;
 uniform float uWarpDrift;
-uniform int   uWarpMode;      // 0=fbm 1=cellular 2=hybrid
+uniform int   uWarpMode;      // 0=fbm 1=cellular 2=hybrid 3=curl (BB.B1 — opt-in Bridson flow warp)
 uniform int   uNoiseOctaves;
 
 // Medium
@@ -142,7 +144,7 @@ float vnoise(vec2 p) {
 // The rotated-octave FBM rotation constant — spliced from the shared chunk (AV.W2
 // — the one FBM_ROT). The loop below stays aurora-local (2.02 lacunarity +
 // uniform-driven octaves, per §3a — only the rotation constant converges).
-${o}
+${a}
 
 float fbm(vec2 p) {
   float v = 0.0;
@@ -156,6 +158,19 @@ float fbm(vec2 p) {
   return v;
 }
 
+// BB.B1 — the shared curl-noise flow chunk (curlFBM), spliced from the single source
+// procedural surfaces share (the AV.W2 shared-chunk precedent, applied to FLOW). The
+// chunk owns ONLY the basis-agnostic curl operator; the host owns the noise basis —
+// aurora wraps its OWN fbm as the scalar potential (so curl rides the same 2.02-
+// lacunarity loop the rest of the warp uses). The curl operator is consumed ONLY on
+// the OPT-IN warpMode == "curl" branch (uWarpMode == 3) — the default warp paths are
+// byte-untouched (no curl call), so every existing aurora gate stays byte-equivalent.
+${ge}
+
+// The scalar fbm potential ψ for curlFBM — aurora's own fbm. The splice above
+// forward-declares potentialFBM; this defines it (ES 3.00 allows the prototype).
+float potentialFBM(vec2 p) { return fbm(p); }
+
 // AX.W12 — the painterly-medium organic noise basis (Jarzynski PCG2D integer-bit hash
 // + 2D simplex gradient noise), SPLICED from the shared procedural-color chunk (the
 // single hash source).
@@ -164,19 +179,19 @@ float fbm(vec2 p) {
 // fbm above (the cost-tiering is preserved per AX.W12 §3a — this is an ADDITIVE basis,
 // not a blanket upgrade). The aurora-local hash21/vnoise/fbm loop stays (legitimately
 // divergent per the chunk's §3a).
-${u}
+${t}
 
 // The sRGB OETF (linearToSrgb) — spliced from the shared procedural-color chunk
 // (AV.W2 — the single OETF source; the AV.W1 local copy is deleted here). Aurora's
 // linear pipeline closes the seam with this transfer before fragColor (the
 // proof:aurora-space-gamma seam).
-${c}
+${r}
 
 // W5 — the four Ottosson OKLab/OKLCh matrices + their space fns, spliced from the
 // SAME shared chunk the goo-blob uses (zero new payload, 1e-6-verified). PI is
 // defined first because the chunk's oklabToOklch folds the hue into [0, 2pi).
 const float PI = 3.141592653589793;
-${n}
+${i}
 
 // ── OKLCh palette ramp (W5 interp · AX.W11 single-sourced twin) ───────────────
 // The muddy-midtone kill: distant-hue stops interpolated by a plain linear mix()
@@ -188,7 +203,7 @@ ${n}
 // the smoothstep-eased samplePaletteRamp dispatcher) is SPLICED from the shared
 // procedural-color chunk (the single ramp source). composition.glsl.ts's
 // samplePalette calls samplePaletteRamp.
-${i}
+${n}
 
 // Interleaved Gradient Noise (Jimenez) — a 1-LSB triangular dither applied in
 // DISPLAY space AFTER the OETF, so 8-bit mid-tone banding on the soft gradient is
@@ -232,6 +247,13 @@ vec2 domainWarp(vec2 p, float t) {
     float c1 = cellular(p * uWarpScale * 1.2);
     float c2 = cellular(p * uWarpScale * 1.2 + vec2(11.0, 7.0));
     warp = mix(r, vec2(c1, c2), 0.5);
+  } else if (uWarpMode == 3) {
+    // BB.B1 — curl (Bridson flow warp). The divergence-free curl of an fbm potential
+    // advects the field along a source-free swirl — folds + stretches like fluid
+    // advection, never the source-y bulge a raw fbm gradient produces. The potential
+    // scrolls on the same warp-drift clock; curlFBM is the shared chunk's operator.
+    vec2 fp = p * uWarpScale + vec2(t * uWarpDrift * K_WARP);
+    warp = curlFBM(fp);
   }
   vec2 warped = p + uWarpAmount * warp;
 
@@ -257,14 +279,35 @@ vec2 domainWarp(vec2 p, float t) {
   return warped;
 }
 
-` + le + A + ue + A + "// ── Color utils (OKLCh — W5; the sRGB YIQ hueShift matrix is DELETED) ──────────\nconst vec3 W_LUMA = vec3(0.2126, 0.7152, 0.0722);\n\n// Linear sRGB → OKLab (the cbrt-LMS path; the goo-blob's srgbToOklab takes gamma,\n// here the color field is already linear so we skip the OETF inverse).\nvec3 linOklab(vec3 lin) {\n  vec3 lms = LINEAR_SRGB_TO_LMS * lin;\n  return LMS_TO_OKLAB * (sign(lms) * pow(abs(lms), vec3(1.0 / 3.0)));\n}\n\n// Broken color (W5) — hue + chroma jitter at fixed PERCEPTUAL lightness, in OKLCh.\n// Broken color is hue variation at constant VALUE, which only OKLCh makes true (the\n// prior YIQ-style sRGB rotation muddied value). This is the seam AW.W4's van-Gogh /\n// oil-pastel per-stroke pigment jitter consumes (the pigment-variation read).\nvec3 brokenColorJitter(vec3 c, float hueSeed, float valueSeed, float strength) {\n  float amt = clamp(uBrokenColor * strength, 0.0, 1.0);\n  if (amt <= 0.001) return c;\n  vec3 lch = oklabToOklch(linOklab(c));\n  // ±~16° hue swing + ±~12% chroma swing; L held so value reads constant.\n  lch.z += (hueSeed - 0.5) * (32.0 * PI / 180.0) * amt;\n  lch.y = max(lch.y * (1.0 + (valueSeed - 0.5) * 0.28 * amt), 0.0);\n  return max(oklabToLinearSrgb(oklchToOklab(lch)), vec3(0.0));\n}\n\n// Saturation (W5) — chroma scale at fixed L/H in OKLCh (the perceptual saturation,\n// replacing the luma-mix sRGB form). amt=1 is identity; amt<1 desaturates toward\n// the achromatic axis with NO lightness shift (the muddy-luma-mix grey is gone).\nvec3 saturate3(vec3 c, float amt) {\n  vec3 lch = oklabToOklch(linOklab(c));\n  lch.y = max(lch.y * amt, 0.0);\n  return max(oklabToLinearSrgb(oklchToOklab(lch)), vec3(0.0));\n}\n\n" + de + A + me + A + fe + A + he + A + "\n// ── Main ──────────────────────────────────────────────────────────────────\nvoid main() {\n  // Normalized 0..1 coordinates for nuclei, domain warp, and medium sampling.\n  vec2 uv = vUv;\n  vec2 pN = uv;\n\n  float t = uTime;\n\n  // Warp in pN space\n  vec2 p_warp = domainWarp(pN, t);\n\n  // Composition\n  float paletteId; float valueMod;\n  nucleiField(p_warp, t, paletteId, valueMod);\n  vec3 col = samplePalette(paletteId);\n  col *= 1.0 + uValueVariance * valueMod;\n\n  // Breath — slow global luminance wobble\n  float breath = sin(t * 6.2831 / max(uBreathPeriod, 1.0));\n  col *= 1.0 + uBreathDepth * breath * 0.5;\n\n  // Medium dispatch (AX.W13 — each is a FIRST-CLASS body, no shared dispatch):\n  //   1 pastel · 2 watercolor · 3 oil · 4 crayon (DRY tooth-multiply) ·\n  //   5 van-Gogh (atomic comma/crescent dabs) · 6 oil-pastel (stroke deposition).\n  // van-Gogh is no longer a mediumOil passthrough; oil-pastel and crayon no longer\n  // share a body — they share the SUBSTRATE (the tooth/placement helpers), not the\n  // dispatch body (slice 8 F0/F1).\n  if (uMedium == 1) col = mediumPastel(col, pN, t);\n  else if (uMedium == 2) col = mediumWatercolor(col, pN, t);\n  else if (uMedium == 3) col = mediumOil(col, pN, t);\n  else if (uMedium == 4) col = mediumCrayon(col, pN, t);\n  else if (uMedium == 5) col = mediumVangogh(col, pN, t);\n  else if (uMedium == 6) col = mediumOilPastel(col, pN, t);\n\n  // Saturation trim\n  col = saturate3(col, uSaturation);\n\n  // Tonemap + film grain\n  col = aces(col);\n  float grain = hash21(gl_FragCoord.xy + t * 17.0);\n  col += (grain - 0.5) * uPaperGrain;\n\n  col = clamp(col * 0.985 + 0.008, 0.0, 1.0);\n\n  // MANDATORY OETF — closes the seam (mirrors metaball.frag.ts:278). The whole\n  // pipeline above (linear palette LUT, ACES tonemap, grain) runs in LINEAR; without\n  // this transfer the aurora ships ~2.2× too dark (linear 0.5 → ~0.215, not ~0.735).\n  col = linearToSrgb(col);\n\n  // 1-LSB triangular IGN dither in DISPLAY space, AFTER the OETF — breaks 8-bit\n  // mid-tone banding on the soft gradient. Texture-free (the canonical Jimenez IGN).\n  col += (1.0 / 255.0) * (ign(gl_FragCoord.xy) - 0.5);\n\n  // Premultiply operates on DISPLAY-space color (matches the blob: OETF, then premul).\n  fragColor = vec4(col * uAlpha, uAlpha);\n}\n", j = "[Aurora]", M = /* @__PURE__ */ "uTime.uPalette.uStopCount.uNucleiCount.uNucleiPos.uNucleiRadius.uNucleiPaletteBias.uNucleiValueBias.uNucleiDriftRadius.uNucleiDriftPhase.uNucleiElong.uNucleiAngle.uSoftmaxBeta.uValueVariance.uWarpAmount.uWarpScale.uWarpDrift.uWarpMode.uNoiseOctaves.uMedium.uHuePath.uFlowPattern.uFlowFocal.uFlowAngle.uFlowCurl.uCursor.uCursorStrength.uCursorRadius.uCursorVelocity.uCursorBurst.uStrokeAmount.uStrokeScale.uStrokeAnisotropy.uStrokeLayers.uStrokeMode.uStrokeOrient.uLightDir.uLightColor.uWetEdge.uGranulation.uImpasto.uBrokenColor.uCanvasGrain.uNucleiDrift.uPaletteDrift.uBreathDepth.uBreathPeriod.uSaturation.uPaperGrain.uAlpha".split(".");
-function N(e, t, n) {
-	let r = l(e, e.VERTEX_SHADER, t, j), i = l(e, e.FRAGMENT_SHADER, n, j), o = a(e, r, i, j);
-	e.useProgram(o);
-	let s = e.createVertexArray();
-	e.bindVertexArray(s);
-	let c = e.createBuffer();
-	e.bindBuffer(e.ARRAY_BUFFER, c), e.bufferData(e.ARRAY_BUFFER, new Float32Array([
+` + _e + k + ve + k + "// ── Color utils (OKLCh — W5; the sRGB YIQ hueShift matrix is DELETED) ──────────\nconst vec3 W_LUMA = vec3(0.2126, 0.7152, 0.0722);\n\n// Linear sRGB → OKLab (the cbrt-LMS path; the goo-blob's srgbToOklab takes gamma,\n// here the color field is already linear so we skip the OETF inverse).\nvec3 linOklab(vec3 lin) {\n  vec3 lms = LINEAR_SRGB_TO_LMS * lin;\n  return LMS_TO_OKLAB * (sign(lms) * pow(abs(lms), vec3(1.0 / 3.0)));\n}\n\n// Broken color (W5) — hue + chroma jitter at fixed PERCEPTUAL lightness, in OKLCh.\n// Broken color is hue variation at constant VALUE, which only OKLCh makes true (the\n// prior YIQ-style sRGB rotation muddied value). This is the seam AW.W4's van-Gogh /\n// oil-pastel per-stroke pigment jitter consumes (the pigment-variation read).\nvec3 brokenColorJitter(vec3 c, float hueSeed, float valueSeed, float strength) {\n  float amt = clamp(uBrokenColor * strength, 0.0, 1.0);\n  if (amt <= 0.001) return c;\n  vec3 lch = oklabToOklch(linOklab(c));\n  // ±~16° hue swing + ±~12% chroma swing; L held so value reads constant.\n  lch.z += (hueSeed - 0.5) * (32.0 * PI / 180.0) * amt;\n  lch.y = max(lch.y * (1.0 + (valueSeed - 0.5) * 0.28 * amt), 0.0);\n  return max(oklabToLinearSrgb(oklchToOklab(lch)), vec3(0.0));\n}\n\n// Saturation (W5) — chroma scale at fixed L/H in OKLCh (the perceptual saturation,\n// replacing the luma-mix sRGB form). amt=1 is identity; amt<1 desaturates toward\n// the achromatic axis with NO lightness shift (the muddy-luma-mix grey is gone).\nvec3 saturate3(vec3 c, float amt) {\n  vec3 lch = oklabToOklch(linOklab(c));\n  lch.y = max(lch.y * amt, 0.0);\n  return max(oklabToLinearSrgb(oklchToOklab(lch)), vec3(0.0));\n}\n\n" + ye + k + Ce + k + be + k + we + k + "\n// ── Main ──────────────────────────────────────────────────────────────────\nvoid main() {\n  // Normalized 0..1 coordinates for nuclei, domain warp, and medium sampling.\n  vec2 uv = vUv;\n  vec2 pN = uv;\n\n  float t = uTime;\n\n  // Warp in pN space\n  vec2 p_warp = domainWarp(pN, t);\n\n  // Composition\n  float paletteId; float valueMod;\n  nucleiField(p_warp, t, paletteId, valueMod);\n  vec3 col = samplePalette(paletteId);\n  col *= 1.0 + uValueVariance * valueMod;\n\n  // Breath — slow global luminance wobble\n  float breath = sin(t * 6.2831 / max(uBreathPeriod, 1.0));\n  col *= 1.0 + uBreathDepth * breath * 0.5;\n\n  // Medium dispatch (AX.W13 — each is a FIRST-CLASS body, no shared dispatch):\n  //   1 pastel · 2 watercolor · 3 oil · 4 crayon (DRY tooth-multiply) ·\n  //   5 van-Gogh (atomic comma/crescent dabs) · 6 oil-pastel (stroke deposition) ·\n  //   7 kuwahara (anisotropic-Kuwahara edge-preserving painterly finish — BB.W-AUR-KUWAHARA).\n  // van-Gogh is no longer a mediumOil passthrough; oil-pastel and crayon no longer\n  // share a body — they share the SUBSTRATE (the tooth/placement helpers), not the\n  // dispatch body (slice 8 F0/F1). Kuwahara is OPT-IN (default-off — no preset selects\n  // it unless a config explicitly sets medium:\"kuwahara\"); the smooth default + the\n  // van-Gogh HERO + the existing painterly bodies are byte-unchanged by its addition.\n  if (uMedium == 1) col = mediumPastel(col, pN, t);\n  else if (uMedium == 2) col = mediumWatercolor(col, pN, t);\n  else if (uMedium == 3) col = mediumOil(col, pN, t);\n  else if (uMedium == 4) col = mediumCrayon(col, pN, t);\n  else if (uMedium == 5) col = mediumVangogh(col, pN, t);\n  else if (uMedium == 6) col = mediumOilPastel(col, pN, t);\n  else if (uMedium == 7) col = mediumKuwahara(col, pN, t);\n\n  // Saturation trim\n  col = saturate3(col, uSaturation);\n\n  // Tonemap + film grain\n  col = aces(col);\n  float grain = hash21(gl_FragCoord.xy + t * 17.0);\n  col += (grain - 0.5) * uPaperGrain;\n\n  col = clamp(col * 0.985 + 0.008, 0.0, 1.0);\n\n  // MANDATORY OETF — closes the seam (mirrors metaball.frag.ts:278). The whole\n  // pipeline above (linear palette LUT, ACES tonemap, grain) runs in LINEAR; without\n  // this transfer the aurora ships ~2.2× too dark (linear 0.5 → ~0.215, not ~0.735).\n  col = linearToSrgb(col);\n\n  // 1-LSB triangular IGN dither in DISPLAY space, AFTER the OETF — breaks 8-bit\n  // mid-tone banding on the soft gradient. Texture-free (the canonical Jimenez IGN).\n  col += (1.0 / 255.0) * (ign(gl_FragCoord.xy) - 0.5);\n\n  // Premultiply operates on DISPLAY-space color (matches the blob: OETF, then premul).\n  fragColor = vec4(col * uAlpha, uAlpha);\n}\n";
+//#endregion
+//#region src/components/custom/aurora/constants/renderMode.ts
+function A() {
+	let e;
+	try {
+		e = f();
+	} catch {
+		return !1;
+	}
+	return e === null ? !1 : e.includes("swiftshader") || e.includes("llvmpipe") || e.includes("software") || e.includes("basic render") || e.includes("microsoft basic");
+}
+function Ee(e, t = {}) {
+	if (e !== "css" && !t.forceWebGLUnderSoftwareRaster && A()) return "css";
+	if (e !== "auto") return e;
+	if (typeof navigator > "u" || typeof window > "u") return "webgl";
+	let n = typeof navigator.hardwareConcurrency == "number" && navigator.hardwareConcurrency <= 4, r = typeof window.matchMedia == "function" && window.matchMedia("(prefers-reduced-motion: reduce)").matches, i = navigator.connection?.saveData === !0;
+	return n || r || i ? "css" : "webgl";
+}
+//#endregion
+//#region src/components/custom/aurora/composables/glSetup.ts
+var j = "[Aurora]", De = /* @__PURE__ */ "uTime.uPalette.uStopCount.uNucleiCount.uNucleiPos.uNucleiRadius.uNucleiPaletteBias.uNucleiValueBias.uNucleiDriftRadius.uNucleiDriftPhase.uNucleiElong.uNucleiAngle.uSoftmaxBeta.uValueVariance.uWarpAmount.uWarpScale.uWarpDrift.uWarpMode.uNoiseOctaves.uMedium.uHuePath.uFlowPattern.uFlowFocal.uFlowAngle.uFlowCurl.uCursor.uCursorStrength.uCursorRadius.uCursorVelocity.uCursorBurst.uStrokeAmount.uStrokeScale.uStrokeAnisotropy.uStrokeLayers.uStrokeMode.uStrokeOrient.uLightDir.uLightColor.uWetEdge.uGranulation.uImpasto.uBrokenColor.uCanvasGrain.uNucleiDrift.uPaletteDrift.uBreathDepth.uBreathPeriod.uSaturation.uPaperGrain.uAlpha".split(".");
+function Oe(e, t, n) {
+	let r = h(e, e.VERTEX_SHADER, t, j), i = h(e, e.FRAGMENT_SHADER, n, j), a = m(e, r, i, j);
+	e.useProgram(a);
+	let o = e.createVertexArray();
+	e.bindVertexArray(o);
+	let s = e.createBuffer();
+	e.bindBuffer(e.ARRAY_BUFFER, s), e.bufferData(e.ARRAY_BUFFER, new Float32Array([
 		-1,
 		-1,
 		3,
@@ -272,134 +315,136 @@ function N(e, t, n) {
 		-1,
 		3
 	]), e.STATIC_DRAW);
-	let u = e.getAttribLocation(o, "aPos");
-	e.enableVertexAttribArray(u), e.vertexAttribPointer(u, 2, e.FLOAT, !1, 0, 0);
-	let d = {};
-	for (let t of M) d[t] = e.getUniformLocation(o, t);
+	let c = e.getAttribLocation(a, "aPos");
+	e.enableVertexAttribArray(c), e.vertexAttribPointer(c, 2, e.FLOAT, !1, 0, 0);
+	let l = {};
+	for (let t of De) l[t] = e.getUniformLocation(a, t);
 	return {
-		program: o,
+		program: a,
 		vs: r,
 		fs: i,
-		uniforms: d,
+		uniforms: l,
 		geometry: {
-			vao: s,
-			buf: c
+			vao: o,
+			buf: s
 		}
 	};
 }
 //#endregion
 //#region src/components/custom/aurora/composables/color.ts
-function P(e, t = 8, n) {
+function M(e, t = 8, n) {
 	let r = n ?? new Float32Array(t * 3), i = Math.min(e.length, t);
 	for (let t = 0; t < i; t++) {
-		let [n, i, a] = h(e[t]);
+		let [n, i, a] = b(e[t]);
 		r[t * 3 + 0] = n, r[t * 3 + 1] = i, r[t * 3 + 2] = a;
 	}
 	for (let e = i; e < t; e++) r[e * 3 + 0] = 0, r[e * 3 + 1] = 0, r[e * 3 + 2] = 0;
 	return r;
 }
-function F(e) {
+function N(e) {
 	if (e.length === 0) return "transparent";
-	if (e.length === 1) return v(e[0]);
+	if (e.length === 1) return x(e[0]);
 	let t = e.length;
 	return `linear-gradient(135deg, ${e.map((e, n) => {
 		let r = Math.round(n / (t - 1) * 100);
-		return `${v(e)} ${r}%`;
+		return `${x(e)} ${r}%`;
 	}).join(", ")})`;
 }
-function I(e) {
-	let t = e.replace("#", ""), n = parseInt(t.slice(0, 2), 16), r = parseInt(t.slice(2, 4), 16), i = parseInt(t.slice(4, 6), 16), [a, o, s] = se(n / 255, r / 255, i / 255), [c, l, u] = oe(a, o, s);
+function ke(e) {
+	let t = e.replace("#", ""), n = parseInt(t.slice(0, 2), 16), r = parseInt(t.slice(2, 4), 16), i = parseInt(t.slice(4, 6), 16), [a, o, s] = me(n / 255, r / 255, i / 255), [c, l, u] = pe(a, o, s);
 	return {
 		L: c,
 		C: l,
 		h: u
 	};
 }
-var L = [.35, .95];
-function R(e, t = {}) {
-	let n = typeof e == "string" ? y(e) : e, { stopCount: r = 4, harmony: i = "analogous", lightnessSpread: a = .32, chromaFalloff: o = .85, hueSpread: s = 28, lightnessEasing: c = "linear", chromaEasing: l = "bell", temperatureShift: u = 0 } = t, d = Math.max(2, Math.min(8, Math.round(r))), [f, p] = L, h = n.L - a / 2, _ = n.L + a / 2;
-	h < f && (_ += f - h, h = f), _ > p && (h -= _ - p, _ = p), h = Math.max(f, h), _ = Math.min(p, _);
-	let v = [];
+var Ae = [.35, .95];
+function P(e, t = {}) {
+	let n = typeof e == "string" ? S(e) : e, { stopCount: r = 4, harmony: i = "analogous", lightnessSpread: a = .32, chromaFalloff: o = .85, hueSpread: s = 28, lightnessEasing: c = "linear", chromaEasing: l = "bell", temperatureShift: u = 0 } = t, d = Math.max(2, Math.min(8, Math.round(r))), [f, p] = Ae, m = n.L - a / 2, h = n.L + a / 2;
+	m < f && (h += f - m, m = f), h > p && (m -= h - p, h = p), m = Math.max(f, m), h = Math.min(p, h);
+	let g = [];
 	for (let e = 0; e < d; e++) {
-		let t = d === 1 ? 0 : e / (d - 1), r = h + (_ - h) * z(t, c), a = n.C * o, f;
-		if (l === "bell") f = a + (n.C - a) * _e(t);
+		let t = d === 1 ? 0 : e / (d - 1), r = m + (h - m) * F(t, c), a = n.C * o, f;
+		if (l === "bell") f = a + (n.C - a) * je(t);
 		else {
-			let e = l === "sine" ? z(t, "sine") : t;
+			let e = l === "sine" ? F(t, "sine") : t;
 			f = n.C - (n.C - a) * e;
 		}
-		let p = g(n.h, i, s, t);
-		u > 0 && (p = xe(p, t, u)), v.push(m({
+		let p = te(n.h, i, s, t);
+		u > 0 && (p = Fe(p, t, u)), g.push(ee({
 			L: r,
 			C: f,
 			h: p
 		}));
 	}
-	return v;
+	return g;
 }
-function z(e, t) {
+function F(e, t) {
 	return t === "sine" ? .5 - .5 * Math.cos(Math.PI * e) : e;
 }
-function _e(e) {
+function je(e) {
 	return Math.sin(Math.PI * e);
 }
-var ve = 70, ye = 250, be = 22;
-function xe(e, t, n) {
+var Me = 70, Ne = 250, Pe = 22;
+function Fe(e, t, n) {
 	let r = (t - .5) * 2;
 	if (r === 0 || n === 0) return e;
-	let i = (ae(e, r >= 0 ? ve : ye, 1, "shorter") - e + 540) % 360 - 180;
-	return ((e + Math.sign(i) * Math.min(Math.abs(r) * n * be, Math.abs(i))) % 360 + 360) % 360;
+	let i = (fe(e, r >= 0 ? Me : Ne, 1, "shorter") - e + 540) % 360 - 180;
+	return ((e + Math.sign(i) * Math.min(Math.abs(r) * n * Pe, Math.abs(i))) % 360 + 360) % 360;
 }
 //#endregion
 //#region src/components/custom/aurora/composables/uniformBridge.ts
-var Se = {
+var I = {
 	smooth: 0,
 	pastel: 1,
 	watercolor: 2,
 	oil: 3,
 	crayon: 4,
 	vangogh: 5,
-	"oil-pastel": 6
-}, Ce = {
+	"oil-pastel": 6,
+	kuwahara: 7
+}, L = {
 	shorter: 0,
 	longer: 1,
 	increasing: 2,
 	decreasing: 3
-}, we = {
+}, Ie = {
 	none: 0,
 	radial: 1,
 	swirl: 2,
 	diagonal: 3,
 	multi: 4,
 	tensor: 5
-}, B = {
+}, R = {
 	flow: 0,
 	tensor: 1
-}, Te = {
+}, z = {
 	fbm: 0,
 	cellular: 1,
-	hybrid: 2
-}, Ee = {
+	hybrid: 2,
+	curl: 3
+}, Le = {
 	oil: 0,
 	knife: 1,
 	chunky: 3
-}, V = {
+}, B = {
 	L: .985,
 	C: .0125,
 	h: 77.5
 };
-function De(e) {
-	return e === void 0 ? _(V.L, V.C, V.h) : Array.isArray(e) ? e : _(e.L, e.C, e.h);
+function Re(e) {
+	return e === void 0 ? ne(B.L, B.C, B.h) : Array.isArray(e) ? e : ne(e.L, e.C, e.h);
 }
-function Oe(e) {
-	return Se[e.medium];
+function ze(e) {
+	return I[e.medium];
 }
-function ke(e) {
-	return Ee[e.strokeMode];
+function Be(e) {
+	return Le[e.strokeMode];
 }
-function H(e) {
-	return e.medium === "vangogh" || e.medium === "oil-pastel" || e.medium === "crayon" ? B.tensor : B[e.strokeOrient ?? "flow"];
+function Ve(e) {
+	return e.medium === "vangogh" || e.medium === "oil-pastel" || e.medium === "crayon" || e.medium === "kuwahara" ? R.tensor : R[e.strokeOrient ?? "flow"];
 }
-function Ae(e, t, n, r) {
+function He(e, t, n, r) {
 	let i = (e) => 1 - e, a = {
 		palette: new Float32Array(24),
 		pos: new Float32Array(12),
@@ -412,7 +457,7 @@ function Ae(e, t, n, r) {
 		angle: new Float32Array(6)
 	};
 	return function(o) {
-		e.useProgram(t), P(o.palette, 8, a.palette), e.uniform3fv(n.uPalette, a.palette), e.uniform1i(n.uStopCount, Math.min(o.palette.length, 8));
+		e.useProgram(t), M(o.palette, 8, a.palette), e.uniform3fv(n.uPalette, a.palette), e.uniform1i(n.uStopCount, Math.min(o.palette.length, 8));
 		let s = Math.min(o.nuclei.length, 6);
 		e.uniform1i(n.uNucleiCount, s);
 		for (let e = 0; e < s; e++) {
@@ -420,13 +465,13 @@ function Ae(e, t, n, r) {
 			a.pos[e * 2 + 0] = t.x, a.pos[e * 2 + 1] = i(t.y), a.rad[e] = t.radius, a.pb[e] = t.paletteBias, a.vb[e] = t.valueBias, a.dr[e] = t.driftRadius, a.dp[e] = t.driftPhase, a.elong[e] = t.elongation ?? 1, a.angle[e] = -(t.angle ?? 0) * Math.PI / 180;
 		}
 		for (let e = s; e < 6; e++) a.pos[e * 2 + 0] = 0, a.pos[e * 2 + 1] = 0, a.rad[e] = 0, a.pb[e] = 0, a.vb[e] = 0, a.dr[e] = 0, a.dp[e] = 0, a.elong[e] = 1, a.angle[e] = 0;
-		e.uniform2fv(n.uNucleiPos, a.pos), e.uniform1fv(n.uNucleiRadius, a.rad), e.uniform1fv(n.uNucleiPaletteBias, a.pb), e.uniform1fv(n.uNucleiValueBias, a.vb), e.uniform1fv(n.uNucleiDriftRadius, a.dr), e.uniform1fv(n.uNucleiDriftPhase, a.dp), e.uniform1fv(n.uNucleiElong, a.elong), e.uniform1fv(n.uNucleiAngle, a.angle), e.uniform1f(n.uSoftmaxBeta, o.softmaxBeta), e.uniform1f(n.uValueVariance, o.valueVariance), e.uniform1f(n.uWarpAmount, o.warpAmount), e.uniform1f(n.uWarpScale, o.warpScale), e.uniform1f(n.uWarpDrift, o.warpDrift), e.uniform1i(n.uWarpMode, Te[o.warpMode]), e.uniform1i(n.uNoiseOctaves, o.noiseOctaves), e.uniform1i(n.uMedium, Oe(o)), e.uniform1i(n.uHuePath, Ce[o.huePath ?? "shorter"]), e.uniform1i(n.uFlowPattern, we[o.flow.pattern]), e.uniform2f(n.uFlowFocal, o.flow.focalX, i(o.flow.focalY)), e.uniform1f(n.uFlowAngle, o.flow.angle), e.uniform1f(n.uFlowCurl, o.flow.curl), e.uniform2f(n.uCursor, r.x, i(r.y)), e.uniform1f(n.uCursorStrength, r.strength), e.uniform1f(n.uCursorRadius, r.radius), e.uniform2f(n.uCursorVelocity, r.velX, -r.velY), e.uniform1f(n.uCursorBurst, r.burst), e.uniform1f(n.uStrokeAmount, o.strokeAmount), e.uniform1f(n.uStrokeScale, o.strokeScale), e.uniform1f(n.uStrokeAnisotropy, o.strokeAnisotropy), e.uniform1i(n.uStrokeLayers, o.strokeLayers), e.uniform1i(n.uStrokeMode, ke(o)), e.uniform1i(n.uStrokeOrient, H(o));
+		e.uniform2fv(n.uNucleiPos, a.pos), e.uniform1fv(n.uNucleiRadius, a.rad), e.uniform1fv(n.uNucleiPaletteBias, a.pb), e.uniform1fv(n.uNucleiValueBias, a.vb), e.uniform1fv(n.uNucleiDriftRadius, a.dr), e.uniform1fv(n.uNucleiDriftPhase, a.dp), e.uniform1fv(n.uNucleiElong, a.elong), e.uniform1fv(n.uNucleiAngle, a.angle), e.uniform1f(n.uSoftmaxBeta, o.softmaxBeta), e.uniform1f(n.uValueVariance, o.valueVariance), e.uniform1f(n.uWarpAmount, o.warpAmount), e.uniform1f(n.uWarpScale, o.warpScale), e.uniform1f(n.uWarpDrift, o.warpDrift), e.uniform1i(n.uWarpMode, z[o.warpMode]), e.uniform1i(n.uNoiseOctaves, o.noiseOctaves), e.uniform1i(n.uMedium, ze(o)), e.uniform1i(n.uHuePath, L[o.huePath ?? "shorter"]), e.uniform1i(n.uFlowPattern, Ie[o.flow.pattern]), e.uniform2f(n.uFlowFocal, o.flow.focalX, i(o.flow.focalY)), e.uniform1f(n.uFlowAngle, o.flow.angle), e.uniform1f(n.uFlowCurl, o.flow.curl), e.uniform2f(n.uCursor, r.x, i(r.y)), e.uniform1f(n.uCursorStrength, r.strength), e.uniform1f(n.uCursorRadius, r.radius), e.uniform2f(n.uCursorVelocity, r.velX, -r.velY), e.uniform1f(n.uCursorBurst, r.burst), e.uniform1f(n.uStrokeAmount, o.strokeAmount), e.uniform1f(n.uStrokeScale, o.strokeScale), e.uniform1f(n.uStrokeAnisotropy, o.strokeAnisotropy), e.uniform1i(n.uStrokeLayers, o.strokeLayers), e.uniform1i(n.uStrokeMode, Be(o)), e.uniform1i(n.uStrokeOrient, Ve(o));
 		{
 			let t = o.lightDir ?? [
 				-.5,
 				.6,
 				.62
-			], r = De(o.lightColor);
+			], r = Re(o.lightColor);
 			e.uniform3f(n.uLightDir, t[0], t[1], t[2]), e.uniform3f(n.uLightColor, r[0], r[1], r[2]);
 		}
 		e.uniform1f(n.uWetEdge, o.wetEdge), e.uniform1f(n.uGranulation, o.granulation), e.uniform1f(n.uImpasto, o.impasto), e.uniform1f(n.uBrokenColor, o.brokenColor), e.uniform1f(n.uCanvasGrain, o.canvasGrain), e.uniform1f(n.uNucleiDrift, o.nucleiDrift), e.uniform1f(n.uPaletteDrift, o.paletteDrift), e.uniform1f(n.uBreathDepth, o.breathDepth), e.uniform1f(n.uBreathPeriod, o.breathPeriod), e.uniform1f(n.uSaturation, o.saturation), e.uniform1f(n.uPaperGrain, o.paperGrain), e.uniform1f(n.uAlpha, o.alpha);
@@ -434,8 +479,8 @@ function Ae(e, t, n, r) {
 }
 //#endregion
 //#region src/components/custom/aurora/composables/cursorModel.ts
-var U = .22, je = .18, Me = .992, Ne = .96, W = .3;
-function Pe() {
+var V = .22, Ue = .18, We = .992, Ge = .96, H = .3;
+function Ke() {
 	return {
 		x: .5,
 		y: .5,
@@ -449,22 +494,22 @@ function Pe() {
 		burst: 0
 	};
 }
-function Fe(e, t = 1) {
+function U(e, t = 1) {
 	let n = Math.max(0, Math.min(1, t));
-	e.x += (e.targetX - e.x) * U, e.y += (e.targetY - e.y) * U, e.strength += (e.targetStrength - e.strength) * je, e.targetStrength *= Me, e.velX *= (1 - W) * n, e.velY *= (1 - W) * n, e.burst *= Ne * n;
+	e.x += (e.targetX - e.x) * V, e.y += (e.targetY - e.y) * V, e.strength += (e.targetStrength - e.strength) * Ue, e.targetStrength *= We, e.velX *= (1 - H) * n, e.velY *= (1 - H) * n, e.burst *= Ge * n;
 }
-function Ie(e, t, n) {
+function qe(e, t, n) {
 	let r = .12, i = Math.max(-r, Math.min(r, t)), a = Math.max(-r, Math.min(r, n));
-	e.velX += (i - e.velX) * W, e.velY += (a - e.velY) * W;
+	e.velX += (i - e.velX) * H, e.velY += (a - e.velY) * H;
 	let o = Math.hypot(i, a);
 	e.burst = Math.min(1, e.burst + o * 4);
 }
-function Le(e) {
+function W(e) {
 	return e.targetStrength > .001 || e.strength > .001 || Math.abs(e.x - e.targetX) > .001 || Math.abs(e.y - e.targetY) > .001;
 }
 //#endregion
 //#region src/components/custom/aurora/composables/frameLoop.ts
-function Re(e) {
+function Je(e) {
 	let { gl: t, prog: n, uniforms: r, cursor: i, getConfig: a, getReducedMotion: o } = e, s = (e) => 1 - e;
 	function c() {
 		return +!o();
@@ -479,47 +524,496 @@ function Re(e) {
 	function u() {
 		if (o()) return !1;
 		let e = a();
-		return e.nucleiDrift !== 0 || e.paletteDrift !== 0 || e.breathDepth !== 0 || e.warpDrift !== 0 || e.interactivity?.light && i.burst > .001 ? !0 : Le(i);
+		return e.nucleiDrift !== 0 || e.paletteDrift !== 0 || e.breathDepth !== 0 || e.warpDrift !== 0 || e.interactivity?.light && i.burst > .001 ? !0 : W(i);
 	}
 	return {
 		frame: (e) => {
-			Fe(i, c()), l(e);
+			U(i, c()), l(e);
 		},
 		needsAnimation: u
 	};
 }
 //#endregion
+//#region src/components/custom/aurora/constants/shaders/aurora.wgsl.ts
+var Ye = `
+const MAX_NUCLEI: i32 = 6;
+const MAX_STOPS: i32 = 8;
+const PI: f32 = 3.141592653589793;
+
+// ── Uniforms (the typed-struct source-of-truth — see aurora.wgsl.uniforms.ts) ──
+// Scalars are packed into vec4 lanes to keep the std140/WGSL 16-byte struct alignment
+// trivial. nuc0[i] = (pos.x, pos.y, radius, paletteBias); nuc1[i] = (valueBias,
+// driftRadius, driftPhase, elong); nuc2[i] = (angle, _, _, _). palette[i].rgb is the
+// CPU-baked linear-sRGB stop.
+struct Uniforms {
+  // scalars0: (uTime, uSoftmaxBeta, uValueVariance, uWarpAmount)
+  scalars0: vec4<f32>,
+  // scalars1: (uWarpScale, uWarpDrift, uSaturation, uPaperGrain)
+  scalars1: vec4<f32>,
+  // scalars2: (uAlpha, uNucleiDrift, uPaletteDrift, uBreathDepth)
+  scalars2: vec4<f32>,
+  // scalars3: (uBreathPeriod, uCursorStrength, uCursorRadius, _pad)
+  scalars3: vec4<f32>,
+  // cursor: (uCursor.x, uCursor.y, _pad, _pad)
+  cursor: vec4<f32>,
+  // ints: (uStopCount, uNucleiCount, uWarpMode, uNoiseOctaves)
+  ints0: vec4<i32>,
+  // ints1: (uMedium, uHuePath, _, _)
+  ints1: vec4<i32>,
+  palette: array<vec4<f32>, 8>,
+  nuc0: array<vec4<f32>, 6>,
+  nuc1: array<vec4<f32>, 6>,
+  nuc2: array<vec4<f32>, 6>,
+};
+
+@group(0) @binding(0) var<uniform> u: Uniforms;
+
+// ── Authoring-scale → rad/sec lifts (byte-identical to aurora.frag.ts) ──
+const K_NUCLEI: f32 = 14.0;
+const K_PAL: f32 = 24.0;
+const K_WARP: f32 = 5.0;
+
+// ── Noise (aurora-local hash21/vnoise/fbm — legitimately divergent per AV.W2 §3a) ──
+fn hash21(p0: vec2<f32>) -> f32 {
+  var p = fract(p0 * vec2<f32>(123.34, 456.21));
+  p = p + dot(p, p + 45.32);
+  return fract(p.x * p.y);
+}
+
+fn hash22(p0: vec2<f32>) -> vec2<f32> {
+  let p = vec2<f32>(dot(p0, vec2<f32>(127.1, 311.7)),
+                    dot(p0, vec2<f32>(269.5, 183.3)));
+  return fract(sin(p) * 43758.5453);
+}
+
+fn vnoise(p: vec2<f32>) -> f32 {
+  let i = floor(p);
+  let f = fract(p);
+  let uu = f * f * (3.0 - 2.0 * f);
+  let a = hash21(i);
+  let b = hash21(i + vec2<f32>(1.0, 0.0));
+  let c = hash21(i + vec2<f32>(0.0, 1.0));
+  let d = hash21(i + vec2<f32>(1.0, 1.0));
+  return mix(mix(a, b, uu.x), mix(c, d, uu.x), uu.y);
+}
+
+${p}
+
+fn fbm(p0: vec2<f32>) -> f32 {
+  var v = 0.0;
+  var a = 0.5;
+  var p = p0;
+  let octaves = u.ints0.w;
+  for (var i = 0; i < 5; i = i + 1) {
+    if (i >= octaves) { break; }
+    v = v + a * vnoise(p);
+    p = FBM_ROT * p * 2.02;
+    a = a * 0.5;
+  }
+  return v;
+}
+
+${u}
+${d}
+${c}
+
+// Interleaved Gradient Noise (Jimenez) — display-space 1-LSB dither.
+fn ign(p: vec2<f32>) -> f32 {
+  return fract(52.9829189 * fract(dot(p, vec2<f32>(0.06711056, 0.00583715))));
+}
+
+// Cellular / Worley f1.
+fn cellular(p: vec2<f32>) -> f32 {
+  let i = floor(p);
+  let f = fract(p);
+  var m = 1e9;
+  for (var y = -1; y <= 1; y = y + 1) {
+    for (var x = -1; x <= 1; x = x + 1) {
+      let g = vec2<f32>(f32(x), f32(y));
+      let o = hash22(i + g);
+      let r = g + o - f;
+      m = min(m, dot(r, r));
+    }
+  }
+  return sqrt(m);
+}
+
+// ── Warp (Quilez double warp + the three warp modes + the cursor swirl) ──
+fn domainWarp(p: vec2<f32>, t: f32) -> vec2<f32> {
+  let warpScale = u.scalars1.x;
+  let warpDrift = u.scalars1.y;
+  let warpAmount = u.scalars0.w;
+  let warpMode = u.ints0.z;
+  let cursorStrength = u.scalars3.y;
+  let cursorRadius = u.scalars3.z;
+  let cursor = u.cursor.xy;
+
+  let q = vec2<f32>(fbm(p * warpScale + vec2<f32>(0.0, 0.0) + t * warpDrift * K_WARP),
+                    fbm(p * warpScale + vec2<f32>(5.2, 1.3) + t * warpDrift * K_WARP));
+  let r = vec2<f32>(fbm(p * warpScale + 4.0 * q + vec2<f32>(1.7, 9.2)),
+                    fbm(p * warpScale + 4.0 * q + vec2<f32>(8.3, 2.8)));
+
+  var warp = r;
+  if (warpMode == 1) {
+    // cellular — chunky territories
+    let c1 = cellular(p * warpScale * 1.5 + vec2<f32>(t * warpDrift * K_WARP * 2.0, 0.0));
+    let c2 = cellular(p * warpScale * 1.5 + vec2<f32>(11.0, 7.0 + t * warpDrift * K_WARP * 2.0));
+    warp = vec2<f32>(c1, c2);
+  } else if (warpMode == 2) {
+    // hybrid — fbm + cellular averaged
+    let c1 = cellular(p * warpScale * 1.2);
+    let c2 = cellular(p * warpScale * 1.2 + vec2<f32>(11.0, 7.0));
+    warp = mix(r, vec2<f32>(c1, c2), 0.5);
+  }
+  var warped = p + warpAmount * warp;
+
+  // Cursor swirl — rotate p around the cursor with radial falloff.
+  if (cursorStrength > 0.001) {
+    let toP = p - cursor;
+    let d = length(toP);
+    let rr = max(cursorRadius, 0.01);
+    let w = exp(-(d * d) / (rr * rr * 0.45));
+    let ang = w * cursorStrength * 2.1;
+    let ca = cos(ang);
+    let sa = sin(ang);
+    var rotated = vec2<f32>(ca * toP.x - sa * toP.y, sa * toP.x + ca * toP.y) + cursor;
+    let pinch = w * cursorStrength * 0.08;
+    rotated = mix(rotated, cursor, pinch);
+    warped = mix(warped, rotated + warpAmount * warp * 0.7, w * cursorStrength);
+  }
+  return warped;
+}
+
+// ── Palette LUT ──
+fn samplePalette(id0: f32) -> vec3<f32> {
+  let id = clamp(id0, 0.0, 1.0);
+  let stopCount = u.ints0.x;
+  let scaled = id * f32(stopCount - 1);
+  let i0 = i32(floor(scaled));
+  let i1 = min(i0 + 1, stopCount - 1);
+  let t = fract(scaled);
+  return samplePaletteRamp(u.palette[i0].rgb, u.palette[i1].rgb, t, u.ints1.y);
+}
+
+// ── Nuclei field (anisotropic Gaussian softmax + drift + palette drift) ──
+struct NucleiResult { paletteId: f32, valueMod: f32 };
+
+fn nucleiField(p: vec2<f32>, t: f32) -> NucleiResult {
+  let nucleiCount = u.ints0.y;
+  let nucleiDrift = u.scalars2.y;
+  let paletteDrift = u.scalars2.z;
+  let softmaxBeta = u.scalars0.y;
+
+  var accumBias = 0.0;
+  var accumValue = 0.0;
+  var accumW = 0.0;
+  for (var i = 0; i < MAX_NUCLEI; i = i + 1) {
+    if (i >= nucleiCount) { break; }
+    let n0 = u.nuc0[i];   // pos.x, pos.y, radius, paletteBias
+    let n1 = u.nuc1[i];   // valueBias, driftRadius, driftPhase, elong
+    let angle = u.nuc2[i].x;
+    let driftRadius = n1.y;
+    let driftPhase = n1.z;
+    let posI = n0.xy + driftRadius * vec2<f32>(
+      cos(t * nucleiDrift * K_NUCLEI + driftPhase),
+      sin(t * nucleiDrift * K_NUCLEI + driftPhase * 1.13)
+    );
+    let diff = p - posI;
+    let ca = cos(angle);
+    let sa = sin(angle);
+    let local = vec2<f32>(ca * diff.x + sa * diff.y, -sa * diff.x + ca * diff.y);
+    let along = local.x / max(n1.w, 0.01);
+    let across = local.y;
+    let d2 = along * along + across * across;
+    let r = max(n0.z, 0.01);
+    let w = exp(-softmaxBeta * d2 / (r * r));
+    accumBias = accumBias + w * n0.w;
+    accumValue = accumValue + w * n1.x;
+    accumW = accumW + w;
+  }
+  var paletteId = accumBias / max(accumW, 1e-4);
+  let valueMod = accumValue / max(accumW, 1e-4);
+
+  let palAmp = clamp(paletteDrift * 6.0, 0.0, 0.06);
+  paletteId = paletteId + palAmp * sin(t * paletteDrift * K_PAL);
+  paletteId = clamp(paletteId, 0.0, 1.0);
+  return NucleiResult(paletteId, valueMod);
+}
+
+// ── Color utils (OKLCh saturation) ──
+fn linOklab(lin: vec3<f32>) -> vec3<f32> {
+  let lms = LINEAR_SRGB_TO_LMS * lin;
+  return LMS_TO_OKLAB * cbrt3(lms);
+}
+
+fn saturate3(c: vec3<f32>, amt: f32) -> vec3<f32> {
+  var lch = oklabToOklch(linOklab(c));
+  lch.y = max(lch.y * amt, 0.0);
+  return max(oklabToLinearSrgb(oklchToOklab(lch)), vec3<f32>(0.0));
+}
+
+// ── PBR-Neutral tonemap (named aces() in the GLSL — the tonemap slot) ──
+fn aces(color0: vec3<f32>) -> vec3<f32> {
+  let startCompression = 0.8 - 0.04;
+  let desaturation = 0.15;
+  var color = color0;
+  let x = min(color.r, min(color.g, color.b));
+  var offset = 0.04;
+  if (x < 0.08) { offset = x - 6.25 * x * x; }
+  color = color - offset;
+  let peak = max(color.r, max(color.g, color.b));
+  if (peak < startCompression) { return clamp(color, vec3<f32>(0.0), vec3<f32>(1.0)); }
+  let d = 1.0 - startCompression;
+  let newPeak = 1.0 - d * d / (peak + d - startCompression);
+  color = color * (newPeak / peak);
+  let g = 1.0 - 1.0 / (desaturation * (peak - newPeak) + 1.0);
+  return clamp(mix(color, newPeak * vec3<f32>(1.0), g), vec3<f32>(0.0), vec3<f32>(1.0));
+}
+
+// ── Full-screen triangle vertex stage (the pilot idiom — no vertex buffer) ──
+struct VSOut {
+  @builtin(position) pos: vec4<f32>,
+  @location(0) uv: vec2<f32>,
+};
+
+@vertex
+fn vs_main(@builtin(vertex_index) vi: u32) -> VSOut {
+  // NDC corners of the covering triangle (-1,-1), (3,-1), (-1,3).
+  var p = vec2<f32>(-1.0, -1.0);
+  if (vi == 1u) { p = vec2<f32>(3.0, -1.0); }
+  else if (vi == 2u) { p = vec2<f32>(-1.0, 3.0); }
+  var out: VSOut;
+  out.pos = vec4<f32>(p, 0.0, 1.0);
+  // vUv = aPos * 0.5 + 0.5 (the aurora.vert convention).
+  out.uv = p * 0.5 + 0.5;
+  return out;
+}
+
+// ── Fragment stage (the main() transcription) ──
+@fragment
+fn fs_main(in: VSOut) -> @location(0) vec4<f32> {
+  let t = u.scalars0.x;
+  let valueVariance = u.scalars0.z;
+  let breathDepth = u.scalars2.w;
+  let breathPeriod = u.scalars3.x;
+  let saturation = u.scalars1.z;
+  let paperGrain = u.scalars1.w;
+  let alpha = u.scalars2.x;
+
+  let uv = in.uv;
+  let pN = uv;
+
+  let p_warp = domainWarp(pN, t);
+
+  let nf = nucleiField(p_warp, t);
+  var col = samplePalette(nf.paletteId);
+  col = col * (1.0 + valueVariance * nf.valueMod);
+
+  // Breath — slow global luminance wobble.
+  let breath = sin(t * 6.2831 / max(breathPeriod, 1.0));
+  col = col * (1.0 + breathDepth * breath * 0.5);
+
+  // The painterly mediums (uMedium 1-6) are the WebGL2 fallback's register; the WGSL
+  // primary renders the smooth core for the default parity capture (the booked
+  // W-AURORA-WGPU-MEDIUMS successor ports the painterly bodies).
+
+  // Saturation trim.
+  col = saturate3(col, saturation);
+
+  // Tonemap + film grain.
+  col = aces(col);
+  // gl_FragCoord.xy → the device-pixel position (the WGSL @builtin(position)).
+  let fragCoord = in.pos.xy;
+  let grain = hash21(fragCoord + t * 17.0);
+  col = col + (grain - 0.5) * paperGrain;
+
+  col = clamp(col * 0.985 + 0.008, vec3<f32>(0.0), vec3<f32>(1.0));
+
+  // MANDATORY OETF — closes the seam (the whole pipeline above runs in LINEAR).
+  col = linearToSrgb(col);
+
+  // 1-LSB triangular IGN dither in DISPLAY space, AFTER the OETF.
+  col = col + (1.0 / 255.0) * (ign(fragCoord) - 0.5);
+
+  // Premultiply on DISPLAY-space color (matches the WebGL2 path: OETF, then premul).
+  return vec4<f32>(col * alpha, alpha);
+}
+`, G = {
+	scalars0: 0,
+	scalars1: 4,
+	scalars2: 8,
+	scalars3: 12,
+	cursor: 16,
+	ints0: 20,
+	ints1: 24,
+	palette: 28,
+	nuc0: 60,
+	nuc1: 84,
+	nuc2: 108
+}, K = (e) => 1 - e;
+function Xe() {
+	let e = /* @__PURE__ */ new ArrayBuffer(528);
+	return {
+		buffer: e,
+		f32: new Float32Array(e),
+		i32: new Int32Array(e),
+		paletteScratch: new Float32Array(24)
+	};
+}
+function Ze(e, t, n, r) {
+	let { f32: i, i32: a, paletteScratch: o } = e;
+	i[G.scalars0 + 0] = r, i[G.scalars0 + 1] = t.softmaxBeta, i[G.scalars0 + 2] = t.valueVariance, i[G.scalars0 + 3] = t.warpAmount, i[G.scalars1 + 0] = t.warpScale, i[G.scalars1 + 1] = t.warpDrift, i[G.scalars1 + 2] = t.saturation, i[G.scalars1 + 3] = t.paperGrain, i[G.scalars2 + 0] = t.alpha, i[G.scalars2 + 1] = t.nucleiDrift, i[G.scalars2 + 2] = t.paletteDrift, i[G.scalars2 + 3] = t.breathDepth, i[G.scalars3 + 0] = t.breathPeriod, i[G.scalars3 + 1] = n.strength, i[G.scalars3 + 2] = n.radius, i[G.scalars3 + 3] = 0, i[G.cursor + 0] = n.x, i[G.cursor + 1] = K(n.y), i[G.cursor + 2] = 0, i[G.cursor + 3] = 0, a[G.ints0 + 0] = Math.min(t.palette.length, 8);
+	let s = Math.min(t.nuclei.length, 6);
+	a[G.ints0 + 1] = s, a[G.ints0 + 2] = z[t.warpMode], a[G.ints0 + 3] = t.noiseOctaves, a[G.ints1 + 0] = I[t.medium], a[G.ints1 + 1] = L[t.huePath ?? "shorter"], a[G.ints1 + 2] = 0, a[G.ints1 + 3] = 0, M(t.palette, 8, o);
+	for (let e = 0; e < 8; e++) {
+		let t = G.palette + e * 4;
+		i[t + 0] = o[e * 3 + 0], i[t + 1] = o[e * 3 + 1], i[t + 2] = o[e * 3 + 2], i[t + 3] = 0;
+	}
+	for (let e = 0; e < 6; e++) {
+		let n = G.nuc0 + e * 4, r = G.nuc1 + e * 4, a = G.nuc2 + e * 4;
+		if (e < s) {
+			let o = t.nuclei[e];
+			i[n + 0] = o.x, i[n + 1] = K(o.y), i[n + 2] = o.radius, i[n + 3] = o.paletteBias, i[r + 0] = o.valueBias, i[r + 1] = o.driftRadius, i[r + 2] = o.driftPhase, i[r + 3] = o.elongation ?? 1, i[a + 0] = -(o.angle ?? 0) * Math.PI / 180, i[a + 1] = 0, i[a + 2] = 0, i[a + 3] = 0;
+		} else i[n + 0] = 0, i[n + 1] = 0, i[n + 2] = 0, i[n + 3] = 0, i[r + 0] = 0, i[r + 1] = 0, i[r + 2] = 0, i[r + 3] = 1, i[a + 0] = 0, i[a + 1] = 0, i[a + 2] = 0, i[a + 3] = 0;
+	}
+	return e;
+}
+//#endregion
+//#region src/components/custom/aurora/composables/wgpuSetup.ts
+var Qe = 64, $e = 8, et = 2;
+function tt(e) {
+	let { canvas: t, cursor: n, getConfig: r, getReducedMotion: i } = e;
+	return function(e, a, o) {
+		let s = e.createShaderModule({
+			label: "[Aurora] aurora.wgsl",
+			code: Ye
+		}), c = e.createBuffer({
+			label: "[Aurora] uniforms",
+			size: 528,
+			usage: Qe | $e
+		}), u = e.createBindGroupLayout({
+			label: "[Aurora] bind-group-0",
+			entries: [{
+				binding: 0,
+				visibility: et,
+				buffer: { type: "uniform" }
+			}]
+		}), d = e.createRenderPipeline({
+			label: "[Aurora] pipeline",
+			layout: e.createPipelineLayout({ bindGroupLayouts: [u] }),
+			vertex: {
+				module: s,
+				entryPoint: "vs_main"
+			},
+			fragment: {
+				module: s,
+				entryPoint: "fs_main",
+				targets: [{
+					format: o,
+					blend: {
+						color: {
+							srcFactor: "one",
+							dstFactor: "one-minus-src-alpha",
+							operation: "add"
+						},
+						alpha: {
+							srcFactor: "one",
+							dstFactor: "one-minus-src-alpha",
+							operation: "add"
+						}
+					}
+				}]
+			},
+			primitive: { topology: "triangle-list" }
+		}), f = e.createBindGroup({
+			label: "[Aurora] bind-group",
+			layout: u,
+			entries: [{
+				binding: 0,
+				resource: { buffer: c }
+			}]
+		}), p = Xe();
+		function m() {
+			let e = l(), n = t.getBoundingClientRect(), r = t.parentElement?.getBoundingClientRect(), i = n.width || r?.width || 1, a = n.height || r?.height || 1, o = Math.max(1, Math.round(i * e)), s = Math.max(1, Math.round(a * e));
+			(t.width !== o || t.height !== s) && (t.width = o, t.height = s);
+		}
+		function h(t) {
+			U(n, +!i()), Ze(p, r(), n, t), e.queue.writeBuffer(c, 0, p.buffer);
+			let o = a.getCurrentTexture().createView(), s = e.createCommandEncoder({ label: "[Aurora] frame" }), l = s.beginRenderPass({ colorAttachments: [{
+				view: o,
+				clearValue: {
+					r: 0,
+					g: 0,
+					b: 0,
+					a: 0
+				},
+				loadOp: "clear",
+				storeOp: "store"
+			}] });
+			l.setPipeline(d), l.setBindGroup(0, f), l.draw(3, 1, 0, 0), l.end(), e.queue.submit([s.finish()]);
+		}
+		function g() {
+			if (i()) return !1;
+			let e = r();
+			return e.nucleiDrift !== 0 || e.paletteDrift !== 0 || e.breathDepth !== 0 || e.warpDrift !== 0 || e.interactivity?.light && n.burst > .001 ? !0 : W(n);
+		}
+		return {
+			frame: h,
+			shouldContinue: g,
+			resize: m,
+			time: (e) => i() ? 3.7 : e,
+			teardown: () => {
+				c.destroy();
+			}
+		};
+	};
+}
+//#endregion
 //#region src/components/custom/aurora/composables/runtime.ts
-function ze(e) {
+function nt(e) {
 	return e.mode === "capture" ? !0 : e.initStrategy === "eager";
 }
-function Be(e) {
+function rt(e) {
 	return typeof e.preserveDrawingBuffer == "boolean" ? e.preserveDrawingBuffer : e.mode === "capture";
 }
-function G(e, n, r = {}) {
-	let i = Be(r), a = n, o = Pe(), c = null, l = s(e, {
-		mode: r.mode === "capture" ? "capture" : "live",
+function q(e, t, n = {}) {
+	let r = rt(n), i = !n.forceWebGLUnderSoftwareRaster && A(), a = t, s = Ke(), c = null, u = i ? {
+		arm: () => {},
+		suspend: () => {},
+		resume: () => {},
+		wake: () => {},
+		renderAt: () => {},
+		dispose: () => {},
+		reducedMotion: !1
+	} : o(e, {
+		mode: n.mode === "capture" ? "capture" : "live",
 		contextAttrs: {
 			antialias: !1,
 			alpha: !0,
 			premultipliedAlpha: !0,
-			preserveDrawingBuffer: i
+			preserveDrawingBuffer: r
 		},
-		setup: (n) => {
-			let { program: r, vs: i, fs: s, uniforms: u, geometry: d } = N(n, ce, ge);
+		setupWGPU: tt({
+			canvas: e,
+			cursor: s,
+			getConfig: () => a,
+			getReducedMotion: () => u.reducedMotion
+		}),
+		setupGL: (t) => {
+			let { program: n, vs: r, fs: i, uniforms: o, geometry: d } = Oe(t, he, Te);
 			function f() {
-				let i = t(), a = e.getBoundingClientRect(), o = e.parentElement?.getBoundingClientRect(), s = a.width || o?.width || 1, c = a.height || o?.height || 1, l = Math.max(1, Math.round(s * i)), u = Math.max(1, Math.round(c * i));
-				(e.width !== l || e.height !== u) && (e.width = l, e.height = u), n.viewport(0, 0, l, u), n.useProgram(r);
+				let r = l(), i = e.getBoundingClientRect(), a = e.parentElement?.getBoundingClientRect(), o = i.width || a?.width || 1, s = i.height || a?.height || 1, c = Math.max(1, Math.round(o * r)), u = Math.max(1, Math.round(s * r));
+				(e.width !== c || e.height !== u) && (e.width = c, e.height = u), t.viewport(0, 0, c, u), t.useProgram(n);
 			}
-			let p = Ae(n, r, u, o), m = Re({
-				gl: n,
-				prog: r,
-				uniforms: u,
-				cursor: o,
+			let p = He(t, n, o, s), m = Je({
+				gl: t,
+				prog: n,
+				uniforms: o,
+				cursor: s,
 				getConfig: () => a,
-				getReducedMotion: () => l.reducedMotion
+				getReducedMotion: () => u.reducedMotion
 			});
-			return n.clearColor(0, 0, 0, 0), n.disable(n.DEPTH_TEST), n.enable(n.BLEND), n.blendFunc(n.ONE, n.ONE_MINUS_SRC_ALPHA), c = (e) => {
+			return t.clearColor(0, 0, 0, 0), t.disable(t.DEPTH_TEST), t.enable(t.BLEND), t.blendFunc(t.ONE, t.ONE_MINUS_SRC_ALPHA), c = (e) => {
 				a = e, p(e);
 			}, c(a), requestAnimationFrame(() => {
 				f(), requestAnimationFrame(f);
@@ -527,55 +1021,60 @@ function G(e, n, r = {}) {
 				frame: m.frame,
 				shouldContinue: m.needsAnimation,
 				resize: f,
-				time: (e) => l.reducedMotion ? 3.7 : e,
+				time: (e) => u.reducedMotion ? 3.7 : e,
 				teardown: () => {
-					n.deleteProgram(r), n.deleteShader(i), n.deleteShader(s), n.deleteBuffer(d.buf), n.deleteVertexArray(d.vao);
+					t.deleteProgram(n), t.deleteShader(r), t.deleteShader(i), t.deleteBuffer(d.buf), t.deleteVertexArray(d.vao);
 				}
 			};
 		}
 	});
-	function u(e, t, n = .8) {
-		o.targetX = e, o.targetY = t, o.targetStrength = n, l.wake();
-	}
-	function d(e, t) {
-		l.reducedMotion || (Ie(o, e, t), l.wake());
+	function d(e, t, n = .8) {
+		s.targetX = e, s.targetY = t, s.targetStrength = n, u.wake();
 	}
-	function f() {
-		o.targetStrength = 0, l.wake();
+	function f(e, t) {
+		u.reducedMotion || (qe(s, e, t), u.wake());
 	}
-	function p(e) {
-		o.radius = e, l.wake();
+	function p() {
+		s.targetStrength = 0, u.wake();
 	}
 	function m(e) {
-		l.wake();
+		s.radius = e, u.wake();
+	}
+	function h(e) {
+		u.wake();
 	}
-	return r.mode !== "capture" && ze(r) && l.arm(), {
-		arm: () => l.arm(),
+	function g() {
+		u.armAsync ? u.armAsync().catch((e) => {
+			n.onInitError?.(e instanceof Error ? e : Error(String(e)));
+		}) : u.arm();
+	}
+	return n.mode !== "capture" && nt(n) && g(), {
+		arm: () => g(),
 		update: (e) => {
-			a = e, c?.(e), l.wake();
+			a = e, c?.(e), u.wake();
 		},
-		setCursor: u,
-		clearCursor: f,
-		setCursorRadius: p,
-		setReducedMotion: m,
-		injectCursorVelocity: d,
+		setCursor: d,
+		clearCursor: p,
+		setCursorRadius: m,
+		setReducedMotion: h,
+		injectCursorVelocity: f,
 		get reducedMotion() {
-			return l.reducedMotion;
+			return u.reducedMotion;
 		},
-		pause: (e = "manual") => l.suspend(e),
-		resume: (e = "manual") => l.resume(e),
-		renderAt: (e) => l.renderAt(e),
-		dispose: () => l.dispose()
+		pause: (e = "manual") => u.suspend(e),
+		resume: (e = "manual") => u.resume(e),
+		renderAt: (e) => u.renderAt(e),
+		dispose: () => u.dispose()
 	};
 }
 //#endregion
 //#region src/components/custom/aurora/composables/configSource.ts
-function K(e) {
+function J(e) {
 	return typeof e == "function" ? e : e && typeof e == "object" && "value" in e ? () => e.value : () => e;
 }
 //#endregion
 //#region src/components/custom/aurora/composables/useAurora.ts
-function Ve(e) {
+function it(e) {
 	if (typeof window > "u") return () => {};
 	if (typeof window.requestIdleCallback == "function") {
 		let t = window.requestIdleCallback(() => e(), { timeout: 2e3 });
@@ -592,8 +1091,8 @@ function Ve(e) {
 		r = !0, cancelAnimationFrame(t), n !== void 0 && clearTimeout(n);
 	};
 }
-function q(e, t, n = {}, r = {}) {
-	let i = K(t), a = r.renderMode === "css", o = null, s = null, c = null, l = null, u = null, d = null, f = null, p = !1, m = D(!1), h = n.mode === "capture" || n.initStrategy === "eager";
+function Y(e, t, n = {}, r = {}) {
+	let i = J(t), a = r.renderMode === "css", o = null, s = null, c = null, l = null, u = null, d = null, f = null, p = !1, m = E(!1), h = n.mode === "capture" || n.initStrategy === "eager";
 	function g() {
 		f && o && o.setReducedMotion(f.matches);
 	}
@@ -614,7 +1113,7 @@ function q(e, t, n = {}, r = {}) {
 				_(e);
 				return;
 			}
-			m.value = !0, f = window.matchMedia("(prefers-reduced-motion: reduce)"), f.addEventListener("change", g), s = k(i, (e) => o?.update(e), { deep: !0 }), i().interactivity?.scroll && e.value && (l = k(te({
+			m.value = !0, f = window.matchMedia("(prefers-reduced-motion: reduce)"), f.addEventListener("change", g), s = O(i, (e) => o?.update(e), { deep: !0 }), i().interactivity?.scroll && e.value && (l = O(ie({
 				target: e,
 				trackExit: !0
 			}), (e) => {
@@ -626,7 +1125,7 @@ function q(e, t, n = {}, r = {}) {
 		let t = e.value;
 		if (!t || a) return;
 		try {
-			o = G(t, i(), {
+			o = q(t, i(), {
 				initStrategy: "deferred",
 				...n
 			});
@@ -638,13 +1137,13 @@ function q(e, t, n = {}, r = {}) {
 			v();
 			return;
 		}
-		d = ee(e, {
+		d = re(e, {
 			pause: () => o?.pause("off-screen"),
 			resume: () => o?.resume("off-screen")
 		}, { pauseWhenHidden: !1 });
 		let r = d.isIntersecting;
-		c = k(r, (e) => {
-			!e || p || u || (u = Ve(() => {
+		c = O(r, (e) => {
+			!e || p || u || (u = it(() => {
 				u = null, r.value && v();
 			}));
 		}, { immediate: !0 });
@@ -658,33 +1157,184 @@ function q(e, t, n = {}, r = {}) {
 		renderAt: (e) => o?.renderAt(e),
 		pause: () => o?.pause(),
 		resume: () => o?.resume(),
-		isArmed: E(m)
+		isArmed: de(m)
 	};
 }
 //#endregion
-//#region src/components/custom/aurora/constants/renderMode.ts
-function J(e) {
-	if (e !== "auto") return e;
-	if (typeof navigator > "u" || typeof window > "u") return "webgl";
-	let t = typeof navigator.hardwareConcurrency == "number" && navigator.hardwareConcurrency <= 4, n = typeof window.matchMedia == "function" && window.matchMedia("(prefers-reduced-motion: reduce)").matches, r = navigator.connection?.saveData === !0;
-	return t || n || r ? "css" : "webgl";
+//#region src/components/custom/aurora/composables/auroraFallbackGround.ts
+function X(e) {
+	let t = Math.max(0, Math.min(1, e));
+	return t <= .0031308 ? t * 12.92 : 1.055 * t ** (1 / 2.4) - .055;
+}
+function at([e, t, n]) {
+	let r = .76, i = Math.min(e, Math.min(t, n)), a = i < .08 ? i - 6.25 * i * i : .04, o = e - a, s = t - a, c = n - a, l = Math.max(o, Math.max(s, c)), u = (e) => Math.max(0, Math.min(1, e));
+	if (l < r) return [
+		u(o),
+		u(s),
+		u(c)
+	];
+	let d = 1 - r, f = 1 - d * d / (l + d - r), p = f / l;
+	o *= p, s *= p, c *= p;
+	let m = 1 - 1 / (.15 * (l - f) + 1);
+	return [
+		u(o + (f - o) * m),
+		u(s + (f - s) * m),
+		u(c + (f - c) * m)
+	];
+}
+function ot([e, t, n]) {
+	return .2126 * e + .7152 * t + .0722 * n;
+}
+function st(e, t) {
+	let n = e.length;
+	if (n === 0) return [
+		0,
+		0,
+		0
+	];
+	if (n === 1) return e[0];
+	let r = Math.max(0, Math.min(1, t)) * (n - 1), i = Math.floor(r), a = Math.min(i + 1, n - 1), o = r - i, s = e[i], c = e[a];
+	return [
+		s[0] + (c[0] - s[0]) * o,
+		s[1] + (c[1] - s[1]) * o,
+		s[2] + (c[2] - s[2]) * o
+	];
+}
+function ct(e, t, n, r) {
+	let i = 0, a = 0, o = 0;
+	for (let s of e) {
+		let e = n - s.x, c = r - s.y, l = (s.angle ?? 0) * Math.PI / 180, u = Math.cos(l), d = Math.sin(l), f = u * e + d * c, p = -d * e + u * c, m = f / Math.max(s.elongation ?? 1, .01), h = p, g = m * m + h * h, _ = Math.max(s.radius, .01), v = Math.exp(-t * g / (_ * _));
+		i += v * s.paletteBias, a += v * s.valueBias, o += v;
+	}
+	let s = Math.max(o, 1e-4);
+	return {
+		paletteId: Math.max(0, Math.min(1, i / s)),
+		valueMod: a / s
+	};
+}
+function lt(e, t, n) {
+	let r = e.palette.map((e) => b(e)), { paletteId: i, valueMod: a } = ct(e.nuclei, e.softmaxBeta, t, n), o = st(r, i), s = 1 + e.valueVariance * a, c = [
+		o[0] * s,
+		o[1] * s,
+		o[2] * s
+	];
+	return c = at(c), c = [
+		Math.max(0, Math.min(1, c[0] * .985 + .008)),
+		Math.max(0, Math.min(1, c[1] * .985 + .008)),
+		Math.max(0, Math.min(1, c[2] * .985 + .008))
+	], [
+		X(c[0]),
+		X(c[1]),
+		X(c[2])
+	];
+}
+function ut([e, t, n]) {
+	let r = (e) => Math.round(Math.max(0, Math.min(1, e)) * 255);
+	return `rgb(${r(e)}, ${r(t)}, ${r(n)})`;
+}
+function dt(e) {
+	return ot(e.map((e) => {
+		let t = Math.max(0, Math.min(1, e));
+		return t <= .04045 ? t / 12.92 : ((t + .055) / 1.055) ** 2.4;
+	}));
+}
+function ft(e, t = {}) {
+	let n = Math.max(2, Math.min(16, Math.round(t.grid ?? 8))), r = [], i = 0, a = [
+		0,
+		0,
+		0,
+		0
+	], o = [
+		0,
+		0,
+		0,
+		0
+	];
+	for (let t = 0; t < n; t++) for (let s = 0; s < n; s++) {
+		let c = (s + .5) / n, l = (t + .5) / n, u = lt(e, c, l);
+		r.push({
+			x: c,
+			y: l,
+			color: u
+		});
+		let d = dt(u);
+		i += d;
+		let f = (l < .5 ? 0 : 2) + (c < .5 ? 0 : 1);
+		a[f] += d, o[f] += 1;
+	}
+	let s = i / r.length, c = [
+		a[0] / Math.max(1, o[0]),
+		a[1] / Math.max(1, o[1]),
+		a[2] / Math.max(1, o[2]),
+		a[3] / Math.max(1, o[3])
+	], l = [
+		0,
+		0,
+		0
+	];
+	for (let e of r) l[0] += e.color[0], l[1] += e.color[1], l[2] += e.color[2];
+	l[0] /= r.length, l[1] /= r.length, l[2] /= r.length;
+	let u = pt(r, n);
+	if (u) return {
+		backgroundImage: `url("${u}")`,
+		backgroundColor: ut(l),
+		metrics: {
+			mean: s,
+			quadrants: c
+		}
+	};
+	let d = 100 / n * 1.4;
+	return {
+		backgroundImage: r.map((e) => {
+			let t = (e.x * 100).toFixed(2), n = (e.y * 100).toFixed(2), r = (e) => Math.round(Math.max(0, Math.min(1, e)) * 255), [i, a, o] = e.color, s = `rgba(${r(i)}, ${r(a)}, ${r(o)}, 0.7)`, c = `rgba(${r(i)}, ${r(a)}, ${r(o)}, 0)`;
+			return `radial-gradient(circle ${d.toFixed(2)}% at ${t}% ${n}%, ${s} 0%, ${c} 100%)`;
+		}).join(", "),
+		backgroundColor: ut(l),
+		metrics: {
+			mean: s,
+			quadrants: c
+		}
+	};
+}
+function pt(e, t) {
+	if (typeof document > "u") return null;
+	let n;
+	try {
+		n = document.createElement("canvas");
+	} catch {
+		return null;
+	}
+	n.width = t, n.height = t;
+	let r = n.getContext("2d");
+	if (!r) return null;
+	let i = r.createImageData(t, t);
+	for (let n of e) {
+		let e = Math.min(t - 1, Math.floor(n.x * t)), r = (Math.min(t - 1, Math.floor(n.y * t)) * t + e) * 4, a = (e) => Math.round(Math.max(0, Math.min(1, e)) * 255);
+		i.data[r] = a(n.color[0]), i.data[r + 1] = a(n.color[1]), i.data[r + 2] = a(n.color[2]), i.data[r + 3] = 255;
+	}
+	r.putImageData(i, 0, 0);
+	try {
+		return n.toDataURL("image/png");
+	} catch {
+		return null;
+	}
 }
 //#endregion
-//#region src/components/custom/aurora/Aurora.vue
-var He = /* @__PURE__ */ e(/* @__PURE__ */ ne({
+//#region src/components/custom/aurora/Aurora.vue?vue&type=script&setup=true&lang.ts
+var mt = ["data-aurora-substrate"], ht = /* @__PURE__ */ e(/* @__PURE__ */ se({
 	__name: "Aurora",
 	props: {
-		config: { default: () => p },
+		config: { default: () => y },
 		runtimeOptions: {},
 		onInitError: {},
 		renderMode: { default: "auto" },
 		opacityCeiling: { default: 1 }
 	},
 	setup(e, { expose: t }) {
-		let n = e, r = J(n.renderMode), i = b(() => Math.max(0, Math.min(1, n.opacityCeiling))), a = D(null), o = q(a, () => n.config, b(() => ({
+		let n = e, r = Ee(n.renderMode, { forceWebGLUnderSoftwareRaster: n.runtimeOptions?.forceWebGLUnderSoftwareRaster ?? !1 }), i = C(() => Math.max(0, Math.min(1, n.opacityCeiling))), a = E(null), o = Y(a, () => n.config, C(() => ({
 			...n.runtimeOptions ?? {},
 			...n.onInitError ? { onInitError: n.onInitError } : {}
-		})).value, { renderMode: r }), s = b(() => F(n.config.palette));
+		})).value, { renderMode: r }), s = C(() => r === "css"), c = C(() => ft(n.config)), l = C(() => s.value ? c.value.backgroundImage : N(n.config.palette)), u = C(() => s.value ? c.value.backgroundColor : "transparent");
 		return t({
 			config: n.config,
 			canvasRef: a,
@@ -696,25 +1346,29 @@ var He = /* @__PURE__ */ e(/* @__PURE__ */ ne({
 			pause: o.pause,
 			resume: o.resume,
 			isArmed: o.isArmed
-		}), (e, t) => (ie(), x("div", {
+		}), (e, t) => (ue(), ae("div", {
 			class: "aurora-root block h-full w-full overflow-hidden",
-			style: C({ "--aurora-opacity-ceiling": i.value })
-		}, [S("div", {
+			"data-aurora-substrate": D(r),
+			style: le({ "--aurora-opacity-ceiling": i.value })
+		}, [oe("div", {
 			class: "aurora-placeholder h-full w-full",
 			"aria-hidden": "true",
-			style: C({ backgroundImage: s.value })
-		}, null, 4), S("canvas", {
+			style: le({
+				backgroundImage: l.value,
+				backgroundColor: u.value
+			})
+		}, null, 4), oe("canvas", {
 			ref_key: "canvasRef",
 			ref: a,
 			"aria-hidden": "true",
-			class: re(["aurora-canvas block h-full w-full", { "aurora-canvas--armed": O(o).isArmed.value }])
-		}, null, 2)], 4));
+			class: ce(["aurora-canvas block h-full w-full", { "aurora-canvas--armed": D(o).isArmed.value }])
+		}, null, 2)], 12, mt));
 	}
-}), [["__scopeId", "data-v-b51a5fc3"]]);
+}), [["__scopeId", "data-v-c05f182a"]]);
 //#endregion
 //#region src/components/custom/aurora/composables/useCursorInteraction.ts
-function Ue(e, t, n) {
-	let r = K(t), i = n.hitFraction ?? .35, a = null, o = null, s = null, c = null;
+function gt(e, t, n) {
+	let r = J(t), i = n.hitFraction ?? .35, a = null, o = null, s = null, c = null;
 	function l(e, t) {
 		let n = t.getBoundingClientRect();
 		return {
@@ -810,12 +1464,11 @@ function Ue(e, t, n) {
 	}), { isDragging: () => a };
 }
 //#endregion
-//#region src/components/custom/aurora/composables/atoms.ts
-var Y = {};
-function X(e, t, n) {
+//#region src/components/custom/aurora/composables/atoms-fields.ts
+function Z(e, t, n) {
 	return e + (t - e) * n;
 }
-var Z = {
+var Q = {
 	saturation: {
 		calm: .85,
 		vivid: 1.2
@@ -832,7 +1485,7 @@ var Z = {
 		calm: .15,
 		vivid: .5
 	}
-}, We = {
+}, _t = {
 	still: {
 		nucleiDrift: 0,
 		paletteDrift: 0,
@@ -840,10 +1493,10 @@ var Z = {
 		breathDepth: 0
 	},
 	breathing: {
-		nucleiDrift: 0,
-		paletteDrift: 0,
-		warpDrift: 0,
-		breathDepth: .05
+		nucleiDrift: .005,
+		paletteDrift: .006,
+		warpDrift: .003,
+		breathDepth: .06
 	},
 	drifting: {
 		nucleiDrift: .015,
@@ -851,34 +1504,74 @@ var Z = {
 		warpDrift: .008,
 		breathDepth: .05
 	}
-}, Ge = .61803398875;
-function Q(e, t = "composed") {
-	let n = Math.max(1, Math.min(6, Math.round(e))), i = [
+};
+function vt(e, t, n) {
+	let r = s(n, 0, 1);
+	switch (t) {
+		case "watercolor":
+			e.wetEdge = r, e.granulation = r * .8;
+			break;
+		case "pastel":
+			e.strokeAmount = r, e.canvasGrain = r * .05;
+			break;
+		case "oil":
+		case "vangogh":
+			e.strokeAmount = r, e.impasto = r, e.canvasGrain = r * .06;
+			break;
+		case "oil-pastel":
+		case "crayon":
+		case "kuwahara":
+			e.strokeAmount = r, e.canvasGrain = r * .05;
+			break;
+	}
+}
+function yt(e, t) {
+	switch (t) {
+		case "watercolor": return s(e.wetEdge, 0, 1);
+		case "oil":
+		case "vangogh":
+		case "pastel":
+		case "oil-pastel":
+		case "crayon":
+		case "kuwahara": return s(e.strokeAmount, 0, 1);
+	}
+}
+function bt(e) {
+	return e.nucleiDrift + e.paletteDrift + e.warpDrift > .001 ? "drifting" : e.breathDepth > .001 ? "breathing" : "still";
+}
+function xt(e, t, n) {
+	return n === t ? 0 : s((e - t) / (n - t), 0, 1);
+}
+//#endregion
+//#region src/components/custom/aurora/composables/atoms.ts
+var $ = {}, St = .61803398875;
+function Ct(e, t = "composed") {
+	let n = Math.max(1, Math.min(6, Math.round(e))), r = [
 		[.33, .33],
 		[.67, .67],
 		[.67, .33],
 		[.33, .67],
 		[.5, .5],
 		[.5, .33]
-	], a = (e) => {
+	], i = (e) => {
 		switch (t) {
 			case "scattered": {
-				let t = e * Ge * Math.PI * 2, i = .12 + .3 * Math.sqrt((e + .5) / n);
-				return [r(.5 + i * Math.cos(t), .1, .9), r(.5 + i * Math.sin(t), .1, .9)];
+				let t = e * St * Math.PI * 2, r = .12 + .3 * Math.sqrt((e + .5) / n);
+				return [s(.5 + r * Math.cos(t), .1, .9), s(.5 + r * Math.sin(t), .1, .9)];
 			}
 			case "centred": {
 				if (n === 1) return [.5, .5];
 				let t = e / n * Math.PI * 2, r = .12;
 				return [.5 + r * Math.cos(t), .5 + r * Math.sin(t)];
 			}
-			default: return i[e % i.length];
+			default: return r[e % r.length];
 		}
 	};
 	return Array.from({ length: n }, (e, t) => {
-		let [r, i] = a(t);
+		let [r, a] = i(t);
 		return {
 			x: r,
-			y: i,
+			y: a,
 			radius: .5,
 			paletteBias: n === 1 ? 0 : t / (n - 1),
 			valueBias: (t % 2 == 0 ? 1 : -1) * .05,
@@ -887,82 +1580,47 @@ function Q(e, t = "composed") {
 		};
 	});
 }
-function Ke(e) {
+function wt(e) {
 	return e < .4 ? "fbm" : e < .75 ? "hybrid" : "cellular";
 }
-function qe(e, t) {
-	let n = r(t, 0, 1);
-	e.warpAmount = r(X(.2, .6, n), 0, .6), e.warpScale = r(X(1, 2.6, n), .5, 3), e.warpMode = Ke(n), e.noiseOctaves = n < .5 ? 3 : n < .85 ? 4 : 5;
+function Tt(e, t) {
+	let n = s(t, 0, 1);
+	e.warpAmount = s(Z(.2, .6, n), 0, .6), e.warpScale = s(Z(1, 2.6, n), .5, 3), e.warpMode = wt(n), e.noiseOctaves = n < .5 ? 3 : n < .85 ? 4 : 5;
 }
-function Je(e, t, n) {
-	let i = r(n, 0, 1);
-	switch (t) {
-		case "watercolor":
-			e.wetEdge = i, e.granulation = i * .8;
-			break;
-		case "pastel":
-			e.strokeAmount = i, e.canvasGrain = i * .05;
-			break;
-		case "oil":
-		case "vangogh":
-			e.strokeAmount = i, e.impasto = i, e.canvasGrain = i * .06;
-			break;
-		case "oil-pastel":
-		case "crayon":
-			e.strokeAmount = i, e.canvasGrain = i * .05;
-			break;
-	}
-}
-function Ye(e = Y, t = p) {
+function Et(e = $, t = y) {
 	let n = {
 		...t,
 		palette: t.palette.map((e) => ({ ...e })),
 		nuclei: t.nuclei.map((e) => ({ ...e })),
 		flow: { ...t.flow }
-	}, i = e.colorEnergy === void 0 ? void 0 : r(e.colorEnergy, 0, 1);
-	if (e.seed !== void 0 && (n.palette = R(e.seed, {
+	}, r = e.colorEnergy === void 0 ? void 0 : s(e.colorEnergy, 0, 1);
+	if (e.seed !== void 0 && (n.palette = P(e.seed, {
 		harmony: e.harmony ?? "analogous",
 		stopCount: 4,
-		temperatureShift: i === void 0 ? 0 : X(Z.temperatureShift.calm, Z.temperatureShift.vivid, i)
-	})), i !== void 0 && (n.saturation = r(X(Z.saturation.calm, Z.saturation.vivid, i), .6, 1.3), n.valueVariance = r(X(Z.valueVariance.calm, Z.valueVariance.vivid, i), 0, .3), n.breathDepth = r(X(Z.breathDepth.calm, Z.breathDepth.vivid, i), 0, .15)), e.zones !== void 0 && (n.nuclei = Q(e.zones.count, e.zones.arrangement ?? "composed")), e.noise !== void 0 && qe(n, e.noise), e.medium !== void 0 && (n.medium = e.medium.kind, (e.medium.kind === "vangogh" || e.medium.kind === "oil-pastel") && (n.strokeOrient = "tensor"), e.medium.kind !== "smooth" && e.medium.amount !== void 0 && Je(n, e.medium.kind, e.medium.amount)), e.motion !== void 0) {
-		let t = We[e.motion];
-		n.nucleiDrift = r(t.nucleiDrift, 0, .05), n.paletteDrift = r(t.paletteDrift, 0, .04), n.warpDrift = r(t.warpDrift, 0, .015), n.breathDepth = r(e.motion === "still" ? 0 : Math.max(n.breathDepth, t.breathDepth), 0, .15);
+		temperatureShift: r === void 0 ? 0 : Z(Q.temperatureShift.calm, Q.temperatureShift.vivid, r)
+	})), r !== void 0 && (n.saturation = s(Z(Q.saturation.calm, Q.saturation.vivid, r), .6, 1.3), n.valueVariance = s(Z(Q.valueVariance.calm, Q.valueVariance.vivid, r), 0, .3), n.breathDepth = s(Z(Q.breathDepth.calm, Q.breathDepth.vivid, r), 0, .15)), e.zones !== void 0 && (n.nuclei = Ct(e.zones.count, e.zones.arrangement ?? "composed")), e.noise !== void 0 && Tt(n, e.noise), e.medium !== void 0 && (n.medium = e.medium.kind, (e.medium.kind === "vangogh" || e.medium.kind === "oil-pastel") && (n.strokeOrient = "tensor"), e.medium.kind !== "smooth" && e.medium.amount !== void 0 && vt(n, e.medium.kind, e.medium.amount)), e.motion !== void 0) {
+		let t = _t[e.motion];
+		n.nucleiDrift = s(t.nucleiDrift, 0, .05), n.paletteDrift = s(t.paletteDrift, 0, .04), n.warpDrift = s(t.warpDrift, 0, .015), n.breathDepth = s(e.motion === "still" ? 0 : Math.max(n.breathDepth, t.breathDepth), 0, .15);
 	}
 	return e.interactivity !== void 0 && (n.interactivity = { ...e.interactivity }), n;
 }
-function $(e, t, n) {
-	return n === t ? 0 : r((e - t) / (n - t), 0, 1);
-}
-function Xe(e) {
-	return e.nucleiDrift + e.paletteDrift + e.warpDrift > .001 ? "drifting" : e.breathDepth > .001 ? "breathing" : "still";
-}
-function Ze(e, t) {
-	switch (t) {
-		case "watercolor": return r(e.wetEdge, 0, 1);
-		case "oil":
-		case "vangogh":
-		case "pastel":
-		case "oil-pastel":
-		case "crayon": return r(e.strokeAmount, 0, 1);
-	}
-}
-function Qe(e) {
+function Dt(e) {
 	let t = e.medium, n = t === "smooth" ? { kind: t } : {
 		kind: t,
-		amount: Ze(e, t)
+		amount: yt(e, t)
 	};
 	return {
-		seed: e.palette.length > 0 ? v(e.palette[0]) : "#3a7bd5",
+		seed: e.palette.length > 0 ? x(e.palette[0]) : "#3a7bd5",
 		harmony: "analogous",
-		colorEnergy: $(e.saturation, Z.saturation.calm, Z.saturation.vivid),
+		colorEnergy: xt(e.saturation, Q.saturation.calm, Q.saturation.vivid),
 		zones: {
 			count: Math.max(1, Math.min(6, e.nuclei.length)),
 			arrangement: "composed"
 		},
-		noise: $(e.warpAmount, .2, .6),
+		noise: xt(e.warpAmount, .2, .6),
 		medium: n,
-		motion: Xe(e)
+		motion: bt(e)
 	};
 }
 //#endregion
-export { He as Aurora, Y as DEFAULT_ATOMS, p as DEFAULT_AURORA_CONFIG, d as MAX_NUCLEI, f as MAX_STOPS, Qe as configToAtoms, G as createAurora, y as cssToOklch, R as deriveAurora, P as flattenPalette, I as hexToOklchStop, Q as nucleiPrior, v as oklchStopToHex, h as oklchToLinear, F as paletteToCssGradient, Ye as resolveAtoms, J as resolveRenderMode, q as useAurora, Ue as useCursorInteraction };
+export { ht as Aurora, $ as DEFAULT_ATOMS, y as DEFAULT_AURORA_CONFIG, _ as MAX_NUCLEI, v as MAX_STOPS, g as PAPER_WASH_GROUND, ft as auroraFallbackGround, Dt as configToAtoms, q as createAurora, S as cssToOklch, P as deriveAurora, M as flattenPalette, ke as hexToOklchStop, A as isSoftwareWebGLRenderer, ct as nucleiFieldStatic, Ct as nucleiPrior, x as oklchStopToHex, b as oklchToLinear, N as paletteToCssGradient, ot as relativeLuminance, Et as resolveAtoms, Ee as resolveRenderMode, lt as sampleAuroraField, Y as useAurora, gt as useCursorInteraction };