@photivra/engine 1.5.0 → 2.0.0

package/dist/gpu/chart-radiance-gpu.js

@@ -0,0 +1,176 @@
+// SPDX-License-Identifier: Apache-2.0
+/**
+ * Host side of the `webgpu-f32` chart-radiance stage (#277): per-material value tables in float64 (the
+ * reference's own operations), preset constants, banded dispatch, readback, and host recomputation of
+ * flagged points with the reference. The float64 reference is `calculateChartRadianceBatch`
+ * (facade/chart-radiance-batch.ts). Internal; the public entry point is the image-kernel accelerator.
+ */
+import { createHostYieldGate } from "../capture/host-yield-cadence.js";
+import { calculateChartRadianceBatchPlan } from "../facade/chart-radiance-batch.js";
+import { calculatePointGridCenters } from "../facade/chart-scenes.js";
+import { ImageKernelExecutionError, webGpuExecution } from "./image-kernel-execution.js";
+import { BUFFER_USAGE, GpuResourceScope, abortedError, linearDispatch, runComputePass } from "./webgpu-runtime.js";
+import { CHART_PRESET_CODE, CHART_RADIANCE_WGSL, CHART_WORKGROUP_SIZE } from "./wgsl-chart-kernel.js";
+/** Frame-unit magnitude limit for GPU points. */
+export const CHART_GPU_MAXIMUM_COORDINATE = 2 ** 20;
+const FLOAT32_SAFE_MAGNITUDE = 2 ** 120;
+function checkAbort(signal) {
+    if (signal.aborted)
+        throw abortedError();
+}
+export async function evaluateChartRadianceOnCpu(input, abort, yieldControl, execution) {
+    const plan = calculateChartRadianceBatchPlan(input), points = input.points.slice(), indices = input.wavelengthIndices?.slice();
+    const values = new Float64Array(plan.queryCount * plan.columns), gate = createHostYieldGate();
+    for (let q = 0; q < plan.queryCount; q++) {
+        for (let c = 0; c < plan.columns; c++)
+            values[q * plan.columns + c] = plan.queryValue(points, indices, q, c);
+        if ((q & 0xffff) === 0xffff) {
+            checkAbort(abort.signal);
+            if (yieldControl && gate.due()) {
+                await yieldControl(abort.signal);
+                gate.restart();
+            }
+        }
+    }
+    checkAbort(abort.signal);
+    return { plan, values, columns: plan.columns, execution, hostRecomputedQueryIndices: new Uint32Array(0) };
+}
+/** Preset constants (uniform words 3–5 and floats 8–15) and the auxiliary table. @internal */
+export function chartGpuConstants(plan) {
+    const chart = plan.chart, p = chart.parameters, aspect = chart.aspect;
+    switch (chart.preset) {
+        case "uniform": return { ints: [0, 0, 0], floats: [], aux: new Float32Array(1) };
+        case "slanted-edge": {
+            // The same float64 constants as resolveChartScene.
+            const t = p.angleDegrees * Math.PI / 180;
+            return { ints: [0, 0, 0], floats: [Math.cos(t), Math.sin(t)], aux: new Float32Array(1) };
+        }
+        case "siemens-star": {
+            const cycles = p.cycles, r = p.radiusFraction * aspect, aux = new Float32Array(2 * (2 * cycles + 1));
+            // Sector boundaries: phase·2·cycles = k  ⇔  θ = π(k/cycles − 1).
+            for (let k = 0; k <= 2 * cycles; k++) {
+                const b = Math.PI * (k / cycles - 1);
+                aux[2 * k] = Math.cos(b);
+                aux[2 * k + 1] = Math.sin(b);
+            }
+            return { ints: [cycles, 0, 0], floats: [r * r], aux };
+        }
+        case "step-wedge": return { ints: [p.steps, 0, 0], floats: [p.steps / 2], aux: new Float32Array(1) };
+        case "grey-ramp": return { ints: [1024, 0, 0], floats: [512], aux: new Float32Array(1) };
+        case "checkerboard": return { ints: [0, 0, 0], floats: [p.squaresAcross / 2, aspect], aux: new Float32Array(1) };
+        case "point-grid": {
+            const columns = p.columns, rows = p.rows, centers = calculatePointGridCenters({ columns, rows, aspect });
+            return { ints: [columns, rows, 0], floats: [columns / 2, aspect, rows / (2 * aspect), p.radius * p.radius],
+                aux: Float32Array.from(centers.flatMap(c => [c.u, c.v])) };
+        }
+        case "patches": {
+            const columns = p.columns, rows = p.rows, count = chart.parameters.patchIds.length;
+            return { ints: [columns, rows, count], floats: [columns / 2, aspect, rows / (2 * aspect), p.gapFraction / 2, 1 - p.gapFraction / 2],
+                aux: new Float32Array(1) };
+        }
+    }
+}
+/** Per-material values in float64, with the reference's operations; a material whose spectrum fails is NaN (recomputed on the host). */
+function materialTable(plan) {
+    const materials = plan.chart.materials.length, nm = plan.wavelengthsNanometers, stride = plan.mode === "radiance" ? nm.length : plan.columns;
+    const table = new Float32Array(materials * stride);
+    for (let m = 0; m < materials; m++) {
+        for (let c = 0; c < stride; c++) {
+            let value;
+            try {
+                if (plan.mode === "radiance")
+                    value = plan.materialRadiance(m, nm[c]);
+                else {
+                    value = 0;
+                    const row = plan.channelWeights[c];
+                    for (let j = 0; j < nm.length; j++)
+                        value += row[j] * plan.materialRadiance(m, nm[j]);
+                }
+            }
+            catch {
+                value = Number.NaN;
+            }
+            if (Math.abs(value) > FLOAT32_SAFE_MAGNITUDE)
+                throw new ImageKernelExecutionError("unsupported-float32-domain", "Chart values exceed the float32 domain.");
+            table[m * stride + c] = value;
+        }
+    }
+    return table;
+}
+export async function evaluateChartRadianceOnGpu(runtime, bandBytes, input, abort, yieldControl, pipelines) {
+    const plan = calculateChartRadianceBatchPlan(input), points = input.points.slice(), indices = input.wavelengthIndices?.slice();
+    for (let i = 0; i < points.length; i++) {
+        if (Math.abs(points[i]) > CHART_GPU_MAXIMUM_COORDINATE)
+            throw new ImageKernelExecutionError("unsupported-float32-domain", "Chart points beyond 2^20 frame units are outside the float32 domain.");
+    }
+    const { ints, floats, aux } = chartGpuConstants(plan), table = materialTable(plan), radianceMode = plan.mode === "radiance";
+    const columns = plan.columns, stride = radianceMode ? plan.wavelengthsNanometers.length : columns, count = plan.queryCount;
+    const bytesPerQuery = 8 + 4 * columns + 4 + 4;
+    const rows = Math.max(1, Math.floor(bandBytes / bytesPerQuery));
+    if (Math.max(table.byteLength, aux.byteLength, Math.min(rows, count) * 4 * columns) > runtime.maximumStorageBindingBytes) {
+        throw new ImageKernelExecutionError("insufficient-limits", "Chart tables or bands exceed the device's storage binding limit.");
+    }
+    const params = new ArrayBuffer(64), words = new Uint32Array(params), f = new Float32Array(params);
+    words[1] = stride;
+    words[2] = CHART_PRESET_CODE[plan.chart.preset];
+    words[3] = ints[0];
+    words[4] = ints[1];
+    words[5] = ints[2];
+    words[6] = radianceMode ? 1 : 0;
+    f.set(floats, 8);
+    const pipeline = await pipelines("photivra-chart-radiance", CHART_RADIANCE_WGSL);
+    const values = new Float32Array(count * columns), recomputed = [];
+    const scope = new GpuResourceScope(runtime);
+    try {
+        const tableBuffer = scope.upload("photivra-chart-table", table, BUFFER_USAGE.STORAGE);
+        const auxBuffer = scope.upload("photivra-chart-aux", aux, BUFFER_USAGE.STORAGE);
+        const gate = createHostYieldGate();
+        for (let first = 0; first < count; first += rows) {
+            checkAbort(abort.signal);
+            const n = Math.min(rows, count - first);
+            words[0] = n;
+            const band = new GpuResourceScope(runtime);
+            try {
+                const pointBuffer = band.upload("photivra-chart-points", new Float32Array(points.subarray(2 * first, 2 * (first + n))), BUFFER_USAGE.STORAGE);
+                const selectors = band.upload("photivra-chart-selectors", indices ? Uint32Array.from(indices.subarray(first, first + n)) : new Uint32Array(1), BUFFER_USAGE.STORAGE);
+                const out = band.buffer("photivra-chart-values", n * 4 * columns, BUFFER_USAGE.STORAGE | BUFFER_USAGE.COPY_SRC);
+                const flagBuffer = band.buffer("photivra-chart-flags", n * 4, BUFFER_USAGE.STORAGE | BUFFER_USAGE.COPY_SRC);
+                const uniform = band.upload("photivra-chart-params", new Uint8Array(params), BUFFER_USAGE.UNIFORM);
+                const bindGroup = runtime.device.createBindGroup({ layout: pipeline.getBindGroupLayout(0), entries: [
+                        { binding: 0, resource: { buffer: uniform } }, { binding: 1, resource: { buffer: pointBuffer } }, { binding: 2, resource: { buffer: tableBuffer } },
+                        { binding: 3, resource: { buffer: auxBuffer } }, { binding: 4, resource: { buffer: selectors } }, { binding: 5, resource: { buffer: out } },
+                        { binding: 6, resource: { buffer: flagBuffer } }
+                    ] });
+                const [valueData, flagData] = await runComputePass(runtime, band, pipeline, bindGroup, linearDispatch(n, CHART_WORKGROUP_SIZE, runtime.maximumWorkgroupsPerDimension), [{ buffer: out, bytes: n * 4 * columns }, { buffer: flagBuffer, bytes: n * 4 }], abort.signal);
+                values.set(new Float32Array(valueData, 0, n * columns), first * columns);
+                const bandFlags = new Uint32Array(flagData, 0, n);
+                for (let i = 0; i < n; i++) {
+                    let recompute = bandFlags[i] !== 0;
+                    for (let c = 0; c < columns && !recompute; c++)
+                        if (!Number.isFinite(values[(first + i) * columns + c]))
+                            recompute = true;
+                    if (recompute)
+                        recomputed.push(first + i);
+                }
+            }
+            finally {
+                band.release();
+            }
+            if (yieldControl && gate.due()) {
+                await yieldControl(abort.signal);
+                gate.restart();
+            }
+        }
+    }
+    finally {
+        scope.release();
+    }
+    checkAbort(abort.signal);
+    // Flagged points (and any whose material spectrum failed): the reference's own value, rounded once to float32.
+    for (const q of recomputed)
+        for (let c = 0; c < columns; c++)
+            values[q * columns + c] = plan.queryValue(points, indices, q, c);
+    const execution = webGpuExecution({ kind: "relative-to-reference", relative: 2 ** -24, absolute: 2 ** -126 }, 0);
+    return { plan, values, columns, execution, hostRecomputedQueryIndices: Uint32Array.from(recomputed) };
+}
+//# sourceMappingURL=chart-radiance-gpu.js.map

