@photivra/engine 2.20.0 → 2.22.0

package/dist/gpu/display-view-image-material-module.js

@@ -0,0 +1,361 @@
+// SPDX-License-Identifier: Apache-2.0
+/**
+ * The image-plane module of the GPU display view (#581): the Speed kernel's `scene_value` (the seam of #478 and #491,
+ * `DisplayViewGpuMaterialModule`) for an image-plane print, so that a photo renders on the device like a chart.
+ *
+ * **What the reference computes** at a query point (expected-signal.ts, the direct path's `signalAt`): the print-local point
+ * cu = u·m/h, cv = v·m/h (outside |cu| ≤ 1, |cv| ≤ aspect the surround: its reflectance times the neutral signal N_c); the
+ * pixel coordinates X = (cu + 1)·W/2, Y = (aspect − cv)·H/(2·aspect); then the bilinear weights of the four nearest pixel
+ * centers (clamped to the edge) times each pixel's channel signal S_j,c = Σ_n A_c(λ_n)·ρ_j(λ_n)·W_n, integrated from the pixel's
+ * reconstructed spectrum ρ_j (its cell's eight node spectra blended, then the paper; `ResolvedImagePlaneScene.pixelSignals`), or
+ * the nearest pixel's. No per-channel values of the lattice nodes are formed.
+ *
+ * **On the device.** The host computes every pixel's S_j,c in float64 and rounds it once to float32 (3 values per pixel). When
+ * they fit one storage binding they are uploaded once per scene and kept on the device while the renderer lives; otherwise each
+ * band of a render builds and uploads only its rows. Either way they are built and written 4 MiB of rows at a time, so the host
+ * never holds more than that of S (owner decision 2026-11-09; the bytes are those of the whole array). The module recovers the exact lattice point from the
+ * kernel's float32 point as the composite module does (`LATTICE_POINT_WGSL`, about 2⁻⁴⁶ relative), forms cu, cv, X and Y
+ * in double-float32 and decides the print's edge in the composite kernel's bound arithmetic (a point within its bound of the
+ * edge is flagged: the display pixel is recomputed on the host). Per axis it splits s = X − ½ into its integer part n and
+ * fraction f exactly (double-float32), rounds the smaller of f and 1 − f once to float32 and forms the other as 1 minus it, so
+ * each weight carries at most 2E relative error (E = 2⁻²³) besides the point's own error. The point's error moves weight
+ * between neighbors: with δ the bound of |X_gpu − X_ref| in pixels (4·(2⁻⁴⁰·(|X| + 1)) + 4·10⁻¹⁴·W/2·k: the recovery, the
+ * double-float32 steps and the reference's own float64 rounding of its query point, generously), the value moves by at most
+ * (δx + δy)·2·max S, and the point is flagged when that exceeds E times the value. So an unflagged value carries a relative
+ * error only: S's rounding 1, the weights 2 + 2, the two products and two sums 4, the flag's budget 1, that is 10 roundings
+ * beyond the table lookup's one (`extraSiteRoundings`), added to the Speed kernel's chains (docs/WEBGPU.md, "Display view (Speed)").
+ * `nearest` takes the pixel under the point, flagging a point within its bound of a pixel edge (`bfloor`). All S are
+ * non-negative (non-negative responses, illuminant and spectra), so the weighted sums need no cancellation term.
+ *
+ * The surround's three channel signals, integrated from its spectrum (its reflectance times the white radiance at each node) like
+ * a pixel's, are the kernel's binding 3 (`spectral_table`, which this module does not otherwise use).
+ */
+import { ImageKernelExecutionError } from "./image-kernel-execution.js";
+import { BUFFER_USAGE, GpuResourceScope } from "./webgpu-runtime.js";
+import { dwSplit } from "./chart-composite-gpu.js";
+import { displayViewModuleCommonWgsl, displayViewLatticeConstants, LATTICE_BINDINGS_WGSL, LATTICE_POINT_WGSL } from "./display-view-composite-material-module.js";
+import { renderDisplayView } from "../facade/display-view.js";
+import { resolveCaptureIlluminant } from "../facade/scene-radiance.js";
+import { resolveMixedIllumination } from "../facade/mixed-illumination.js";
+import { imagePlaneSupport } from "../facade/image-plane-support.js";
+import { CONFORMANCE_FRAME, conformanceChannelAt, conformanceResponses, conformanceSetup, conformanceVerdict, renderDisplayViewSpeedOnGpu, withConformanceAbort } from "./display-view-speed-gpu.js";
+import { renderDisplayViewQualityOnGpu } from "./display-view-quality-gpu.js";
+import { imagePlaneQualityPasses } from "./display-view-image-quality.js";
+import { ChartCompositeProbeError, probeChartCompositeDevice } from "./chart-composite-gpu.js";
+import { createComputePipeline } from "./webgpu-runtime.js";
+/** The module's id per product path. */
+export function imagePlaneMaterialModuleId(path) { return path === "dekker" ? "image-plane" : "image-plane-fma"; }
+/** Roundings of the module's value beyond the table lookup's one (see the module comment). */
+export const IMAGE_PLANE_EXTRA_SITE_ROUNDINGS = 10;
+/** Layout of the photo uniform (float32 / u32 words, 64 bytes). */
+export const IMAGE_PLANE_PARAMS = /* @__PURE__ */ Object.freeze({ width: 0, height: 1, nearest: 2, row0: 3, halfWHi: 4, halfWLo: 5, halfHHi: 6, halfHLo: 7, aspectHi: 8, aspectLo: 9,
+    slackX: 10, slackY: 11, words: 16 });
+/** Table values above this are outside the float32 domain. */
+const FLOAT32_SAFE_MAGNITUDE = 2 ** 120;
+/** Host bytes of pixel signals built per write (whole uploads and bands). */
+const UPLOAD_CHUNK_BYTES = 4 << 20;
+/** Scenes whose pixel signals stay uploaded per renderer. */
+const KEPT_SCENES = 2;
+/** The module's WGSL with one product path (see the module comment). */
+export function imagePlaneMaterialWgsl(path = "dekker") {
+    return `${displayViewModuleCommonWgsl(path)}
+@group(1) @binding(0) var<storage, read> photo_signal: array<f32>;
+struct PhotoParams {
+  width: u32, height: u32, nearest: u32, row0: u32,
+  half_w_hi: f32, half_w_lo: f32, half_h_hi: f32, half_h_lo: f32,
+  aspect_hi: f32, aspect_lo: f32, slack_x: f32, slack_y: f32,
+  pad1: f32, pad2: f32, pad3: f32, pad4: f32,
+};
+@group(1) @binding(1) var<uniform> photo: PhotoParams;
+${LATTICE_BINDINGS_WGSL}
+${LATTICE_POINT_WGSL}
+// The point's error bound in pixels along one axis (see the module comment).
+fn pixel_slack(x: D, slack: f32) -> f32 { return 4.0 * (DW * (abs(x.h) + 1.0)) + 4.0 * slack; }
+// One axis of the bilinear weights: the pixel pair (clamped to the edge) and the weights of s = X − 1/2.
+struct Axis { i0: u32, i1: u32, w0: f32, w1: f32 }
+fn bilinear_axis(s: D, n: u32) -> Axis {
+  var f0 = floor(s.h);
+  if (f0 == s.h && s.l < 0.0) { f0 = f0 - 1.0; }
+  let frac = d_sub(s, D(f0, 0.0));
+  let comp = d_sub(D(ob(f0 + 1.0), 0.0), s);
+  var w0 = 0.0; var w1 = 0.0;
+  if (frac.h <= 0.5) { w1 = ob(frac.h + frac.l); w0 = ob(1.0 - w1); } else { w0 = ob(comp.h + comp.l); w1 = ob(1.0 - w0); }
+  let last = i32(n) - 1;
+  return Axis(u32(clamp(i32(f0), 0, last)), u32(clamp(i32(f0) + 1, 0, last)), w0, w1);
+}
+fn photo_at(i: u32, j: u32, c: u32) -> f32 { return photo_signal[((j - photo.row0) * photo.width + i) * 3u + c]; }
+fn scene_value(p: vec2<f32>, c: u32) -> f32 {
+  let k = D(lattice.chart_k_hi, lattice.chart_k_lo);
+  let u = d_mul(lattice_point(p.x, D(lattice.scale_u_hi, lattice.scale_u_lo), 0u), k);
+  let v = d_mul(lattice_point(p.y, D(lattice.scale_v_hi, lattice.scale_v_lo), 1u), k);
+  let aspect = D(photo.aspect_hi, photo.aspect_lo);
+  if (!(le(babs(bd(u)), bx(1.0)) && le(babs(bd(v)), bd(aspect)))) { return spectral_table[c]; }
+  let x = d_mul(d_add(u, D(1.0, 0.0)), D(photo.half_w_hi, photo.half_w_lo));
+  let y = d_mul(d_sub(aspect, v), D(photo.half_h_hi, photo.half_h_lo));
+  if (photo.nearest != 0u) {
+    let last_x = f32(photo.width - 1u);
+    let last_y = f32(photo.height - 1u);
+    let i = u32(clamp(bfloor(bd(x), 1.0, last_x), 0.0, last_x));
+    let j = u32(clamp(bfloor(bd(y), 1.0, last_y), 0.0, last_y));
+    return photo_at(i, j, c);
+  }
+  let ax = bilinear_axis(d_sub(x, D(0.5, 0.0)), photo.width);
+  let ay = bilinear_axis(d_sub(y, D(0.5, 0.0)), photo.height);
+  let s00 = photo_at(ax.i0, ay.i0, c);
+  let s10 = photo_at(ax.i1, ay.i0, c);
+  let s01 = photo_at(ax.i0, ay.i1, c);
+  let s11 = photo_at(ax.i1, ay.i1, c);
+  let value = ob(ob(ay.w0 * ob(ob(ax.w0 * s00) + ob(ax.w1 * s10))) + ob(ay.w1 * ob(ob(ax.w0 * s01) + ob(ax.w1 * s11))));
+  let largest = max(max(s00, s10), max(s01, s11));
+  let moved = (pixel_slack(x, photo.slack_x) + pixel_slack(y, photo.slack_y)) * 2.0 * largest;
+  if (!(moved <= E * value)) { material_uncertain = true; }
+  return value;
+}
+`;
+}
+/** The image-plane print behind a GPU setup's expected signal (`ImagePlaneSignalScene.image`), or undefined. */
+export function imagePlaneOf(gpu) {
+    return gpu.signalInput.scene.image;
+}
+/**
+ * The image-plane modules of a renderer; `probe` is the composite module's device probe (single flight), whose double-word product
+ * path the module's double-float32 takes.
+ */
+export function createDisplayViewImageModules(runtime, probe) {
+    const uploads = [], wgsl = new Map();
+    // Pixel rows [j0, j1)'s channel signals, built on demand so that the host holds one chunk or band (#581): each pixel's signals
+    // integrated from its reconstructed spectrum (`pixelSignals`, the reference's own float64 values), rounded once.
+    const pixel = new Float64Array(3);
+    const rowsOf = (image, j0, j1) => {
+        const { scene, nodes } = image, values = new Float32Array((j1 - j0) * scene.width * 3);
+        for (let j = j0 * scene.width, o = 0; j < j1 * scene.width; j++) {
+            scene.pixelSignals(j, nodes, pixel, 0);
+            for (let c = 0; c < 3; c++) {
+                const value = pixel[c];
+                if (!(value >= 0 && value <= FLOAT32_SAFE_MAGNITUDE))
+                    throw new ImageKernelExecutionError("unsupported-float32-domain", "A pixel's channel signal is outside the float32 domain.");
+                values[o++] = value;
+            }
+        }
+        return values;
+    };
+    // Rows [a, b)'s signals written to `buffer` (whose first row is `base`), built UPLOAD_CHUNK_BYTES of rows at a time.
+    const fill = (buffer, image, a, b, base) => {
+        const rowBytes = image.scene.width * 12, chunk = Math.max(1, Math.floor(UPLOAD_CHUNK_BYTES / rowBytes));
+        for (let j = a; j < b; j += chunk) {
+            const rows = rowsOf(image, j, Math.min(b, j + chunk));
+            runtime.device.queue.writeBuffer(buffer, (j - base) * rowBytes, rows.buffer, 0, rows.byteLength);
+        }
+    };
+    const rowsBuffer = (scope, label, image, j0, j1) => {
+        const buffer = scope.buffer(label, (j1 - j0) * image.scene.width * 12, BUFFER_USAGE.STORAGE | BUFFER_USAGE.COPY_DST);
+        fill(buffer, image, j0, j1, j0);
+        return buffer;
+    };
+    // The latest band's rows stay on the device until the next band: the rows the two share are copied there on the device, and
+    // only the others are built on the host (the same bytes; the host still holds one chunk at most).
+    let lastBand = null;
+    const bandBuffer = (image, j0, j1) => {
+        const rowBytes = image.scene.width * 12, scope = new GpuResourceScope(runtime), last = lastBand?.image === image ? lastBand : null;
+        try {
+            const buffer = scope.buffer("photivra-display-view-image-signal-band", (j1 - j0) * rowBytes, BUFFER_USAGE.STORAGE | BUFFER_USAGE.COPY_DST | BUFFER_USAGE.COPY_SRC);
+            const o0 = last ? Math.max(j0, last.j0) : j0, o1 = last ? Math.min(j1, last.j1) : j0;
+            if (o1 > o0) {
+                const encoder = runtime.device.createCommandEncoder({ label: "photivra-display-view-image-band-rows" });
+                encoder.copyBufferToBuffer(last.buffer, (o0 - last.j0) * rowBytes, buffer, (o0 - j0) * rowBytes, (o1 - o0) * rowBytes);
+                runtime.device.queue.submit([encoder.finish()]);
+                fill(buffer, image, j0, o0, j0);
+                fill(buffer, image, o1, j1, j0);
+            }
+            else
+                fill(buffer, image, j0, j1, j0);
+            lastBand?.scope.release();
+            lastBand = { image, j0, j1, scope, buffer };
+            return buffer;
+        }
+        catch (error) {
+            scope.release();
+            throw error;
+        }
+    };
+    const upload = (image) => {
+        const key = image, index = uploads.findIndex((e) => e.key === key);
+        if (index >= 0) {
+            const [entry] = uploads.splice(index, 1);
+            uploads.push(entry);
+            return entry;
+        }
+        const scope = new GpuResourceScope(runtime);
+        let buffer;
+        try {
+            buffer = rowsBuffer(scope, "photivra-display-view-image-signal", image, 0, image.scene.height);
+        }
+        catch (error) {
+            scope.release();
+            throw error;
+        }
+        const entry = { key, scope, buffer, surround: Float32Array.from(image.surround) };
+        uploads.push(entry);
+        while (uploads.length > KEPT_SCENES)
+            uploads.shift().scope.release();
+        return entry;
+    };
+    let checked = null;
+    const build = (gpu, path) => {
+        const image = imagePlaneOf(gpu);
+        if (!image || !gpu.placement)
+            throw new ImageKernelExecutionError("unsupported-on-backend", "The setup carries no image-plane print and placement; use cpu-reference.");
+        const scene = image.scene, lattice = displayViewLatticeConstants(gpu.signalInput, gpu.placement, 0), rowBytes = scene.width * 12, banded = scene.height * rowBytes > runtime.maximumStorageBindingBytes;
+        const params = new ArrayBuffer(4 * IMAGE_PLANE_PARAMS.words), words = new Uint32Array(params), floats = new Float32Array(params), P = IMAGE_PLANE_PARAMS;
+        words[P.width] = scene.width;
+        words[P.height] = scene.height;
+        words[P.nearest] = scene.interpolation === "nearest" ? 1 : 0;
+        const halfW = scene.width / 2, halfH = scene.height / (2 * scene.aspect);
+        [floats[P.halfWHi], floats[P.halfWLo]] = dwSplit(halfW);
+        [floats[P.halfHHi], floats[P.halfHLo]] = dwSplit(halfH);
+        [floats[P.aspectHi], floats[P.aspectLo]] = dwSplit(scene.aspect);
+        // The reference's float64 query point carries up to 10⁻¹⁴ frame units (the composite kernel's slack) times the pixels per frame unit.
+        floats[P.slackX] = 1e-14 * gpu.placement.chartUnitsPerFrameUnit * halfW;
+        floats[P.slackY] = 1e-14 * gpu.placement.chartUnitsPerFrameUnit * halfH;
+        let code = wgsl.get(path);
+        if (code === undefined) {
+            code = imagePlaneMaterialWgsl(path);
+            wgsl.set(path, code);
+        }
+        const base = { id: imagePlaneMaterialModuleId(path), wgsl: code, extraSiteRoundings: IMAGE_PLANE_EXTRA_SITE_ROUNDINGS, siteTable: Float32Array.from(image.surround) };
+        if (!banded) {
+            const uploaded = upload(image);
+            return { ...base, bindings: [{ binding: 0, kind: "storage", buffer: uploaded.buffer }, { binding: 1, kind: "uniform", data: new Uint8Array(params) },
+                    { binding: 2, kind: "uniform", data: new Uint8Array(lattice.params) }, { binding: 3, kind: "storage", data: lattice.offsets }] };
+        }
+        // A photo larger than one binding: per band the pixel rows its native site rows reach. The print-local v of site row y and
+        // query offset o is (y + c_o)·scale_v·k; the photo row is (aspect − v)·H/(2·aspect), plus two rows of margin either side.
+        const offsetsV = Array.from({ length: lattice.offsets.length / 4 }, (_, o) => lattice.offsets[4 * o + 2] + lattice.offsets[4 * o + 3]);
+        const rowOf = (y, c) => (scene.aspect - (y + c) * lattice.scaleV * gpu.placement.chartUnitsPerFrameUnit) * halfH;
+        const rate = Math.abs(lattice.scaleV) * gpu.placement.chartUnitsPerFrameUnit * halfH, maximumRows = Math.floor(runtime.maximumStorageBindingBytes / rowBytes);
+        return { ...base, bindings: [],
+            bandRows: (sitesPerDisplayRow) => Math.max(1, Math.floor(((maximumRows - 6) / rate - 6) / sitesPerDisplayRow)),
+            band(scope, layout, y0, y1) {
+                let low = Infinity, high = -Infinity;
+                for (const y of [y0, y1])
+                    for (const c of offsetsV) {
+                        const r = rowOf(y, c);
+                        low = Math.min(low, r);
+                        high = Math.max(high, r);
+                    }
+                const j0 = Math.max(0, Math.min(scene.height - 1, Math.floor(low) - 2)), j1 = Math.max(j0, Math.min(scene.height - 1, Math.floor(high) + 2));
+                if ((j1 - j0 + 1) * rowBytes > runtime.maximumStorageBindingBytes)
+                    throw new ImageKernelExecutionError("insufficient-limits", "A band's pixel signals exceed the device's storage binding limit.");
+                const bandParams = params.slice(0);
+                new Uint32Array(bandParams)[P.row0] = j0;
+                return runtime.device.createBindGroup({ layout, entries: [
+                        { binding: 0, resource: { buffer: bandBuffer(image, j0, j1 + 1) } },
+                        { binding: 1, resource: { buffer: scope.upload("photivra-display-view-image-params", new Uint8Array(bandParams), BUFFER_USAGE.UNIFORM) } },
+                        { binding: 2, resource: { buffer: scope.upload("photivra-display-view-image-lattice", new Uint8Array(lattice.params), BUFFER_USAGE.UNIFORM) } },
+                        { binding: 3, resource: { buffer: scope.upload("photivra-display-view-image-offsets", lattice.offsets, BUFFER_USAGE.STORAGE) } }
+                    ] });
+            } };
+    };
+    // The module's self-check, once per runtime (single flight; a failure is retried on the next photo).
+    const check = async () => {
+        checked ??= runImagePlaneConformance(runtime, undefined).catch((error) => { checked = null; throw error; });
+        const verdict = await checked;
+        if (verdict !== true) {
+            checked = null;
+            throw new ImageKernelExecutionError("unsupported-on-backend", "The image-plane module failed its self-check on this device: " + verdict);
+        }
+    };
+    return {
+        async moduleFor(gpu, signal) { const path = (await probe(signal)).productPath ?? "dekker"; await check(); return build(gpu, path); },
+        async qualityPasses(gpu, signal) { const path = (await probe(signal)).productPath ?? "dekker"; await check(); return imagePlaneQualityPasses(runtime, gpu, path); },
+        build,
+        release() { for (const entry of uploads.splice(0))
+            entry.scope.release(); lastBand?.scope.release(); lastBand = null; }
+    };
+}
+/**
+ * #581's run-time self-check of the module, run on a device's first photo: a 12 × 8 photo (a color ramp with hard edges, every
+ * lattice cell's corner colors among its pixels) printed at 0.8 of the frame's width on the 18 % surround of a 38 × 26 frame shown
+ * at 13 × 9, the GPU against `renderDisplayView` within the a-priori bounds, every flagged pixel recomputed: the Speed kernel with
+ * the 2 × 2 box, then the Quality passes with the four-spot filter's 16 offsets and synthetic radius-1 kernels at three nodes, and
+ * (#597) the same passes under mixed light (a three-source gradient with a hard edge and a point source).
+ * A device that fails the double-float32 probe refuses photos (unsupported-on-backend, as composite charts), so the check passes
+ * there without keeping the device's charts off the GPU.
+ */
+export async function runImagePlaneConformance(handle, signal) {
+    const runtime = handle, support = imagePlaneSupport();
+    if (!support)
+        return true;
+    const { nativeWidth, nativeHeight, frameWidth, nodes } = CONFORMANCE_FRAME, g = 2 / frameWidth;
+    const width = 12, height = 8, data = new Uint8Array(width * height * 3);
+    for (let y = 0; y < height; y++)
+        for (let x = 0; x < width; x++)
+            for (let c = 0; c < 3; c++)
+                data[(y * width + x) * 3 + c] = [(x * 23) % 256, (y * 37 + x * 5) % 256, x < 6 ? 255 - y * 30 : 10][c];
+    const printWidth = 1.6, aspect = height / width, k = 2 / printWidth, surround = 0.18, placement = { chartUnitsPerFrameUnit: k, aspect };
+    const scene = await support.resolveImagePlaneScene({ kind: "image-plane", sceneId: "conformance", image: { width, height, channels: 3, encoding: "srgb-8", data } }, aspect);
+    const white = Float64Array.of(0.5, 0.9, 0.7), A = conformanceResponses(1);
+    // The print's radiance as the provider forms it: the print-local point, the surround outside, times the white radiance at the node.
+    const radianceAt = (u, v, nm) => {
+        const cu = u * k, cv = v * k;
+        return (Math.abs(cu) <= 1 && Math.abs(cv) <= aspect ? scene.reflectanceAt(cu, cv, nm) : surround) * white[nodes.indexOf(nm)];
+    };
+    const printScene = scene.signalScene(nodes, A, white, radianceAt, k, 1, surround), readout = { perElectron: 1 / 60000, darkElectrons: 0.07, constantVariance: 1.5e-9 };
+    const speedOffsets = [[-0.25, -0.25], [0.25, -0.25], [-0.25, 0.25], [0.25, 0.25]].map(([x, y]) => ({ x: x, y: y, weight: 0.25 }));
+    const speedInput = { nativeWidth, nativeHeight, outputWidth: frameWidth, wavelengthsNm: nodes, electronsPerRadiance: A, offsets: speedOffsets, scene: printScene,
+        channelAt: conformanceChannelAt, queryPoint: (x, y, o, out) => { out[0] = (x + 0.5 + speedOffsets[o].x - nativeWidth / 2) * g; out[1] = -(y + 0.5 + speedOffsets[o].y - nativeHeight / 2) * g; } };
+    const speed = conformanceSetup(speedInput, readout, 581, 2, signal, { gpu: { signalInput: speedInput, chart: null, placement } });
+    const qualityOffsets = [];
+    for (const [nx, ny] of [[-0.25, -0.25], [0.25, -0.25], [-0.25, 0.25], [0.25, 0.25]])
+        for (const [ax, ay] of [[-0.5, -0.5], [0.5, -0.5], [-0.5, 0.5], [0.5, 0.5]])
+            qualityOffsets.push({ x: nx - ax, y: ny - ay, weight: 1 / 16 });
+    const kernel = (sigma) => { const w = Float64Array.from({ length: 9 }, (_, t) => Math.exp(-0.5 * (((t % 3) - 1) ** 2 + (Math.floor(t / 3) - 1) ** 2) / (sigma * sigma))); const sum = w.reduce((a, b) => a + b, 0); return w.map((v) => v / sum); };
+    const qualityInput = { nativeWidth, nativeHeight, outputWidth: frameWidth, wavelengthsNm: nodes, electronsPerRadiance: A, offsets: qualityOffsets,
+        kernels: { radius: 1, perNode: [kernel(0.6), kernel(0.7), kernel(0.8)] }, scene: printScene, channelAt: conformanceChannelAt };
+    const quality = conformanceSetup(qualityInput, readout, 582, 5, signal, { gpu: { signalInput: qualityInput, chart: null, placement } });
+    // #597: the same Quality case under mixed light: a gradient of three sources with a hard edge (equal stops) across the print, and a
+    // point source over its upper left, the reference's irradiance from the same parts in float64.
+    const day = { kind: "planckian", temperatureKelvin: 6500 }, lamp = { kind: "cie-a" };
+    const light = resolveMixedIllumination({ sources: [{ illuminant: day, illuminanceLux: 1.2e5 }, { illuminant: lamp, illuminanceLux: 8e4 }, { illuminant: day, illuminanceLux: 3e4 }],
+        blend: { kind: "linear-gradient", from: { u: -0.9, v: 0.1 }, to: { u: 0.8, v: -0.2 }, stops: [0, 0.45, 0.45] },
+        points: [{ illuminant: lamp, position: { xM: -0.3, yM: 0.2, heightM: 0.35 }, intensity: { illuminanceLux: 6e4, distanceM: 0.35 } }] }, resolveCaptureIlluminant);
+    // The print's half width is 1 m here (print-local units are metres), so the point sources read the print-local point.
+    const litIrradiance = (cu, cv, nm) => light.irradianceAt(cu, cv, cu, cv, nm);
+    const litRadianceAt = (u, v, nm) => {
+        const cu = u * k, cv = v * k;
+        return (Math.abs(cu) <= 1 && Math.abs(cv) <= aspect ? scene.reflectanceAt(cu, cv, nm) : surround) * litIrradiance(cu, cv, nm) / Math.PI;
+    };
+    const litInput = { ...qualityInput, scene: scene.signalScene(nodes, A, white, litRadianceAt, k, 1, surround, light) };
+    const lit = conformanceSetup(litInput, readout, 597, 5, signal, { gpu: { signalInput: litInput, chart: null, placement } });
+    const [speedReference, qualityReference, litReference] = [await renderDisplayView(speed), await renderDisplayView(quality), await renderDisplayView(lit)];
+    const images = createDisplayViewImageModules(runtime, (s) => probeChartCompositeDevice(runtime, (label, code) => createComputePipeline(runtime, label, code, s), s));
+    try {
+        return await withConformanceAbort(runtime, signal, async (abort, pipelines) => {
+            let path;
+            try {
+                path = (await probeChartCompositeDevice(runtime, pipelines, abort.signal)).productPath ?? "dekker";
+            }
+            catch (error) {
+                if (error instanceof ChartCompositeProbeError)
+                    return true;
+                throw error;
+            }
+            const verdict = conformanceVerdict(speedReference, await renderDisplayViewSpeedOnGpu(runtime, 1 << 20, { setup: speed, signalInput: speedInput, material: images.build(speed.gpu, path) }, abort, undefined, pipelines));
+            if (verdict !== true)
+                return "Speed: " + verdict;
+            const gpu = await renderDisplayViewQualityOnGpu(runtime, 1 << 20, { setup: quality, signalInput: qualityInput, material: { id: "image-plane", wgsl: "", bindings: [] },
+                image: imagePlaneQualityPasses(runtime, quality.gpu, path) }, abort, undefined, pipelines);
+            const qualityVerdict = conformanceVerdict(qualityReference, gpu);
+            if (qualityVerdict !== true)
+                return "Quality: " + qualityVerdict;
+            const litGpu = await renderDisplayViewQualityOnGpu(runtime, 1 << 20, { setup: lit, signalInput: litInput, material: { id: "image-plane", wgsl: "", bindings: [] },
+                image: imagePlaneQualityPasses(runtime, lit.gpu, path) }, abort, undefined, pipelines);
+            const litVerdict = conformanceVerdict(litReference, litGpu);
+            return litVerdict === true ? true : "Quality under mixed light: " + litVerdict;
+        });
+    }
+    finally {
+        images.release();
+    }
+}
+//# sourceMappingURL=display-view-image-material-module.js.map

