@woosh/meep-engine 3.59.0 → 3.60.1

package/src/shade/renderer/global_illumination/probe_clipmap/probe/shader_probe_clipmap_light.js

@@ -1,34 +1,50 @@
 import { MESH_SDF_BRICK_VOXELS } from "../../../geometry/sdf/MESH_SDF.js";
 import { chunk_light_cluster_counts } from "../../../light/cluster/chunk_light_cluster_counts.js";
 import { LIGHT_DATA_RESOURCE } from "../../../light/LIGHT_DATA_RESOURCE.js";
+import { LIGHT_DATABASE_SPEC } from "../../../light/LIGHT_DATABASE_SPEC.js";
+import { LIGHT_FLAGS } from "../../../light/model/io/LIGHT_FLAGS.js";
 import {
     chunk_world_light_cluster_lookup_at,
     WORLD_CLUSTER_DATA_RESOURCE,
     WORLD_CLUSTER_LOOKUP_RESOURCE
 } from "../../../light/world_cluster/chunk_world_light_cluster_lookup.js";
+import { chunk_rgb_to_luminance } from "../../../shader/chunk/color/chunk_rgb_to_luminance.js";
+import { chunk_rgbe9995_decode } from "../../../shader/chunk/color/rgbe9995/chunk_rgbe9995_decode.js";
+import { chunk_rgbe9995_encode } from "../../../shader/chunk/color/rgbe9995/chunk_rgbe9995_encode.js";
 import { chunk_uv_octahedral_unit_decode } from "../../../shader/chunk/geometry/octahedron/chunk_uv_octahedral_unit_decode.js";
+import { chunk_sphere_sample_volume } from "../../../shader/chunk/geometry/sphere/chunk_sphere_sample_volume.js";
 import { chunk_folded_workgroup_index } from "../../../shader/chunk/grid/chunk_folded_workgroup_index.js";
 import { RECIPROCAL_PI_DECLARATION } from "../../../shader/chunk/math/RECIPROCAL_PI_DECLARATION.js";
+import { chunk_hash_to_float01 } from "../../../shader/chunk/random/chunk_hash_to_float01.js";
+import { chunk_pcg4d } from "../../../shader/chunk/random/pcg/chunk_pcg4d.js";
 import {
     chunk_resolve_irradiance_sh_to_diffuse
 } from "../../../shader/chunk/spherical_harmonics/l1i/chunk_resolve_irradiance_sh_to_diffuse.js";
 import { CodeChunk } from "../../../shader/compiler/CodeChunk.js";
 import { ComputeShader } from "../../../shader/ComputeShader.js";
 import { ShaderResourceSetDescriptor } from "../../../shader/resource/ShaderResourceSetDescriptor.js";
-import { SHADOWMAP_ATLAS_RESOURCE } from "../../../shadow/map/SHADOWMAP_RESOURCES.js";
 import { STBN_VEC2_RESOURCE } from "../../../STATIC_GRAPHICS_ENGINE_ASSET_RESOURCES.js";
-import { chunk_volumetrics_local_light_shadow } from "../../../volumetrics/chunk_volumetrics_local_light_shadow.js";
 import { chunk_irradiance_clipmap_local_light } from "../../../volumetrics/irradiance_clipmap/chunk_irradiance_clipmap_local_light.js";
 import { IRRADIANCE_CLIPMAP_LIGHT_BUDGET } from "../../../volumetrics/irradiance_clipmap/IRRADIANCE_CLIPMAP.js";
 import { BRICK_TRACE_HIT } from "../brick/BRICK_CLIPMAP.js";
