three-realtime-rt 0.3.0 → 0.3.2

package/src/RTLightingPass.js

@@ -1,680 +1,730 @@
-import * as THREE from "three";
-import { shaderStructs, shaderIntersectFunction } from "three-mesh-bvh";
-import { MAX_LIGHTS } from "./SceneCompiler.js";
-import { SKY_GLSL } from "./sky.glsl.js";
-import { BVH_ANY_HIT_GLSL } from "./bvhAnyHit.glsl.js";
-
-const fullscreenVert = /* glsl */ `
-out vec2 vUv;
-void main() {
-  vUv = uv;
-  gl_Position = vec4(position.xy, 0.0, 1.0);
-}
-`;
-
-const rtLightingFrag = /* glsl */ `
-precision highp float;
-precision highp isampler2D;
-precision highp usampler2D;
-
-${shaderStructs}
-${shaderIntersectFunction}
-${BVH_ANY_HIT_GLSL}
-${SKY_GLSL}
-
-#define MAX_LIGHTS ${MAX_LIGHTS}
-#define PI 3.14159265358979
-
-layout(location = 0) out vec4 outIrradiance;
-
-in vec2 vUv;
-
-// Two-level BVH: static (uploaded once) + dynamic (small, refit each frame).
-uniform BVH bvhStatic;
-uniform BVH bvhDynamic;
-uniform bool uHasDynamic;
-// One packed per-vertex texture per level: normal.xyz + materialIndex.w.
-// (Two BVH structs already use 8 samplers; WebGL2 guarantees only 16 total.)
-uniform sampler2D uAttrStatic;
-uniform sampler2D uAttrDynamic;
-uniform sampler2D uMaterialsTex;        // 2 texels per material (shared)
-
-uniform sampler2D uGWorldPos;
-uniform sampler2D uGNormalMetal;
-
-// temporal reprojection (stage 2). Validation is plane-distance only — the
-// normal test was dropped to free a sampler for the ReSTIR reservoir (same
-// simplification the TAA resolve already made, no observed regressions).
-uniform sampler2D uPrevAccum;        // rgb = irradiance history, a = sample count
-uniform sampler2D uPrevGWorldPos;    // previous frame's G-buffer, for validation
-uniform sampler2D uReservoir;        // ReSTIR winner per pixel (see RestirPass)
-uniform mat4 uPrevViewProj;
-uniform mat4 uViewProj;
-uniform vec3 uCameraPos;
-uniform float uMaxHistory;
-uniform bool uTemporalReprojection;
-uniform float uFireflyClamp;
-
-uniform vec4 uLightPosType[MAX_LIGHTS];     // xyz pos|dir, w: 0 point, 1 directional
-uniform vec4 uLightColorRadius[MAX_LIGHTS]; // rgb color*intensity, w radius
-uniform int uLightCount;
-uniform int uEmissiveCount; // NEE area-light triangles in row 1 of uMaterialsTex
-uniform bool uReflEnabled;  // traced reflections on metallic surfaces
-uniform bool uRefrEnabled;  // traced refraction on transmissive surfaces
-uniform float uIor;         // index of refraction for transmissive materials
-uniform bool uLightStochastic; // 1 direct shadow ray/pixel/frame instead of 1/light
-uniform bool uRestirEnabled;   // shade the reservoir winner instead of sampling
-
-uniform vec3 uEnvColor;
-uniform float uEnvIntensity;
-uniform float uFrame;
-uniform float uEps;
-uniform bool uGIEnabled;
-
-// Procedural sky (when enabled, replaces the flat env colour as the "miss" term
-// for GI rays — this is what gives natural outdoor bounce light).
-uniform bool uSkyEnabled;
-uniform vec3 uSunDir;      // direction toward the sun
-uniform vec3 uSunColor;
-uniform vec3 uSkyZenith;
-uniform vec3 uSkyHorizon;
-uniform float uSkyIntensity;
-
-// ---------- RNG ----------
-// The FIRST four random numbers each frame come from a 64x64 blue-noise tile
-// (rows 2..65 of the scene-data texture), rotated over time with an R2
-// low-discrepancy sequence. Those dimensions drive direct lighting — light
-// pick + area-sample position — where noise is most visible; blue noise turns
-// the residual error high-frequency, which temporal accumulation and the
-// denoiser remove far better than white-noise clumps. Later dimensions fall
-// back to PCG white noise (correlating many dimensions hurts).
-uint gSeed;
-int gBnDim;
-vec4 gBlueNoise;
-uint pcgHash(uint s) {
-  uint state = s * 747796405u + 2891336453u;
-  uint word = ((state >> ((state >> 28u) + 4u)) ^ state) * 277803737u;
-  return (word >> 22u) ^ word;
-}
-float rand() {
-  if (gBnDim < 4) {
-    float v = gBlueNoise[gBnDim];
-    gBnDim++;
-    return v;
-  }
-  gSeed = pcgHash(gSeed);
-  return float(gSeed) * (1.0 / 4294967296.0);
-}
-vec2 rand2() { return vec2(rand(), rand()); }
-
-vec4 fetchBlueNoise() {
-  ivec2 p = ivec2(gl_FragCoord.xy) & 63;
-  vec4 bn = texelFetch(uMaterialsTex, ivec2(p.x, 2 + p.y), 0);
-  // R2 sequence: per-frame toroidal shift, decorrelated per channel.
-  vec4 shift = fract(float(uFrame) * vec4(0.6180340, 0.7548777, 0.5698403, 0.8191725));
-  return fract(bn + shift);
-}
-
-// Branchless orthonormal basis (Duff et al. 2017) — cheaper and stable for
-// every n, including the poles the old cross-product picker handled branchily.
-void orthoBasis(vec3 n, out vec3 t, out vec3 b) {
-  float s = n.z >= 0.0 ? 1.0 : -1.0;
-  float a = -1.0 / (s + n.z);
-  float m = n.x * n.y * a;
-  t = vec3(1.0 + s * n.x * n.x * a, s * m, -s * n.x);
-  b = vec3(m, s + n.y * n.y * a, -n.y);
-}
-
-vec3 cosineSampleHemisphere(vec3 n, vec2 u) {
-  float a = 2.0 * PI * u.x;
-  float r = sqrt(u.y);
-  vec3 t, b;
-  orthoBasis(n, t, b);
-  return normalize(t * (r * cos(a)) + b * (r * sin(a)) + n * sqrt(max(0.0, 1.0 - u.y)));
-}
-
-vec3 randUnitVector() {
-  vec2 u = rand2();
-  float z = u.x * 2.0 - 1.0;
-  float a = u.y * 2.0 * PI;
-  float r = sqrt(max(0.0, 1.0 - z * z));
-  return vec3(r * cos(a), r * sin(a), z);
-}
-
-// ---------- two-level BVH helpers ----------
-// Closest hit across both levels; isDyn says which one so the caller samples
-// the matching vertex-attribute textures. (No backticks in these GLSL comments —
-// they would terminate the enclosing JS template literal.)
-bool traceBoth(vec3 ro, vec3 rd, out uvec4 fi, out vec3 bary, out float dist, out bool isDyn) {
-  uvec4 fiS; vec3 fnS; vec3 bcS; float sideS; float distS;
-  bool hitS = bvhIntersectFirstHit(bvhStatic, ro, rd, fiS, fnS, bcS, sideS, distS);
-  uvec4 fiD; vec3 fnD; vec3 bcD; float sideD; float distD;
-  bool hitD = uHasDynamic && bvhIntersectFirstHit(bvhDynamic, ro, rd, fiD, fnD, bcD, sideD, distD);
-  if (hitS && (!hitD || distS <= distD)) { fi = fiS; bary = bcS; dist = distS; isDyn = false; return true; }
-  if (hitD) { fi = fiD; bary = bcD; dist = distD; isDyn = true; return true; }
-  return false;
-}
-
-// Shadow rays only need to know IF something blocks, not what's closest —
-// the unordered any-hit traversal early-outs on the first blocker.
-bool occluded(vec3 ro, vec3 rd, float maxDist) {
-  if (bvhIntersectAnyHit(bvhStatic, ro, rd, maxDist - 2.0 * uEps)) return true;
-  if (uHasDynamic && bvhIntersectAnyHit(bvhDynamic, ro, rd, maxDist - 2.0 * uEps)) return true;
-  return false;
-}
-
-void fetchMaterial(float matIndex, out vec3 albedo, out float roughness,
-                   out vec3 emissive, out float metalness) {
-  int mi = int(round(matIndex)) * 2;
-  vec4 t0 = texelFetch(uMaterialsTex, ivec2(mi, 0), 0);
-  vec4 t1 = texelFetch(uMaterialsTex, ivec2(mi + 1, 0), 0);
-  albedo = t0.rgb;
-  roughness = t0.a;
-  emissive = t1.rgb;
-  metalness = t1.a;
-}
-
-// ---------- lighting ----------
-// Direct irradiance (demodulated: no albedo) at point P with normal N,
-// from light i, with one shadow ray. Area-samples point lights for soft shadows.
-vec3 lightContribution(int i, vec3 P, vec3 N) {
-  vec4 posType = uLightPosType[i];
-  vec4 colRad = uLightColorRadius[i];
-
-  vec3 L;
-  float dist2 = 1.0;
-  float maxDist = 1e7;
-
-  if (posType.w < 0.5) {
-    // point light: sample a point on its sphere for soft shadows
-    vec3 lp = posType.xyz + randUnitVector() * colRad.w;
-    vec3 d = lp - P;
-    float dl = length(d);
-    if (dl < 1e-5) return vec3(0.0);
-    L = d / dl;
-    dist2 = dl * dl;
-    maxDist = dl;
-  } else {
-    // directional light: jitter within a small cone
-    L = normalize(-posType.xyz + randUnitVector() * colRad.w);
-    dist2 = 1.0;
-  }
-
-  float NdotL = dot(N, L);
-  if (NdotL <= 0.0) return vec3(0.0);
-
-  if (occluded(P + N * uEps, L, maxDist)) return vec3(0.0);
-  return colRad.rgb * (NdotL / dist2);
-}
-
-// Direct light at a GI bounce hit: sample ONE random light (weighted by count).
-vec3 sampleOneLight(vec3 P, vec3 N) {
-  if (uLightCount == 0) return vec3(0.0);
-  int i = min(int(rand() * float(uLightCount)), uLightCount - 1);
-  return lightContribution(i, P, N) * float(uLightCount);
-}
-
-// Next-event estimation on emissive-mesh triangles (row 1 of uMaterialsTex):
-// pick one triangle, sample a point on it, cast one shadow ray, convert the
-// area pdf to solid angle. Turns emitters into proper soft area lights instead
-// of surfaces a GI ray has to hit by luck.
-vec3 sampleEmissiveTri(vec3 P, vec3 N) {
-  if (uEmissiveCount == 0) return vec3(0.0);
-  int i = min(int(rand() * float(uEmissiveCount)), uEmissiveCount - 1) * 4;
-  vec4 t0 = texelFetch(uMaterialsTex, ivec2(i, 1), 0);     // v0 | area
-  vec4 t1 = texelFetch(uMaterialsTex, ivec2(i + 1, 1), 0); // e1 | emit.r
-  vec4 t2 = texelFetch(uMaterialsTex, ivec2(i + 2, 1), 0); // e2 | emit.g
-  vec4 t3 = texelFetch(uMaterialsTex, ivec2(i + 3, 1), 0); // n  | emit.b
-
-  vec2 u = rand2();
-  if (u.x + u.y > 1.0) u = 1.0 - u; // uniform over the triangle
-  vec3 lp = t0.xyz + t1.xyz * u.x + t2.xyz * u.y;
-
-  vec3 d = lp - P;
-  float d2 = dot(d, d);
-  float dist = sqrt(d2);
-  if (dist < 1e-4) return vec3(0.0);
-  vec3 wi = d / dist;
-
-  float cosS = dot(N, wi);
-  // abs(): double-sided emission, matching what a GI ray hitting either face sees.
-  float cosL = abs(dot(t3.xyz, wi));
-  if (cosS <= 0.0 || cosL < 1e-4) return vec3(0.0);
-  if (occluded(P + N * uEps, wi, dist)) return vec3(0.0);
-
-  // Uniform pick of 1-of-count tris + uniform point: pdf_area = 1/(count·area).
-  // Solid-angle conversion gives irradiance Le · cosS · cosL / (d² · pdf_area).
-  vec3 e = vec3(t1.w, t2.w, t3.w) * (cosS * cosL * float(uEmissiveCount) * t0.w / max(d2, 1e-6));
-
-  // Uniform-area sampling has huge single-sample variance for receivers close
-  // to a big emitter (sampled point can land almost on top of P, d² → 0);
-  // those 100× spikes read as speckles because the EMA decays them only as
-  // 1/count. Clamp at 2× the indirect firefly limit — slight bias right next
-  // to the emitter, stable everywhere.
-  float eLum = dot(e, vec3(0.299, 0.587, 0.114));
-  float eCap = uFireflyClamp * 2.0;
-  if (eLum > eCap) e *= eCap / eLum;
-  return e;
-}
-
-// Shade this pixel's ReSTIR reservoir winner: recompute the (unshadowed)
-// contribution — MUST match RestirPass.candidateContribution — then pay the
-// one visibility ray and weight by W = wSum / (M · p̂). Analytic lights
-// re-draw their soft-radius jitter here (the reservoir stores which light,
-// not the jitter). The estimator inherently tames near-emitter spikes: a huge
-// contribution comes with a proportionally huge p̂, and W divides it out.
-vec3 shadeReservoir(vec3 P, vec3 N) {
-  // Spatial-stage encoding: r = id, a = precomputed W (vs. centroid score).
-  vec4 res = texture(uReservoir, vUv);
-  if (res.a <= 0.0) return vec3(0.0);
-  float id = res.r;
-
-  vec3 C;
-  vec3 wi;
-  float maxDist;
-  if (id < float(MAX_LIGHTS)) {
-    int i = int(id);
-    vec4 posType = uLightPosType[i];
-    vec4 colRad = uLightColorRadius[i];
-    if (posType.w < 0.5) {
-      vec3 d = posType.xyz - P;
-      float dl = length(d);
-      if (dl < 1e-5) return vec3(0.0);
-      float NdotL = dot(N, d / dl);
-      if (NdotL <= 0.0) return vec3(0.0);
-      C = colRad.rgb * (NdotL / (dl * dl));
-      vec3 lp = posType.xyz + randUnitVector() * colRad.w; // soft shadows
-      vec3 dj = lp - P;
-      maxDist = length(dj);
-      if (maxDist < 1e-5) return vec3(0.0);
-      wi = dj / maxDist;
-    } else {
-      float NdotL = dot(N, -posType.xyz);
-      if (NdotL <= 0.0) return vec3(0.0);
-      C = colRad.rgb * NdotL;
-      wi = normalize(-posType.xyz + randUnitVector() * colRad.w);
-      maxDist = 1e7;
-    }
-  } else {
-    int t = (int(id) - MAX_LIGHTS) * 4;
-    vec4 t0 = texelFetch(uMaterialsTex, ivec2(t, 1), 0);
-    vec4 t1 = texelFetch(uMaterialsTex, ivec2(t + 1, 1), 0);
-    vec4 t2 = texelFetch(uMaterialsTex, ivec2(t + 2, 1), 0);
-    vec4 t3 = texelFetch(uMaterialsTex, ivec2(t + 3, 1), 0);
-    // v3: the reservoir chose the TRIANGLE; draw a FRESH point on it every
-    // frame so the area light keeps averaging (no frozen-point noise). W was
-    // normalized against the centroid score, and E[point sample] = the
-    // triangle's true contribution, so the estimator stays consistent.
-    vec2 uv = rand2();
-    if (uv.x + uv.y > 1.0) uv = 1.0 - uv;
-    vec3 lp = t0.xyz + t1.xyz * uv.x + t2.xyz * uv.y;
-    vec3 d = lp - P;
-    float d2 = dot(d, d);
-    maxDist = sqrt(d2);
-    if (maxDist < 1e-4) return vec3(0.0);
-    wi = d / maxDist;
-    float cosS = dot(N, wi);
-    float cosL = abs(dot(t3.xyz, wi));
-    if (cosS <= 0.0 || cosL < 1e-4) return vec3(0.0);
-    C = vec3(t1.w, t2.w, t3.w) * (cosS * cosL * t0.w / max(d2, 1e-6));
-  }
-
-  if (occluded(P + N * uEps, wi, maxDist)) return vec3(0.0);
-  vec3 e = C * res.a;
-  // Safety clamp, same budget as the emissive direct clamp elsewhere.
-  float l = dot(e, vec3(0.299, 0.587, 0.114));
-  float cap = uFireflyClamp * 2.0;
-  if (l > cap) e *= cap / l;
-  return e;
-}
-
-// ONE light sample for secondary path vertices: stochastically pick either the
-// analytic lights or the emissive set (weighted 1/p). Costs a single shadow
-// ray — same ray budget the GI bounce had before emissive NEE existed —
-// instead of two; the estimator stays unbiased and temporal accumulation
-// averages out the extra variance.
-vec3 sampleOneAny(vec3 P, vec3 N) {
-  bool hasL = uLightCount > 0;
-  bool hasE = uEmissiveCount > 0;
-  if (hasL && hasE) {
-    return rand() < 0.5
-      ? sampleOneLight(P, N) * 2.0
-      : sampleEmissiveTri(P, N) * 2.0;
-  }
-  if (hasL) return sampleOneLight(P, N);
-  if (hasE) return sampleEmissiveTri(P, N);
-  return vec3(0.0);
-}
-
-// Incoming radiance along rd: trace, shade the hit with direct + NEE lighting,
-// sky/env on a miss. Specular rays keep emitter emission on hit (NEE at the ray
-// origin cannot cover a specular path); diffuse GI rays drop it for NEE-listed
-// (static) emitters so that light isn't counted twice.
-vec3 traceRadiance(vec3 ro, vec3 rd, bool specular) {
-  uvec4 fi; vec3 bary; float dist; bool isDyn;
-  if (!traceBoth(ro, rd, fi, bary, dist, isDyn)) {
-    return uSkyEnabled
-      ? skyColor(rd, uSunDir, uSunColor, uSkyZenith, uSkyHorizon, uSkyIntensity)
-      : uEnvColor * uEnvIntensity;
-  }
-  vec4 attr = isDyn
-    ? textureSampleBarycoord(uAttrDynamic, bary, fi.xyz)
-    : textureSampleBarycoord(uAttrStatic, bary, fi.xyz);
-  vec3 hAlbedo; float hRough; vec3 hEmissive; float hMetal;
-  fetchMaterial(attr.w, hAlbedo, hRough, hEmissive, hMetal);
-  vec3 hN = normalize(attr.xyz);
-  if (dot(hN, rd) > 0.0) hN = -hN;
-  vec3 hP = ro + rd * dist;
-  vec3 Ld = sampleOneAny(hP + hN * uEps, hN);
-  vec3 hLe = (!specular && uEmissiveCount > 0 && !isDyn) ? vec3(0.0) : hEmissive;
-  return hLe + hAlbedo * Ld * (1.0 / PI);
-}
-
-float schlick(float cosT, float eta) {
-  float r0 = (1.0 - eta) / (1.0 + eta);
-  r0 *= r0;
-  return r0 + (1.0 - r0) * pow(1.0 - cosT, 5.0);
-}
-
-// Roughness-jittered mirror direction (glossy cone approximation).
-vec3 glossyReflect(vec3 V, vec3 N, float rough) {
-  vec3 refl = reflect(V, N);
-  if (rough > 0.0) {
-    refl = normalize(mix(refl, cosineSampleHemisphere(N, rand2()), rough * rough));
-  }
-  return refl;
-}
-
-// Glass: Fresnel-weighted blend of a surface reflection and a two-interface
-// refraction (enter at P, march to the exit surface, refract again).
-vec3 glassRadiance(vec3 P, vec3 N, vec3 V, float rough) {
-  vec3 refl = glossyReflect(V, N, rough);
-  vec3 reflRad = dot(refl, N) > 0.0 ? traceRadiance(P + N * uEps, refl, true) : vec3(0.0);
-
-  float eta = 1.0 / uIor;
-  vec3 rd = refract(V, N, eta);
-  if (rd == vec3(0.0)) return reflRad; // total internal reflection at entry
-  float fres = schlick(clamp(-dot(V, N), 0.0, 1.0), eta);
-
-  vec3 ro = P - N * (2.0 * uEps);
-  vec3 refrRad;
-  uvec4 fi; vec3 bary; float dist; bool isDyn;
-  if (traceBoth(ro, rd, fi, bary, dist, isDyn)) {
-    // Exit interface: refract back out (or bounce once on internal reflection).
-    vec4 attr = isDyn
-      ? textureSampleBarycoord(uAttrDynamic, bary, fi.xyz)
-      : textureSampleBarycoord(uAttrStatic, bary, fi.xyz);
-    vec3 xN = normalize(attr.xyz);
-    if (dot(xN, rd) > 0.0) xN = -xN;
-    vec3 xP = ro + rd * dist;
-    vec3 rd2 = refract(rd, xN, uIor);
-    if (rd2 == vec3(0.0)) rd2 = reflect(rd, xN);
-    refrRad = traceRadiance(xP - xN * uEps, rd2, true);
-  } else {
-    refrRad = uSkyEnabled
-      ? skyColor(rd, uSunDir, uSunColor, uSkyZenith, uSkyHorizon, uSkyIntensity)
-      : uEnvColor * uEnvIntensity;
-  }
-  return mix(refrRad, reflRad, fres);
-}
-
-void main() {
-  vec4 wp = texture(uGWorldPos, vUv);
-  if (wp.w < 0.5) {
-    outIrradiance = vec4(0.0);
-    return;
-  }
-
-  ivec2 px = ivec2(gl_FragCoord.xy);
-  gSeed = uint(px.x) * 1973u + uint(px.y) * 9277u + uint(uFrame) * 26699u;
-  gSeed = pcgHash(gSeed);
-  gBlueNoise = fetchBlueNoise();
-  gBnDim = 0;
-
-  vec3 P = wp.xyz;
-  vec4 nmSample = texture(uGNormalMetal, vUv);
-  vec3 N = normalize(nmSample.xyz);
-  // Decode the packed material word (see GBufferPass): >= 2 → glass, else metal.
-  float transmission = nmSample.w >= 2.0 ? clamp(nmSample.w - 2.0, 0.0, 1.0) : 0.0;
-  float metal = nmSample.w >= 2.0 ? 0.0 : nmSample.w;
-  float rough = clamp(wp.w - 1.0, 0.0, 1.0);
-
-  // --- direct lighting ---
-  // ReSTIR: shade the reservoir's winner with one visibility ray (flat cost in
-  // light count). Stochastic: one blind random sample. Full: one shadow ray
-  // per light + one for the emissive set.
-  vec3 direct = vec3(0.0);
-  if (uRestirEnabled) {
-    direct = shadeReservoir(P, N);
-  } else if (uLightStochastic) {
-    direct = sampleOneAny(P, N);
-  } else {
-    for (int i = 0; i < MAX_LIGHTS; i++) {
-      if (i >= uLightCount) break;
-      direct += lightContribution(i, P, N);
-    }
-    // Emissive meshes as area lights (next-event estimation, one shadow ray).
-    direct += sampleEmissiveTri(P, N);
-  }
-
-  // --- 1-bounce indirect (cosine-weighted; pdf cancels the NdotL/PI).
-  // traceRadiance shades the hit with direct + NEE light, or returns the
-  // sky/env colour when the ray escapes (the natural ambient bounce).
-  vec3 indirect = vec3(0.0);
-  if (uGIEnabled) {
-    indirect = traceRadiance(P + N * uEps, cosineSampleHemisphere(N, rand2()), false);
-  }
-
-  // Firefly clamp: suppress rare huge GI samples (big perceived-noise win,
-  // slightly biased). Applied to indirect only; direct is analytic.
-  float lum = dot(indirect, vec3(0.299, 0.587, 0.114));
-  if (lum > uFireflyClamp) indirect *= uFireflyClamp / lum;
-
-  vec3 sampleIrr = direct + indirect;
-
-  // --- traced specular: mirror/glossy reflections on metals ---
-  if (uReflEnabled && metal > 0.001) {
-    vec3 V = normalize(P - uCameraPos);
-    vec3 refl = glossyReflect(V, N, rough);
-    if (dot(refl, N) > 0.0) {
-      // Metals have no diffuse term: replace by metalness. The composite's
-      // albedo multiply then tints the reflection (F0 = albedo for metals).
-      sampleIrr = mix(sampleIrr, traceRadiance(P + N * uEps, refl, true), metal);
-    }
-  }
-
-  // --- traced glass: Fresnel reflection + two-interface refraction ---
-  if (uRefrEnabled && transmission > 0.001) {
-    vec3 V = normalize(P - uCameraPos);
-    sampleIrr = mix(sampleIrr, glassRadiance(P, N, V, rough), transmission);
-  }
-
-  // A single NaN/Inf sample would poison the EMA history for good (mix() with
-  // NaN stays NaN until a disocclusion resets the pixel) — sanitize first.
-  if (any(isnan(sampleIrr)) || any(isinf(sampleIrr))) sampleIrr = vec3(0.0);
-
-  // --- temporal reprojection: pull validated history from last frame ---
-  float count = 1.0;
-  vec3 history = vec3(0.0);
-  if (uTemporalReprojection) {
-    vec4 clip = uPrevViewProj * vec4(P, 1.0);
-    vec4 clipC = uViewProj * vec4(P, 1.0);
-    if (clip.w > 0.0 && clipC.w > 0.0) {
-      vec2 prevUv = (clip.xy / clip.w) * 0.5 + 0.5;
-      // P comes from a full-res G-buffer texel, which sits sub-pixel off this
-      // half-res fragment's center. That constant offset would bias bilinear
-      // history reads every frame (content drifts/smears at renderScale < 1).
-      // Cancel it: measure P's offset in the CURRENT frame and subtract.
-      vec2 currUv = (clipC.xy / clipC.w) * 0.5 + 0.5;
-      prevUv -= currUv - vUv;
-      if (prevUv.x >= 0.0 && prevUv.x <= 1.0 && prevUv.y >= 0.0 && prevUv.y <= 1.0) {
-        vec4 prevPos = texture(uPrevGWorldPos, prevUv);
-        // Plane-distance test: robust at grazing angles (position error from
-        // texel quantization lies along the surface, not along the normal).
-        float distToCam = distance(P, uCameraPos);
-        float tol = 0.005 * distToCam + 20.0 * uEps;
-        bool valid = prevPos.w > 0.5
-          && abs(dot(P - prevPos.xyz, N)) < tol;
-        if (valid) {
-          vec4 h = texture(uPrevAccum, prevUv); // bilinear
-          // Mirror-like pixels keep a SHORT history: their reflected content
-          // moves differently from the surface, so long history smears the
-          // reflection under camera motion — and specular rays are nearly
-          // deterministic, so they don't need the accumulation anyway.
-          float specHist = max(metal, transmission) * (1.0 - rough);
-          float histCap = mix(uMaxHistory, min(uMaxHistory, 10.0), specHist);
-          count = clamp(h.a, 0.0, histCap) + 1.0;
-          history = h.rgb;
-        }
-      }
-    }
-  }
-
-  // Exponential moving average; count=1 (disocclusion / first frame) means
-  // the fresh sample is used as-is.
-  vec3 blended = mix(history, sampleIrr, 1.0 / count);
-  outIrradiance = vec4(blended, count);
-}
-`;
-
-/**
- * Fullscreen pass: for every G-buffer pixel, trace shadow rays to every light and
- * one cosine-weighted GI bounce against the BVH. Outputs demodulated irradiance,
- * progressively accumulated into a ping-pong float target while the camera is still.
- */
-export class RTLightingPass {
-  constructor(width, height) {
-    this.targetA = this._makeTarget(width, height);
-    this.targetB = this._makeTarget(width, height);
-
-    this.material = new THREE.ShaderMaterial({
-      glslVersion: THREE.GLSL3,
-      vertexShader: fullscreenVert,
-      fragmentShader: rtLightingFrag,
-      uniforms: {
-        bvhStatic: { value: null },
-        bvhDynamic: { value: null },
-        uHasDynamic: { value: false },
-        uAttrStatic: { value: null },
-        uAttrDynamic: { value: null },
-        uMaterialsTex: { value: null },
-        uGWorldPos: { value: null },
-        uGNormalMetal: { value: null },
-        uPrevAccum: { value: null },
-        uPrevGWorldPos: { value: null },
-        uReservoir: { value: null },
-        uRestirEnabled: { value: false },
-        uPrevViewProj: { value: new THREE.Matrix4() },
-        uViewProj: { value: new THREE.Matrix4() },
-        uCameraPos: { value: new THREE.Vector3() },
-        uMaxHistory: { value: 128 },
-        uTemporalReprojection: { value: true },
-        uFireflyClamp: { value: 4.0 },
-        uLightPosType: { value: [] },
-        uLightColorRadius: { value: [] },
-        uLightCount: { value: 0 },
-        uEmissiveCount: { value: 0 },
-        uReflEnabled: { value: true },
-        uRefrEnabled: { value: true },
-        uIor: { value: 1.5 },
-        uLightStochastic: { value: false },
-        uEnvColor: { value: new THREE.Color(0.03, 0.04, 0.06) },
-        uEnvIntensity: { value: 1.0 },
-        uFrame: { value: 0 },
-        uEps: { value: 1e-3 },
-        uGIEnabled: { value: true },
-        uSkyEnabled: { value: false },
-        uSunDir: { value: new THREE.Vector3(0.4, 0.8, 0.45).normalize() },
-        uSunColor: { value: new THREE.Color(1.0, 0.9, 0.75) },
-        uSkyZenith: { value: new THREE.Color(0.18, 0.34, 0.62) },
-        uSkyHorizon: { value: new THREE.Color(0.7, 0.8, 0.9) },
-        uSkyIntensity: { value: 1.0 },
-      },
-      depthTest: false,
-      depthWrite: false,
-    });
-
-    this.scene = new THREE.Scene();
-    this.camera = new THREE.OrthographicCamera(-1, 1, 1, -1, 0, 1);
-    this.quad = new THREE.Mesh(new THREE.PlaneGeometry(2, 2), this.material);
-    this.quad.frustumCulled = false;
-    this.scene.add(this.quad);
-  }
-
-  _makeTarget(width, height) {
-    // Half-float + linear: history is sampled bilinearly at reprojected UVs,
-    // and fp16 halves the bandwidth (EMA blending never accumulates a raw sum,
-    // so fp16 precision is sufficient).
-    const t = new THREE.WebGLRenderTarget(width, height, {
-      minFilter: THREE.LinearFilter,
-      magFilter: THREE.LinearFilter,
-      format: THREE.RGBAFormat,
-      type: THREE.HalfFloatType,
-      depthBuffer: false,
-      stencilBuffer: false,
-    });
-    t.texture.generateMipmaps = false;
-    return t;
-  }
-
-  /** Clear both history buffers (e.g. after scene recompile or resize). */
-  clearHistory(renderer) {
-    const prevTarget = renderer.getRenderTarget();
-    const prevColor = new THREE.Color();
-    renderer.getClearColor(prevColor);
-    const prevAlpha = renderer.getClearAlpha();
-    renderer.setClearColor(0x000000, 0);
-    for (const t of [this.targetA, this.targetB]) {
-      renderer.setRenderTarget(t);
-      renderer.clear(true, false, false);
-    }
-    renderer.setRenderTarget(prevTarget);
-    renderer.setClearColor(prevColor, prevAlpha);
-  }
-
-  setSize(width, height) {
-    this.targetA.setSize(width, height);
-    this.targetB.setSize(width, height);
-  }
-
-  setCompiledScene(compiled) {
-    const u = this.material.uniforms;
-    u.bvhStatic.value = compiled.staticBvhUniform;
-    u.bvhDynamic.value = compiled.dynamicBvhUniform;
-    u.uHasDynamic.value = compiled.hasDynamic;
-    u.uAttrStatic.value = compiled.staticAttrTex;
-    u.uAttrDynamic.value = compiled.dynamicAttrTex;
-    u.uMaterialsTex.value = compiled.materialsTex;
-    u.uLightPosType.value = compiled.lightPosType;
-    u.uLightColorRadius.value = compiled.lightColorRadius;
-    u.uLightCount.value = compiled.lightCount;
-    u.uEmissiveCount.value = compiled.emissiveTriCount;
-  }
-
-  /** Renders into targetA (reading targetB as history), then swaps. Returns the fresh accum texture. */
-  render(renderer, gbuffer, frame, reservoirTexture = null) {
-    const u = this.material.uniforms;
-    u.uGWorldPos.value = gbuffer.worldPos;
-    u.uGNormalMetal.value = gbuffer.normalMetal;
-    u.uPrevGWorldPos.value = gbuffer.prevWorldPos;
-    u.uPrevAccum.value = this.targetB.texture;
-    u.uReservoir.value = reservoirTexture;
-    u.uRestirEnabled.value = reservoirTexture !== null;
-    u.uFrame.value = frame;
-
-    renderer.setRenderTarget(this.targetA);
-    renderer.render(this.scene, this.camera);
-    renderer.setRenderTarget(null);
-
-    const out = this.targetA;
-    [this.targetA, this.targetB] = [this.targetB, this.targetA];
-    return out.texture;
-  }
-
-  dispose() {
-    this.targetA.dispose();
-    this.targetB.dispose();
-    this.material.dispose();
-    this.quad.geometry.dispose();
-  }
-}
+import * as THREE from "three";
+import { shaderStructs, shaderIntersectFunction } from "three-mesh-bvh";
+import { MAX_LIGHTS } from "./SceneCompiler.js";
+import { SKY_GLSL } from "./sky.glsl.js";
+import { BVH_ANY_HIT_GLSL } from "./bvhAnyHit.glsl.js";
+
+const fullscreenVert = /* glsl */ `
+out vec2 vUv;
+void main() {
+  vUv = uv;
+  gl_Position = vec4(position.xy, 0.0, 1.0);
+}
+`;
+
+const rtLightingFrag = /* glsl */ `
+precision highp float;
+precision highp isampler2D;
+precision highp usampler2D;
+
+${shaderStructs}
+${shaderIntersectFunction}
+${BVH_ANY_HIT_GLSL}
+${SKY_GLSL}
+
+#define MAX_LIGHTS ${MAX_LIGHTS}
+#define PI 3.14159265358979
+
+layout(location = 0) out vec4 outIrradiance;
+
+in vec2 vUv;
+
+// Two-level BVH: static (uploaded once) + dynamic (small, refit each frame).
+uniform BVH bvhStatic;
+uniform BVH bvhDynamic;
+uniform bool uHasDynamic;
+// One packed per-vertex texture per level: normal.xyz + materialIndex.w.
+// (Two BVH structs already use 8 samplers; WebGL2 guarantees only 16 total.)
+uniform sampler2D uAttrStatic;
+uniform sampler2D uAttrDynamic;
+uniform sampler2D uMaterialsTex;        // 2 texels per material (shared)
+
+uniform sampler2D uGWorldPos;
+uniform sampler2D uGNormalMetal;
+
+// temporal reprojection (stage 2). Validation is plane-distance only — the
+// normal test was dropped to free a sampler for the ReSTIR reservoir (same
+// simplification the TAA resolve already made, no observed regressions).
+uniform sampler2D uPrevAccum;        // rgb = irradiance history, a = sample count
+uniform sampler2D uPrevGWorldPos;    // previous frame's G-buffer, for validation
+uniform sampler2D uReservoir;        // ReSTIR winner per pixel (see RestirPass)
+uniform mat4 uPrevViewProj;
+uniform mat4 uViewProj;
+uniform vec3 uCameraPos;
+uniform float uMaxHistory;
+uniform bool uTemporalReprojection;
+uniform float uFireflyClamp;
+
+uniform vec4 uLightPosType[MAX_LIGHTS];     // xyz pos|dir, w: 0 point, 1 directional, >=2 spot (w-2 = cosInner)
+uniform vec4 uLightColorRadius[MAX_LIGHTS]; // rgb color*intensity, w radius
+uniform vec4 uLightDirCone[MAX_LIGHTS];     // spot: direction.xyz + cos(outer angle)
+uniform int uLightCount;
+uniform int uEmissiveCount; // NEE area-light triangles in row 1 of uMaterialsTex
+uniform bool uReflEnabled;  // traced reflections on metallic surfaces
+uniform bool uRefrEnabled;  // traced refraction on transmissive surfaces
+uniform float uIor;         // index of refraction for transmissive materials
+uniform bool uLightStochastic; // 1 direct shadow ray/pixel/frame instead of 1/light
+uniform bool uRestirEnabled;   // shade the reservoir winner instead of sampling
+uniform bool uGIHalfRate;      // GI ray on alternating checkerboard, doubled
+
+uniform vec3 uEnvColor;
+uniform float uEnvIntensity;
+uniform float uFrame;
+uniform float uEps;
+uniform bool uGIEnabled;
+
+// Procedural sky (when enabled, replaces the flat env colour as the "miss" term
+// for GI rays — this is what gives natural outdoor bounce light).
+uniform bool uSkyEnabled;
+uniform vec3 uSunDir;      // direction toward the sun
+uniform vec3 uSunColor;
+uniform vec3 uSkyZenith;
+uniform vec3 uSkyHorizon;
+uniform float uSkyIntensity;
+
+// ---------- RNG ----------
+// The FIRST four random numbers each frame come from a 64x64 blue-noise tile
+// (rows 2..65 of the scene-data texture), rotated over time with an R2
+// low-discrepancy sequence. Those dimensions drive direct lighting — light
+// pick + area-sample position — where noise is most visible; blue noise turns
+// the residual error high-frequency, which temporal accumulation and the
+// denoiser remove far better than white-noise clumps. Later dimensions fall
+// back to PCG white noise (correlating many dimensions hurts).
+uint gSeed;
+int gBnDim;
+vec4 gBlueNoise;
+uint pcgHash(uint s) {
+  uint state = s * 747796405u + 2891336453u;
+  uint word = ((state >> ((state >> 28u) + 4u)) ^ state) * 277803737u;
+  return (word >> 22u) ^ word;
+}
+float rand() {
+  if (gBnDim < 4) {
+    float v = gBlueNoise[gBnDim];
+    gBnDim++;
+    return v;
+  }
+  gSeed = pcgHash(gSeed);
+  return float(gSeed) * (1.0 / 4294967296.0);
+}
+vec2 rand2() { return vec2(rand(), rand()); }
+
+vec4 fetchBlueNoise() {
+  ivec2 p = ivec2(gl_FragCoord.xy) & 63;
+  vec4 bn = texelFetch(uMaterialsTex, ivec2(p.x, 2 + p.y), 0);
+  // R2 sequence: per-frame toroidal shift, decorrelated per channel.
+  vec4 shift = fract(float(uFrame) * vec4(0.6180340, 0.7548777, 0.5698403, 0.8191725));
+  return fract(bn + shift);
+}
+
+// Branchless orthonormal basis (Duff et al. 2017) — cheaper and stable for
+// every n, including the poles the old cross-product picker handled branchily.
+void orthoBasis(vec3 n, out vec3 t, out vec3 b) {
+  float s = n.z >= 0.0 ? 1.0 : -1.0;
+  float a = -1.0 / (s + n.z);
+  float m = n.x * n.y * a;
+  t = vec3(1.0 + s * n.x * n.x * a, s * m, -s * n.x);
+  b = vec3(m, s + n.y * n.y * a, -n.y);
+}
+
+vec3 cosineSampleHemisphere(vec3 n, vec2 u) {
+  float a = 2.0 * PI * u.x;
+  float r = sqrt(u.y);
+  vec3 t, b;
+  orthoBasis(n, t, b);
+  return normalize(t * (r * cos(a)) + b * (r * sin(a)) + n * sqrt(max(0.0, 1.0 - u.y)));
+}
+
+vec3 randUnitVector() {
+  vec2 u = rand2();
+  float z = u.x * 2.0 - 1.0;
+  float a = u.y * 2.0 * PI;
+  float r = sqrt(max(0.0, 1.0 - z * z));
+  return vec3(r * cos(a), r * sin(a), z);
+}
+
+// ---------- two-level BVH helpers ----------
+// Closest hit across both levels; isDyn says which one so the caller samples
+// the matching vertex-attribute textures. (No backticks in these GLSL comments —
+// they would terminate the enclosing JS template literal.)
+bool traceBoth(vec3 ro, vec3 rd, out uvec4 fi, out vec3 bary, out float dist, out bool isDyn) {
+  uvec4 fiS; vec3 fnS; vec3 bcS; float sideS; float distS;
+  bool hitS = bvhIntersectFirstHit(bvhStatic, ro, rd, fiS, fnS, bcS, sideS, distS);
+  uvec4 fiD; vec3 fnD; vec3 bcD; float sideD; float distD;
+  bool hitD = uHasDynamic && bvhIntersectFirstHit(bvhDynamic, ro, rd, fiD, fnD, bcD, sideD, distD);
+  if (hitS && (!hitD || distS <= distD)) { fi = fiS; bary = bcS; dist = distS; isDyn = false; return true; }
+  if (hitD) { fi = fiD; bary = bcD; dist = distD; isDyn = true; return true; }
+  return false;
+}
+
+// Shadow rays only need to know IF something blocks, not what's closest —
+// the unordered any-hit traversal early-outs on the first blocker.
+bool occluded(vec3 ro, vec3 rd, float maxDist) {
+  if (bvhIntersectAnyHit(bvhStatic, ro, rd, maxDist - 2.0 * uEps)) return true;
+  if (uHasDynamic && bvhIntersectAnyHit(bvhDynamic, ro, rd, maxDist - 2.0 * uEps)) return true;
+  return false;
+}
+
+void fetchMaterial(float matIndex, out vec3 albedo, out float roughness,
+                   out vec3 emissive, out float metalness) {
+  int mi = int(round(matIndex)) * 2;
+  vec4 t0 = texelFetch(uMaterialsTex, ivec2(mi, 0), 0);
+  vec4 t1 = texelFetch(uMaterialsTex, ivec2(mi + 1, 0), 0);
+  albedo = t0.rgb;
+  roughness = t0.a;
+  emissive = t1.rgb;
+  metalness = t1.a;
+}
+
+// ---------- lighting ----------
+// Direct irradiance (demodulated: no albedo) at point P with normal N,
+// from light i, with one shadow ray. Area-samples point lights for soft shadows.
+// Spot cone falloff: smooth between the outer and inner cone cosines
+// (posType.w = 2 + cosInner; dirCone.w = cosOuter).
+float spotFalloff(int i, vec3 lightToP) {
+  vec4 posType = uLightPosType[i];
+  if (posType.w < 1.5) return 1.0;
+  vec4 dc = uLightDirCone[i];
+  return smoothstep(dc.w, posType.w - 2.0, dot(dc.xyz, lightToP));
+}
+
+vec3 lightContribution(int i, vec3 P, vec3 N) {
+  vec4 posType = uLightPosType[i];
+  vec4 colRad = uLightColorRadius[i];
+
+  vec3 L;
+  float dist2 = 1.0;
+  float maxDist = 1e7;
+  float cone = 1.0;
+
+  if (posType.w < 0.5 || posType.w >= 1.5) {
+    // point/spot light: sample a point on its sphere for soft shadows
+    vec3 lp = posType.xyz + randUnitVector() * colRad.w;
+    vec3 d = lp - P;
+    float dl = length(d);
+    if (dl < 1e-5) return vec3(0.0);
+    L = d / dl;
+    dist2 = dl * dl;
+    maxDist = dl;
+    cone = spotFalloff(i, -L);
+    if (cone <= 0.0) return vec3(0.0);
+  } else {
+    // directional light: jitter within a small cone
+    L = normalize(-posType.xyz + randUnitVector() * colRad.w);
+    dist2 = 1.0;
+  }
+
+  float NdotL = dot(N, L);
+  if (NdotL <= 0.0) return vec3(0.0);
+
+  if (occluded(P + N * uEps, L, maxDist)) return vec3(0.0);
+  return colRad.rgb * (cone * NdotL / dist2);
+}
+
+// Direct light at a GI bounce hit: sample ONE random light (weighted by count).
+vec3 sampleOneLight(vec3 P, vec3 N) {
+  if (uLightCount == 0) return vec3(0.0);
+  int i = min(int(rand() * float(uLightCount)), uLightCount - 1);
+  return lightContribution(i, P, N) * float(uLightCount);
+}
+
+// Next-event estimation on emissive-mesh triangles (row 1 of uMaterialsTex):
+// pick one triangle, sample a point on it, cast one shadow ray, convert the
+// area pdf to solid angle. Turns emitters into proper soft area lights instead
+// of surfaces a GI ray has to hit by luck.
+vec3 sampleEmissiveTri(vec3 P, vec3 N) {
+  if (uEmissiveCount == 0) return vec3(0.0);
+  int i = min(int(rand() * float(uEmissiveCount)), uEmissiveCount - 1) * 4;
+  vec4 t0 = texelFetch(uMaterialsTex, ivec2(i, 1), 0);     // v0 | area
+  vec4 t1 = texelFetch(uMaterialsTex, ivec2(i + 1, 1), 0); // e1 | emit.r
+  vec4 t2 = texelFetch(uMaterialsTex, ivec2(i + 2, 1), 0); // e2 | emit.g
+  vec4 t3 = texelFetch(uMaterialsTex, ivec2(i + 3, 1), 0); // n  | emit.b
+
+  vec2 u = rand2();
+  if (u.x + u.y > 1.0) u = 1.0 - u; // uniform over the triangle
+  vec3 lp = t0.xyz + t1.xyz * u.x + t2.xyz * u.y;
+
+  vec3 d = lp - P;
+  float d2 = dot(d, d);
+  float dist = sqrt(d2);
+  if (dist < 1e-4) return vec3(0.0);
+  vec3 wi = d / dist;
+
+  float cosS = dot(N, wi);
+  // abs(): double-sided emission, matching what a GI ray hitting either face sees.
+  float cosL = abs(dot(t3.xyz, wi));
+  if (cosS <= 0.0 || cosL < 1e-4) return vec3(0.0);
+  if (occluded(P + N * uEps, wi, dist)) return vec3(0.0);
+
+  // Uniform pick of 1-of-count tris + uniform point: pdf_area = 1/(count·area).
+  // Solid-angle conversion gives irradiance Le · cosS · cosL / (d² · pdf_area).
+  vec3 e = vec3(t1.w, t2.w, t3.w) * (cosS * cosL * float(uEmissiveCount) * t0.w / max(d2, 1e-6));
+
+  // Uniform-area sampling has huge single-sample variance for receivers close
+  // to a big emitter (sampled point can land almost on top of P, d² → 0);
+  // those 100× spikes read as speckles because the EMA decays them only as
+  // 1/count. Clamp at 2× the indirect firefly limit — slight bias right next
+  // to the emitter, stable everywhere.
+  float eLum = dot(e, vec3(0.299, 0.587, 0.114));
+  float eCap = uFireflyClamp * 2.0;
+  if (eLum > eCap) e *= eCap / eLum;
+  return e;
+}
+
+// Shade this pixel's ReSTIR reservoir winner: recompute the (unshadowed)
+// contribution — MUST match RestirPass.candidateContribution — then pay the
+// one visibility ray and weight by W = wSum / (M · p̂). Analytic lights
+// re-draw their soft-radius jitter here (the reservoir stores which light,
+// not the jitter). The estimator inherently tames near-emitter spikes: a huge
+// contribution comes with a proportionally huge p̂, and W divides it out.
+vec3 shadeReservoir(vec3 P, vec3 N) {
+  // Spatial-stage encoding: r = id, a = precomputed W (vs. centroid score).
+  vec4 res = texture(uReservoir, vUv);
+  if (res.a <= 0.0) return vec3(0.0);
+  float id = res.r;
+
+  vec3 C;
+  vec3 wi;
+  float maxDist;
+  if (id < float(MAX_LIGHTS)) {
+    int i = int(id);
+    vec4 posType = uLightPosType[i];
+    vec4 colRad = uLightColorRadius[i];
+    if (posType.w < 0.5 || posType.w >= 1.5) {
+      vec3 d = posType.xyz - P;
+      float dl = length(d);
+      if (dl < 1e-5) return vec3(0.0);
+      float NdotL = dot(N, d / dl);
+      if (NdotL <= 0.0) return vec3(0.0);
+      float cone = spotFalloff(i, -d / dl);
+      if (cone <= 0.0) return vec3(0.0);
+      C = colRad.rgb * (cone * NdotL / (dl * dl));
+      vec3 lp = posType.xyz + randUnitVector() * colRad.w; // soft shadows
+      vec3 dj = lp - P;
+      maxDist = length(dj);
+      if (maxDist < 1e-5) return vec3(0.0);
+      wi = dj / maxDist;
+    } else {
+      float NdotL = dot(N, -posType.xyz);
+      if (NdotL <= 0.0) return vec3(0.0);
+      C = colRad.rgb * NdotL;
+      wi = normalize(-posType.xyz + randUnitVector() * colRad.w);
+      maxDist = 1e7;
+    }
+  } else {
+    int t = (int(id) - MAX_LIGHTS) * 4;
+    vec4 t0 = texelFetch(uMaterialsTex, ivec2(t, 1), 0);
+    vec4 t1 = texelFetch(uMaterialsTex, ivec2(t + 1, 1), 0);
+    vec4 t2 = texelFetch(uMaterialsTex, ivec2(t + 2, 1), 0);
+    vec4 t3 = texelFetch(uMaterialsTex, ivec2(t + 3, 1), 0);
+    // v3: the reservoir chose the TRIANGLE; draw a FRESH point on it every
+    // frame so the area light keeps averaging (no frozen-point noise). W was
+    // normalized against the centroid score, and E[point sample] = the
+    // triangle's true contribution, so the estimator stays consistent.
+    vec2 uv = rand2();
+    if (uv.x + uv.y > 1.0) uv = 1.0 - uv;
+    vec3 lp = t0.xyz + t1.xyz * uv.x + t2.xyz * uv.y;
+    vec3 d = lp - P;
+    float d2 = dot(d, d);
+    maxDist = sqrt(d2);
+    if (maxDist < 1e-4) return vec3(0.0);
+    wi = d / maxDist;
+    float cosS = dot(N, wi);
+    float cosL = abs(dot(t3.xyz, wi));
+    if (cosS <= 0.0 || cosL < 1e-4) return vec3(0.0);
+    C = vec3(t1.w, t2.w, t3.w) * (cosS * cosL * t0.w / max(d2, 1e-6));
+  }
+
+  if (occluded(P + N * uEps, wi, maxDist)) return vec3(0.0);
+  vec3 e = C * res.a;
+  // Safety clamp, same budget as the emissive direct clamp elsewhere.
+  float l = dot(e, vec3(0.299, 0.587, 0.114));
+  float cap = uFireflyClamp * 2.0;
+  if (l > cap) e *= cap / l;
+  return e;
+}
+
+// ONE light sample for secondary path vertices: stochastically pick either the
+// analytic lights or the emissive set (weighted 1/p). Costs a single shadow
+// ray — same ray budget the GI bounce had before emissive NEE existed —
+// instead of two; the estimator stays unbiased and temporal accumulation
+// averages out the extra variance.
+vec3 sampleOneAny(vec3 P, vec3 N) {
+  bool hasL = uLightCount > 0;
+  bool hasE = uEmissiveCount > 0;
+  if (hasL && hasE) {
+    return rand() < 0.5
+      ? sampleOneLight(P, N) * 2.0
+      : sampleEmissiveTri(P, N) * 2.0;
+  }
+  if (hasL) return sampleOneLight(P, N);
+  if (hasE) return sampleEmissiveTri(P, N);
+  return vec3(0.0);
+}
+
+// Incoming radiance along rd: trace, shade the hit with direct + NEE lighting,
+// sky/env on a miss. Specular rays keep emitter emission on hit (NEE at the ray
+// origin cannot cover a specular path); diffuse GI rays drop it for NEE-listed
+// (static) emitters so that light isn't counted twice.
+vec3 traceRadiance(vec3 ro, vec3 rd, bool specular) {
+  uvec4 fi; vec3 bary; float dist; bool isDyn;
+  if (!traceBoth(ro, rd, fi, bary, dist, isDyn)) {
+    return uSkyEnabled
+      ? skyColor(rd, uSunDir, uSunColor, uSkyZenith, uSkyHorizon, uSkyIntensity)
+      : uEnvColor * uEnvIntensity;
+  }
+  vec4 attr = isDyn
+    ? textureSampleBarycoord(uAttrDynamic, bary, fi.xyz)
+    : textureSampleBarycoord(uAttrStatic, bary, fi.xyz);
+  vec3 hAlbedo; float hRough; vec3 hEmissive; float hMetal;
+  fetchMaterial(attr.w, hAlbedo, hRough, hEmissive, hMetal);
+  vec3 hN = normalize(attr.xyz);
+  if (dot(hN, rd) > 0.0) hN = -hN;
+  vec3 hP = ro + rd * dist;
+  vec3 Ld = sampleOneAny(hP + hN * uEps, hN);
+  vec3 hLe = (!specular && uEmissiveCount > 0 && !isDyn) ? vec3(0.0) : hEmissive;
+  return hLe + hAlbedo * Ld * (1.0 / PI);
+}
+
+float schlick(float cosT, float eta) {
+  float r0 = (1.0 - eta) / (1.0 + eta);
+  r0 *= r0;
+  return r0 + (1.0 - r0) * pow(1.0 - cosT, 5.0);
+}
+
+// Roughness-jittered mirror direction (glossy cone approximation).
+vec3 glossyReflect(vec3 V, vec3 N, float rough) {
+  vec3 refl = reflect(V, N);
+  if (rough > 0.0) {
+    refl = normalize(mix(refl, cosineSampleHemisphere(N, rand2()), rough * rough));
+  }
+  return refl;
+}
+
+// Glass: Fresnel-weighted blend of a surface reflection and a two-interface
+// refraction (enter at P, march to the exit surface, refract again).
+vec3 glassRadiance(vec3 P, vec3 N, vec3 V, float rough) {
+  vec3 refl = glossyReflect(V, N, rough);
+  vec3 reflRad = dot(refl, N) > 0.0 ? traceRadiance(P + N * uEps, refl, true) : vec3(0.0);
+
+  float eta = 1.0 / uIor;
+  vec3 rd = refract(V, N, eta);
+  if (rd == vec3(0.0)) return reflRad; // total internal reflection at entry
+  float fres = schlick(clamp(-dot(V, N), 0.0, 1.0), eta);
+
+  vec3 ro = P - N * (2.0 * uEps);
+  vec3 refrRad;
+  uvec4 fi; vec3 bary; float dist; bool isDyn;
+  if (traceBoth(ro, rd, fi, bary, dist, isDyn)) {
+    // Exit interface: refract back out (or bounce once on internal reflection).
+    vec4 attr = isDyn
+      ? textureSampleBarycoord(uAttrDynamic, bary, fi.xyz)
+      : textureSampleBarycoord(uAttrStatic, bary, fi.xyz);
+    vec3 xN = normalize(attr.xyz);
+    if (dot(xN, rd) > 0.0) xN = -xN;
+    vec3 xP = ro + rd * dist;
+    vec3 rd2 = refract(rd, xN, uIor);
+    if (rd2 == vec3(0.0)) rd2 = reflect(rd, xN);
+    refrRad = traceRadiance(xP - xN * uEps, rd2, true);
+  } else {
+    refrRad = uSkyEnabled
+      ? skyColor(rd, uSunDir, uSunColor, uSkyZenith, uSkyHorizon, uSkyIntensity)
+      : uEnvColor * uEnvIntensity;
+  }
+  return mix(refrRad, reflRad, fres);
+}
+
+void main() {
+  vec4 wp = texture(uGWorldPos, vUv);
+  if (wp.w < 0.5) {
+    outIrradiance = vec4(0.0);
+    return;
+  }
+
+  ivec2 px = ivec2(gl_FragCoord.xy);
+  gSeed = uint(px.x) * 1973u + uint(px.y) * 9277u + uint(uFrame) * 26699u;
+  gSeed = pcgHash(gSeed);
+  gBlueNoise = fetchBlueNoise();
+  gBnDim = 0;
+
+  vec3 P = wp.xyz;
+  vec4 nmSample = texture(uGNormalMetal, vUv);
+  vec3 N = normalize(nmSample.xyz);
+  // Decode the packed material word (see GBufferPass): >= 2 → glass, else metal.
+  float transmission = nmSample.w >= 2.0 ? clamp(nmSample.w - 2.0, 0.0, 1.0) : 0.0;
+  float metal = nmSample.w >= 2.0 ? 0.0 : nmSample.w;
+  float rough = clamp(wp.w - 1.0, 0.0, 1.0);
+
+  // --- direct lighting ---
+  // ReSTIR: shade the reservoir's winner with one visibility ray (flat cost in
+  // light count). Stochastic: one blind random sample. Full: one shadow ray
+  // per light + one for the emissive set.
+  vec3 direct = vec3(0.0);
+  if (uRestirEnabled) {
+    direct = shadeReservoir(P, N);
+  } else if (uLightStochastic) {
+    direct = sampleOneAny(P, N);
+  } else {
+    for (int i = 0; i < MAX_LIGHTS; i++) {
+      if (i >= uLightCount) break;
+      direct += lightContribution(i, P, N);
+    }
+    // Emissive meshes as area lights (next-event estimation, one shadow ray).
+    direct += sampleEmissiveTri(P, N);
+  }
+
+  // --- 1-bounce indirect (cosine-weighted; pdf cancels the NdotL/PI).
+  // traceRadiance shades the hit with direct + NEE light, or returns the
+  // sky/env colour when the ray escapes (the natural ambient bounce).
+  // Half-rate mode traces on alternating checkerboard parity each frame,
+  // DOUBLED — the temporal average converges to the same brightness
+  // (unbiased) while GI's ray cost halves; accumulation + denoise absorb
+  // the alternation.
+  vec3 indirect = vec3(0.0);
+  if (uGIEnabled) {
+    bool trace = !uGIHalfRate || (((px.x + px.y + int(uFrame)) & 1) == 0);
+    if (trace) {
+      indirect = traceRadiance(P + N * uEps, cosineSampleHemisphere(N, rand2()), false);
+      if (uGIHalfRate) indirect *= 2.0;
+    }
+  }
+
+  // Firefly clamp: suppress rare huge GI samples (big perceived-noise win,
+  // slightly biased). Applied to indirect only; direct is analytic.
+  float lum = dot(indirect, vec3(0.299, 0.587, 0.114));
+  if (lum > uFireflyClamp) indirect *= uFireflyClamp / lum;
+
+  vec3 sampleIrr = direct + indirect;
+
+  // --- traced specular: mirror/glossy reflections on metals ---
+  if (uReflEnabled && metal > 0.001) {
+    vec3 V = normalize(P - uCameraPos);
+    vec3 refl = glossyReflect(V, N, rough);
+    if (dot(refl, N) > 0.0) {
+      // Metals have no diffuse term: replace by metalness. The composite's
+      // albedo multiply then tints the reflection (F0 = albedo for metals).
+      sampleIrr = mix(sampleIrr, traceRadiance(P + N * uEps, refl, true), metal);
+    }
+  }
+
+  // --- traced glass: Fresnel reflection + two-interface refraction ---
+  if (uRefrEnabled && transmission > 0.001) {
+    vec3 V = normalize(P - uCameraPos);
+    sampleIrr = mix(sampleIrr, glassRadiance(P, N, V, rough), transmission);
+  }
+
+  // A single NaN/Inf sample would poison the EMA history for good (mix() with
+  // NaN stays NaN until a disocclusion resets the pixel) — sanitize first.
+  if (any(isnan(sampleIrr)) || any(isinf(sampleIrr))) sampleIrr = vec3(0.0);
+
+  // --- temporal reprojection: pull validated history from last frame ---
+  float count = 1.0;
+  vec3 history = vec3(0.0);
+  if (uTemporalReprojection) {
+    vec4 clip = uPrevViewProj * vec4(P, 1.0);
+    vec4 clipC = uViewProj * vec4(P, 1.0);
+    if (clip.w > 0.0 && clipC.w > 0.0) {
+      vec2 prevUv = (clip.xy / clip.w) * 0.5 + 0.5;
+      // P comes from a full-res G-buffer texel, which sits sub-pixel off this
+      // half-res fragment's center. That constant offset would bias bilinear
+      // history reads every frame (content drifts/smears at renderScale < 1).
+      // Cancel it: measure P's offset in the CURRENT frame and subtract.
+      vec2 currUv = (clipC.xy / clipC.w) * 0.5 + 0.5;
+      prevUv -= currUv - vUv;
+      if (prevUv.x >= 0.0 && prevUv.x <= 1.0 && prevUv.y >= 0.0 && prevUv.y <= 1.0) {
+        vec4 prevPos = texture(uPrevGWorldPos, prevUv);
+        // Plane-distance test: robust at grazing angles (position error from
+        // texel quantization lies along the surface, not along the normal).
+        float distToCam = distance(P, uCameraPos);
+        float tol = 0.005 * distToCam + 20.0 * uEps;
+        bool valid = prevPos.w > 0.5
+          && abs(dot(P - prevPos.xyz, N)) < tol;
+        if (valid) {
+          vec4 h = texture(uPrevAccum, prevUv); // bilinear
+          // Mirror-like pixels keep a SHORT history: their reflected content
+          // moves differently from the surface, so long history smears the
+          // reflection under camera motion — and specular rays are nearly
+          // deterministic, so they don't need the accumulation anyway.
+          float specHist = max(metal, transmission) * (1.0 - rough);
+          float histCap = mix(uMaxHistory, min(uMaxHistory, 10.0), specHist);
+          count = clamp(h.a, 0.0, histCap) + 1.0;
+          history = h.rgb;
+        }
+      }
+    }
+  }
+
+  // Exponential moving average; count=1 (disocclusion / first frame) means
+  // the fresh sample is used as-is.
+  vec3 blended = mix(history, sampleIrr, 1.0 / count);
+  outIrradiance = vec4(blended, count);
+}
+`;
+
+/**
+ * Fullscreen pass: for every G-buffer pixel, trace shadow rays to every light and
+ * one cosine-weighted GI bounce against the BVH. Outputs demodulated irradiance,
+ * progressively accumulated into a ping-pong float target while the camera is still.
+ */
+export class RTLightingPass {
+  constructor(width, height) {
+    this.targetA = this._makeTarget(width, height);
+    this.targetB = this._makeTarget(width, height);
+
+    this.material = new THREE.ShaderMaterial({
+      glslVersion: THREE.GLSL3,
+      vertexShader: fullscreenVert,
+      fragmentShader: rtLightingFrag,
+      uniforms: {
+        bvhStatic: { value: null },
+        bvhDynamic: { value: null },
+        uHasDynamic: { value: false },
+        uAttrStatic: { value: null },
+        uAttrDynamic: { value: null },
+        uMaterialsTex: { value: null },
+        uGWorldPos: { value: null },
+        uGNormalMetal: { value: null },
+        uPrevAccum: { value: null },
+        uPrevGWorldPos: { value: null },
+        uReservoir: { value: null },
+        uRestirEnabled: { value: false },
+        uPrevViewProj: { value: new THREE.Matrix4() },
+        uViewProj: { value: new THREE.Matrix4() },
+        uCameraPos: { value: new THREE.Vector3() },
+        uMaxHistory: { value: 128 },
+        uTemporalReprojection: { value: true },
+        uFireflyClamp: { value: 4.0 },
+        uLightPosType: { value: [] },
+        uLightColorRadius: { value: [] },
+        uLightDirCone: { value: [] },
+        uLightCount: { value: 0 },
+        uEmissiveCount: { value: 0 },
+        uReflEnabled: { value: true },
+        uRefrEnabled: { value: true },
+        uIor: { value: 1.5 },
+        uLightStochastic: { value: false },
+        uGIHalfRate: { value: false },
+        uEnvColor: { value: new THREE.Color(0.03, 0.04, 0.06) },
+        uEnvIntensity: { value: 1.0 },
+        uFrame: { value: 0 },
+        uEps: { value: 1e-3 },
+        uGIEnabled: { value: true },
+        uSkyEnabled: { value: false },
+        uSunDir: { value: new THREE.Vector3(0.4, 0.8, 0.45).normalize() },
+        uSunColor: { value: new THREE.Color(1.0, 0.9, 0.75) },
+        uSkyZenith: { value: new THREE.Color(0.18, 0.34, 0.62) },
+        uSkyHorizon: { value: new THREE.Color(0.7, 0.8, 0.9) },
+        uSkyIntensity: { value: 1.0 },
+      },
+      depthTest: false,
+      depthWrite: false,
+    });
+
+    this.scene = new THREE.Scene();
+    this.camera = new THREE.OrthographicCamera(-1, 1, 1, -1, 0, 1);
+    this.quad = new THREE.Mesh(new THREE.PlaneGeometry(2, 2), this.material);
+    this.quad.frustumCulled = false;
+    this.scene.add(this.quad);
+  }
+
+  _makeTarget(width, height) {
+    // Half-float + linear: history is sampled bilinearly at reprojected UVs,
+    // and fp16 halves the bandwidth (EMA blending never accumulates a raw sum,
+    // so fp16 precision is sufficient).
+    const t = new THREE.WebGLRenderTarget(width, height, {
+      minFilter: THREE.LinearFilter,
+      magFilter: THREE.LinearFilter,
+      format: THREE.RGBAFormat,
+      type: THREE.HalfFloatType,
+      depthBuffer: false,
+      stencilBuffer: false,
+    });
+    t.texture.generateMipmaps = false;
+    return t;
+  }
+
+  /** Clear both history buffers (e.g. after scene recompile or resize). */
+  clearHistory(renderer) {
+    const prevTarget = renderer.getRenderTarget();
+    const prevColor = new THREE.Color();
+    renderer.getClearColor(prevColor);
+    const prevAlpha = renderer.getClearAlpha();
+    renderer.setClearColor(0x000000, 0);
+    for (const t of [this.targetA, this.targetB]) {
+      renderer.setRenderTarget(t);
+      renderer.clear(true, false, false);
+    }
+    renderer.setRenderTarget(prevTarget);
+    renderer.setClearColor(prevColor, prevAlpha);
+  }
+
+  setSize(width, height) {
+    this.targetA.setSize(width, height);
+    this.targetB.setSize(width, height);
+  }
+
+  /**
+   * Reallocate the history targets to a new size while PRESERVING the
+   * accumulated irradiance. The plain setSize + clearHistory path dumps every
+   * temporal sample, which strobes the image back to 1-spp noise on every
+   * governor renderScale step — this carries the history over instead.
+   *
+   * The freshest history is targetB (last frame's output — see the swap in
+   * render()); it is resampled through copyPass into the new targetB, its
+   * per-pixel sample count (alpha) clamped to `carryFrames` so the EMA
+   * reconverges smoothly at the new resolution rather than freezing on stale
+   * values. targetA is overwritten on the next render, so it only needs the
+   * fresh allocation, not a copy.
+   */
+  resizeCarry(renderer, copyPass, width, height, carryFrames) {
+    const newA = this._makeTarget(width, height);
+    const newB = this._makeTarget(width, height);
+    copyPass.blit(renderer, this.targetB.texture, newB, carryFrames);
+    this.targetA.dispose();
+    this.targetB.dispose();
+    this.targetA = newA;
+    this.targetB = newB;
+  }
+
+  setCompiledScene(compiled) {
+    const u = this.material.uniforms;
+    u.bvhStatic.value = compiled.staticBvhUniform;
+    u.bvhDynamic.value = compiled.dynamicBvhUniform;
+    u.uHasDynamic.value = compiled.hasDynamic;
+    u.uAttrStatic.value = compiled.staticAttrTex;
+    u.uAttrDynamic.value = compiled.dynamicAttrTex;
+    u.uMaterialsTex.value = compiled.materialsTex;
+    u.uLightPosType.value = compiled.lightPosType;
+    u.uLightColorRadius.value = compiled.lightColorRadius;
+    u.uLightDirCone.value = compiled.lightDirCone;
+    u.uLightCount.value = compiled.lightCount;
+    u.uEmissiveCount.value = compiled.emissiveTriCount;
+  }
+
+  /** Renders into targetA (reading targetB as history), then swaps. Returns the fresh accum texture. */
+  render(renderer, gbuffer, frame, reservoirTexture = null) {
+    const u = this.material.uniforms;
+    u.uGWorldPos.value = gbuffer.worldPos;
+    u.uGNormalMetal.value = gbuffer.normalMetal;
+    u.uPrevGWorldPos.value = gbuffer.prevWorldPos;
+    u.uPrevAccum.value = this.targetB.texture;
+    u.uReservoir.value = reservoirTexture;
+    u.uRestirEnabled.value = reservoirTexture !== null;
+    u.uFrame.value = frame;
+
+    renderer.setRenderTarget(this.targetA);
+    renderer.render(this.scene, this.camera);
+    renderer.setRenderTarget(null);
+
+    const out = this.targetA;
+    [this.targetA, this.targetB] = [this.targetB, this.targetA];
+    return out.texture;
+  }
+
+  dispose() {
+    this.targetA.dispose();
+    this.targetB.dispose();
+    this.material.dispose();
+    this.quad.geometry.dispose();
+  }
+}