package/dist/gpu/chart-radiance-gpu.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"chart-radiance-gpu.js","sourceRoot":"","sources":["../../src/gpu/chart-radiance-gpu.ts"],"names":[],"mappings":"AAAA,sCAAsC;AAEtC;;;;;GAKG;AACH,OAAO,EAAE,mBAAmB,EAAE,MAAM,kCAAkC,CAAC;AACvE,OAAO,EAAE,+BAA+B,EAA6D,MAAM,mCAAmC,CAAC;AAC/I,OAAO,EAAE,yBAAyB,EAAE,MAAM,2BAA2B,CAAC;AACtE,OAAO,EAAE,yBAAyB,EAAE,eAAe,EAA6B,MAAM,6BAA6B,CAAC;AAEpH,OAAO,EAAE,YAAY,EAAE,gBAAgB,EAAE,YAAY,EAAE,cAAc,EAAE,cAAc,EAAsB,MAAM,qBAAqB,CAAC;AACvI,OAAO,EAAE,iBAAiB,EAAE,mBAAmB,EAAE,oBAAoB,EAAE,MAAM,wBAAwB,CAAC;AAKtG,iDAAiD;AACjD,MAAM,CAAC,MAAM,4BAA4B,GAAG,CAAC,IAAI,EAAE,CAAC;AACpD,MAAM,sBAAsB,GAAG,CAAC,IAAI,GAAG,CAAC;AAGxC,SAAS,UAAU,CAAC,MAAmB;IACrC,IAAI,MAAM,CAAC,OAAO;QAAE,MAAM,YAAY,EAAE,CAAC;AAC3C,CAAC;AAED,MAAM,CAAC,KAAK,UAAU,0BAA0B,CAAC,KAA8B,EAAE,KAAsB,EAAE,YAA0B,EACjI,SAA+B;IAC/B,MAAM,IAAI,GAAG,+BAA+B,CAAC,KAAK,CAAC,EAAE,MAAM,GAAG,KAAK,CAAC,MAAM,CAAC,KAAK,EAAE,EAAE,OAAO,GAAG,KAAK,CAAC,iBAAiB,EAAE,KAAK,EAAE,CAAC;IAC/H,MAAM,MAAM,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,UAAU,GAAG,IAAI,CAAC,OAAO,CAAC,EAAE,IAAI,GAAG,mBAAmB,EAAE,CAAC;IAC9F,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,UAAU,EAAE,CAAC,EAAE,EAAE,CAAC;QACzC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,OAAO,EAAE,CAAC,EAAE;YAAE,MAAM,CAAC,CAAC,GAAG,IAAI,CAAC,OAAO,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,UAAU,CAAC,MAAM,EAAE,OAAO,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAC7G,IAAI,CAAC,CAAC,GAAG,MAAM,CAAC,KAAK,MAAM,EAAE,CAAC;YAC5B,UAAU,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;YACzB,IAAI,YAAY,IAAI,IAAI,CAAC,GAAG,EAAE,EAAE,CAAC;gBAAC,MAAM,YAAY,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;gBAAC,IAAI,CAAC,OAAO,EAAE,CAAC;YAAC,CAAC;QACvF,CAAC;IACH,CAAC;IACD,UAAU,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;IACzB,OAAO,EAAE,IAAI,EAAE,MAAM,EAAE,OAAO,EAAE,IAAI,CAAC,OAAO,EAAE,SAAS,EAAE,0BAA0B,EAAE,IAAI,WAAW,CAAC,CAAC,CAAC,EAAE,CAAC;AAC5G,CAAC;AAED,8FAA8F;AAC9F,MAAM,UAAU,iBAAiB,CAAC,IAA4B;IAC5D,MAAM,KAAK,GAAG,IAAI,CAAC,KAAK,EAAE,CAAC,GAAG,KAAK,CAAC,UAAoC,EAAE,MAAM,GAAG,KAAK,CAAC,MAAM,CAAC;IAChG,QAAQ,KAAK,CAAC,MAAM,EAAE,CAAC;QACrB,KAAK,SAAS,CAAC,CAAC,OAAO,EAAE,IAAI,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,GAAG,EAAE,IAAI,YAAY,CAAC,CAAC,CAAC,EAAE,CAAC;QACjF,KAAK,cAAc,CAAC,CAAC,CAAC;YACpB,mDAAmD;YACnD,MAAM,CAAC,GAAG,CAAC,CAAC,YAAa,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC;YAC1C,OAAO,EAAE,IAAI,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,IAAI,YAAY,CAAC,CAAC,CAAC,EAAE,CAAC;QAC3F,CAAC;QACD,KAAK,cAAc,CAAC,CAAC,CAAC;YACpB,MAAM,MAAM,GAAG,CAAC,CAAC,MAAO,EAAE,CAAC,GAAG,CAAC,CAAC,cAAe,GAAG,MAAM,EAAE,GAAG,GAAG,IAAI,YAAY,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC;YACvG,iEAAiE;YACjE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,GAAG,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;gBAAC,MAAM,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,MAAM,GAAG,CAAC,CAAC,CAAC;gBAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;gBAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YAAC,CAAC;YACvI,OAAO,EAAE,IAAI,EAAE,CAAC,MAAM,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC;QACxD,CAAC;QACD,KAAK,YAAY,CAAC,CAAC,OAAO,EAAE,IAAI,EAAE,CAAC,CAAC,CAAC,KAAM,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,CAAC,CAAC,KAAM,GAAG,CAAC,CAAC,EAAE,GAAG,EAAE,IAAI,YAAY,CAAC,CAAC,CAAC,EAAE,CAAC;QACvG,KAAK,WAAW,CAAC,CAAC,OAAO,EAAE,IAAI,EAAE,CAAC,IAAI,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,IAAI,YAAY,CAAC,CAAC,CAAC,EAAE,CAAC;QACzF,KAAK,cAAc,CAAC,CAAC,OAAO,EAAE,IAAI,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,CAAC,CAAC,aAAc,GAAG,CAAC,EAAE,MAAM,CAAC,EAAE,GAAG,EAAE,IAAI,YAAY,CAAC,CAAC,CAAC,EAAE,CAAC;QAClH,KAAK,YAAY,CAAC,CAAC,CAAC;YAClB,MAAM,OAAO,GAAG,CAAC,CAAC,OAAQ,EAAE,IAAI,GAAG,CAAC,CAAC,IAAK,EAAE,OAAO,GAAG,yBAAyB,CAAC,EAAE,OAAO,EAAE,IAAI,EAAE,MAAM,EAAE,CAAC,CAAC;YAC3G,OAAO,EAAE,IAAI,EAAE,CAAC,OAAO,EAAE,IAAI,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,OAAO,GAAG,CAAC,EAAE,MAAM,EAAE,IAAI,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,EAAE,CAAC,CAAC,MAAO,GAAG,CAAC,CAAC,MAAO,CAAC;gBAC1G,GAAG,EAAE,YAAY,CAAC,IAAI,CAAC,OAAO,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;QAC/D,CAAC;QACD,KAAK,SAAS,CAAC,CAAC,CAAC;YACf,MAAM,OAAO,GAAG,CAAC,CAAC,OAAQ,EAAE,IAAI,GAAG,CAAC,CAAC,IAAK,EAAE,KAAK,GAAI,KAAK,CAAC,UAAU,CAAC,QAA8B,CAAC,MAAM,CAAC;YAC5G,OAAO,EAAE,IAAI,EAAE,CAAC,OAAO,EAAE,IAAI,EAAE,KAAK,CAAC,EAAE,MAAM,EAAE,CAAC,OAAO,GAAG,CAAC,EAAE,MAAM,EAAE,IAAI,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,EAAE,CAAC,CAAC,WAAY,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,WAAY,GAAG,CAAC,CAAC;gBACnI,GAAG,EAAE,IAAI,YAAY,CAAC,CAAC,CAAC,EAAE,CAAC;QAC/B,CAAC;IACH,CAAC;AACH,CAAC;AAED,wIAAwI;AACxI,SAAS,aAAa,CAAC,IAA4B;IACjD,MAAM,SAAS,GAAG,IAAI,CAAC,KAAK,CAAC,SAAS,CAAC,MAAM,EAAE,EAAE,GAAG,IAAI,CAAC,qBAAqB,EAAE,MAAM,GAAG,IAAI,CAAC,IAAI,KAAK,UAAU,CAAC,CAAC,CAAC,EAAE,CAAC,MAAM,CAAC,CAAC,CAAC,IAAI,CAAC,OAAO,CAAC;IAC7I,MAAM,KAAK,GAAG,IAAI,YAAY,CAAC,SAAS,GAAG,MAAM,CAAC,CAAC;IACnD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,SAAS,EAAE,CAAC,EAAE,EAAE,CAAC;QACnC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;YAChC,IAAI,KAAa,CAAC;YAClB,IAAI,CAAC;gBACH,IAAI,IAAI,CAAC,IAAI,KAAK,UAAU;oBAAE,KAAK,GAAG,IAAI,CAAC,gBAAgB,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAE,CAAC,CAAC;qBAClE,CAAC;oBACJ,KAAK,GAAG,CAAC,CAAC;oBACV,MAAM,GAAG,GAAG,IAAI,CAAC,cAAe,CAAC,CAAC,CAAE,CAAC;oBACrC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,MAAM,EAAE,CAAC,EAAE;wBAAE,KAAK,IAAI,GAAG,CAAC,CAAC,CAAE,GAAG,IAAI,CAAC,gBAAgB,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAE,CAAC,CAAC;gBAC1F,CAAC;YACH,CAAC;YAAC,MAAM,CAAC;gBACP,KAAK,GAAG,MAAM,CAAC,GAAG,CAAC;YACrB,CAAC;YACD,IAAI,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,GAAG,sBAAsB;gBAAE,MAAM,IAAI,yBAAyB,CAAC,4BAA4B,EAAE,yCAAyC,CAAC,CAAC;YAC3J,KAAK,CAAC,CAAC,GAAG,MAAM,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC;QAChC,CAAC;IACH,CAAC;IACD,OAAO,KAAK,CAAC;AACf,CAAC;AAED,MAAM,CAAC,KAAK,UAAU,0BAA0B,CAAC,OAAsB,EAAE,SAAiB,EAAE,KAA8B,EAAE,KAAsB,EAChJ,YAA0B,EAAE,SAAyB;IACrD,MAAM,IAAI,GAAG,+BAA+B,CAAC,KAAK,CAAC,EAAE,MAAM,GAAG,KAAK,CAAC,MAAM,CAAC,KAAK,EAAE,EAAE,OAAO,GAAG,KAAK,CAAC,iBAAiB,EAAE,KAAK,EAAE,CAAC;IAC/H,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACvC,IAAI,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAE,CAAC,GAAG,4BAA4B;YAAE,MAAM,IAAI,yBAAyB,CAAC,4BAA4B,EAAE,sEAAsE,CAAC,CAAC;IACrM,CAAC;IACD,MAAM,EAAE,IAAI,EAAE,MAAM,EAAE,GAAG,EAAE,GAAG,iBAAiB,CAAC,IAAI,CAAC,EAAE,KAAK,GAAG,aAAa,CAAC,IAAI,CAAC,EAAE,YAAY,GAAG,IAAI,CAAC,IAAI,KAAK,UAAU,CAAC;IAC5H,MAAM,OAAO,GAAG,IAAI,CAAC,OAAO,EAAE,MAAM,GAAG,YAAY,CAAC,CAAC,CAAC,IAAI,CAAC,qBAAqB,CAAC,MAAM,CAAC,CAAC,CAAC,OAAO,EAAE,KAAK,GAAG,IAAI,CAAC,UAAU,CAAC;IAC3H,MAAM,aAAa,GAAG,CAAC,GAAG,CAAC,GAAG,OAAO,GAAG,CAAC,GAAG,CAAC,CAAC;IAC9C,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,SAAS,GAAG,aAAa,CAAC,CAAC,CAAC;IAChE,IAAI,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,UAAU,EAAE,GAAG,CAAC,UAAU,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,KAAK,CAAC,GAAG,CAAC,GAAG,OAAO,CAAC,GAAG,OAAO,CAAC,0BAA0B,EAAE,CAAC;QACzH,MAAM,IAAI,yBAAyB,CAAC,qBAAqB,EAAE,kEAAkE,CAAC,CAAC;IACjI,CAAC;IACD,MAAM,MAAM,GAAG,IAAI,WAAW,CAAC,EAAE,CAAC,EAAE,KAAK,GAAG,IAAI,WAAW,CAAC,MAAM,CAAC,EAAE,CAAC,GAAG,IAAI,YAAY,CAAC,MAAM,CAAC,CAAC;IAClG,KAAK,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC;IAAC,KAAK,CAAC,CAAC,CAAC,GAAG,iBAAiB,CAAC,IAAI,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;IAAC,KAAK,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;IAAC,KAAK,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;IAAC,KAAK,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;IAAC,KAAK,CAAC,CAAC,CAAC,GAAG,YAAY,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAChK,CAAC,CAAC,GAAG,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC;IACjB,MAAM,QAAQ,GAAG,MAAM,SAAS,CAAC,yBAAyB,EAAE,mBAAmB,CAAC,CAAC;IACjF,MAAM,MAAM,GAAG,IAAI,YAAY,CAAC,KAAK,GAAG,OAAO,CAAC,EAAE,UAAU,GAAa,EAAE,CAAC;IAC5E,MAAM,KAAK,GAAG,IAAI,gBAAgB,CAAC,OAAO,CAAC,CAAC;IAC5C,IAAI,CAAC;QACH,MAAM,WAAW,GAAG,KAAK,CAAC,MAAM,CAAC,sBAAsB,EAAE,KAAK,EAAE,YAAY,CAAC,OAAO,CAAC,CAAC;QACtF,MAAM,SAAS,GAAG,KAAK,CAAC,MAAM,CAAC,oBAAoB,EAAE,GAAG,EAAE,YAAY,CAAC,OAAO,CAAC,CAAC;QAChF,MAAM,IAAI,GAAG,mBAAmB,EAAE,CAAC;QACnC,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,KAAK,EAAE,KAAK,IAAI,IAAI,EAAE,CAAC;YACjD,UAAU,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;YACzB,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,KAAK,GAAG,KAAK,CAAC,CAAC;YACxC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YACb,MAAM,IAAI,GAAG,IAAI,gBAAgB,CAAC,OAAO,CAAC,CAAC;YAC3C,IAAI,CAAC;gBACH,MAAM,WAAW,GAAG,IAAI,CAAC,MAAM,CAAC,uBAAuB,EAAE,IAAI,YAAY,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,GAAG,KAAK,EAAE,CAAC,GAAG,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,YAAY,CAAC,OAAO,CAAC,CAAC;gBAC9I,MAAM,SAAS,GAAG,IAAI,CAAC,MAAM,CAAC,0BAA0B,EAAE,OAAO,CAAC,CAAC,CAAC,WAAW,CAAC,IAAI,CAAC,OAAO,CAAC,QAAQ,CAAC,KAAK,EAAE,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,WAAW,CAAC,CAAC,CAAC,EAAE,YAAY,CAAC,OAAO,CAAC,CAAC;gBACrK,MAAM,GAAG,GAAG,IAAI,CAAC,MAAM,CAAC,uBAAuB,EAAE,CAAC,GAAG,CAAC,GAAG,OAAO,EAAE,YAAY,CAAC,OAAO,GAAG,YAAY,CAAC,QAAQ,CAAC,CAAC;gBAChH,MAAM,UAAU,GAAG,IAAI,CAAC,MAAM,CAAC,sBAAsB,EAAE,CAAC,GAAG,CAAC,EAAE,YAAY,CAAC,OAAO,GAAG,YAAY,CAAC,QAAQ,CAAC,CAAC;gBAC5G,MAAM,OAAO,GAAG,IAAI,CAAC,MAAM,CAAC,uBAAuB,EAAE,IAAI,UAAU,CAAC,MAAM,CAAC,EAAE,YAAY,CAAC,OAAO,CAAC,CAAC;gBACnG,MAAM,SAAS,GAAG,OAAO,CAAC,MAAM,CAAC,eAAe,CAAC,EAAE,MAAM,EAAE,QAAQ,CAAC,kBAAkB,CAAC,CAAC,CAAC,EAAE,OAAO,EAAE;wBAClG,EAAE,OAAO,EAAE,CAAC,EAAE,QAAQ,EAAE,EAAE,MAAM,EAAE,OAAO,EAAE,EAAE,EAAE,EAAE,OAAO,EAAE,CAAC,EAAE,QAAQ,EAAE,EAAE,MAAM,EAAE,WAAW,EAAE,EAAE,EAAE,EAAE,OAAO,EAAE,CAAC,EAAE,QAAQ,EAAE,EAAE,MAAM,EAAE,WAAW,EAAE,EAAE;wBACnJ,EAAE,OAAO,EAAE,CAAC,EAAE,QAAQ,EAAE,EAAE,MAAM,EAAE,SAAS,EAAE,EAAE,EAAE,EAAE,OAAO,EAAE,CAAC,EAAE,QAAQ,EAAE,EAAE,MAAM,EAAE,SAAS,EAAE,EAAE,EAAE,EAAE,OAAO,EAAE,CAAC,EAAE,QAAQ,EAAE,EAAE,MAAM,EAAE,GAAG,EAAE,EAAE;wBAC3I,EAAE,OAAO,EAAE,CAAC,EAAE,QAAQ,EAAE,EAAE,MAAM,EAAE,UAAU,EAAE,EAAE;qBAAC,EAAE,CAAC,CAAC;gBACvD,MAAM,CAAC,SAAS,EAAE,QAAQ,CAAC,GAAG,MAAM,cAAc,CAAC,OAAO,EAAE,IAAI,EAAE,QAAQ,EAAE,SAAS,EAAE,cAAc,CAAC,CAAC,EAAE,oBAAoB,EAAE,OAAO,CAAC,6BAA6B,CAAC,EACnK,CAAC,EAAE,MAAM,EAAE,GAAG,EAAE,KAAK,EAAE,CAAC,GAAG,CAAC,GAAG,OAAO,EAAE,EAAE,EAAE,MAAM,EAAE,UAAU,EAAE,KAAK,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,KAAK,CAAC,MAAM,CAAC,CAAC;gBACjG,MAAM,CAAC,GAAG,CAAC,IAAI,YAAY,CAAC,SAAU,EAAE,CAAC,EAAE,CAAC,GAAG,OAAO,CAAC,EAAE,KAAK,GAAG,OAAO,CAAC,CAAC;gBAC1E,MAAM,SAAS,GAAG,IAAI,WAAW,CAAC,QAAS,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;gBACnD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC3B,IAAI,SAAS,GAAG,SAAS,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC;oBACnC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,IAAI,CAAC,SAAS,EAAE,CAAC,EAAE;wBAAE,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC,GAAG,OAAO,GAAG,CAAC,CAAE,CAAC;4BAAE,SAAS,GAAG,IAAI,CAAC;oBAC3H,IAAI,SAAS;wBAAE,UAAU,CAAC,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC;gBAC5C,CAAC;YACH,CAAC;oBAAS,CAAC;gBACT,IAAI,CAAC,OAAO,EAAE,CAAC;YACjB,CAAC;YACD,IAAI,YAAY,IAAI,IAAI,CAAC,GAAG,EAAE,EAAE,CAAC;gBAAC,MAAM,YAAY,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;gBAAC,IAAI,CAAC,OAAO,EAAE,CAAC;YAAC,CAAC;QACvF,CAAC;IACH,CAAC;YAAS,CAAC;QACT,KAAK,CAAC,OAAO,EAAE,CAAC;IAClB,CAAC;IACD,UAAU,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;IACzB,+GAA+G;IAC/G,KAAK,MAAM,CAAC,IAAI,UAAU;QAAE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,EAAE,CAAC,EAAE;YAAE,MAAM,CAAC,CAAC,GAAG,OAAO,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,UAAU,CAAC,MAAM,EAAE,OAAO,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAC/H,MAAM,SAAS,GAAG,eAAe,CAAC,EAAE,IAAI,EAAE,uBAAuB,EAAE,QAAQ,EAAE,CAAC,IAAI,CAAC,EAAE,EAAE,QAAQ,EAAE,CAAC,IAAI,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC,CAAC;IACjH,OAAO,EAAE,IAAI,EAAE,MAAM,EAAE,OAAO,EAAE,SAAS,EAAE,0BAA0B,EAAE,WAAW,CAAC,IAAI,CAAC,UAAU,CAAC,EAAE,CAAC;AACxG,CAAC","sourcesContent":["// SPDX-License-Identifier: Apache-2.0\n\n/**\n * Host side of the `webgpu-f32` chart-radiance stage (#277): per-material value tables in float64 (the\n * reference's own operations), preset constants, banded dispatch, readback, and host recomputation of\n * flagged points with the reference. The float64 reference is `calculateChartRadianceBatch`\n * (facade/chart-radiance-batch.ts). Internal; the public entry point is the image-kernel accelerator.\n */\nimport { createHostYieldGate } from \"../capture/host-yield-cadence.js\";\nimport { calculateChartRadianceBatchPlan, type ChartRadianceBatchInput, type ChartRadianceBatchPlan } from \"../facade/chart-radiance-batch.js\";\nimport { calculatePointGridCenters } from \"../facade/chart-scenes.js\";\nimport { ImageKernelExecutionError, webGpuExecution, type ImageKernelExecution } from \"./image-kernel-execution.js\";\nimport type { ChartRadianceResult } from \"./image-kernel-accelerator.js\";\nimport { BUFFER_USAGE, GpuResourceScope, abortedError, linearDispatch, runComputePass, type WebGpuRuntime } from \"./webgpu-runtime.js\";\nimport { CHART_PRESET_CODE, CHART_RADIANCE_WGSL, CHART_WORKGROUP_SIZE } from \"./wgsl-chart-kernel.js\";\n\ntype PipelineSource = (label: string, code: string) => Promise<GPUComputePipeline>;\ntype YieldControl = ((signal: AbortSignal) => Promise<void>) | undefined;\n\n/** Frame-unit magnitude limit for GPU points. */\nexport const CHART_GPU_MAXIMUM_COORDINATE = 2 ** 20;\nconst FLOAT32_SAFE_MAGNITUDE = 2 ** 120;\n\n\nfunction checkAbort(signal: AbortSignal): void {\n  if (signal.aborted) throw abortedError();\n}\n\nexport async function evaluateChartRadianceOnCpu(input: ChartRadianceBatchInput, abort: AbortController, yieldControl: YieldControl,\n  execution: ImageKernelExecution): Promise<ChartRadianceResult> {\n  const plan = calculateChartRadianceBatchPlan(input), points = input.points.slice(), indices = input.wavelengthIndices?.slice();\n  const values = new Float64Array(plan.queryCount * plan.columns), gate = createHostYieldGate();\n  for (let q = 0; q < plan.queryCount; q++) {\n    for (let c = 0; c < plan.columns; c++) values[q * plan.columns + c] = plan.queryValue(points, indices, q, c);\n    if ((q & 0xffff) === 0xffff) {\n      checkAbort(abort.signal);\n      if (yieldControl && gate.due()) { await yieldControl(abort.signal); gate.restart(); }\n    }\n  }\n  checkAbort(abort.signal);\n  return { plan, values, columns: plan.columns, execution, hostRecomputedQueryIndices: new Uint32Array(0) };\n}\n\n/** Preset constants (uniform words 3–5 and floats 8–15) and the auxiliary table. @internal */\nexport function chartGpuConstants(plan: ChartRadianceBatchPlan): { ints: [number, number, number]; floats: number[]; aux: Float32Array<ArrayBuffer> } {\n  const chart = plan.chart, p = chart.parameters as Record<string, number>, aspect = chart.aspect;\n  switch (chart.preset) {\n    case \"uniform\": return { ints: [0, 0, 0], floats: [], aux: new Float32Array(1) };\n    case \"slanted-edge\": {\n      // The same float64 constants as resolveChartScene.\n      const t = p.angleDegrees! * Math.PI / 180;\n      return { ints: [0, 0, 0], floats: [Math.cos(t), Math.sin(t)], aux: new Float32Array(1) };\n    }\n    case \"siemens-star\": {\n      const cycles = p.cycles!, r = p.radiusFraction! * aspect, aux = new Float32Array(2 * (2 * cycles + 1));\n      // Sector boundaries: phase·2·cycles = k  ⇔  θ = π(k/cycles − 1).\n      for (let k = 0; k <= 2 * cycles; k++) { const b = Math.PI * (k / cycles - 1); aux[2 * k] = Math.cos(b); aux[2 * k + 1] = Math.sin(b); }\n      return { ints: [cycles, 0, 0], floats: [r * r], aux };\n    }\n    case \"step-wedge\": return { ints: [p.steps!, 0, 0], floats: [p.steps! / 2], aux: new Float32Array(1) };\n    case \"grey-ramp\": return { ints: [1024, 0, 0], floats: [512], aux: new Float32Array(1) };\n    case \"checkerboard\": return { ints: [0, 0, 0], floats: [p.squaresAcross! / 2, aspect], aux: new Float32Array(1) };\n    case \"point-grid\": {\n      const columns = p.columns!, rows = p.rows!, centers = calculatePointGridCenters({ columns, rows, aspect });\n      return { ints: [columns, rows, 0], floats: [columns / 2, aspect, rows / (2 * aspect), p.radius! * p.radius!],\n        aux: Float32Array.from(centers.flatMap(c => [c.u, c.v])) };\n    }\n    case \"patches\": {\n      const columns = p.columns!, rows = p.rows!, count = (chart.parameters.patchIds as readonly string[]).length;\n      return { ints: [columns, rows, count], floats: [columns / 2, aspect, rows / (2 * aspect), p.gapFraction! / 2, 1 - p.gapFraction! / 2],\n        aux: new Float32Array(1) };\n    }\n  }\n}\n\n/** Per-material values in float64, with the reference's operations; a material whose spectrum fails is NaN (recomputed on the host). */\nfunction materialTable(plan: ChartRadianceBatchPlan): Float32Array<ArrayBuffer> {\n  const materials = plan.chart.materials.length, nm = plan.wavelengthsNanometers, stride = plan.mode === \"radiance\" ? nm.length : plan.columns;\n  const table = new Float32Array(materials * stride);\n  for (let m = 0; m < materials; m++) {\n    for (let c = 0; c < stride; c++) {\n      let value: number;\n      try {\n        if (plan.mode === \"radiance\") value = plan.materialRadiance(m, nm[c]!);\n        else {\n          value = 0;\n          const row = plan.channelWeights![c]!;\n          for (let j = 0; j < nm.length; j++) value += row[j]! * plan.materialRadiance(m, nm[j]!);\n        }\n      } catch {\n        value = Number.NaN;\n      }\n      if (Math.abs(value) > FLOAT32_SAFE_MAGNITUDE) throw new ImageKernelExecutionError(\"unsupported-float32-domain\", \"Chart values exceed the float32 domain.\");\n      table[m * stride + c] = value;\n    }\n  }\n  return table;\n}\n\nexport async function evaluateChartRadianceOnGpu(runtime: WebGpuRuntime, bandBytes: number, input: ChartRadianceBatchInput, abort: AbortController,\n  yieldControl: YieldControl, pipelines: PipelineSource): Promise<ChartRadianceResult> {\n  const plan = calculateChartRadianceBatchPlan(input), points = input.points.slice(), indices = input.wavelengthIndices?.slice();\n  for (let i = 0; i < points.length; i++) {\n    if (Math.abs(points[i]!) > CHART_GPU_MAXIMUM_COORDINATE) throw new ImageKernelExecutionError(\"unsupported-float32-domain\", \"Chart points beyond 2^20 frame units are outside the float32 domain.\");\n  }\n  const { ints, floats, aux } = chartGpuConstants(plan), table = materialTable(plan), radianceMode = plan.mode === \"radiance\";\n  const columns = plan.columns, stride = radianceMode ? plan.wavelengthsNanometers.length : columns, count = plan.queryCount;\n  const bytesPerQuery = 8 + 4 * columns + 4 + 4;\n  const rows = Math.max(1, Math.floor(bandBytes / bytesPerQuery));\n  if (Math.max(table.byteLength, aux.byteLength, Math.min(rows, count) * 4 * columns) > runtime.maximumStorageBindingBytes) {\n    throw new ImageKernelExecutionError(\"insufficient-limits\", \"Chart tables or bands exceed the device's storage binding limit.\");\n  }\n  const params = new ArrayBuffer(64), words = new Uint32Array(params), f = new Float32Array(params);\n  words[1] = stride; words[2] = CHART_PRESET_CODE[plan.chart.preset]; words[3] = ints[0]; words[4] = ints[1]; words[5] = ints[2]; words[6] = radianceMode ? 1 : 0;\n  f.set(floats, 8);\n  const pipeline = await pipelines(\"photivra-chart-radiance\", CHART_RADIANCE_WGSL);\n  const values = new Float32Array(count * columns), recomputed: number[] = [];\n  const scope = new GpuResourceScope(runtime);\n  try {\n    const tableBuffer = scope.upload(\"photivra-chart-table\", table, BUFFER_USAGE.STORAGE);\n    const auxBuffer = scope.upload(\"photivra-chart-aux\", aux, BUFFER_USAGE.STORAGE);\n    const gate = createHostYieldGate();\n    for (let first = 0; first < count; first += rows) {\n      checkAbort(abort.signal);\n      const n = Math.min(rows, count - first);\n      words[0] = n;\n      const band = new GpuResourceScope(runtime);\n      try {\n        const pointBuffer = band.upload(\"photivra-chart-points\", new Float32Array(points.subarray(2 * first, 2 * (first + n))), BUFFER_USAGE.STORAGE);\n        const selectors = band.upload(\"photivra-chart-selectors\", indices ? Uint32Array.from(indices.subarray(first, first + n)) : new Uint32Array(1), BUFFER_USAGE.STORAGE);\n        const out = band.buffer(\"photivra-chart-values\", n * 4 * columns, BUFFER_USAGE.STORAGE | BUFFER_USAGE.COPY_SRC);\n        const flagBuffer = band.buffer(\"photivra-chart-flags\", n * 4, BUFFER_USAGE.STORAGE | BUFFER_USAGE.COPY_SRC);\n        const uniform = band.upload(\"photivra-chart-params\", new Uint8Array(params), BUFFER_USAGE.UNIFORM);\n        const bindGroup = runtime.device.createBindGroup({ layout: pipeline.getBindGroupLayout(0), entries: [\n          { binding: 0, resource: { buffer: uniform } }, { binding: 1, resource: { buffer: pointBuffer } }, { binding: 2, resource: { buffer: tableBuffer } },\n          { binding: 3, resource: { buffer: auxBuffer } }, { binding: 4, resource: { buffer: selectors } }, { binding: 5, resource: { buffer: out } },\n          { binding: 6, resource: { buffer: flagBuffer } }] });\n        const [valueData, flagData] = await runComputePass(runtime, band, pipeline, bindGroup, linearDispatch(n, CHART_WORKGROUP_SIZE, runtime.maximumWorkgroupsPerDimension),\n          [{ buffer: out, bytes: n * 4 * columns }, { buffer: flagBuffer, bytes: n * 4 }], abort.signal);\n        values.set(new Float32Array(valueData!, 0, n * columns), first * columns);\n        const bandFlags = new Uint32Array(flagData!, 0, n);\n        for (let i = 0; i < n; i++) {\n          let recompute = bandFlags[i] !== 0;\n          for (let c = 0; c < columns && !recompute; c++) if (!Number.isFinite(values[(first + i) * columns + c]!)) recompute = true;\n          if (recompute) recomputed.push(first + i);\n        }\n      } finally {\n        band.release();\n      }\n      if (yieldControl && gate.due()) { await yieldControl(abort.signal); gate.restart(); }\n    }\n  } finally {\n    scope.release();\n  }\n  checkAbort(abort.signal);\n  // Flagged points (and any whose material spectrum failed): the reference's own value, rounded once to float32.\n  for (const q of recomputed) for (let c = 0; c < columns; c++) values[q * columns + c] = plan.queryValue(points, indices, q, c);\n  const execution = webGpuExecution({ kind: \"relative-to-reference\", relative: 2 ** -24, absolute: 2 ** -126 }, 0);\n  return { plan, values, columns, execution, hostRecomputedQueryIndices: Uint32Array.from(recomputed) };\n}\n"]}