package/dist/gpu/display-view-image-material-module.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"display-view-image-material-module.js","sourceRoot":"","sources":["../../src/gpu/display-view-image-material-module.ts"],"names":[],"mappings":"AAAA,sCAAsC;AAEtC;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA8BG;AAEH,OAAO,EAAE,yBAAyB,EAAkC,MAAM,6BAA6B,CAAC;AACxG,OAAO,EAAE,YAAY,EAAE,gBAAgB,EAAsB,MAAM,qBAAqB,CAAC;AACzF,OAAO,EAAE,OAAO,EAAE,MAAM,0BAA0B,CAAC;AAEnD,OAAO,EAAE,2BAA2B,EAAE,2BAA2B,EAAE,qBAAqB,EAAE,kBAAkB,EAAE,MAAM,6CAA6C,CAAC;AAElK,OAAO,EAAE,iBAAiB,EAA4B,MAAM,2BAA2B,CAAC;AAExF,OAAO,EAAE,wBAAwB,EAA8B,MAAM,6BAA6B,CAAC;AACnG,OAAO,EAAE,wBAAwB,EAAE,MAAM,iCAAiC,CAAC;AAC3E,OAAO,EAAE,iBAAiB,EAAE,MAAM,kCAAkC,CAAC;AACrE,OAAO,EAAE,iBAAiB,EAAE,oBAAoB,EAAE,oBAAoB,EAAE,gBAAgB,EAAE,kBAAkB,EAAE,2BAA2B,EAAE,oBAAoB,EAAE,MAAM,6BAA6B,CAAC;AACrM,OAAO,EAAE,6BAA6B,EAAsC,MAAM,+BAA+B,CAAC;AAClH,OAAO,EAAE,uBAAuB,EAAE,MAAM,iCAAiC,CAAC;AAC1E,OAAO,EAAE,wBAAwB,EAAE,yBAAyB,EAAE,MAAM,0BAA0B,CAAC;AAC/F,OAAO,EAAE,qBAAqB,EAAE,MAAM,qBAAqB,CAAC;AAE5D,wCAAwC;AACxC,MAAM,UAAU,0BAA0B,CAAC,IAA+B,IAAY,OAAO,IAAI,KAAK,QAAQ,CAAC,CAAC,CAAC,aAAa,CAAC,CAAC,CAAC,iBAAiB,CAAC,CAAC,CAAC;AACrJ,8FAA8F;AAC9F,MAAM,CAAC,MAAM,gCAAgC,GAAG,EAAE,CAAC;AACnD,mEAAmE;AACnE,MAAM,CAAC,MAAM,kBAAkB,GAAG,eAAe,CAAC,MAAM,CAAC,MAAM,CAAC,EAAE,KAAK,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,IAAI,EAAE,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,QAAQ,EAAE,CAAC,EAAE,QAAQ,EAAE,CAAC;IAClL,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,KAAK,EAAE,EAAE,EAAE,CAAC,CAAC;AACvC,8DAA8D;AAC9D,MAAM,sBAAsB,GAAG,CAAC,IAAI,GAAG,CAAC;AACxC,6EAA6E;AAC7E,MAAM,kBAAkB,GAAG,CAAC,IAAI,EAAE,CAAC;AACnC,6DAA6D;AAC7D,MAAM,WAAW,GAAG,CAAC,CAAC;AAEtB,wEAAwE;AACxE,MAAM,UAAU,sBAAsB,CAAC,OAAkC,QAAQ;IAC/E,OAAO,GAAG,2BAA2B,CAAC,IAAI,CAAC;;;;;;;;;EAS3C,qBAAqB;EACrB,kBAAkB;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;CA2CnB,CAAC;AACF,CAAC;AAED,gHAAgH;AAChH,MAAM,UAAU,YAAY,CAAC,GAAwB;IACnD,OAAQ,GAAG,CAAC,WAAW,CAAC,KAAoD,CAAC,KAAK,CAAC;AACrF,CAAC;AAeD;;;GAGG;AACH,MAAM,UAAU,6BAA6B,CAAC,OAAsB,EAAE,KAAkE;IACtI,MAAM,OAAO,GAAoB,EAAE,EAAE,IAAI,GAAG,IAAI,GAAG,EAAqC,CAAC;IACzF,+HAA+H;IAC/H,iHAAiH;IACjH,MAAM,KAAK,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;IAClC,MAAM,MAAM,GAAG,CAAC,KAAqC,EAAE,EAAU,EAAE,EAAU,EAA6B,EAAE;QAC1G,MAAM,EAAE,KAAK,EAAE,KAAK,EAAE,GAAG,KAAK,EAAE,MAAM,GAAG,IAAI,YAAY,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,KAAK,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC;QACvF,KAAK,IAAI,CAAC,GAAG,EAAE,GAAG,KAAK,CAAC,KAAK,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,GAAG,KAAK,CAAC,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;YAChE,KAAK,CAAC,YAAY,CAAC,CAAC,EAAE,KAAK,EAAE,KAAK,EAAE,CAAC,CAAC,CAAC;YACvC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,MAAM,KAAK,GAAG,KAAK,CAAC,CAAC,CAAE,CAAC;gBACxB,IAAI,CAAC,CAAC,KAAK,IAAI,CAAC,IAAI,KAAK,IAAI,sBAAsB,CAAC;oBAAE,MAAM,IAAI,yBAAyB,CAAC,4BAA4B,EAAE,yDAAyD,CAAC,CAAC;gBACnL,MAAM,CAAC,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC;YACtB,CAAC;QACH,CAAC;QACD,OAAO,MAAM,CAAC;IAChB,CAAC,CAAC;IACF,qHAAqH;IACrH,MAAM,IAAI,GAAG,CAAC,MAAiB,EAAE,KAAqC,EAAE,CAAS,EAAE,CAAS,EAAE,IAAY,EAAQ,EAAE;QAClH,MAAM,QAAQ,GAAG,KAAK,CAAC,KAAK,CAAC,KAAK,GAAG,EAAE,EAAE,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,kBAAkB,GAAG,QAAQ,CAAC,CAAC,CAAC;QACxG,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,KAAK,EAAE,CAAC;YAAC,MAAM,IAAI,GAAG,MAAM,CAAC,KAAK,EAAE,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;YAAC,OAAO,CAAC,MAAM,CAAC,KAAK,CAAC,WAAW,CAAC,MAAM,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,QAAQ,EAAE,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE,IAAI,CAAC,UAAU,CAAC,CAAC;QAAC,CAAC;IACjM,CAAC,CAAC;IACF,MAAM,UAAU,GAAG,CAAC,KAAuB,EAAE,KAAa,EAAE,KAAqC,EAAE,EAAU,EAAE,EAAU,EAAa,EAAE;QACtI,MAAM,MAAM,GAAG,KAAK,CAAC,MAAM,CAAC,KAAK,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,KAAK,CAAC,KAAK,CAAC,KAAK,GAAG,EAAE,EAAE,YAAY,CAAC,OAAO,GAAG,YAAY,CAAC,QAAQ,CAAC,CAAC;QACrH,IAAI,CAAC,MAAM,EAAE,KAAK,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;QAChC,OAAO,MAAM,CAAC;IAChB,CAAC,CAAC;IACF,4HAA4H;IAC5H,kGAAkG;IAClG,IAAI,QAAQ,GAAiG,IAAI,CAAC;IAClH,MAAM,UAAU,GAAG,CAAC,KAAqC,EAAE,EAAU,EAAE,EAAU,EAAa,EAAE;QAC9F,MAAM,QAAQ,GAAG,KAAK,CAAC,KAAK,CAAC,KAAK,GAAG,EAAE,EAAE,KAAK,GAAG,IAAI,gBAAgB,CAAC,OAAO,CAAC,EAAE,IAAI,GAAG,QAAQ,EAAE,KAAK,KAAK,KAAK,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC;QACnI,IAAI,CAAC;YACH,MAAM,MAAM,GAAG,KAAK,CAAC,MAAM,CAAC,yCAAyC,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,QAAQ,EAAE,YAAY,CAAC,OAAO,GAAG,YAAY,CAAC,QAAQ,GAAG,YAAY,CAAC,QAAQ,CAAC,CAAC;YACnK,MAAM,EAAE,GAAG,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,EAAE,GAAG,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;YACrF,IAAI,EAAE,GAAG,EAAE,EAAE,CAAC;gBACZ,MAAM,OAAO,GAAG,OAAO,CAAC,MAAM,CAAC,oBAAoB,CAAC,EAAE,KAAK,EAAE,uCAAuC,EAAE,CAAC,CAAC;gBACxG,OAAO,CAAC,kBAAkB,CAAC,IAAK,CAAC,MAAM,EAAE,CAAC,EAAE,GAAG,IAAK,CAAC,EAAE,CAAC,GAAG,QAAQ,EAAE,MAAM,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,QAAQ,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,QAAQ,CAAC,CAAC;gBACzH,OAAO,CAAC,MAAM,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC,OAAO,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC;gBAChD,IAAI,CAAC,MAAM,EAAE,KAAK,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;gBAAC,IAAI,CAAC,MAAM,EAAE,KAAK,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;YACnE,CAAC;;gBAAM,IAAI,CAAC,MAAM,EAAE,KAAK,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;YACvC,QAAQ,EAAE,KAAK,CAAC,OAAO,EAAE,CAAC;YAC1B,QAAQ,GAAG,EAAE,KAAK,EAAE,EAAE,EAAE,EAAE,EAAE,KAAK,EAAE,MAAM,EAAE,CAAC;YAC5C,OAAO,MAAM,CAAC;QAChB,CAAC;QAAC,OAAO,KAAK,EAAE,CAAC;YAAC,KAAK,CAAC,OAAO,EAAE,CAAC;YAAC,MAAM,KAAK,CAAC;QAAC,CAAC;IACnD,CAAC,CAAC;IACF,MAAM,MAAM,GAAG,CAAC,KAAqC,EAAiB,EAAE;QACtE,MAAM,GAAG,GAAG,KAAe,EAAE,KAAK,GAAG,OAAO,CAAC,SAAS,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC;QAC7E,IAAI,KAAK,IAAI,CAAC,EAAE,CAAC;YAAC,MAAM,CAAC,KAAK,CAAC,GAAG,OAAO,CAAC,MAAM,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC;YAAC,OAAO,CAAC,IAAI,CAAC,KAAM,CAAC,CAAC;YAAC,OAAO,KAAM,CAAC;QAAC,CAAC;QAClG,MAAM,KAAK,GAAG,IAAI,gBAAgB,CAAC,OAAO,CAAC,CAAC;QAC5C,IAAI,MAAiB,CAAC;QACtB,IAAI,CAAC;YAAC,MAAM,GAAG,UAAU,CAAC,KAAK,EAAE,oCAAoC,EAAE,KAAK,EAAE,CAAC,EAAE,KAAK,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;QAAC,CAAC;QAAC,OAAO,KAAK,EAAE,CAAC;YAAC,KAAK,CAAC,OAAO,EAAE,CAAC;YAAC,MAAM,KAAK,CAAC;QAAC,CAAC;QACvJ,MAAM,KAAK,GAAkB,EAAE,GAAG,EAAE,KAAK,EAAE,MAAM,EAAE,QAAQ,EAAE,YAAY,CAAC,IAAI,CAAC,KAAK,CAAC,QAAQ,CAAC,EAAE,CAAC;QACjG,OAAO,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QACpB,OAAO,OAAO,CAAC,MAAM,GAAG,WAAW;YAAE,OAAO,CAAC,KAAK,EAAG,CAAC,KAAK,CAAC,OAAO,EAAE,CAAC;QACtE,OAAO,KAAK,CAAC;IACf,CAAC,CAAC;IACF,IAAI,OAAO,GAAkC,IAAI,CAAC;IAClD,MAAM,KAAK,GAAG,CAAC,GAAwB,EAAE,IAA+B,EAAgC,EAAE;QACtG,MAAM,KAAK,GAAG,YAAY,CAAC,GAAG,CAAC,CAAC;QAChC,IAAI,CAAC,KAAK,IAAI,CAAC,GAAG,CAAC,SAAS;YAAE,MAAM,IAAI,yBAAyB,CAAC,wBAAwB,EAAE,0EAA0E,CAAC,CAAC;QACxK,MAAM,KAAK,GAAG,KAAK,CAAC,KAAK,EAAE,OAAO,GAAG,2BAA2B,CAAC,GAAG,CAAC,WAAW,EAAE,GAAG,CAAC,SAAS,EAAE,CAAC,CAAC,EAAE,QAAQ,GAAG,KAAK,CAAC,KAAK,GAAG,EAAE,EAAE,MAAM,GAAG,KAAK,CAAC,MAAM,GAAG,QAAQ,GAAG,OAAO,CAAC,0BAA0B,CAAC;QACxM,MAAM,MAAM,GAAG,IAAI,WAAW,CAAC,CAAC,GAAG,kBAAkB,CAAC,KAAK,CAAC,EAAE,KAAK,GAAG,IAAI,WAAW,CAAC,MAAM,CAAC,EAAE,MAAM,GAAG,IAAI,YAAY,CAAC,MAAM,CAAC,EAAE,CAAC,GAAG,kBAAkB,CAAC;QACzJ,KAAK,CAAC,CAAC,CAAC,KAAK,CAAC,GAAG,KAAK,CAAC,KAAK,CAAC;QAAC,KAAK,CAAC,CAAC,CAAC,MAAM,CAAC,GAAG,KAAK,CAAC,MAAM,CAAC;QAAC,KAAK,CAAC,CAAC,CAAC,OAAO,CAAC,GAAG,KAAK,CAAC,aAAa,KAAK,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC3H,MAAM,KAAK,GAAG,KAAK,CAAC,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,KAAK,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,MAAM,CAAC,CAAC;QACzE,CAAC,MAAM,CAAC,CAAC,CAAC,OAAO,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC;QAAC,CAAC,MAAM,CAAC,CAAC,CAAC,OAAO,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC;QAAC,CAAC,MAAM,CAAC,CAAC,CAAC,QAAQ,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;QACnL,sIAAsI;QACtI,MAAM,CAAC,CAAC,CAAC,MAAM,CAAC,GAAG,KAAK,GAAG,GAAG,CAAC,SAAS,CAAC,sBAAsB,GAAG,KAAK,CAAC;QAAC,MAAM,CAAC,CAAC,CAAC,MAAM,CAAC,GAAG,KAAK,GAAG,GAAG,CAAC,SAAS,CAAC,sBAAsB,GAAG,KAAK,CAAC;QACjJ,IAAI,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;QAC1B,IAAI,IAAI,KAAK,SAAS,EAAE,CAAC;YAAC,IAAI,GAAG,sBAAsB,CAAC,IAAI,CAAC,CAAC;YAAC,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC;QAAC,CAAC;QACtF,MAAM,IAAI,GAAG,EAAE,EAAE,EAAE,0BAA0B,CAAC,IAAI,CAAC,EAAE,IAAI,EAAE,IAAI,EAAE,kBAAkB,EAAE,gCAAgC,EAAE,SAAS,EAAE,YAAY,CAAC,IAAI,CAAC,KAAK,CAAC,QAAQ,CAAC,EAAE,CAAC;QACtK,IAAI,CAAC,MAAM,EAAE,CAAC;YACZ,MAAM,QAAQ,GAAG,MAAM,CAAC,KAAK,CAAC,CAAC;YAC/B,OAAO,EAAE,GAAG,IAAI,EAAE,QAAQ,EAAE,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,IAAI,EAAE,SAAS,EAAE,MAAM,EAAE,QAAQ,CAAC,MAAM,EAAE,EAAE,EAAE,OAAO,EAAE,CAAC,EAAE,IAAI,EAAE,SAAS,EAAE,IAAI,EAAE,IAAI,UAAU,CAAC,MAAM,CAAC,EAAE;oBAClJ,EAAE,OAAO,EAAE,CAAC,EAAE,IAAI,EAAE,SAAS,EAAE,IAAI,EAAE,IAAI,UAAU,CAAC,OAAO,CAAC,MAAM,CAAC,EAAE,EAAE,EAAE,OAAO,EAAE,CAAC,EAAE,IAAI,EAAE,SAAS,EAAE,IAAI,EAAE,OAAO,CAAC,OAAO,EAAE,CAAC,EAAE,CAAC;QACrI,CAAC;QACD,2HAA2H;QAC3H,0HAA0H;QAC1H,MAAM,QAAQ,GAAG,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,OAAO,CAAC,OAAO,CAAC,MAAM,GAAG,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,OAAO,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAE,GAAG,OAAO,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAE,CAAC,CAAC;QACzI,MAAM,KAAK,GAAG,CAAC,CAAS,EAAE,CAAS,EAAU,EAAE,CAAC,CAAC,KAAK,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,MAAM,GAAG,GAAG,CAAC,SAAU,CAAC,sBAAsB,CAAC,GAAG,KAAK,CAAC;QAC1I,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,MAAM,CAAC,GAAG,GAAG,CAAC,SAAS,CAAC,sBAAsB,GAAG,KAAK,EAAE,WAAW,GAAG,IAAI,CAAC,KAAK,CAAC,OAAO,CAAC,0BAA0B,GAAG,QAAQ,CAAC,CAAC;QAC9J,OAAO,EAAE,GAAG,IAAI,EAAE,QAAQ,EAAE,EAAE;YAC5B,QAAQ,EAAE,CAAC,kBAA0B,EAAU,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,WAAW,GAAG,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC,GAAG,kBAAkB,CAAC,CAAC;YAC9H,IAAI,CAAC,KAAK,EAAE,MAAM,EAAE,EAAE,EAAE,EAAE;gBACxB,IAAI,GAAG,GAAG,QAAQ,EAAE,IAAI,GAAG,CAAC,QAAQ,CAAC;gBACrC,KAAK,MAAM,CAAC,IAAI,CAAC,EAAE,EAAE,EAAE,CAAC;oBAAE,KAAK,MAAM,CAAC,IAAI,QAAQ,EAAE,CAAC;wBAAC,MAAM,CAAC,GAAG,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;wBAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;wBAAC,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;oBAAC,CAAC;gBAChI,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,MAAM,GAAG,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,MAAM,GAAG,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;gBAC7I,IAAI,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,QAAQ,GAAG,OAAO,CAAC,0BAA0B;oBAAE,MAAM,IAAI,yBAAyB,CAAC,qBAAqB,EAAE,mEAAmE,CAAC,CAAC;gBACnM,MAAM,UAAU,GAAG,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;gBACnC,IAAI,WAAW,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,EAAE,CAAC;gBACzC,OAAO,OAAO,CAAC,MAAM,CAAC,eAAe,CAAC,EAAE,MAAM,EAAE,OAAO,EAAE;wBACvD,EAAE,OAAO,EAAE,CAAC,EAAE,QAAQ,EAAE,EAAE,MAAM,EAAE,UAAU,CAAC,KAAK,EAAE,EAAE,EAAE,EAAE,GAAG,CAAC,CAAC,EAAE,EAAE;wBACnE,EAAE,OAAO,EAAE,CAAC,EAAE,QAAQ,EAAE,EAAE,MAAM,EAAE,KAAK,CAAC,MAAM,CAAC,oCAAoC,EAAE,IAAI,UAAU,CAAC,UAAU,CAAC,EAAE,YAAY,CAAC,OAAO,CAAC,EAAE,EAAE;wBAC1I,EAAE,OAAO,EAAE,CAAC,EAAE,QAAQ,EAAE,EAAE,MAAM,EAAE,KAAK,CAAC,MAAM,CAAC,qCAAqC,EAAE,IAAI,UAAU,CAAC,OAAO,CAAC,MAAM,CAAC,EAAE,YAAY,CAAC,OAAO,CAAC,EAAE,EAAE;wBAC/I,EAAE,OAAO,EAAE,CAAC,EAAE,QAAQ,EAAE,EAAE,MAAM,EAAE,KAAK,CAAC,MAAM,CAAC,qCAAqC,EAAE,OAAO,CAAC,OAAO,EAAE,YAAY,CAAC,OAAO,CAAC,EAAE,EAAE;qBAAC,EAAE,CAAC,CAAC;YACzI,CAAC,EAAE,CAAC;IACV,CAAC,CAAC;IACF,qGAAqG;IACrG,MAAM,KAAK,GAAG,KAAK,IAAmB,EAAE;QACtC,OAAO,KAAK,wBAAwB,CAAC,OAAO,EAAE,SAAS,CAAC,CAAC,KAAK,CAAC,CAAC,KAAc,EAAE,EAAE,GAAG,OAAO,GAAG,IAAI,CAAC,CAAC,MAAM,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;QACrH,MAAM,OAAO,GAAG,MAAM,OAAO,CAAC;QAC9B,IAAI,OAAO,KAAK,IAAI,EAAE,CAAC;YAAC,OAAO,GAAG,IAAI,CAAC;YAAC,MAAM,IAAI,yBAAyB,CAAC,wBAAwB,EAAE,+DAA+D,GAAG,OAAO,CAAC,CAAC;QAAC,CAAC;IACrL,CAAC,CAAC;IACF,OAAO;QACL,KAAK,CAAC,SAAS,CAAC,GAAG,EAAE,MAAM,IAA2C,MAAM,IAAI,GAAG,CAAC,MAAM,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,WAAW,IAAI,QAAQ,CAAC,CAAC,MAAM,KAAK,EAAE,CAAC,CAAC,OAAO,KAAK,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;QAC3K,KAAK,CAAC,aAAa,CAAC,GAAG,EAAE,MAAM,IAA4C,MAAM,IAAI,GAAG,CAAC,MAAM,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,WAAW,IAAI,QAAQ,CAAC,CAAC,MAAM,KAAK,EAAE,CAAC,CAAC,OAAO,uBAAuB,CAAC,OAAO,EAAE,GAAG,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;QAC3M,KAAK;QACL,OAAO,KAAW,KAAK,MAAM,KAAK,IAAI,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC;YAAE,KAAK,CAAC,KAAK,CAAC,OAAO,EAAE,CAAC,CAAC,QAAQ,EAAE,KAAK,CAAC,OAAO,EAAE,CAAC,CAAC,QAAQ,GAAG,IAAI,CAAC,CAAC,CAAC;KAC9H,CAAC;AACJ,CAAC;AAED;;;;;;;;GAQG;AACH,MAAM,CAAC,KAAK,UAAU,wBAAwB,CAAC,MAAe,EAAE,MAA+B;IAC7F,MAAM,OAAO,GAAG,MAAuB,EAAE,OAAO,GAAG,iBAAiB,EAAE,CAAC;IACvE,IAAI,CAAC,OAAO;QAAE,OAAO,IAAI,CAAC;IAC1B,MAAM,EAAE,WAAW,EAAE,YAAY,EAAE,UAAU,EAAE,KAAK,EAAE,GAAG,iBAAiB,EAAE,CAAC,GAAG,CAAC,GAAG,UAAU,CAAC;IAC/F,MAAM,KAAK,GAAG,EAAE,EAAE,MAAM,GAAG,CAAC,EAAE,IAAI,GAAG,IAAI,UAAU,CAAC,KAAK,GAAG,MAAM,GAAG,CAAC,CAAC,CAAC;IACxE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,EAAE,CAAC,EAAE;QAAE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE;YAAE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;gBAAE,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAE,CAAC;IACrM,MAAM,UAAU,GAAG,GAAG,EAAE,MAAM,GAAG,MAAM,GAAG,KAAK,EAAE,CAAC,GAAG,CAAC,GAAG,UAAU,EAAE,QAAQ,GAAG,IAAI,EAAE,SAAS,GAAG,EAAE,sBAAsB,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC;IACxI,MAAM,KAAK,GAAG,MAAM,OAAO,CAAC,sBAAsB,CAAC,EAAE,IAAI,EAAE,aAAa,EAAE,OAAO,EAAE,aAAa,EAAE,KAAK,EAAE,EAAE,KAAK,EAAE,MAAM,EAAE,QAAQ,EAAE,CAAC,EAAE,QAAQ,EAAE,QAAQ,EAAE,IAAI,EAAE,EAAE,EAAE,MAAM,CAAC,CAAC;IAC7K,MAAM,KAAK,GAAG,YAAY,CAAC,EAAE,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,EAAE,CAAC,GAAG,oBAAoB,CAAC,CAAC,CAAC,CAAC;IAC1E,oIAAoI;IACpI,MAAM,UAAU,GAAG,CAAC,CAAS,EAAE,CAAS,EAAE,EAAU,EAAU,EAAE;QAC9D,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;QAC7B,OAAO,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,MAAM,CAAC,CAAC,CAAC,KAAK,CAAC,aAAa,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,GAAG,KAAK,CAAC,KAAK,CAAC,OAAO,CAAC,EAAE,CAAC,CAAE,CAAC;IAChI,CAAC,CAAC;IACF,MAAM,UAAU,GAAG,KAAK,CAAC,WAAW,CAAC,KAAK,EAAE,CAAC,EAAE,KAAK,EAAE,UAAU,EAAE,CAAC,EAAE,CAAC,EAAE,QAAQ,CAAC,EAAE,OAAO,GAAG,EAAE,WAAW,EAAE,CAAC,GAAG,KAAK,EAAE,aAAa,EAAE,IAAI,EAAE,gBAAgB,EAAE,MAAM,EAAE,CAAC;IACvK,MAAM,YAAY,GAAG,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,IAAI,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,IAAI,EAAE,IAAI,CAAC,EAAE,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAE,EAAE,CAAC,EAAE,CAAE,EAAE,MAAM,EAAE,IAAI,EAAE,CAAC,CAAC,CAAC;IACpI,MAAM,UAAU,GAAmC,EAAE,WAAW,EAAE,YAAY,EAAE,WAAW,EAAE,UAAU,EAAE,aAAa,EAAE,KAAK,EAAE,oBAAoB,EAAE,CAAC,EAAE,OAAO,EAAE,YAAY,EAAE,KAAK,EAAE,UAAU;QAC9L,SAAS,EAAE,oBAAoB,EAAE,UAAU,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,GAAG,YAAY,CAAC,CAAC,CAAE,CAAC,CAAC,GAAG,WAAW,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,YAAY,CAAC,CAAC,CAAE,CAAC,CAAC,GAAG,YAAY,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;IACvM,MAAM,KAAK,GAAG,gBAAgB,CAAC,UAAU,EAAE,OAAO,EAAE,GAAG,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,GAAG,EAAE,EAAE,WAAW,EAAE,UAAU,EAAE,KAAK,EAAE,IAAI,EAAE,SAAS,EAAE,EAAE,CAAC,CAAC;IAClI,MAAM,cAAc,GAA+C,EAAE,CAAC;IACtE,KAAK,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,IAAI,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,IAAI,EAAE,IAAI,CAAC,EAAE,CAAC,IAAI,EAAE,IAAI,CAAC,CAAU;QAAE,KAAK,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,GAAG,CAAC,CAAU;YAAE,cAAc,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,EAAE,MAAM,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC;IACpP,MAAM,MAAM,GAAG,CAAC,KAAa,EAAgB,EAAE,GAAG,MAAM,CAAC,GAAG,YAAY,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,KAAK,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,GAAG,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IAC1Q,MAAM,YAAY,GAAmC,EAAE,WAAW,EAAE,YAAY,EAAE,WAAW,EAAE,UAAU,EAAE,aAAa,EAAE,KAAK,EAAE,oBAAoB,EAAE,CAAC,EAAE,OAAO,EAAE,cAAc;QAC/K,OAAO,EAAE,EAAE,MAAM,EAAE,CAAC,EAAE,OAAO,EAAE,CAAC,MAAM,CAAC,GAAG,CAAC,EAAE,MAAM,CAAC,GAAG,CAAC,EAAE,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,KAAK,EAAE,UAAU,EAAE,SAAS,EAAE,oBAAoB,EAAE,CAAC;IACjI,MAAM,OAAO,GAAG,gBAAgB,CAAC,YAAY,EAAE,OAAO,EAAE,GAAG,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,GAAG,EAAE,EAAE,WAAW,EAAE,YAAY,EAAE,KAAK,EAAE,IAAI,EAAE,SAAS,EAAE,EAAE,CAAC,CAAC;IACxI,oIAAoI;IACpI,+FAA+F;IAC/F,MAAM,GAAG,GAAG,EAAE,IAAI,EAAE,WAAoB,EAAE,iBAAiB,EAAE,IAAI,EAAE,EAAE,IAAI,GAAG,EAAE,IAAI,EAAE,OAAgB,EAAE,CAAC;IACvG,MAAM,KAAK,GAAG,wBAAwB,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,UAAU,EAAE,GAAG,EAAE,cAAc,EAAE,KAAK,EAAE,EAAE,EAAE,UAAU,EAAE,IAAI,EAAE,cAAc,EAAE,GAAG,EAAE,EAAE,EAAE,UAAU,EAAE,GAAG,EAAE,cAAc,EAAE,GAAG,EAAE,CAAC;QACjL,KAAK,EAAE,EAAE,IAAI,EAAE,iBAAiB,EAAE,IAAI,EAAE,EAAE,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,KAAK,EAAE,CAAC,CAAC,EAAE,IAAI,EAAE,IAAI,CAAC,EAAE;QAC9G,MAAM,EAAE,CAAC,EAAE,UAAU,EAAE,IAAI,EAAE,QAAQ,EAAE,EAAE,EAAE,EAAE,CAAC,GAAG,EAAE,EAAE,EAAE,GAAG,EAAE,OAAO,EAAE,IAAI,EAAE,EAAE,SAAS,EAAE,EAAE,cAAc,EAAE,GAAG,EAAE,SAAS,EAAE,IAAI,EAAE,EAAE,CAAC,EAAE,EAAE,wBAAwB,CAAC,CAAC;IACnK,sHAAsH;IACtH,MAAM,aAAa,GAAG,CAAC,EAAU,EAAE,EAAU,EAAE,EAAU,EAAU,EAAE,CAAC,KAAK,CAAC,YAAY,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;IAC7G,MAAM,aAAa,GAAG,CAAC,CAAS,EAAE,CAAS,EAAE,EAAU,EAAU,EAAE;QACjE,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;QAC7B,OAAO,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,MAAM,CAAC,CAAC,CAAC,KAAK,CAAC,aAAa,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,GAAG,aAAa,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;IAC1I,CAAC,CAAC;IACF,MAAM,QAAQ,GAAmC,EAAE,GAAG,YAAY,EAAE,KAAK,EAAE,KAAK,CAAC,WAAW,CAAC,KAAK,EAAE,CAAC,EAAE,KAAK,EAAE,aAAa,EAAE,CAAC,EAAE,CAAC,EAAE,QAAQ,EAAE,KAAK,CAAC,EAAE,CAAC;IACtJ,MAAM,GAAG,GAAG,gBAAgB,CAAC,QAAQ,EAAE,OAAO,EAAE,GAAG,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,GAAG,EAAE,EAAE,WAAW,EAAE,QAAQ,EAAE,KAAK,EAAE,IAAI,EAAE,SAAS,EAAE,EAAE,CAAC,CAAC;IAC5H,MAAM,CAAC,cAAc,EAAE,gBAAgB,EAAE,YAAY,CAAC,GAAG,CAAC,MAAM,iBAAiB,CAAC,KAAK,CAAC,EAAE,MAAM,iBAAiB,CAAC,OAAO,CAAC,EAAE,MAAM,iBAAiB,CAAC,GAAG,CAAC,CAAC,CAAC;IAC1J,MAAM,MAAM,GAAG,6BAA6B,CAAC,OAAO,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,yBAAyB,CAAC,OAAO,EAAE,CAAC,KAAK,EAAE,IAAI,EAAE,EAAE,CAAC,qBAAqB,CAAC,OAAO,EAAE,KAAK,EAAE,IAAI,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IACrK,IAAI,CAAC;QACH,OAAO,MAAM,oBAAoB,CAAC,OAAO,EAAE,MAAM,EAAE,KAAK,EAAE,KAAK,EAAE,SAAS,EAAE,EAAE;YAC5E,IAAI,IAAsB,CAAC;YAC3B,IAAI,CAAC;gBAAC,IAAI,GAAG,CAAC,MAAM,yBAAyB,CAAC,OAAO,EAAE,SAAS,EAAE,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,WAAW,IAAI,QAAQ,CAAC;YAAC,CAAC;YAAC,OAAO,KAAK,EAAE,CAAC;gBACzH,IAAI,KAAK,YAAY,wBAAwB;oBAAE,OAAO,IAAI,CAAC;gBAC3D,MAAM,KAAK,CAAC;YACd,CAAC;YACD,MAAM,OAAO,GAAG,kBAAkB,CAAC,cAAc,EAAE,MAAM,2BAA2B,CAAC,OAAO,EAAE,CAAC,IAAI,EAAE,EAAE,EAAE,KAAK,EAAE,KAAK,EAAE,WAAW,EAAE,UAAU,EAAE,QAAQ,EAAE,MAAM,CAAC,KAAK,CAAC,KAAK,CAAC,GAAI,EAAE,IAAI,CAAC,EAAE,EAAE,KAAK,EAAE,SAAS,EAAE,SAAS,CAAC,CAAC,CAAC;YAC1N,IAAI,OAAO,KAAK,IAAI;gBAAE,OAAO,SAAS,GAAG,OAAO,CAAC;YACjD,MAAM,GAAG,GAAG,MAAM,6BAA6B,CAAC,OAAO,EAAE,CAAC,IAAI,EAAE,EAAE,EAAE,KAAK,EAAE,OAAO,EAAE,WAAW,EAAE,YAAY,EAAE,QAAQ,EAAE,EAAE,EAAE,EAAE,aAAa,EAAE,IAAI,EAAE,EAAE,EAAE,QAAQ,EAAE,EAAE,EAAE;gBACpK,KAAK,EAAE,uBAAuB,CAAC,OAAO,EAAE,OAAO,CAAC,GAAI,EAAE,IAAI,CAAC,EAAE,EAAE,KAAK,EAAE,SAAS,EAAE,SAAS,CAAC,CAAC;YAC9F,MAAM,cAAc,GAAG,kBAAkB,CAAC,gBAAgB,EAAE,GAAG,CAAC,CAAC;YACjE,IAAI,cAAc,KAAK,IAAI;gBAAE,OAAO,WAAW,GAAG,cAAc,CAAC;YACjE,MAAM,MAAM,GAAG,MAAM,6BAA6B,CAAC,OAAO,EAAE,CAAC,IAAI,EAAE,EAAE,EAAE,KAAK,EAAE,GAAG,EAAE,WAAW,EAAE,QAAQ,EAAE,QAAQ,EAAE,EAAE,EAAE,EAAE,aAAa,EAAE,IAAI,EAAE,EAAE,EAAE,QAAQ,EAAE,EAAE,EAAE;gBAC/J,KAAK,EAAE,uBAAuB,CAAC,OAAO,EAAE,GAAG,CAAC,GAAI,EAAE,IAAI,CAAC,EAAE,EAAE,KAAK,EAAE,SAAS,EAAE,SAAS,CAAC,CAAC;YAC1F,MAAM,UAAU,GAAG,kBAAkB,CAAC,YAAY,EAAE,MAAM,CAAC,CAAC;YAC5D,OAAO,UAAU,KAAK,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,6BAA6B,GAAG,UAAU,CAAC;QACjF,CAAC,CAAC,CAAC;IACL,CAAC;YAAS,CAAC;QAAC,MAAM,CAAC,OAAO,EAAE,CAAC;IAAC,CAAC;AACjC,CAAC","sourcesContent":["// SPDX-License-Identifier: Apache-2.0\n\n/**\n * The image-plane module of the GPU display view (#581): the Speed kernel's `scene_value` (the seam of #478 and #491,\n * `DisplayViewGpuMaterialModule`) for an image-plane print, so that a photo renders on the device like a chart.\n *\n * **What the reference computes** at a query point (expected-signal.ts, the direct path's `signalAt`): the print-local point\n * cu = u·m/h, cv = v·m/h (outside |cu| ≤ 1, |cv| ≤ aspect the surround: its reflectance times the neutral signal N_c); the\n * pixel coordinates X = (cu + 1)·W/2, Y = (aspect − cv)·H/(2·aspect); then the bilinear weights of the four nearest pixel\n * centers (clamped to the edge) times each pixel's channel signal S_j,c = Σ_n A_c(λ_n)·ρ_j(λ_n)·W_n, integrated from the pixel's\n * reconstructed spectrum ρ_j (its cell's eight node spectra blended, then the paper; `ResolvedImagePlaneScene.pixelSignals`), or\n * the nearest pixel's. No per-channel values of the lattice nodes are formed.\n *\n * **On the device.** The host computes every pixel's S_j,c in float64 and rounds it once to float32 (3 values per pixel). When\n * they fit one storage binding they are uploaded once per scene and kept on the device while the renderer lives; otherwise each\n * band of a render builds and uploads only its rows. Either way they are built and written 4 MiB of rows at a time, so the host\n * never holds more than that of S (owner decision 2026-11-09; the bytes are those of the whole array). The module recovers the exact lattice point from the\n * kernel's float32 point as the composite module does (`LATTICE_POINT_WGSL`, about 2⁻⁴⁶ relative), forms cu, cv, X and Y\n * in double-float32 and decides the print's edge in the composite kernel's bound arithmetic (a point within its bound of the\n * edge is flagged: the display pixel is recomputed on the host). Per axis it splits s = X − ½ into its integer part n and\n * fraction f exactly (double-float32), rounds the smaller of f and 1 − f once to float32 and forms the other as 1 minus it, so\n * each weight carries at most 2E relative error (E = 2⁻²³) besides the point's own error. The point's error moves weight\n * between neighbors: with δ the bound of |X_gpu − X_ref| in pixels (4·(2⁻⁴⁰·(|X| + 1)) + 4·10⁻¹⁴·W/2·k: the recovery, the\n * double-float32 steps and the reference's own float64 rounding of its query point, generously), the value moves by at most\n * (δx + δy)·2·max S, and the point is flagged when that exceeds E times the value. So an unflagged value carries a relative\n * error only: S's rounding 1, the weights 2 + 2, the two products and two sums 4, the flag's budget 1, that is 10 roundings\n * beyond the table lookup's one (`extraSiteRoundings`), added to the Speed kernel's chains (docs/WEBGPU.md, \"Display view (Speed)\").\n * `nearest` takes the pixel under the point, flagging a point within its bound of a pixel edge (`bfloor`). All S are\n * non-negative (non-negative responses, illuminant and spectra), so the weighted sums need no cancellation term.\n *\n * The surround's three channel signals, integrated from its spectrum (its reflectance times the white radiance at each node) like\n * a pixel's, are the kernel's binding 3 (`spectral_table`, which this module does not otherwise use).\n */\n\nimport { ImageKernelExecutionError, type ChartCompositeProbeRecord } from \"./image-kernel-execution.js\";\nimport { BUFFER_USAGE, GpuResourceScope, type WebGpuRuntime } from \"./webgpu-runtime.js\";\nimport { dwSplit } from \"./chart-composite-gpu.js\";\nimport type { ChartCompositeProductPath } from \"./wgsl-chart-composite-kernel.js\";\nimport { displayViewModuleCommonWgsl, displayViewLatticeConstants, LATTICE_BINDINGS_WGSL, LATTICE_POINT_WGSL } from \"./display-view-composite-material-module.js\";\nimport type { DisplayViewGpuMaterialModule } from \"./display-view-speed-gpu.js\";\nimport { renderDisplayView, type DisplayViewGpuSetup } from \"../facade/display-view.js\";\nimport type { DisplayViewExpectedSignalInput } from \"../facade/expected-signal.js\";\nimport { resolveCaptureIlluminant, type ImagePlaneSignalScene } from \"../facade/scene-radiance.js\";\nimport { resolveMixedIllumination } from \"../facade/mixed-illumination.js\";\nimport { imagePlaneSupport } from \"../facade/image-plane-support.js\";\nimport { CONFORMANCE_FRAME, conformanceChannelAt, conformanceResponses, conformanceSetup, conformanceVerdict, renderDisplayViewSpeedOnGpu, withConformanceAbort } from \"./display-view-speed-gpu.js\";\nimport { renderDisplayViewQualityOnGpu, type DisplayViewQualityImagePasses } from \"./display-view-quality-gpu.js\";\nimport { imagePlaneQualityPasses } from \"./display-view-image-quality.js\";\nimport { ChartCompositeProbeError, probeChartCompositeDevice } from \"./chart-composite-gpu.js\";\nimport { createComputePipeline } from \"./webgpu-runtime.js\";\n\n/** The module's id per product path. */\nexport function imagePlaneMaterialModuleId(path: ChartCompositeProductPath): string { return path === \"dekker\" ? \"image-plane\" : \"image-plane-fma\"; }\n/** Roundings of the module's value beyond the table lookup's one (see the module comment). */\nexport const IMAGE_PLANE_EXTRA_SITE_ROUNDINGS = 10;\n/** Layout of the photo uniform (float32 / u32 words, 64 bytes). */\nexport const IMAGE_PLANE_PARAMS = /* @__PURE__ */ Object.freeze({ width: 0, height: 1, nearest: 2, row0: 3, halfWHi: 4, halfWLo: 5, halfHHi: 6, halfHLo: 7, aspectHi: 8, aspectLo: 9,\n  slackX: 10, slackY: 11, words: 16 });\n/** Table values above this are outside the float32 domain. */\nconst FLOAT32_SAFE_MAGNITUDE = 2 ** 120;\n/** Host bytes of pixel signals built per write (whole uploads and bands). */\nconst UPLOAD_CHUNK_BYTES = 4 << 20;\n/** Scenes whose pixel signals stay uploaded per renderer. */\nconst KEPT_SCENES = 2;\n\n/** The module's WGSL with one product path (see the module comment). */\nexport function imagePlaneMaterialWgsl(path: ChartCompositeProductPath = \"dekker\"): string {\n  return `${displayViewModuleCommonWgsl(path)}\n@group(1) @binding(0) var<storage, read> photo_signal: array<f32>;\nstruct PhotoParams {\n  width: u32, height: u32, nearest: u32, row0: u32,\n  half_w_hi: f32, half_w_lo: f32, half_h_hi: f32, half_h_lo: f32,\n  aspect_hi: f32, aspect_lo: f32, slack_x: f32, slack_y: f32,\n  pad1: f32, pad2: f32, pad3: f32, pad4: f32,\n};\n@group(1) @binding(1) var<uniform> photo: PhotoParams;\n${LATTICE_BINDINGS_WGSL}\n${LATTICE_POINT_WGSL}\n// The point's error bound in pixels along one axis (see the module comment).\nfn pixel_slack(x: D, slack: f32) -> f32 { return 4.0 * (DW * (abs(x.h) + 1.0)) + 4.0 * slack; }\n// One axis of the bilinear weights: the pixel pair (clamped to the edge) and the weights of s = X − 1/2.\nstruct Axis { i0: u32, i1: u32, w0: f32, w1: f32 }\nfn bilinear_axis(s: D, n: u32) -> Axis {\n  var f0 = floor(s.h);\n  if (f0 == s.h && s.l < 0.0) { f0 = f0 - 1.0; }\n  let frac = d_sub(s, D(f0, 0.0));\n  let comp = d_sub(D(ob(f0 + 1.0), 0.0), s);\n  var w0 = 0.0; var w1 = 0.0;\n  if (frac.h <= 0.5) { w1 = ob(frac.h + frac.l); w0 = ob(1.0 - w1); } else { w0 = ob(comp.h + comp.l); w1 = ob(1.0 - w0); }\n  let last = i32(n) - 1;\n  return Axis(u32(clamp(i32(f0), 0, last)), u32(clamp(i32(f0) + 1, 0, last)), w0, w1);\n}\nfn photo_at(i: u32, j: u32, c: u32) -> f32 { return photo_signal[((j - photo.row0) * photo.width + i) * 3u + c]; }\nfn scene_value(p: vec2<f32>, c: u32) -> f32 {\n  let k = D(lattice.chart_k_hi, lattice.chart_k_lo);\n  let u = d_mul(lattice_point(p.x, D(lattice.scale_u_hi, lattice.scale_u_lo), 0u), k);\n  let v = d_mul(lattice_point(p.y, D(lattice.scale_v_hi, lattice.scale_v_lo), 1u), k);\n  let aspect = D(photo.aspect_hi, photo.aspect_lo);\n  if (!(le(babs(bd(u)), bx(1.0)) && le(babs(bd(v)), bd(aspect)))) { return spectral_table[c]; }\n  let x = d_mul(d_add(u, D(1.0, 0.0)), D(photo.half_w_hi, photo.half_w_lo));\n  let y = d_mul(d_sub(aspect, v), D(photo.half_h_hi, photo.half_h_lo));\n  if (photo.nearest != 0u) {\n    let last_x = f32(photo.width - 1u);\n    let last_y = f32(photo.height - 1u);\n    let i = u32(clamp(bfloor(bd(x), 1.0, last_x), 0.0, last_x));\n    let j = u32(clamp(bfloor(bd(y), 1.0, last_y), 0.0, last_y));\n    return photo_at(i, j, c);\n  }\n  let ax = bilinear_axis(d_sub(x, D(0.5, 0.0)), photo.width);\n  let ay = bilinear_axis(d_sub(y, D(0.5, 0.0)), photo.height);\n  let s00 = photo_at(ax.i0, ay.i0, c);\n  let s10 = photo_at(ax.i1, ay.i0, c);\n  let s01 = photo_at(ax.i0, ay.i1, c);\n  let s11 = photo_at(ax.i1, ay.i1, c);\n  let value = ob(ob(ay.w0 * ob(ob(ax.w0 * s00) + ob(ax.w1 * s10))) + ob(ay.w1 * ob(ob(ax.w0 * s01) + ob(ax.w1 * s11))));\n  let largest = max(max(s00, s10), max(s01, s11));\n  let moved = (pixel_slack(x, photo.slack_x) + pixel_slack(y, photo.slack_y)) * 2.0 * largest;\n  if (!(moved <= E * value)) { material_uncertain = true; }\n  return value;\n}\n`;\n}\n\n/** The image-plane print behind a GPU setup's expected signal (`ImagePlaneSignalScene.image`), or undefined. */\nexport function imagePlaneOf(gpu: DisplayViewGpuSetup): ImagePlaneSignalScene[\"image\"] | undefined {\n  return (gpu.signalInput.scene as { image?: ImagePlaneSignalScene[\"image\"] }).image;\n}\n/** The pixel signals of a scene on the device (per renderer), keyed by the signal scene's image record. */\ninterface UploadedImage { key: object; scope: GpuResourceScope; buffer: GPUBuffer; surround: Float32Array<ArrayBuffer> }\n\n/** The image-plane modules of one display-view renderer: `moduleFor` builds the setup's module (its uploads kept), `release` frees them. */\nexport interface DisplayViewImageModules {\n  /** The setup's module; the first call on a runtime runs the module's self-check (`runImagePlaneConformance`) and refuses the device if it fails. */\n  moduleFor(gpu: DisplayViewGpuSetup, signal: AbortSignal): Promise<DisplayViewGpuMaterialModule>;\n  /** The setup's module without the self-check (the self-check itself uses it). */\n  build(gpu: DisplayViewGpuSetup, path: ChartCompositeProductPath): DisplayViewGpuMaterialModule;\n  /** #581: the print's Quality passes 1 and 2 (src/gpu/display-view-image-quality.ts), after the same self-check. */\n  qualityPasses(gpu: DisplayViewGpuSetup, signal: AbortSignal): Promise<DisplayViewQualityImagePasses>;\n  release(): void;\n}\n\n/**\n * The image-plane modules of a renderer; `probe` is the composite module's device probe (single flight), whose double-word product\n * path the module's double-float32 takes.\n */\nexport function createDisplayViewImageModules(runtime: WebGpuRuntime, probe: (signal: AbortSignal) => Promise<ChartCompositeProbeRecord>): DisplayViewImageModules {\n  const uploads: UploadedImage[] = [], wgsl = new Map<ChartCompositeProductPath, string>();\n  // Pixel rows [j0, j1)'s channel signals, built on demand so that the host holds one chunk or band (#581): each pixel's signals\n  // integrated from its reconstructed spectrum (`pixelSignals`, the reference's own float64 values), rounded once.\n  const pixel = new Float64Array(3);\n  const rowsOf = (image: ImagePlaneSignalScene[\"image\"], j0: number, j1: number): Float32Array<ArrayBuffer> => {\n    const { scene, nodes } = image, values = new Float32Array((j1 - j0) * scene.width * 3);\n    for (let j = j0 * scene.width, o = 0; j < j1 * scene.width; j++) {\n      scene.pixelSignals(j, nodes, pixel, 0);\n      for (let c = 0; c < 3; c++) {\n        const value = pixel[c]!;\n        if (!(value >= 0 && value <= FLOAT32_SAFE_MAGNITUDE)) throw new ImageKernelExecutionError(\"unsupported-float32-domain\", \"A pixel's channel signal is outside the float32 domain.\");\n        values[o++] = value;\n      }\n    }\n    return values;\n  };\n  // Rows [a, b)'s signals written to `buffer` (whose first row is `base`), built UPLOAD_CHUNK_BYTES of rows at a time.\n  const fill = (buffer: GPUBuffer, image: ImagePlaneSignalScene[\"image\"], a: number, b: number, base: number): void => {\n    const rowBytes = image.scene.width * 12, chunk = Math.max(1, Math.floor(UPLOAD_CHUNK_BYTES / rowBytes));\n    for (let j = a; j < b; j += chunk) { const rows = rowsOf(image, j, Math.min(b, j + chunk)); runtime.device.queue.writeBuffer(buffer, (j - base) * rowBytes, rows.buffer, 0, rows.byteLength); }\n  };\n  const rowsBuffer = (scope: GpuResourceScope, label: string, image: ImagePlaneSignalScene[\"image\"], j0: number, j1: number): GPUBuffer => {\n    const buffer = scope.buffer(label, (j1 - j0) * image.scene.width * 12, BUFFER_USAGE.STORAGE | BUFFER_USAGE.COPY_DST);\n    fill(buffer, image, j0, j1, j0);\n    return buffer;\n  };\n  // The latest band's rows stay on the device until the next band: the rows the two share are copied there on the device, and\n  // only the others are built on the host (the same bytes; the host still holds one chunk at most).\n  let lastBand: { image: object; j0: number; j1: number; scope: GpuResourceScope; buffer: GPUBuffer } | null = null;\n  const bandBuffer = (image: ImagePlaneSignalScene[\"image\"], j0: number, j1: number): GPUBuffer => {\n    const rowBytes = image.scene.width * 12, scope = new GpuResourceScope(runtime), last = lastBand?.image === image ? lastBand : null;\n    try {\n      const buffer = scope.buffer(\"photivra-display-view-image-signal-band\", (j1 - j0) * rowBytes, BUFFER_USAGE.STORAGE | BUFFER_USAGE.COPY_DST | BUFFER_USAGE.COPY_SRC);\n      const o0 = last ? Math.max(j0, last.j0) : j0, o1 = last ? Math.min(j1, last.j1) : j0;\n      if (o1 > o0) {\n        const encoder = runtime.device.createCommandEncoder({ label: \"photivra-display-view-image-band-rows\" });\n        encoder.copyBufferToBuffer(last!.buffer, (o0 - last!.j0) * rowBytes, buffer, (o0 - j0) * rowBytes, (o1 - o0) * rowBytes);\n        runtime.device.queue.submit([encoder.finish()]);\n        fill(buffer, image, j0, o0, j0); fill(buffer, image, o1, j1, j0);\n      } else fill(buffer, image, j0, j1, j0);\n      lastBand?.scope.release();\n      lastBand = { image, j0, j1, scope, buffer };\n      return buffer;\n    } catch (error) { scope.release(); throw error; }\n  };\n  const upload = (image: ImagePlaneSignalScene[\"image\"]): UploadedImage => {\n    const key = image as object, index = uploads.findIndex((e) => e.key === key);\n    if (index >= 0) { const [entry] = uploads.splice(index, 1); uploads.push(entry!); return entry!; }\n    const scope = new GpuResourceScope(runtime);\n    let buffer: GPUBuffer;\n    try { buffer = rowsBuffer(scope, \"photivra-display-view-image-signal\", image, 0, image.scene.height); } catch (error) { scope.release(); throw error; }\n    const entry: UploadedImage = { key, scope, buffer, surround: Float32Array.from(image.surround) };\n    uploads.push(entry);\n    while (uploads.length > KEPT_SCENES) uploads.shift()!.scope.release();\n    return entry;\n  };\n  let checked: Promise<true | string> | null = null;\n  const build = (gpu: DisplayViewGpuSetup, path: ChartCompositeProductPath): DisplayViewGpuMaterialModule => {\n      const image = imagePlaneOf(gpu);\n      if (!image || !gpu.placement) throw new ImageKernelExecutionError(\"unsupported-on-backend\", \"The setup carries no image-plane print and placement; use cpu-reference.\");\n      const scene = image.scene, lattice = displayViewLatticeConstants(gpu.signalInput, gpu.placement, 0), rowBytes = scene.width * 12, banded = scene.height * rowBytes > runtime.maximumStorageBindingBytes;\n      const params = new ArrayBuffer(4 * IMAGE_PLANE_PARAMS.words), words = new Uint32Array(params), floats = new Float32Array(params), P = IMAGE_PLANE_PARAMS;\n      words[P.width] = scene.width; words[P.height] = scene.height; words[P.nearest] = scene.interpolation === \"nearest\" ? 1 : 0;\n      const halfW = scene.width / 2, halfH = scene.height / (2 * scene.aspect);\n      [floats[P.halfWHi], floats[P.halfWLo]] = dwSplit(halfW); [floats[P.halfHHi], floats[P.halfHLo]] = dwSplit(halfH); [floats[P.aspectHi], floats[P.aspectLo]] = dwSplit(scene.aspect);\n      // The reference's float64 query point carries up to 10⁻¹⁴ frame units (the composite kernel's slack) times the pixels per frame unit.\n      floats[P.slackX] = 1e-14 * gpu.placement.chartUnitsPerFrameUnit * halfW; floats[P.slackY] = 1e-14 * gpu.placement.chartUnitsPerFrameUnit * halfH;\n      let code = wgsl.get(path);\n      if (code === undefined) { code = imagePlaneMaterialWgsl(path); wgsl.set(path, code); }\n      const base = { id: imagePlaneMaterialModuleId(path), wgsl: code, extraSiteRoundings: IMAGE_PLANE_EXTRA_SITE_ROUNDINGS, siteTable: Float32Array.from(image.surround) };\n      if (!banded) {\n        const uploaded = upload(image);\n        return { ...base, bindings: [{ binding: 0, kind: \"storage\", buffer: uploaded.buffer }, { binding: 1, kind: \"uniform\", data: new Uint8Array(params) },\n          { binding: 2, kind: \"uniform\", data: new Uint8Array(lattice.params) }, { binding: 3, kind: \"storage\", data: lattice.offsets }] };\n      }\n      // A photo larger than one binding: per band the pixel rows its native site rows reach. The print-local v of site row y and\n      // query offset o is (y + c_o)·scale_v·k; the photo row is (aspect − v)·H/(2·aspect), plus two rows of margin either side.\n      const offsetsV = Array.from({ length: lattice.offsets.length / 4 }, (_, o) => lattice.offsets[4 * o + 2]! + lattice.offsets[4 * o + 3]!);\n      const rowOf = (y: number, c: number): number => (scene.aspect - (y + c) * lattice.scaleV * gpu.placement!.chartUnitsPerFrameUnit) * halfH;\n      const rate = Math.abs(lattice.scaleV) * gpu.placement.chartUnitsPerFrameUnit * halfH, maximumRows = Math.floor(runtime.maximumStorageBindingBytes / rowBytes);\n      return { ...base, bindings: [],\n        bandRows: (sitesPerDisplayRow: number): number => Math.max(1, Math.floor(((maximumRows - 6) / rate - 6) / sitesPerDisplayRow)),\n        band(scope, layout, y0, y1): GPUBindGroup {\n          let low = Infinity, high = -Infinity;\n          for (const y of [y0, y1]) for (const c of offsetsV) { const r = rowOf(y, c); low = Math.min(low, r); high = Math.max(high, r); }\n          const j0 = Math.max(0, Math.min(scene.height - 1, Math.floor(low) - 2)), j1 = Math.max(j0, Math.min(scene.height - 1, Math.floor(high) + 2));\n          if ((j1 - j0 + 1) * rowBytes > runtime.maximumStorageBindingBytes) throw new ImageKernelExecutionError(\"insufficient-limits\", \"A band's pixel signals exceed the device's storage binding limit.\");\n          const bandParams = params.slice(0);\n          new Uint32Array(bandParams)[P.row0] = j0;\n          return runtime.device.createBindGroup({ layout, entries: [\n            { binding: 0, resource: { buffer: bandBuffer(image, j0, j1 + 1) } },\n            { binding: 1, resource: { buffer: scope.upload(\"photivra-display-view-image-params\", new Uint8Array(bandParams), BUFFER_USAGE.UNIFORM) } },\n            { binding: 2, resource: { buffer: scope.upload(\"photivra-display-view-image-lattice\", new Uint8Array(lattice.params), BUFFER_USAGE.UNIFORM) } },\n            { binding: 3, resource: { buffer: scope.upload(\"photivra-display-view-image-offsets\", lattice.offsets, BUFFER_USAGE.STORAGE) } }] });\n        } };\n  };\n  // The module's self-check, once per runtime (single flight; a failure is retried on the next photo).\n  const check = async (): Promise<void> => {\n    checked ??= runImagePlaneConformance(runtime, undefined).catch((error: unknown) => { checked = null; throw error; });\n    const verdict = await checked;\n    if (verdict !== true) { checked = null; throw new ImageKernelExecutionError(\"unsupported-on-backend\", \"The image-plane module failed its self-check on this device: \" + verdict); }\n  };\n  return {\n    async moduleFor(gpu, signal): Promise<DisplayViewGpuMaterialModule> { const path = (await probe(signal)).productPath ?? \"dekker\"; await check(); return build(gpu, path); },\n    async qualityPasses(gpu, signal): Promise<DisplayViewQualityImagePasses> { const path = (await probe(signal)).productPath ?? \"dekker\"; await check(); return imagePlaneQualityPasses(runtime, gpu, path); },\n    build,\n    release(): void { for (const entry of uploads.splice(0)) entry.scope.release(); lastBand?.scope.release(); lastBand = null; }\n  };\n}\n\n/**\n * #581's run-time self-check of the module, run on a device's first photo: a 12 × 8 photo (a color ramp with hard edges, every\n * lattice cell's corner colors among its pixels) printed at 0.8 of the frame's width on the 18 % surround of a 38 × 26 frame shown\n * at 13 × 9, the GPU against `renderDisplayView` within the a-priori bounds, every flagged pixel recomputed: the Speed kernel with\n * the 2 × 2 box, then the Quality passes with the four-spot filter's 16 offsets and synthetic radius-1 kernels at three nodes, and\n * (#597) the same passes under mixed light (a three-source gradient with a hard edge and a point source).\n * A device that fails the double-float32 probe refuses photos (unsupported-on-backend, as composite charts), so the check passes\n * there without keeping the device's charts off the GPU.\n */\nexport async function runImagePlaneConformance(handle: unknown, signal: AbortSignal | undefined): Promise<true | string> {\n  const runtime = handle as WebGpuRuntime, support = imagePlaneSupport();\n  if (!support) return true;\n  const { nativeWidth, nativeHeight, frameWidth, nodes } = CONFORMANCE_FRAME, g = 2 / frameWidth;\n  const width = 12, height = 8, data = new Uint8Array(width * height * 3);\n  for (let y = 0; y < height; y++) for (let x = 0; x < width; x++) for (let c = 0; c < 3; c++) data[(y * width + x) * 3 + c] = [(x * 23) % 256, (y * 37 + x * 5) % 256, x < 6 ? 255 - y * 30 : 10][c]!;\n  const printWidth = 1.6, aspect = height / width, k = 2 / printWidth, surround = 0.18, placement = { chartUnitsPerFrameUnit: k, aspect };\n  const scene = await support.resolveImagePlaneScene({ kind: \"image-plane\", sceneId: \"conformance\", image: { width, height, channels: 3, encoding: \"srgb-8\", data } }, aspect);\n  const white = Float64Array.of(0.5, 0.9, 0.7), A = conformanceResponses(1);\n  // The print's radiance as the provider forms it: the print-local point, the surround outside, times the white radiance at the node.\n  const radianceAt = (u: number, v: number, nm: number): number => {\n    const cu = u * k, cv = v * k;\n    return (Math.abs(cu) <= 1 && Math.abs(cv) <= aspect ? scene.reflectanceAt(cu, cv, nm) : surround) * white[nodes.indexOf(nm)]!;\n  };\n  const printScene = scene.signalScene(nodes, A, white, radianceAt, k, 1, surround), readout = { perElectron: 1 / 60000, darkElectrons: 0.07, constantVariance: 1.5e-9 };\n  const speedOffsets = [[-0.25, -0.25], [0.25, -0.25], [-0.25, 0.25], [0.25, 0.25]].map(([x, y]) => ({ x: x!, y: y!, weight: 0.25 }));\n  const speedInput: DisplayViewExpectedSignalInput = { nativeWidth, nativeHeight, outputWidth: frameWidth, wavelengthsNm: nodes, electronsPerRadiance: A, offsets: speedOffsets, scene: printScene,\n    channelAt: conformanceChannelAt, queryPoint: (x, y, o, out) => { out[0] = (x + 0.5 + speedOffsets[o]!.x - nativeWidth / 2) * g; out[1] = -(y + 0.5 + speedOffsets[o]!.y - nativeHeight / 2) * g; } };\n  const speed = conformanceSetup(speedInput, readout, 581, 2, signal, { gpu: { signalInput: speedInput, chart: null, placement } });\n  const qualityOffsets: { x: number; y: number; weight: number }[] = [];\n  for (const [nx, ny] of [[-0.25, -0.25], [0.25, -0.25], [-0.25, 0.25], [0.25, 0.25]] as const) for (const [ax, ay] of [[-0.5, -0.5], [0.5, -0.5], [-0.5, 0.5], [0.5, 0.5]] as const) qualityOffsets.push({ x: nx - ax, y: ny - ay, weight: 1 / 16 });\n  const kernel = (sigma: number): Float64Array => { const w = Float64Array.from({ length: 9 }, (_, t) => Math.exp(-0.5 * (((t % 3) - 1) ** 2 + (Math.floor(t / 3) - 1) ** 2) / (sigma * sigma))); const sum = w.reduce((a, b) => a + b, 0); return w.map((v) => v / sum); };\n  const qualityInput: DisplayViewExpectedSignalInput = { nativeWidth, nativeHeight, outputWidth: frameWidth, wavelengthsNm: nodes, electronsPerRadiance: A, offsets: qualityOffsets,\n    kernels: { radius: 1, perNode: [kernel(0.6), kernel(0.7), kernel(0.8)] }, scene: printScene, channelAt: conformanceChannelAt };\n  const quality = conformanceSetup(qualityInput, readout, 582, 5, signal, { gpu: { signalInput: qualityInput, chart: null, placement } });\n  // #597: the same Quality case under mixed light: a gradient of three sources with a hard edge (equal stops) across the print, and a\n  // point source over its upper left, the reference's irradiance from the same parts in float64.\n  const day = { kind: \"planckian\" as const, temperatureKelvin: 6500 }, lamp = { kind: \"cie-a\" as const };\n  const light = resolveMixedIllumination({ sources: [{ illuminant: day, illuminanceLux: 1.2e5 }, { illuminant: lamp, illuminanceLux: 8e4 }, { illuminant: day, illuminanceLux: 3e4 }],\n    blend: { kind: \"linear-gradient\", from: { u: -0.9, v: 0.1 }, to: { u: 0.8, v: -0.2 }, stops: [0, 0.45, 0.45] },\n    points: [{ illuminant: lamp, position: { xM: -0.3, yM: 0.2, heightM: 0.35 }, intensity: { illuminanceLux: 6e4, distanceM: 0.35 } }] }, resolveCaptureIlluminant);\n  // The print's half width is 1 m here (print-local units are metres), so the point sources read the print-local point.\n  const litIrradiance = (cu: number, cv: number, nm: number): number => light.irradianceAt(cu, cv, cu, cv, nm);\n  const litRadianceAt = (u: number, v: number, nm: number): number => {\n    const cu = u * k, cv = v * k;\n    return (Math.abs(cu) <= 1 && Math.abs(cv) <= aspect ? scene.reflectanceAt(cu, cv, nm) : surround) * litIrradiance(cu, cv, nm) / Math.PI;\n  };\n  const litInput: DisplayViewExpectedSignalInput = { ...qualityInput, scene: scene.signalScene(nodes, A, white, litRadianceAt, k, 1, surround, light) };\n  const lit = conformanceSetup(litInput, readout, 597, 5, signal, { gpu: { signalInput: litInput, chart: null, placement } });\n  const [speedReference, qualityReference, litReference] = [await renderDisplayView(speed), await renderDisplayView(quality), await renderDisplayView(lit)];\n  const images = createDisplayViewImageModules(runtime, (s) => probeChartCompositeDevice(runtime, (label, code) => createComputePipeline(runtime, label, code, s), s));\n  try {\n    return await withConformanceAbort(runtime, signal, async (abort, pipelines) => {\n      let path: \"dekker\" | \"fma\";\n      try { path = (await probeChartCompositeDevice(runtime, pipelines, abort.signal)).productPath ?? \"dekker\"; } catch (error) {\n        if (error instanceof ChartCompositeProbeError) return true;\n        throw error;\n      }\n      const verdict = conformanceVerdict(speedReference, await renderDisplayViewSpeedOnGpu(runtime, 1 << 20, { setup: speed, signalInput: speedInput, material: images.build(speed.gpu!, path) }, abort, undefined, pipelines));\n      if (verdict !== true) return \"Speed: \" + verdict;\n      const gpu = await renderDisplayViewQualityOnGpu(runtime, 1 << 20, { setup: quality, signalInput: qualityInput, material: { id: \"image-plane\", wgsl: \"\", bindings: [] },\n        image: imagePlaneQualityPasses(runtime, quality.gpu!, path) }, abort, undefined, pipelines);\n      const qualityVerdict = conformanceVerdict(qualityReference, gpu);\n      if (qualityVerdict !== true) return \"Quality: \" + qualityVerdict;\n      const litGpu = await renderDisplayViewQualityOnGpu(runtime, 1 << 20, { setup: lit, signalInput: litInput, material: { id: \"image-plane\", wgsl: \"\", bindings: [] },\n        image: imagePlaneQualityPasses(runtime, lit.gpu!, path) }, abort, undefined, pipelines);\n      const litVerdict = conformanceVerdict(litReference, litGpu);\n      return litVerdict === true ? true : \"Quality under mixed light: \" + litVerdict;\n    });\n  } finally { images.release(); }\n}\n"]}