-import { BRICK_LEVELS_RESOURCE } from "../brick/chunk_brick_clipmap_trace.js";
+import {
+    BRICK_ATLAS_TEXTURE_RESOURCE,
+    BRICK_LEVELS_RESOURCE,
+    BRICK_MAP_TEXTURE_RESOURCE,
+    BRICK_OCCUPANCY_TEXTURE_RESOURCE,
+    BRICK_SAMPLER_RESOURCE,
+    chunk_brick_clipmap_trace
+} from "../brick/chunk_brick_clipmap_trace.js";
 import {
     chunk_probe_clipmap,
     chunk_probe_clipmap_ray_jitter,
+    chunk_probe_clipmap_ray_layout,
+    chunk_probe_clipmap_ray_origin,
     chunk_probe_clipmap_ray_texel,
     PROBE_DEPTH_HISTORY_RESOURCE,
     PROBE_LEVELS_RESOURCE,
     PROBE_ORIGINS_READ_RESOURCE,
+    PROBE_RAY_SUMS_READ_RESOURCE,
     PROBE_RAYS_WRITE_RESOURCE,
     PROBE_RECORDS_READ_RESOURCE,
     PROBE_STATES_READ_RESOURCE,
@@ -43,7 +59,8 @@ import {
 import { PROBE_THREAD_WORKGROUP_SIZE } from "./PROBE_CLIPMAP.js";
 import { PROBE_TRACE_SETTINGS } from "./shader_probe_clipmap_trace.js";
 
-const HALF_MAX = 65504;
+const POINT = LIGHT_DATABASE_SPEC.get('light_point');
+const SPOT = LIGHT_DATABASE_SPEC.get('light_spot');
 
 const resources = new ShaderResourceSetDescriptor();
 
@@ -53,6 +70,7 @@ resources.createGroup()
     .add(PROBE_TRACE_READ_RESOURCE)
     .add(PROBE_ORIGINS_READ_RESOURCE)
     .add(PROBE_RAYS_WRITE_RESOURCE)
+    .add(PROBE_RAY_SUMS_READ_RESOURCE)
     .add(PROBE_STATES_READ_RESOURCE)
     .add(STBN_VEC2_RESOURCE);
 
@@ -66,21 +84,35 @@ resources.createGroup()
 
 resources.createGroup()
     .add(BRICK_LEVELS_RESOURCE)
+    .add(BRICK_MAP_TEXTURE_RESOURCE)
+    .add(BRICK_OCCUPANCY_TEXTURE_RESOURCE)
+    .add(BRICK_ATLAS_TEXTURE_RESOURCE)
+    .add(BRICK_SAMPLER_RESOURCE)
     .add(LIGHT_DATA_RESOURCE)
     .add(WORLD_CLUSTER_LOOKUP_RESOURCE)
-    .add(WORLD_CLUSTER_DATA_RESOURCE)
-    .add(SHADOWMAP_ATLAS_RESOURCE);
+    .add(WORLD_CLUSTER_DATA_RESOURCE);
 
 /**
  * The light a front-face hit takes from beyond the trace (PROBE_CLIPMAP_PROPOSAL_2026_10_27.md §8.7): a
- * thread a ray, after shade through the same launch, `PROBE_DISPATCH_RAYS`. It adds ρ_d·(E_loc + E_ind)/π
- * to the radiance shade left, ρ_d being the albedo shade left in the ray's \`row\`.
+ * thread a ray, after shade through the same launch, `PROBE_DISPATCH_RAYS`, its probe the one trace wrote
+ * into its flags. It adds ρ_d·(E_loc + E_ind)/π to the radiance shade left, ρ_d being the albedo shade
+ * left in the ray's \`row\`.
  *
- * - **E_loc, the local lights**: the point and spot lights of the hit's cell of the world light clusters,
- *   built over the bricks' cascades so that every hit has one, at most `IRRADIANCE_CLIPMAP_LIGHT_BUDGET` of
- *   each; past them lights are dropped. Each is the irradiance clipmap's, at the footprint of a voxel of
- *   the hit's cascade, its cosine at the hit's normal, and its shadow map's nearest texel from a voxel off
- *   the surface, past the field's half-voxel expand. A light whose map is cleared lights hits unshadowed.
+ * - **E_loc, the local lights**: one of the point and spot lights of the hit's cell of the world light
+ *   clusters, built over the bricks' cascades so that every hit has one, at most
+ *   `IRRADIANCE_CLIPMAP_LIGHT_BUDGET` of each; past them lights are dropped. Each light's unshadowed
+ *   irradiance is the irradiance clipmap's, at the footprint of a voxel of the hit's cascade, by its
+ *   cosine at the hit's normal.
+ *   - One is picked by weighted reservoir sampling over the cell's lights, its weight the luminance of
+ *     that irradiance, with one number stretched over the candidates as `sample_discrete_wrs` does; it
+ *     lights the hit by its irradiance over the chance of its pick, w/W. The probe's blend of its rays
+ *     averages the picks to the cell's sum.
+ *   - The picked light's shadow is a ray traced through the brick clipmap, as the hit's own ray was:
+ *     from a voxel off the surface, past the field's half-voxel expand, toward a point of the light's
+ *     sphere drawn as the path tracer draws it, stopping short of that point by the light's near clip or
+ *     a voxel, whichever is more, for the expand at that end. It sees the light if it misses; a hit, or a
+ *     march that spends its steps, is no light. A light that casts no shadow traces none.
+ *   - The pick and the point are drawn from the ray's slot and direction, fresh every update.
  * - **E_ind, the next bounce**: this mode's own diffuse at the hit, `probe_clipmap_sample_sh` from the
  *   hit along its normal, the view back along the ray, with the ray's own jitter as the pick's noise.
  *   The projection writes the records later in the frame, so they hold last frame's SH: each frame takes
@@ -100,47 +132,112 @@ export const shader_probe_clipmap_light = ComputeShader.from({
     body: CodeChunk.from(
         //language=WGSL
         `
-// E_loc at a front-face hit in brick cascade \`cascade\`: its world cluster cell's point and spot lights
-fn probe_clipmap_local_irradiance(cascade: u32, position: vec3<f32>, normal: vec3<f32>) -> vec3<f32> {
+// E_loc at a front-face hit in brick cascade \`cascade\`: one of its world cluster cell's point and spot lights,
+// picked by \`random.x\` in proportion to its light at the hit and over the chance of its pick, shadowed by a ray
+// toward the point of the light's sphere \`random.yzw\` draws
+fn probe_clipmap_local_irradiance(cascade: u32, position: vec3<f32>, normal: vec3<f32>, random: vec4<f32>) -> vec3<f32> {
     let level = brick_levels[cascade];
     let voxel = level.cell_size / ${MESH_SDF_BRICK_VOXELS}.0;
     let cell = world_light_cluster_metadata_at(cascade, level, 0.0, position);
     let point_count = min(light_cluster_count_point(cell.counts), ${IRRADIANCE_CLIPMAP_LIGHT_BUDGET}u);
     let spot_count = min(light_cluster_count_spot(cell.counts), ${IRRADIANCE_CLIPMAP_LIGHT_BUDGET}u);
 
-    // the maps' receiver: a voxel off the surface, past the half voxel the field's surface may sit off it
-    let lifted = position + normal * voxel;
+    // the cell's spot lights follow its point lights
+    let spots = cell.offset + light_cluster_count_point(cell.counts);
+
+    // the reservoir: the weights so far, and the light it holds
+    var u = random.x;
+    var total = 0.0;
+    var picked = vec3<f32>(0.0);
+    var picked_weight = 0.0;
+    var picked_index = 0u;
+    var picked_spot = false;
 
-    var irradiance = vec3<f32>(0.0);
+    for (var j = 0u; j < point_count + spot_count; j++) {
+        let spot = j >= point_count;
+        let index = world_cluster_data[select(cell.offset + j, spots + (j - point_count), spot)];
 
-    for (var j = 0u; j < point_count; j++) {
-        let index = world_cluster_data[cell.offset + j];
-        let light = irradiance_clipmap_get_point_light_info_by_index(&light_data, index, position, voxel);
-        let cosine = max(dot(normal, light.direction), 0.0);
+        var irradiance: vec3<f32>;
+
+        if (spot) {
+            let light = irradiance_clipmap_get_spot_light_info_by_index(&light_data, index, position, voxel);
+
+            irradiance = light.color * max(dot(normal, light.direction), 0.0);
+        } else {
+            let light = irradiance_clipmap_get_point_light_info_by_index(&light_data, index, position, voxel);
+
+            irradiance = light.color * max(dot(normal, light.direction), 0.0);
+        }
 
-        if (cosine == 0.0 || all(light.color == vec3<f32>(0.0))) {
+        let weight = rgb_to_luminance(irradiance);
+
+        if (weight <= 0.0) {
             continue;
         }
 
-        irradiance += light.color * (cosine * volumetrics_point_light_shadow(&light_data, index, lifted));
+        // the light takes the reservoir by its share of the weights so far, and the number is stretched back
+        // over [0, 1) from whichever side of the share it fell
+        total += weight;
+
+        let share = weight / total;
+
+        if (u < share) {
+            u = u / share;
+            picked = irradiance;
+            picked_weight = weight;
+            picked_index = index;
+            picked_spot = spot;
+        } else {
+            u = (u - share) / (1.0 - share);
+        }
     }
 
-    // the cell's spot lights follow its point lights
-    let spots = cell.offset + light_cluster_count_point(cell.counts);
+    if (total <= 0.0) {
+        return vec3<f32>(0.0);
+    }
 
-    for (var k = 0u; k < spot_count; k++) {
-        let index = world_cluster_data[spots + k];
-        let light = irradiance_clipmap_get_spot_light_info_by_index(&light_data, index, position, voxel);
-        let cosine = max(dot(normal, light.direction), 0.0);
+    // the pick's light over the chance of the pick, w / W
+    let estimate = picked * (total / picked_weight);
 
-        if (cosine == 0.0 || all(light.color == vec3<f32>(0.0))) {
-            continue;
-        }
+    var source: vec3<f32>;
+    var radius: f32;
+    var near_clip: f32;
+    var flags: u32;
+
+    if (picked_spot) {
+        let light = ${SPOT.marshalling_method_read}(&light_data, picked_index);
+
+        source = light.position;
+        radius = light.radius;
+        near_clip = light.near_clip_distance;
+        flags = light.flags;
+    } else {
+        let light = ${POINT.marshalling_method_read}(&light_data, picked_index);
+
+        source = light.position.xyz;
+        radius = light.radius;
+        near_clip = light.near_clip_distance;
+        flags = light.flags;
+    }
+
+    if ((flags & ${LIGHT_FLAGS.CASTS_SHADOW_BIT}u) == 0u) {
+        return estimate;
+    }
+
+    // from a voxel off the surface, past the half voxel the field's surface may sit off it, to a point of the
+    // light, short of it by the light's near clip or by a voxel for the expand at that end
+    let lifted = position + normal * voxel;
+    let to_light = source + radius * sphere_sample_volume(random.yzw) - lifted;
+    let span = length(to_light);
+    let reach = span - max(near_clip, voxel);
 
-        irradiance += light.color * (cosine * volumetrics_spot_light_shadow(&light_data, index, lifted));
+    if (reach <= 0.0) {
+        return estimate;
     }
 
-    return irradiance;
+    let shadow = brick_clipmap_trace(lifted, to_light / span, 0.0, reach, settings.brick_cascade_count);
+
+    return select(vec3<f32>(0.0), estimate, shadow.kind == BRICK_TRACE_MISS);
 }
 
 @compute @workgroup_size(${PROBE_THREAD_WORKGROUP_SIZE})
@@ -149,12 +246,6 @@ fn main(
     @builtin(local_invocation_index) lane: u32,
 ) {
     let ray = folded_workgroup_index(workgroup_id) * PROBE_THREAD_WORKGROUP_SIZE + lane;
-    let rays = settings.rays_per_probe;
-
-    if (ray >= min(probe_trace.count, settings.trace_capacity) * rays) {
-        return;
-    }
-
     let record = probe_rays[ray];
     let albedo = unpack4x8unorm(record.row).rgb;
     let front_hit = (record.flags & PROBE_RAY_KIND_MASK) == ${BRICK_TRACE_HIT}u && (record.flags & PROBE_RAY_BACK_FACE) == 0u;
@@ -164,18 +255,27 @@ fn main(
         return;
     }
 
-    let slot = probe_trace.elements[ray / rays];
-    let i = ray % rays;
+    // the probe's update's rays, and this one's place among them
+    let k = record.flags >> PROBE_RAY_PROBE_SHIFT;
+    let span = probe_clipmap_probe_rays(k);
+    let rays = span.y;
+    let i = ray - span.x;
+
+    let slot = probe_trace.elements[k];
     let n = settings.radiance_resolution;
     let updates = probe_state_updates(probe_states[slot].x);
     let texel = probe_clipmap_ray_texel(i, updates, rays, n * n);
 
     let direction = uv_octahedral_unit_decode(unpack2x16unorm(record.direction));
-    let position = probe_origins[ray / rays].xyz + record.t * direction;
+    let position = probe_clipmap_ray_origin(probe_origins[k], slot, record.direction) + record.t * direction;
     let normal = uv_octahedral_unit_decode(unpack2x16unorm(record.normal));
 
+    // the local lights' pick and their shadow ray's point, apart from the origin's draw
+    let hash = pcg4d(vec4<u32>(slot, record.direction, 0xA54FF53Au, 0u));
+    let random = vec4<f32>(hash_to_float01(hash.x), hash_to_float01(hash.y), hash_to_float01(hash.z), hash_to_float01(hash.w));
+
     // E/π: the local lights', then the next bounce's, which the resolve gives as E/π
-    var e = probe_clipmap_local_irradiance(record.flags >> PROBE_RAY_CASCADE_SHIFT, position, normal) * RECIPROCAL_PI;
+    var e = probe_clipmap_local_irradiance((record.flags >> PROBE_RAY_CASCADE_SHIFT) & PROBE_RAY_CASCADE_MASK, position, normal, random) * RECIPROCAL_PI;
 
     let probes = probe_clipmap_sample_sh(position, normal, -direction, probe_clipmap_ray_jitter(slot, texel, i, updates, rays, n));
 
@@ -187,16 +287,14 @@ fn main(
         return;
     }
 
-    let light = vec3<f32>(unpack2x16float(record.light_rg), unpack2x16float(record.light_b).x) + albedo * e;
-    let held = min(light, vec3<f32>(${HALF_MAX}.0));
-
-    probe_rays[ray].light_rg = pack2x16float(held.rg);
-    probe_rays[ray].light_b = pack2x16float(vec2<f32>(held.b, 0.0));
+    probe_rays[ray].light = rgbe9995_encode(rgbe9995_decode(record.light) + albedo * e);
 }
 `, [
             chunk_probe_clipmap,
+            chunk_probe_clipmap_ray_layout,
             chunk_probe_clipmap_ray_texel,
             chunk_probe_clipmap_ray_jitter,
+            chunk_probe_clipmap_ray_origin,
             chunk_folded_workgroup_index,
             chunk_uv_octahedral_unit_decode,
             chunk_probe_clipmap_pick,
@@ -205,6 +303,14 @@ fn main(
             chunk_world_light_cluster_lookup_at,
             chunk_light_cluster_counts,
             chunk_irradiance_clipmap_local_light,
-            chunk_volumetrics_local_light_shadow,
+            chunk_rgb_to_luminance,
+            chunk_brick_clipmap_trace,
+            chunk_pcg4d,
+            chunk_hash_to_float01,
+            chunk_sphere_sample_volume,
+            POINT.chunk_read,
+            SPOT.chunk_read,
+            chunk_rgbe9995_decode,
+            chunk_rgbe9995_encode,
         ]),
 });

package/src/shade/renderer/global_illumination/probe_clipmap/probe/shader_probe_clipmap_order.d.ts

@@ -0,0 +1,23 @@
+/**
+ * @type {WebGPUStruct}
+ */
+export const PROBE_ORDER_SETTINGS: WebGPUStruct;
+/**
+ * The trace list in the probes' order (§8.6): a workgroup a run of `PROBE_ORDER_RUN` listed
+ * probes, dispatched through `PROBE_DISPATCH_ORDER` once the due passes have listed the frame's updates,
+ * sorts its run in place by each probe's cascade and the Morton code of its lattice point in the
+ * cascade's window, the cascade above. Every pass after it reads a probe's place in the list from the
+ * sorted list, the rays' layout among them, so the rays of neighbouring probes run in neighbouring
+ * waves. Past one run the list is ordered within each run, not across them.
+ *
+ * Measured (PROBE_RAYS_VARIABLE_AND_SORTED_SCOPE_2026_10_30.md): the rays' stretch 1–4% shorter on
+ * the iGPU at 131k rays a frame, level on the 4090. The order is the only one of the rays' orders
+ * measured that pays: every sort of the rays themselves, global or within chunks, by direction, by
+ * place or by both, costs more than it buys, because the rays of one probe already share their waves.
+ *
+ * @type {ComputeShader}
+ */
+export const shader_probe_clipmap_order: ComputeShader;
+import { WebGPUStruct } from "../../../shader/type/WebGPUStruct.js";
+import { ComputeShader } from "../../../shader/ComputeShader.js";
+//# sourceMappingURL=shader_probe_clipmap_order.d.ts.map

package/src/shade/renderer/global_illumination/probe_clipmap/probe/shader_probe_clipmap_order.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"shader_probe_clipmap_order.d.ts","sourceRoot":"","sources":["../../../../../../../src/shade/renderer/global_illumination/probe_clipmap/probe/shader_probe_clipmap_order.js"],"names":[],"mappings":"AAeA;;GAEG;AACH,mCAFU,YAAY,CASG;AAWzB;;;;;;;;;;;;;;GAcG;AACH,yCAFU,aAAa,CA4DpB;6BArG0B,sCAAsC;8BAHrC,kCAAkC"}

package/src/shade/renderer/global_illumination/probe_clipmap/probe/shader_probe_clipmap_order.js

@@ -0,0 +1,110 @@
+import { chunk_sort_bitonic_workgroup } from "../../../gpu_primitive/sort/chunk_sort_bitonic_workgroup.js";
+import { SORT_KEYS_PER_THREAD, SORT_TILE_SIZE, SORT_WORKGROUP_SIZE } from "../../../gpu_primitive/sort/defs.js";
+import { chunk_v3_morton_encode } from "../../../shader/chunk/bvh/morton/chunk_v3_morton_encode.js";
+import { chunk_folded_workgroup_index } from "../../../shader/chunk/grid/chunk_folded_workgroup_index.js";
+import { CodeChunk } from "../../../shader/compiler/CodeChunk.js";
+import { ComputeShader } from "../../../shader/ComputeShader.js";
+import { ShaderResourceSetDescriptor } from "../../../shader/resource/ShaderResourceSetDescriptor.js";
+import { u32 } from "../../../shader/type/WEBGPU_WGSL_PRIMITIVE_TYPES.js";
+import { WebGPUStruct } from "../../../shader/type/WebGPUStruct.js";
+import { assert } from "../../../../../core/assert.js";
+import { chunk_probe_clipmap, PROBE_LEVELS_RESOURCE, PROBE_RECORDS_READ_RESOURCE, PROBE_TRACE_ORDER_RESOURCE } from "./chunk_probe_clipmap.js";
+import { PROBE_ORDER_RUN } from "./PROBE_CLIPMAP.js";
+
+assert.equal(PROBE_ORDER_RUN, SORT_TILE_SIZE, 'PROBE_ORDER_RUN is a tile of the bitonic sort');
+
+/**
+ * @type {WebGPUStruct}
+ */
+export const PROBE_ORDER_SETTINGS = WebGPUStruct.from({
+    /** probes a frame updates at most: the trace list's bound */
+    trace_capacity: u32,
+    /** lattice points a side of a cascade's window */
+    resolution: u32,
+    /** bits an axis of the Morton code, ⌈log2 resolution⌉ */
+    axis_bits: u32,
+}, "ProbeOrderSettings");
+
+const resources = new ShaderResourceSetDescriptor();
+
+resources.setImmediate("settings", PROBE_ORDER_SETTINGS);
+
+resources.createGroup()
+    .add(PROBE_TRACE_ORDER_RESOURCE)
+    .add(PROBE_RECORDS_READ_RESOURCE)
+    .add(PROBE_LEVELS_RESOURCE);
+
+/**
+ * The trace list in the probes' order (§8.6): a workgroup a run of `PROBE_ORDER_RUN` listed
+ * probes, dispatched through `PROBE_DISPATCH_ORDER` once the due passes have listed the frame's updates,
+ * sorts its run in place by each probe's cascade and the Morton code of its lattice point in the
+ * cascade's window, the cascade above. Every pass after it reads a probe's place in the list from the
+ * sorted list, the rays' layout among them, so the rays of neighbouring probes run in neighbouring
+ * waves. Past one run the list is ordered within each run, not across them.
+ *
+ * Measured (PROBE_RAYS_VARIABLE_AND_SORTED_SCOPE_2026_10_30.md): the rays' stretch 1–4% shorter on
+ * the iGPU at 131k rays a frame, level on the 4090. The order is the only one of the rays' orders
+ * measured that pays: every sort of the rays themselves, global or within chunks, by direction, by
+ * place or by both, costs more than it buys, because the rays of one probe already share their waves.
+ *
+ * @type {ComputeShader}
+ */
+export const shader_probe_clipmap_order = ComputeShader.from({
+    label: "Probe Clipmap Order",
+    resources,
+    body: CodeChunk.from(
+        //language=WGSL
+        `
+// the run's keys and the slots they carry, which the bitonic sort orders by name
+var<workgroup> keys: array<u32, ${SORT_TILE_SIZE}>;
+var<workgroup> values: array<u32, ${SORT_TILE_SIZE}>;
+
+var<workgroup> wg_count: u32;
+
+@compute @workgroup_size(${SORT_WORKGROUP_SIZE}, 1, 1)
+fn main(
+    @builtin(workgroup_id) workgroup_id: vec3<u32>,
+    @builtin(local_invocation_index) thread_id: u32,
+) {
+    let base = folded_workgroup_index(workgroup_id) * ${SORT_TILE_SIZE}u;
+
+    if (thread_id == 0u) {
+        let listed = min(probe_trace.count, settings.trace_capacity);
+
+        wg_count = min(select(0u, listed - base, listed > base), ${SORT_TILE_SIZE}u);
+    }
+
+    // the sort's barriers sit under loops its count bounds, so the count is the workgroup's
+    let count = workgroupUniformLoad(&wg_count);
+    let last = i32(settings.resolution) - 1;
+
+    for (var k = 0u; k < ${SORT_KEYS_PER_THREAD}u; k++) {
+        let at = k * ${SORT_WORKGROUP_SIZE}u + thread_id;
+
+        if (at < count) {
+            let slot = probe_trace.elements[base + at];
+            let cascade = probe_record_cascade(probe_records[slot].flags);
+            let local = vec3<u32>(clamp(probe_records[slot].cell - probe_levels[cascade].corner, vec3<i32>(0), vec3<i32>(last)));
+
+            keys[at] = (cascade << (3u * settings.axis_bits)) | v3_morton_encode(local.x, local.y, local.z);
+            values[at] = slot;
+        }
+    }
+
+    sort_bitonic_workgroup(thread_id, count);
+
+    for (var k = 0u; k < ${SORT_KEYS_PER_THREAD}u; k++) {
+        let at = k * ${SORT_WORKGROUP_SIZE}u + thread_id;
+
+        if (at < count) {
+            probe_trace.elements[base + at] = values[at];
+        }
+    }
+}
+`, [
+            chunk_probe_clipmap,
+            chunk_folded_workgroup_index,
+            chunk_v3_morton_encode,
+            chunk_sort_bitonic_workgroup,
+        ]),
+});

package/src/shade/renderer/global_illumination/probe_clipmap/probe/shader_probe_clipmap_project.d.ts

@@ -6,8 +6,9 @@ export const PROBE_PROJECT_SETTINGS: WebGPUStruct;
  * Projection, rows (PROBE_CLIPMAP_PROPOSAL_2026_10_27.md §8.8): a thread a row of each listed probe's
  * radiance map, dispatched through `PROBE_DISPATCH_ROWS`. It sums the row's texels, new and kept, as
  * irradiance SH, `convert_to_irradiance_sh(L, ω)·dω/(2π)` over each texel's solid angle, so that the
- * resolve gives E/π; and the row's rays this update, resolved and back faces. One row of the probe's
- * projection each, which the probe's pass sums.
+ * resolve gives E/π; the row's rays this update, resolved and back faces; and how far the update's
+ * samples fell from their means, with that sum's variance. One row of the probe's projection each, which
+ * the probe's pass sums.
  *
  * @type {ComputeShader}
  */
@@ -18,6 +19,27 @@ export const shader_probe_clipmap_project_rows: ComputeShader;
  * by the update's resolved rays: dead past a quarter of back faces, else alive, its update counted and
  * its time the frame's. An update with no resolved ray keeps the probe's kind.
  *
+ * Only uniform rays classify. Those are a uniform update's, and a guided update's whose guide weighs
+ * nothing (a dark history), which spreads its rays evenly. A guided update's rays shun back faces, whose
+ * radiance is nothing, so their share would read a probe inside a solid as alive. Such an update keeps
+ * the probe's kind and its back-face share.
+ *
+ * **The probe's ray request** (§8.6): the rays its next guided update asks for, in `probe_requests`.
+ * - After a probe's first update, the most: Enshrouded's new probe that asks for all it can have.
+ * - When its light has moved, the most again. The update's change score sums its blended texels' falls
+ *   from their means, luminance, each in its own deviation and held to `PROBE_CHANGE_SIGNIFICANCE` by
+ *   accumulation, over the root of their count. A test texel by texel sees no change under a one-ray
+ *   texel's own noise (measured: the estimator's inconsistency barely moves on a doubling); a sum of the
+ *   falls unbounded is one bright sample's.
+ * - The score is held to the probe's own: moved means past `PROBE_CHANGE_SIGNIFICANCE` of its running
+ *   deviation, at least one, from its running mean, both over the estimator's variance window (2τ). A
+ *   probe whose means sit off its samples (the estimator's own bias on texels lit by few rays) scores
+ *   steadily off zero.
+ * - Measured, Sponza on the 4090, the camera and the light still, 3,600 frames settled: of the probes
+ *   past 64 updates, 4.2% asked for the most against a sum of the raw falls held to zero, 3.2% with each
+ *   texel's fall bounded, and 1.1–1.3% held to their own statistics.
+ * - Otherwise half what it asked for last, down to the fewest: a converged probe sips.
+ *
  * @type {ComputeShader}
  */
 export const shader_probe_clipmap_project: ComputeShader;

package/src/shade/renderer/global_illumination/probe_clipmap/probe/shader_probe_clipmap_project.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"shader_probe_clipmap_project.d.ts","sourceRoot":"","sources":["../../../../../../../src/shade/renderer/global_illumination/probe_clipmap/probe/shader_probe_clipmap_project.js"],"names":[],"mappings":"AAqBA;;GAEG;AACH,qCAFU,YAAY,CAWK;AAW3B;;;;;;;;GAQG;AACH,gDAFU,aAAa,CAoDpB;AAYH;;;;;;;GAOG;AACH,2CAFU,aAAa,CAuDpB;6BAxK0B,sCAAsC;8BAHrC,kCAAkC"}
+{"version":3,"file":"shader_probe_clipmap_project.d.ts","sourceRoot":"","sources":["../../../../../../../src/shade/renderer/global_illumination/probe_clipmap/probe/shader_probe_clipmap_project.js"],"names":[],"mappings":"AAyBA;;GAEG;AACH,qCAFU,YAAY,CAqBK;AAW3B;;;;;;;;;GASG;AACH,gDAFU,aAAa,CAuDpB;AAcH;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA4BG;AACH,2CAFU,aAAa,CA8FpB;6BAxP0B,sCAAsC;8BAHrC,kCAAkC"}