package/dist/gpu/image-kernel-accelerator.d.ts

@@ -0,0 +1,132 @@
+import { type NativeRawDevelopmentInput, type NativeRawDevelopmentOutput, type NativeRawDevelopmentProvider } from "../capture/native-raw-development.js";
+import type { SdrRenderingProfile } from "../output/sdr-rendering.js";
+import { type NativeCaptureSdrInput, type NativeCaptureSdrOutput, type NativeCaptureSdrProvider } from "../output/native-capture-sdr.js";
+import { type IrradianceConvolutionInput, type IrradianceConvolutionPlan } from "../optics/irradiance-convolution.js";
+import { type ImageKernelBackend, type ImageKernelExecution, type ImageKernelFallbackReason, type ImageKernelRelativeErrorBound } from "./image-kernel-execution.js";
+import type { PackedSensorRawInput, PackedSensorRawResult } from "../capture/packed-sensor-raw.js";
+import { type TemporalExposureSumInput, type TemporalExposureSumPlan } from "../motion/temporal-exposure-sum.js";
+import type { ChartRadianceBatchInput, ChartRadianceBatchPlan } from "../facade/chart-radiance-batch.js";
+export interface ChartRadianceResult {
+    plan: ChartRadianceBatchPlan;
+    /** Row-major (query, column): float64 for `cpu-reference`, float32 for `webgpu-f32`. */
+    values: Float64Array | Float32Array;
+    columns: number;
+    execution: ImageKernelExecution;
+    /**
+     * Queries whose float32 classification fell inside its bound, recomputed on the host with the
+     * reference (then rounded to float32). Empty for `cpu-reference`.
+     */
+    hostRecomputedQueryIndices: Uint32Array;
+}
+export interface AcceleratedPackedSensorRawResult extends PackedSensorRawResult {
+    execution: ImageKernelExecution;
+    /**
+     * Region-relative indices of the flagged sites, ascending (empty for `cpu-reference`). Every other site
+     * equals the reference; recomputing these sites with the reference makes the whole region bit-identical.
+     */
+    flaggedSiteIndices: Uint32Array;
+}
+export interface TemporalExposureSumResult {
+    plan: TemporalExposureSumPlan;
+    /** Float64 for `cpu-reference`, float32 for `webgpu-f32`; row-major, same grid as the input. */
+    values: Float64Array | Float32Array;
+    execution: ImageKernelExecution;
+}
+export interface ImageKernelAcceleratorOptions {
+    /** Requested backend. `cpu-reference` is the reference and never touches WebGPU. */
+    backend: ImageKernelBackend;
+    /** Required with `webgpu-f32`: use the CPU reference (labelled) or fail when WebGPU cannot be acquired. */
+    fallback?: "cpu-reference" | "fail";
+    /** WebGPU entry point; defaults to `navigator.gpu`. Supply a worker's `navigator.gpu` explicitly if needed. */
+    gpu?: unknown;
+    powerPreference?: "low-power" | "high-performance";
+    /** Largest single GPU band buffer in bytes (4 KiB up to the device limit); default min(64 MiB, device limit). Affects memory only, never results. */
+    maximumBandBytes?: number;
+    /** Aborts adapter/device acquisition. */
+    signal?: AbortSignal;
+}
+export interface ImageKernelRunOptions {
+    signal?: AbortSignal;
+    /** Yields to the host event loop between CPU reference chunks and GPU bands. Convolution, sensor and temporal stages; tasks keep their own yields. */
+    yieldControl?: (signal: AbortSignal) => Promise<void>;
+}
+export type ImageKernelAcceleratorState = "ready" | "lost" | "disposed";
+export interface AcceleratedNativeCaptureSdrOutput extends NativeCaptureSdrOutput {
+    execution: ImageKernelExecution;
+}
+export interface AcceleratedNativeRawDevelopmentOutput extends NativeRawDevelopmentOutput {
+    execution: ImageKernelExecution;
+}
+export interface IrradianceConvolutionResult {
+    plan: IrradianceConvolutionPlan;
+    /** Float64 for `cpu-reference`, float32 for `webgpu-f32`; row-major, same grid as the input. */
+    values: Float64Array | Float32Array;
+    /** Float64 sums of the input and output samples. Their difference is the amount lost across the grid edge. */
+    inputSum: number;
+    outputSum: number;
+    execution: ImageKernelExecution;
+}
+/** Runs one operation at a time. dispose() releases the device and every buffer. */
+export interface ImageKernelAccelerator {
+    readonly backend: ImageKernelBackend;
+    readonly requestedBackend: ImageKernelBackend;
+    readonly fallbackReason: ImageKernelFallbackReason | null;
+    readonly state: ImageKernelAcceleratorState;
+    renderNativeCaptureSdr(input: NativeCaptureSdrInput, provider: NativeCaptureSdrProvider, options?: ImageKernelRunOptions): Promise<AcceleratedNativeCaptureSdrOutput>;
+    /** Takes ownership of `input.raw` exactly like `createNativeRawDevelopmentTask` (the caller's arrays are detached). */
+    developNativeRaw(input: NativeRawDevelopmentInput, provider: NativeRawDevelopmentProvider | ((signal: AbortSignal) => Promise<void>), options?: ImageKernelRunOptions): Promise<AcceleratedNativeRawDevelopmentOutput>;
+    convolveIrradiance(input: IrradianceConvolutionInput, options?: ImageKernelRunOptions): Promise<IrradianceConvolutionResult>;
+    /** Sensor stage on typed arrays: seeded charge, read noise and ADC to RAW codes (see `simulatePackedSensorRaw`). */
+    simulateSensorRaw(input: PackedSensorRawInput, options?: ImageKernelRunOptions): Promise<AcceleratedPackedSensorRawResult>;
+    /** Exposure and motion as a weighted sum of shifted frames (see `calculateTemporalExposureSum`). */
+    sumTemporalExposure(input: TemporalExposureSumInput, options?: ImageKernelRunOptions): Promise<TemporalExposureSumResult>;
+    /** Chart scene radiance or its spectral integral at many frame points (see `calculateChartRadianceBatch`). */
+    evaluateChartRadiance(input: ChartRadianceBatchInput, options?: ImageKernelRunOptions): Promise<ChartRadianceResult>;
+    dispose(): void;
+}
+/** Requests a backend. Never swaps backends silently: a fallback must be set in the options and is shown in every result label. */
+export declare function createImageKernelAccelerator(options: ImageKernelAcceleratorOptions): Promise<ImageKernelAccelerator>;
+/**
+ * A-priori float32 bound for every convolution output (all terms are non-negative):
+ * |v_f32 − v_ref| ≤ (γ_{T+3} + 2^-40)·v_ref + 2T·2^-126 + 2·S_max·Σ_t max_p w_{p,t}·[0 < w_{p,t} < 2^-126].
+ * The last term covers weights below the float32 normal range: the host conversion may round them to zero or
+ * a subnormal (|f32(w) − w| ≤ w) and the GPU may flush them, so their whole product w·s (s ≤ S_max) is at risk,
+ * with a factor 2 for the conversion and the later roundings. It is zero when no such weight exists.
+ */
+export declare function irradianceConvolutionErrorBound(plan: IrradianceConvolutionPlan, kernels: Float32Array | Float64Array): ImageKernelRelativeErrorBound;
+/** Uniform values of a tone curve for the float32 kernels (see wgsl-kernels.ts, LUMINANCE_FUNCTIONS). */
+interface ToneUniform {
+    mode: number;
+    tone: [number, number, number, number];
+    lum: [number, number, number, number];
+}
+/**
+ * The tone curve's kernel constants. Luminance curves run on the GPU within this float32 domain:
+ * extended Reinhard with Lw ≥ 2^-20; Naka–Rushton with exponent n ≥ 1 (below 1 its ratio f(L)/L is
+ * unbounded as L → 0, so no useful bound exists near L = 0), 1/g(Lw) ≤ 2^60 and σⁿ ≥ 2^-60.
+ * @internal
+ */
+export declare function toneUniform(profile: Pick<SdrRenderingProfile, "toneCurve" | "toneCurveParameters">): ToneUniform;
+/** Inputs of the a-priori float32 forward-error bound for packed RAW development (see docs/WEBGPU.md). */
+export interface RawDevelopmentBoundInput {
+    maximumTaps: number;
+    maximumWeightSum: number;
+    maximumSiteMagnitude: number;
+    maximumGain: number;
+    reconstruct: boolean;
+    cameraToXyzNorm: number;
+    xyzToRgbNorm: number;
+    scale: number;
+}
+/**
+ * Forward-error bounds (Higham γ_k, binary32) through the development chain: channel sums → gains →
+ * highlight clip/blend → camera→XYZ → XYZ→linear sRGB → rendering exposure. Returns the bound on the
+ * exposed value, and the classification bands for the white-balance clip and reconstruction-start
+ * decisions. Internal; exported only for unit tests.
+ */
+export declare function rawDevelopmentErrorBounds(b: RawDevelopmentBoundInput): {
+    exposed: number;
+    gain: number;
+    channel: number;
+};
+export {};