package/dist/gpu/display-view-image-quality.d.ts

@@ -0,0 +1,87 @@
+import { type WebGpuRuntime } from "./webgpu-runtime.js";
+import type { ChartCompositeProductPath } from "./wgsl-chart-composite-kernel.js";
+import { type DisplayViewQualityImagePasses, type SignalGather } from "./display-view-quality-gpu.js";
+import { type DisplayViewExpectedSignalInput } from "../facade/expected-signal.js";
+import type { DisplayViewGpuSetup } from "../facade/display-view.js";
+import { type ImagePlaneSignalScene } from "../facade/scene-radiance.js";
+/** Relative roundings of a lattice point's radiance (see the module comment). */
+export declare const IMAGE_PLANE_LATTICE_U = 32;
+/**
+ * #597: relative roundings of the mixed light at a lattice point (a point whose computed bound exceeds IMAGE_PLANE_LIGHT_U − 1 is
+ * flagged; see the module comment), added to IMAGE_PLANE_LATTICE_U for a lit print.
+ */
+export declare const IMAGE_PLANE_LIGHT_U = 40;
+/** #597: layout of the light table (float32 / u32 words; see `imagePlaneLightTable`). */
+export declare const IMAGE_PLANE_LIGHT: Readonly<{
+    kind: 0;
+    sources: 1;
+    points: 2;
+    regions: 3;
+    outside: 4;
+    regionsAt: 6;
+    pointsAt: 7;
+    tablesAt: 8;
+    slack: 9;
+    surround: 10;
+    maximumInverse: 11;
+    halfHi: 12;
+    halfLo: 13;
+    aU: 16;
+    aV: 18;
+    dU: 20;
+    dV: 22;
+    inverse: 24;
+    stops: 26;
+    regionWords: 9;
+    pointWords: 6;
+    maximumPoints: 8;
+}>;
+/** Layout of the photo uniform (u32 / float32 words, 96 bytes). */
+export declare const IMAGE_PLANE_QUALITY_PARAMS: Readonly<{
+    width: 0;
+    height: 1;
+    nearest: 2;
+    nodes: 3;
+    row0: 4;
+    zero: 5;
+    lit: 6;
+    suHi: 8;
+    suLo: 9;
+    svHi: 10;
+    svLo: 11;
+    kHi: 12;
+    kLo: 13;
+    halfWHi: 14;
+    halfWLo: 15;
+    halfHHi: 16;
+    halfHLo: 17;
+    aspectHi: 18;
+    aspectLo: 19;
+    slackX: 20;
+    slackY: 21;
+    black: 22;
+    span: 23;
+    words: 24;
+}>;
+/**
+ * #597: the light table of a lit print at the spectral nodes (the layout `IMAGE_PLANE_LIGHT`): the blend (kind, sources, its geometry
+ * and stops with each segment's inverse width, 0 at a hard edge), the regions, the point sources (x, y, h in metres), then per node
+ * T_i,n = k_i·S_i(λ_n)/π for each blend source and c_j·S_j(λ_n)/π for each point source: spectra only, no channel value. `slack` is the
+ * reference's own rounding of a query point in print-local units.
+ */
+export declare function imagePlaneLightTable(light: NonNullable<ImagePlaneSignalScene["image"]["light"]>, wavelengths: readonly number[], surroundReflectance: number, slack: number): {
+    table: Float32Array<ArrayBuffer>;
+    largestIrradiance: number;
+};
+/** Pass 1's WGSL (see the module comment). */
+export declare function imagePlaneQualityLatticeWgsl(path: ChartCompositeProductPath): string;
+/** Pass 2's WGSL (see the module comment): the blocked sum over the nodes, then the Quality site tail. */
+export declare const IMAGE_PLANE_QUALITY_SITES_WGSL: string;
+/**
+ * The print's gather: per group its dense tap rectangle (with site kernels the reference's combined window; without them the
+ * members' span), the lattice extent beyond the site window and the point constants, and the table C_n,c(t) at (n·3 + c)·T + t,
+ * the reference's float64 values rounded once.
+ */
+export declare function imagePlaneQualityGather(signalInput: DisplayViewExpectedSignalInput, largestLitRadiance?: number): SignalGather;
+/** The print's Quality passes for a setup (see `DisplayViewQualityImagePasses`), with the double-word product path the device's probe chose. */
+export declare function imagePlaneQualityPasses(runtime: WebGpuRuntime, gpu: DisplayViewGpuSetup, path: ChartCompositeProductPath): DisplayViewQualityImagePasses;