reze-engine 0.11.1 → 0.11.3

package/src/shaders/{metal.ts → materials/metal.ts}

@@ -2,6 +2,12 @@
 // + NPR toon/AO emission stack (Strength=8.1), MixShader Fac=0.6967.
 // Base color uses a Voronoi pattern sampled in reflection-coord space (Blender 纹理坐标.Reflection)
 // to add subtle metallic sparkle variation. No Normal link in the graph.
+//
+// Graph's base color chain is: 纹理坐标.Reflection → 矢量运算.007(CROSS, Vec2=(0,1,0)) →
+// 沃罗诺伊纹理(F1, Color out) → 颜色渐变(linear) → 混合.005. The dumper did not capture the
+// VectorMath operation — CROSS is the assumed op based on the hardcoded (0,1,0) Vector_001
+// constant (MULTIPLY would zero X/Z producing 1D bands; CROSS produces horizontal ring
+// patterns consistent with metallic anisotropy).
 
 import { NODES_WGSL } from "./nodes"
 
@@ -129,40 +135,30 @@ struct FSOut {
   let out_alpha = material.alpha * tex_s.a;
   if (out_alpha < 0.001) { discard; }
 
-  // ═══ NPR toon stack (图像 → HSV.007 Val=0.8 → 转接点.001) ═══
   let tex_tint = hue_sat_id(1.0, 0.800000011920929, 1.0, tex_rgb);
   let lum_shade = shader_to_rgb_diffuse(n, l, sun, amb, shadow);
   let ramp008 = ramp_constant_edge_aa(lum_shade, METAL_TOON_EDGE, vec4f(0,0,0,1), vec4f(1,1,1,1));
   let mix04_fac = math_multiply(ramp008.r, METAL_MIX04_MUL);
 
-  // 混合.004: A=HSV.002(Val=0.2 dark), B=tex_tint
   let dark_tex = hue_sat_id(1.0, 0.19999998807907104, 1.0, tex_tint);
   let mix04 = mix_blend(mix04_fac, dark_tex, tex_tint);
 
-  // AO white/black ramp → 混合.002 factor
-  let ao = 1.0; // ao_fake(n, v) — no SSAO yet; inline 1.0 so the ramp/mix chain folds at compile time.
-  let ao_ramp_c = ramp_linear(ao, 0.0, vec4f(1,1,1,1), 0.8808, vec4f(0,0,0,1));
-  let overlay_fac = mix(1.0, 0.0, ao_ramp_c.r);
-
-  // 混合.002 OVERLAY: A=HSV.008(Val=1.0 identity) ← mix04, B=HSV.004(Val=2.0 bright) ← mix04
-  let hue008 = mix04; // identity HSV
   let hue004 = hue_sat_id(1.0, 2.0, 1.0, mix04);
-  let npr_rgb = mix_overlay(overlay_fac, hue008, hue004);
+  let npr_rgb = mix_overlay(1.0, mix04, hue004);
   let npr_emission = npr_rgb * METAL_EMIT_STR;
 
-  // ═══ Metallic Principled base color ═══
-  // Reflection-coord Voronoi for metallic sparkle:
-  //   纹理坐标.Reflection → 矢量运算 → 沃罗诺伊(Scale=4.3) → 颜色渐变 → 混合.005
+  // Reflection-coord Voronoi produces the metallic sparkle variation.
+  // VALTORGB takes Color → Fac via Blender's BT.601 implicit color_to_value.
   let refl_dir = reflect(-v, n);
-  let voro = tex_voronoi_f1(refl_dir, METAL_VORONOI_SCALE);
-  let voro_ramp = ramp_linear(voro, 0.0, vec4f(0,0,0,1), 1.0, vec4f(1,1,1,1)).r;
-  // 混合.005: Fac=voro_ramp, A=voro_color(grayscale), B=HSV.006(Hue=0.5 Sat=1.5 Val=1.3)
+  let voro_input = cross(refl_dir, vec3f(0.0, 1.0, 0.0));
+  let voro_rgb = tex_voronoi_color(voro_input, METAL_VORONOI_SCALE);
+  let voro_scalar = color_to_value(voro_rgb);
+  let voro_ramp = ramp_linear(voro_scalar, 0.0, vec4f(0,0,0,1), 1.0, vec4f(1,1,1,1)).r;
   let hue006 = hue_sat_id(1.5, 1.2999999523162842, 1.0, tex_tint);
   let albedo = mix_blend(voro_ramp, vec3f(voro_ramp), hue006);
 
-  // 原理化BSDF (EEVEE port): metallic=1.0, specular=1.0, roughness=0.3.
-  // Per Blender principled wiring: f0 = mix((0.08*spec)*dielectric, albedo, metallic) → with
-  // metallic=1 this is just albedo (specular_tint is dielectric-only and ignored here).
+  // Principled BSDF (EEVEE port): metallic=1 collapses f0 = mix(dielectric, albedo, 1) = albedo;
+  // specular_tint is dielectric-only and ignored here.
   let f0 = albedo;
   let f90 = mix(f0, vec3f(1.0), sqrt(METAL_SPECULAR));
   let NL = max(dot(n, l), 0.0);

package/src/shaders/materials/nodes.ts

@@ -0,0 +1,483 @@
+// Shared WGSL primitives for Blender-style NPR material nodes.
+// Every function here maps 1:1 to a Blender shader node type used in the preset JSONs.
+// Hand-ported material shaders concatenate this block before their own code.
+
+export const NODES_WGSL = /* wgsl */ `
+
+// Baked 64×64 rgba8unorm combined BRDF LUT — created once at engine init by dfg_lut.ts.
+//   .rg = split-sum DFG (Karis: tint = f0·x + f90·y)  → F_brdf_*_scatter
+//   .ba = Heitz 2016 LTC magnitude (ltc_mag_ggx)       → ltc_brdf_scale_from_lut
+// Paired with group(0) binding(2) diffuseSampler (linear filter). Sample once per
+// fragment via brdf_lut_sample() — callers feed .rg and the whole vec4 into the
+// helpers below, halving LUT taps on the default Principled path.
+@group(0) @binding(9) var brdfLut: texture_2d<f32>;
+
+// ─── RGB ↔ HSV ──────────────────────────────────────────────────────
+
+fn rgb_to_hsv(rgb: vec3f) -> vec3f {
+  let c_max = max(rgb.r, max(rgb.g, rgb.b));
+  let c_min = min(rgb.r, min(rgb.g, rgb.b));
+  let delta = c_max - c_min;
+
+  var h = 0.0;
+  if (delta > 1e-6) {
+    if (c_max == rgb.r) {
+      h = (rgb.g - rgb.b) / delta;
+      if (h < 0.0) { h += 6.0; }
+    } else if (c_max == rgb.g) {
+      h = 2.0 + (rgb.b - rgb.r) / delta;
+    } else {
+      h = 4.0 + (rgb.r - rgb.g) / delta;
+    }
+    h /= 6.0;
+  }
+  let s = select(0.0, delta / c_max, c_max > 1e-6);
+  return vec3f(h, s, c_max);
+}
+
+fn hsv_to_rgb(hsv: vec3f) -> vec3f {
+  let h = hsv.x;
+  let s = hsv.y;
+  let v = hsv.z;
+  if (s < 1e-6) { return vec3f(v); }
+
+  let hh = fract(h) * 6.0;
+  let sector = u32(hh);
+  let f = hh - f32(sector);
+  let p = v * (1.0 - s);
+  let q = v * (1.0 - s * f);
+  let t = v * (1.0 - s * (1.0 - f));
+
+  switch (sector) {
+    case 0u: { return vec3f(v, t, p); }
+    case 1u: { return vec3f(q, v, p); }
+    case 2u: { return vec3f(p, v, t); }
+    case 3u: { return vec3f(p, q, v); }
+    case 4u: { return vec3f(t, p, v); }
+    default: { return vec3f(v, p, q); }
+  }
+}
+
+// ─── HUE_SAT node ───────────────────────────────────────────────────
+
+fn hue_sat(hue: f32, saturation: f32, value: f32, fac: f32, color: vec3f) -> vec3f {
+  var hsv = rgb_to_hsv(color);
+  hsv.x = fract(hsv.x + hue - 0.5);
+  hsv.y = clamp(hsv.y * saturation, 0.0, 1.0);
+  hsv.z *= value;
+  return mix(color, hsv_to_rgb(hsv), fac);
+}
+
+// hue_sat specialization for hue=0.5 (identity hue shift — fract(h + 0.5 - 0.5) = h).
+// Branchless equivalent that skips the rgb_to_hsv → hsv_to_rgb roundtrip: WebKit's
+// Metal backend serializes the 3-way if chain in rgb_to_hsv and the 6-way switch in
+// hsv_to_rgb, where this form compiles to linear SIMD ops + a single select.
+fn hue_sat_id(saturation: f32, value: f32, fac: f32, color: vec3f) -> vec3f {
+  let m = max(max(color.r, color.g), color.b);
+  let n = min(min(color.r, color.g), color.b);
+  // Unclamped (sat*old_s ≤ 1): reproj = mix(vec3f(m), color, saturation).
+  // Clamped (saturated to 1):   reproj = (color - n) * m / (m - n).
+  let range = max(m - n, 1e-6);
+  let unclamped = mix(vec3f(m), color, saturation);
+  let clamped = (color - vec3f(n)) * m / range;
+  let needs_clamp = (m - n) * saturation >= m;
+  let reproj = select(unclamped, clamped, needs_clamp);
+  return mix(color, reproj * value, fac);
+}
+
+// ─── BRIGHTCONTRAST node ────────────────────────────────────────────
+
+fn bright_contrast(color: vec3f, bright: f32, contrast: f32) -> vec3f {
+  let a = 1.0 + contrast;
+  let b = bright - contrast * 0.5;
+  return max(vec3f(0.0), color * a + vec3f(b));
+}
+
+// ─── INVERT node ────────────────────────────────────────────────────
+
+fn invert(fac: f32, color: vec3f) -> vec3f {
+  return mix(color, vec3f(1.0) - color, fac);
+}
+
+fn invert_f(fac: f32, val: f32) -> f32 {
+  return mix(val, 1.0 - val, fac);
+}
+
+// ─── Color ramp (VALTORGB) — 2-stop variants ───────────────────────
+// All 7 presets use exclusively 2-stop ramps.
+
+fn ramp_constant(f: f32, p0: f32, c0: vec4f, p1: f32, c1: vec4f) -> vec4f {
+  return select(c0, c1, f >= p1);
+}
+
+// CONSTANT ramp with screen-space edge AA — kills sparkle where fwidth(f) straddles a hard step (NPR terminator)
+fn ramp_constant_edge_aa(f: f32, edge: f32, c0: vec4f, c1: vec4f) -> vec4f {
+  let w = max(fwidth(f) * 1.75, 6e-6);
+  let t = smoothstep(edge - w, edge + w, f);
+  return mix(c0, c1, t);
+}
+
+fn ramp_linear(f: f32, p0: f32, c0: vec4f, p1: f32, c1: vec4f) -> vec4f {
+  let t = saturate((f - p0) / max(p1 - p0, 1e-6));
+  return mix(c0, c1, t);
+}
+
+fn ramp_cardinal(f: f32, p0: f32, c0: vec4f, p1: f32, c1: vec4f) -> vec4f {
+  // cardinal spline with 2 stops degrades to smoothstep
+  let t = saturate((f - p0) / max(p1 - p0, 1e-6));
+  let ss = t * t * (3.0 - 2.0 * t);
+  return mix(c0, c1, ss);
+}
+
+// ─── MATH node operations ───────────────────────────────────────────
+
+fn math_add(a: f32, b: f32) -> f32 { return a + b; }
+fn math_multiply(a: f32, b: f32) -> f32 { return a * b; }
+fn math_power(a: f32, b: f32) -> f32 { return pow(max(a, 0.0), b); }
+fn math_greater_than(a: f32, b: f32) -> f32 { return select(0.0, 1.0, a > b); }
+
+// Blender's implicit Color → Float socket conversion uses BT.601 grayscale
+// (rgb_to_grayscale in blenkernel/intern/node.cc). When a material graph plugs a
+// Color output into a Math node's Value input, this is the scalar it actually sees.
+fn color_to_value(c: vec3f) -> f32 {
+  return 0.299 * c.r + 0.587 * c.g + 0.114 * c.b;
+}
+
+// ─── MIX node (blend_type variants) ────────────────────────────────
+
+fn mix_blend(fac: f32, a: vec3f, b: vec3f) -> vec3f {
+  return mix(a, b, fac);
+}
+
+fn mix_overlay(fac: f32, a: vec3f, b: vec3f) -> vec3f {
+  let lo = 2.0 * a * b;
+  let hi = vec3f(1.0) - 2.0 * (vec3f(1.0) - a) * (vec3f(1.0) - b);
+  let overlay = select(hi, lo, a < vec3f(0.5));
+  return mix(a, overlay, fac);
+}
+
+fn mix_multiply(fac: f32, a: vec3f, b: vec3f) -> vec3f {
+  return mix(a, a * b, fac);
+}
+
+fn mix_lighten(fac: f32, a: vec3f, b: vec3f) -> vec3f {
+  return mix(a, max(a, b), fac);
+}
+
+// Blender Mix (Color) blend LINEAR_LIGHT: result = mix(A, A + 2*B - 1, Fac)
+fn mix_linear_light(fac: f32, a: vec3f, b: vec3f) -> vec3f {
+  return mix(a, a + 2.0 * b - vec3f(1.0), fac);
+}
+
+// Luminance for Shader→RGB scalar gates (linear RGB, Rec.709 weights)
+fn luminance_rec709_linear(c: vec3f) -> f32 {
+  return dot(max(c, vec3f(0.0)), vec3f(0.2126, 0.7152, 0.0722));
+}
+
+// ─── FRESNEL node ───────────────────────────────────────────────────
+// Schlick approximation matching Blender's Fresnel node
+
+fn fresnel(ior: f32, n: vec3f, v: vec3f) -> f32 {
+  let r = (ior - 1.0) / (ior + 1.0);
+  let f0 = r * r;
+  let cos_theta = clamp(dot(n, v), 0.0, 1.0);
+  let m = 1.0 - cos_theta;
+  let m2 = m * m;
+  let m5 = m2 * m2 * m;
+  return f0 + (1.0 - f0) * m5;
+}
+
+// ─── LAYER_WEIGHT node ──────────────────────────────────────────────
+
+fn layer_weight_fresnel(blend: f32, n: vec3f, v: vec3f) -> f32 {
+  let eta = max(1.0 - blend, 1e-4);
+  let r = (1.0 - eta) / (1.0 + eta);
+  let f0 = r * r;
+  let cos_theta = clamp(abs(dot(n, v)), 0.0, 1.0);
+  let m = 1.0 - cos_theta;
+  let m2 = m * m;
+  let m5 = m2 * m2 * m;
+  return f0 + (1.0 - f0) * m5;
+}
+
+fn layer_weight_facing(blend: f32, n: vec3f, v: vec3f) -> f32 {
+  var facing = abs(dot(n, v));
+  let b = clamp(blend, 0.0, 0.99999);
+  if (b != 0.5) {
+    let exponent = select(2.0 * b, 0.5 / (1.0 - b), b >= 0.5);
+    facing = pow(facing, exponent);
+  }
+  return 1.0 - facing;
+}
+
+// ─── SHADER_TO_RGB (white DiffuseBSDF) ──────────────────────────────
+// Eevee captures lit diffuse: (albedo/π)*sun*N·L*shadow + ambient (linear). Albedo=1.
+// Matches default.ts direct term scale so VALTORGB thresholds from Blender JSON stay valid.
+
+fn shader_to_rgb_diffuse(n: vec3f, l: vec3f, sun_rgb: vec3f, ambient_rgb: vec3f, shadow: f32) -> f32 {
+  const PI_S: f32 = 3.141592653589793;
+  let ndotl = max(dot(n, l), 0.0);
+  let rgb = sun_rgb * (ndotl * shadow / PI_S) + ambient_rgb;
+  return luminance_rec709_linear(rgb);
+}
+
+// ─── BUMP node ──────────────────────────────────────────────────────
+// Screen-space bump from a scalar height field. Needs dFdx/dFdy which
+// WGSL provides as dpdx/dpdy.
+
+fn bump(strength: f32, height: f32, normal: vec3f, world_pos: vec3f) -> vec3f {
+  let dhdx = dpdx(height);
+  let dhdy = dpdy(height);
+  let dpdx_pos = dpdx(world_pos);
+  let dpdy_pos = dpdy(world_pos);
+  let perturbed = normalize(normal) - strength * (dhdx * normalize(cross(dpdy_pos, normal)) + dhdy * normalize(cross(normal, dpdx_pos)));
+  return normalize(perturbed);
+}
+
+// LH engine + WebGPU fragment Y: flip dhdy contribution so height peaks read as outward bumps vs Blender reference
+fn bump_lh(strength: f32, height: f32, normal: vec3f, world_pos: vec3f) -> vec3f {
+  let dhdx = dpdx(height);
+  let dhdy = dpdy(height);
+  let dpdx_pos = dpdx(world_pos);
+  let dpdy_pos = dpdy(world_pos);
+  let perturbed = normalize(normal) - strength * (dhdx * normalize(cross(dpdy_pos, normal)) - dhdy * normalize(cross(normal, dpdx_pos)));
+  return normalize(perturbed);
+}
+
+// ─── NOISE texture (Perlin-style) ───────────────────────────────────
+// Simplified gradient noise matching Blender's default noise output.
+
+// PCG-style integer hash. Replaces the classic 'fract(sin(q) * LARGE)' trick because
+// WebKit's Metal backend compiles 'sin' to a full transcendental op (slow), while
+// Safari's Apple-GPU scalar ALU handles int muls/xors near free. Inputs arrive as
+// integer-valued floats (floor(p) + unit offsets) from _noise3, so vec3i cast is exact.
+fn _hash33(p: vec3f) -> vec3f {
+  var h = vec3u(vec3i(p) + vec3i(32768));
+  h = h * vec3u(1664525u, 1013904223u, 2654435761u);
+  h = (h.yzx ^ h) * vec3u(2246822519u, 3266489917u, 668265263u);
+  h = h ^ (h >> vec3u(16u));
+  // Mask to 24 bits — above that f32 loses precision on the u32→f32 convert.
+  let hm = h & vec3u(16777215u);
+  return vec3f(hm) * (2.0 / 16777216.0) - 1.0;
+}
+
+fn _noise3(p: vec3f) -> f32 {
+  let i = floor(p);
+  let f = fract(p);
+  let u = f * f * (3.0 - 2.0 * f);
+
+  return mix(
+    mix(
+      mix(dot(_hash33(i + vec3f(0,0,0)), f - vec3f(0,0,0)),
+          dot(_hash33(i + vec3f(1,0,0)), f - vec3f(1,0,0)), u.x),
+      mix(dot(_hash33(i + vec3f(0,1,0)), f - vec3f(0,1,0)),
+          dot(_hash33(i + vec3f(1,1,0)), f - vec3f(1,1,0)), u.x), u.y),
+    mix(
+      mix(dot(_hash33(i + vec3f(0,0,1)), f - vec3f(0,0,1)),
+          dot(_hash33(i + vec3f(1,0,1)), f - vec3f(1,0,1)), u.x),
+      mix(dot(_hash33(i + vec3f(0,1,1)), f - vec3f(0,1,1)),
+          dot(_hash33(i + vec3f(1,1,1)), f - vec3f(1,1,1)), u.x), u.y),
+    u.z);
+}
+
+fn tex_noise(p: vec3f, scale: f32, detail: f32, roughness: f32, distortion: f32) -> f32 {
+  var q = p;
+  if (abs(distortion) > 1e-6) {
+    let w = _noise3(p * scale * 1.37 + vec3f(2.31, 5.17, 8.09));
+    q = p + (w * 2.0 - 1.0) * distortion;
+  }
+  let coords = q * scale;
+  var value = 0.0;
+  var amplitude = 1.0;
+  var frequency = 1.0;
+  var total_amp = 0.0;
+  let octaves = i32(clamp(detail, 0.0, 15.0)) + 1;
+  for (var i = 0; i < octaves; i++) {
+    value += amplitude * _noise3(coords * frequency);
+    total_amp += amplitude;
+    amplitude *= roughness;
+    frequency *= 2.0;
+  }
+  return value / max(total_amp, 1e-6) * 0.5 + 0.5;
+}
+
+// tex_noise specialization: detail=2.0 (3 octaves), roughness=0.5, distortion=0.
+// WebKit can't unroll tex_noise's for-loop because 'octaves' is a runtime value;
+// this variant is fully unrolled with constants folded (total_amp = 1.75).
+fn tex_noise_d2(p: vec3f, scale: f32) -> f32 {
+  let c = p * scale;
+  let v = _noise3(c) + 0.5 * _noise3(c * 2.0) + 0.25 * _noise3(c * 4.0);
+  return v * (1.0 / 1.75) * 0.5 + 0.5;
+}
+
+// ─── TEX_GRADIENT (linear) ──────────────────────────────────────────
+// Used by Stockings preset. Maps the input vector's X to a 0–1 gradient.
+
+fn tex_gradient_linear(uv: vec3f) -> f32 {
+  return clamp(uv.x, 0.0, 1.0);
+}
+
+// ─── TEX_VORONOI ────────────────────────────────────────────────────
+// 3D F1 voronoi. _f1 returns Distance; _color returns per-cell hash color
+// (matches Blender voronoi.cc: outColor = hash_int3_to_float3(cell + targetOffset)).
+
+fn tex_voronoi_f1(p: vec3f, scale: f32) -> f32 {
+  let coords = p * scale;
+  let i = floor(coords);
+  let f = fract(coords);
+  var min_dist = 1e10;
+  for (var z = -1; z <= 1; z++) {
+    for (var y = -1; y <= 1; y++) {
+      for (var x = -1; x <= 1; x++) {
+        let neighbor = vec3f(f32(x), f32(y), f32(z));
+        let point = _hash33(i + neighbor) * 0.5 + 0.5;
+        let diff = neighbor + point - f;
+        min_dist = min(min_dist, dot(diff, diff));
+      }
+    }
+  }
+  return sqrt(min_dist);
+}
+
+// The per-cell jitter hash IS the Color output in Blender — reuse the same hash
+// tap for jitter + color instead of computing two.
+fn tex_voronoi_color(p: vec3f, scale: f32) -> vec3f {
+  let coords = p * scale;
+  let i = floor(coords);
+  let f = fract(coords);
+  var min_dist = 1e10;
+  var min_hash = vec3f(0.5);
+  for (var z = -1; z <= 1; z++) {
+    for (var y = -1; y <= 1; y++) {
+      for (var x = -1; x <= 1; x++) {
+        let neighbor = vec3f(f32(x), f32(y), f32(z));
+        let h = _hash33(i + neighbor) * 0.5 + 0.5;
+        let diff = neighbor + h - f;
+        let d = dot(diff, diff);
+        if (d < min_dist) {
+          min_dist = d;
+          min_hash = h;
+        }
+      }
+    }
+  }
+  return min_hash;
+}
+
+// ─── SEPXYZ node ────────────────────────────────────────────────────
+
+fn separate_xyz(v: vec3f) -> vec3f { return v; }
+
+// ─── VECT_MATH (cross product) ──────────────────────────────────────
+
+fn vect_math_cross(a: vec3f, b: vec3f) -> vec3f { return cross(a, b); }
+
+// ─── MAPPING node ───────────────────────────────────────────────────
+// Point-type mapping: scale, rotate (euler XYZ), translate.
+
+fn mapping_point(v: vec3f, loc: vec3f, rot: vec3f, scl: vec3f) -> vec3f {
+  var p = v * scl;
+  // simplified: skip rotation when all angles are zero (common case)
+  if (abs(rot.x) + abs(rot.y) + abs(rot.z) > 1e-6) {
+    let cx = cos(rot.x); let sx = sin(rot.x);
+    let cy = cos(rot.y); let sy = sin(rot.y);
+    let cz = cos(rot.z); let sz = sin(rot.z);
+    let rx = vec3f(p.x, cx*p.y - sx*p.z, sx*p.y + cx*p.z);
+    let ry = vec3f(cy*rx.x + sy*rx.z, rx.y, -sy*rx.x + cy*rx.z);
+    p = vec3f(cz*ry.x - sz*ry.y, sz*ry.x + cz*ry.y, ry.z);
+  }
+  return p + loc;
+}
+
+// ─── NORMAL_MAP node (tangent-space) ────────────────────────────────
+// Applies a tangent-space normal map. Requires TBN from vertex stage.
+
+fn normal_map(strength: f32, map_color: vec3f, normal: vec3f, tangent: vec3f, bitangent: vec3f) -> vec3f {
+  let ts = map_color * 2.0 - 1.0;
+  let perturbed = normalize(tangent * ts.x + bitangent * ts.y + normal * ts.z);
+  return normalize(mix(normal, perturbed, strength));
+}
+
+// ─── EEVEE Principled BSDF primitives ───────────────────────────────
+// Ports from Blender 3.6 source/blender/draw/engines/eevee/shaders/
+//   bsdf_common_lib.glsl + gpu_shader_material_principled.glsl.
+// Usage pattern (see material shaders): direct spec = bsdf_ggx × sun × shadow
+// (NL baked in, no F yet); ambient spec = probe_radiance; tint both with
+// reflection_color = F_brdf_multi_scatter(f0, f90, split_sum) AFTER summing.
+
+const EEVEE_PI: f32 = 3.141592653589793;
+
+// Fused analytic GGX specular (direct lights). Returns BRDF × NL.
+// 4·NL·NV is cancelled via G1_Smith reciprocal form — see bsdf_common_lib.glsl:115.
+// Caller passes NL, NV (already computed for diffuse + brdf_lut_sample) so WebKit
+// can reuse them instead of recomputing dot products across the function boundary.
+fn bsdf_ggx(N: vec3f, L: vec3f, V: vec3f, NL_in: f32, NV_in: f32, roughness: f32) -> f32 {
+  let a = max(roughness, 1e-4);
+  let a2 = a * a;
+  let H = normalize(L + V);
+  let NH = max(dot(N, H), 1e-8);
+  let NL = max(NL_in, 1e-8);
+  let NV = max(NV_in, 1e-8);
+  // G1_Smith_GGX_opti reciprocal form — denominator piece only.
+  let G1L = NL + sqrt(NL * (NL - NL * a2) + a2);
+  let G1V = NV + sqrt(NV * (NV - NV * a2) + a2);
+  let G = G1L * G1V;
+  // D_ggx_opti = pi * denom² — reciprocal D × a².
+  let tmp = (NH * a2 - NH) * NH + 1.0;
+  let D_opti = EEVEE_PI * tmp * tmp;
+  return NL * a2 / (D_opti * G);
+}
+
+// Split-sum DFG LUT — Karis 2013 curve fit stand-in for the 64×64 baked LUT.
+// Returns (lut.x, lut.y) in Blender convention: tint = f0·lut.x + f90·lut.y.
+fn brdf_lut_approx(NV: f32, roughness: f32) -> vec2f {
+  let c0 = vec4f(-1.0, -0.0275, -0.572, 0.022);
+  let c1 = vec4f(1.0, 0.0425, 1.04, -0.04);
+  let r = roughness * c0 + c1;
+  let a004 = min(r.x * r.x, exp2(-9.28 * NV)) * r.x + r.y;
+  return vec2f(-1.04, 1.04) * a004 + r.zw;
+}
+
+// Baked combined BRDF LUT — exact port of Blender bsdf_lut_frag.glsl packed with
+// ltc_mag_ggx from eevee_lut.c. Single sample returns DFG (.rg) and LTC mag (.ba).
+// Addressed as Blender's common_utiltex_lib.glsl:lut_coords:
+//   coords = (roughness, sqrt(1 - NV)), then half-texel bias for filtering.
+// Requires group(0) binding(9) brdfLut + binding(2) diffuseSampler in the host shader.
+fn brdf_lut_sample(NV: f32, roughness: f32) -> vec4f {
+  let LUT_SIZE: f32 = 64.0;
+  var uv = vec2f(saturate(roughness), sqrt(saturate(1.0 - NV)));
+  uv = uv * ((LUT_SIZE - 1.0) / LUT_SIZE) + 0.5 / LUT_SIZE;
+  return textureSampleLevel(brdfLut, diffuseSampler, uv, 0.0);
+}
+
+fn F_brdf_single_scatter(f0: vec3f, f90: vec3f, lut: vec2f) -> vec3f {
+  return lut.y * f90 + lut.x * f0;
+}
+
+// Fdez-Agüera 2019 multi-scatter compensation (EEVEE do_multiscatter=1).
+fn F_brdf_multi_scatter(f0: vec3f, f90: vec3f, lut: vec2f) -> vec3f {
+  let FssEss = lut.y * f90 + lut.x * f0;
+  let Ess = lut.x + lut.y;
+  let Ems = 1.0 - Ess;
+  let Favg = f0 + (1.0 - f0) / 21.0;
+  let Fms = FssEss * Favg / (1.0 - (1.0 - Ess) * Favg);
+  return FssEss + Fms * Ems;
+}
+
+// EEVEE direct-specular energy compensation factor — closure_eval_glossy_lib.glsl:79-81:
+//   ltc_brdf_scale = (ltc.x + ltc.y) / (split_sum.x + split_sum.y)
+// Blender evaluates direct lights via LTC (Heitz 2016) but indirect via split-sum;
+// direct radiance is rescaled so total-energy matches the split-sum LUT.
+// Takes a pre-sampled vec4f from brdf_lut_sample() to share the fetch with
+// F_brdf_multi_scatter on the same fragment.
+fn ltc_brdf_scale_from_lut(lut: vec4f) -> f32 {
+  return (lut.z + lut.w) / max(lut.x + lut.y, 1e-6);
+}
+
+// Luminance-normalized hue extraction — Blender tint_from_color (isolates hue+sat).
+fn tint_from_color(color: vec3f) -> vec3f {
+  let lum = dot(color, vec3f(0.3, 0.6, 0.1));
+  return select(vec3f(1.0), color / lum, lum > 0.0);
+}
+
+`;

package/src/shaders/passes/bloom.ts

@@ -0,0 +1,121 @@
+// EEVEE 3.6 bloom pyramid: blit (Karis prefilter) → 13-tap downsamples → 9-tap tent upsamples.
+// Mirrors source/blender/draw/engines/eevee/shaders/effect_bloom_frag.glsl. Firefly suppression
+// lives in the blit (Karis luminance-weighted 4-tap average). A single-pass Gaussian cannot
+// reproduce this — hot pixels dominate and produce the sparkle halo.
+
+const FULLSCREEN_VS = /* wgsl */ `
+@vertex fn vs(@builtin(vertex_index) vi: u32) -> @builtin(position) vec4f {
+  let x = f32((vi & 1u) << 2u) - 1.0;
+  let y = f32((vi & 2u) << 1u) - 1.0;
+  return vec4f(x, y, 0.0, 1.0);
+}
+`
+
+// Full-res HDR → half-res. Karis 4-tap firefly average + EEVEE quadratic knee threshold + clamp.
+export const BLOOM_BLIT_SHADER_WGSL = `${FULLSCREEN_VS}
+@group(0) @binding(0) var hdrTex: texture_2d<f32>;
+@group(0) @binding(1) var<uniform> prefilter: vec4<f32>; // threshold, knee, clamp, _unused
+@group(0) @binding(2) var maskTex: texture_2d<f32>;
+
+fn luminance(c: vec3f) -> f32 {
+  return dot(max(c, vec3f(0.0)), vec3f(0.2126, 0.7152, 0.0722));
+}
+fn fetch(c: vec2<i32>, clampV: f32) -> vec3f {
+  let d = vec2<i32>(textureDimensions(hdrTex));
+  let cc = clamp(c, vec2<i32>(0), d - vec2<i32>(1));
+  let s = textureLoad(hdrTex, cc, 0);
+  // Scene pass uses src-alpha blend with clear alpha 0 → premultiplied. Unpremultiply.
+  let rgb = max(s.rgb / max(s.a, 1e-6), vec3f(0.0));
+  // Bloom mask: MRT r8unorm written by material shaders (1.0 = bloom, 0.0 = skip).
+  let mask = textureLoad(maskTex, cc, 0).r;
+  let masked = rgb * mask;
+  // Blender clamps each tap BEFORE Karis average (eevee_bloom: color = min(clampIntensity, color)).
+  return select(masked, min(masked, vec3f(clampV)), clampV > 0.0);
+}
+
+@fragment fn fs(@builtin(position) p: vec4f) -> @location(0) vec4f {
+  let dst = vec2<i32>(p.xy - vec2f(0.5));
+  let base = dst * 2;
+  let clampV = prefilter.z;
+  let a = fetch(base + vec2<i32>(0, 0), clampV);
+  let b = fetch(base + vec2<i32>(1, 0), clampV);
+  let c = fetch(base + vec2<i32>(0, 1), clampV);
+  let d = fetch(base + vec2<i32>(1, 1), clampV);
+  // Karis partial average: weight each tap by 1/(1+luma) — suppresses fireflies.
+  let wa = 1.0 / (1.0 + luminance(a));
+  let wb = 1.0 / (1.0 + luminance(b));
+  let wc = 1.0 / (1.0 + luminance(c));
+  let wd = 1.0 / (1.0 + luminance(d));
+  let avg = (a * wa + b * wb + c * wc + d * wd) / max(wa + wb + wc + wd, 1e-6);
+  // EEVEE quadratic threshold (brightness = max-channel, then soft-knee curve).
+  let bright = max(avg.r, max(avg.g, avg.b));
+  let soft = clamp(bright - prefilter.x + prefilter.y, 0.0, 2.0 * prefilter.y);
+  let q = (soft * soft) / (4.0 * max(prefilter.y, 1e-4) + 1e-6);
+  let contrib = max(q, bright - prefilter.x) / max(bright, 1e-4);
+  return vec4f(max(avg * contrib, vec3f(0.0)), 1.0);
+}
+`
+
+// Jimenez/COD 13-tap dual-box — 5 weighted 2×2 averages, rejects nyquist ringing.
+export const BLOOM_DOWNSAMPLE_SHADER_WGSL = `${FULLSCREEN_VS}
+@group(0) @binding(0) var srcTex: texture_2d<f32>;
+@group(0) @binding(1) var srcSamp: sampler;
+
+fn samp(uv: vec2f, off: vec2f) -> vec3f {
+  return textureSampleLevel(srcTex, srcSamp, uv + off, 0.0).rgb;
+}
+
+@fragment fn fs(@builtin(position) p: vec4f) -> @location(0) vec4f {
+  let srcDims = vec2f(textureDimensions(srcTex));
+  let t = 1.0 / srcDims;
+  // fragCoord.xy reports pixel centers (e.g. 0.5,0.5 for first pixel) — divide by dst dims directly.
+  let dstDims = srcDims * 0.5;
+  let uv = p.xy / max(dstDims, vec2f(1.0));
+  let A = samp(uv, t * vec2f(-2.0, -2.0));
+  let B = samp(uv, t * vec2f( 0.0, -2.0));
+  let C = samp(uv, t * vec2f( 2.0, -2.0));
+  let D = samp(uv, t * vec2f(-1.0, -1.0));
+  let E = samp(uv, t * vec2f( 1.0, -1.0));
+  let F = samp(uv, t * vec2f(-2.0,  0.0));
+  let G = samp(uv, t * vec2f( 0.0,  0.0));
+  let H = samp(uv, t * vec2f( 2.0,  0.0));
+  let I = samp(uv, t * vec2f(-1.0,  1.0));
+  let J = samp(uv, t * vec2f( 1.0,  1.0));
+  let K = samp(uv, t * vec2f(-2.0,  2.0));
+  let L = samp(uv, t * vec2f( 0.0,  2.0));
+  let M = samp(uv, t * vec2f( 2.0,  2.0));
+  var o = (D + E + I + J) * (0.5 / 4.0);
+  o = o + (A + B + G + F) * (0.125 / 4.0);
+  o = o + (B + C + H + G) * (0.125 / 4.0);
+  o = o + (F + G + L + K) * (0.125 / 4.0);
+  o = o + (G + H + M + L) * (0.125 / 4.0);
+  return vec4f(o, 1.0);
+}
+`
+
+// 9-tap tent, progressively added to matching downsample mip. Blender radius = sample scale.
+export const BLOOM_UPSAMPLE_SHADER_WGSL = `${FULLSCREEN_VS}
+@group(0) @binding(0) var srcTex: texture_2d<f32>;   // coarser accumulator
+@group(0) @binding(1) var baseTex: texture_2d<f32>;  // matching downsample mip
+@group(0) @binding(2) var srcSamp: sampler;
+@group(0) @binding(3) var<uniform> upU: vec4<f32>;   // sampleScale, _, _, _
+
+@fragment fn fs(@builtin(position) p: vec4f) -> @location(0) vec4f {
+  let srcDims = vec2f(textureDimensions(srcTex));
+  let baseDims = vec2f(textureDimensions(baseTex));
+  let uv = p.xy / max(baseDims, vec2f(1.0));
+  let t = upU.x / srcDims;
+  var o = textureSampleLevel(srcTex, srcSamp, uv + t * vec2f(-1.0, -1.0), 0.0).rgb * 1.0;
+  o = o + textureSampleLevel(srcTex, srcSamp, uv + t * vec2f( 0.0, -1.0), 0.0).rgb * 2.0;
+  o = o + textureSampleLevel(srcTex, srcSamp, uv + t * vec2f( 1.0, -1.0), 0.0).rgb * 1.0;
+  o = o + textureSampleLevel(srcTex, srcSamp, uv + t * vec2f(-1.0,  0.0), 0.0).rgb * 2.0;
+  o = o + textureSampleLevel(srcTex, srcSamp, uv + t * vec2f( 0.0,  0.0), 0.0).rgb * 4.0;
+  o = o + textureSampleLevel(srcTex, srcSamp, uv + t * vec2f( 1.0,  0.0), 0.0).rgb * 2.0;
+  o = o + textureSampleLevel(srcTex, srcSamp, uv + t * vec2f(-1.0,  1.0), 0.0).rgb * 1.0;
+  o = o + textureSampleLevel(srcTex, srcSamp, uv + t * vec2f( 0.0,  1.0), 0.0).rgb * 2.0;
+  o = o + textureSampleLevel(srcTex, srcSamp, uv + t * vec2f( 1.0,  1.0), 0.0).rgb * 1.0;
+  o = o * (1.0 / 16.0);
+  let base = textureSampleLevel(baseTex, srcSamp, uv, 0.0).rgb;
+  return vec4f(o + base, 1.0);
+}
+`