@photivra/engine 2.21.0 → 2.23.0

package/dist/color/neighbor-highlight-reconstruction.js

@@ -0,0 +1,257 @@
+// SPDX-License-Identifier: Apache-2.0
+/**
+ * #602 (2.22.0): highlight reconstruction from neighboring pixels, `highlightMode: "reconstruct-neighbors"` (opt-in; an estimate,
+ * not a recovery of the clipped scene signal). The engine's own method, acting on the developed, white-balanced camera values
+ * u = g·c of the whole image at use (no tables):
+ *
+ * 1. A pixel is clipped when one of its normalized camera channels c_k reaches `NEIGHBOR_HIGHLIGHT_THRESHOLD` of white; clipped
+ *    pixels form regions (8-connected).
+ * 2. Each unclipped pixel has the color ratios ρ = u/(u₀ + u₁ + u₂) (its chromaticity in white-balanced camera space). Pixels
+ *    within `NEIGHBOR_HIGHLIGHT_MARGIN` pixels of a clipped one give no ratios: their demosaic reads clipped sites, which lowers
+ *    the clipped channel (and noise selects them by the threshold); the regions below are the clipped pixels widened by it.
+ * 3. The clipped region takes its ratios from its unclipped boundary. **Fast:** the mean of ρ over the unclipped pixels within
+ *    `NEIGHBOR_HIGHLIGHT_RING` pixels (Chebyshev distance) of the region. **Full:** the discrete harmonic (Laplace) interpolation
+ *    of ρ inward: ρ at each clipped pixel is the mean of its 4-neighbors', with the unclipped neighbors' ratios as Dirichlet
+ *    boundary values (pixels outside the image are left out), solved by successive over-relaxation from the fast values. A region
+ *    with no unclipped neighbor keeps the neutral ratios.
+ * 4. A clipped pixel keeps its unclipped channels; its scale s = Σ_unclipped u_k / Σ_unclipped ρ_k, and each clipped channel
+ *    becomes max(u_k, s·ρ_k) (never below what was recorded). Where every channel is clipped, s = max_k u_k/ρ_k.
+ * The result is not clipped at the white level: the tone curve and gamut handling then render it.
+ *
+ * Mathematics only, general and free: the discrete Laplace equation and its Dirichlet problem, and successive over-relaxation
+ * (e.g. G. Strang, MIT OpenCourseWare 18.086, "Mathematical Methods for Engineers II", finite differences for Laplace's equation
+ * and iterative solvers, free at https://ocw.mit.edu/courses/18-086-mathematical-methods-for-engineers-ii-spring-2006/).
+ */
+/** Normalized channel level at which a channel counts as clipped. */
+export const NEIGHBOR_HIGHLIGHT_THRESHOLD = 0.98;
+/** Pixels this close (Chebyshev distance) to a clipped pixel give no ratios and are filled with the region. */
+export const NEIGHBOR_HIGHLIGHT_MARGIN = 2;
+/** The fast method's ring width in pixels around a (widened) clipped region. */
+export const NEIGHBOR_HIGHLIGHT_RING = 2;
+/** The full method's tolerance (the largest change of a sweep) and iteration cap; the relaxation factor is `neighborHighlightOmega`. */
+export const NEIGHBOR_HIGHLIGHT_SOR = /* @__PURE__ */ Object.freeze({ tolerance: 1e-9, maximumIterations: 20000 });
+/**
+ * #605 (2.23.0): the relaxation factor for a region whose members span a × b pixels: the optimal ω = 2/(1 + √(1 − ρ²)) of SOR
+ * for the discrete Laplace equation on the a × b rectangle, whose Jacobi iteration has spectral radius
+ * ρ = (cos(π/(a + 1)) + cos(π/(b + 1)))/2 (Young's theory of successive over-relaxation; e.g. G. Strang, MIT OpenCourseWare
+ * 18.086, and R. J. LeVeque, "Finite Difference Methods for Ordinary and Partial Differential Equations", SIAM 2007, §4.2).
+ * A region inside the rectangle has a Jacobi radius at most that of the rectangle, so this ω is at or above its optimum (where
+ * SOR converges at the rate ω − 1). The system (each member's neighbor count times its ratio minus its neighbors') is
+ * symmetric positive definite, so SOR converges for every ω in (0, 2) (Ostrowski–Reich). 2.22.0 used ω = 1.8 for every region:
+ * 5,728 sweeps on a 260 × 260 region against about 400 (docs/DISPLAY_VIEW.md).
+ */
+export function neighborHighlightOmega(a, b) {
+    const rho = (Math.cos(Math.PI / (a + 1)) + Math.cos(Math.PI / (b + 1))) / 2;
+    return 2 / (1 + Math.sqrt(Math.max(0, 1 - rho * rho)));
+}
+/**
+ * The region ratios of an image: per pixel the ratios ρ its rebuild uses (`ratios`, 3 per pixel; for an unclipped pixel its own),
+ * and which pixels are clipped. `clippedAt(i)` tells whether pixel i is clipped; `valueAt(i, out)` writes pixel i's u.
+ * Pixels are i = y·width + x.
+ */
+export function neighborHighlightRatios(width, height, clippedAt, valueAt, fast) {
+    const n = width * height, clipped = new Uint8Array(n), unknown = new Uint8Array(n), ratios = new Float64Array(3 * n), u = new Float64Array(3), D = NEIGHBOR_HIGHLIGHT_MARGIN;
+    for (let i = 0; i < n; i++)
+        if (clippedAt(i))
+            clipped[i] = 1;
+    // The clipped pixels widened by the margin (a separable dilation).
+    const rows = new Uint8Array(n);
+    for (let y = 0; y < height; y++)
+        for (let x = 0; x < width; x++)
+            if (clipped[y * width + x] === 1)
+                rows.fill(1, y * width + Math.max(0, x - D), y * width + Math.min(width, x + D + 1));
+    for (let x = 0; x < width; x++)
+        for (let y = 0; y < height; y++)
+            if (rows[y * width + x] === 1)
+                for (let yy = Math.max(0, y - D); yy <= Math.min(height - 1, y + D); yy++)
+                    unknown[yy * width + x] = 1;
+    for (let i = 0; i < n; i++) {
+        if (unknown[i] === 1)
+            continue;
+        valueAt(i, u);
+        const sum = u[0] + u[1] + u[2];
+        if (sum > 0) {
+            ratios[3 * i] = u[0] / sum;
+            ratios[3 * i + 1] = u[1] / sum;
+            ratios[3 * i + 2] = u[2] / sum;
+        }
+        else {
+            ratios[3 * i] = 1 / 3;
+            ratios[3 * i + 1] = 1 / 3;
+            ratios[3 * i + 2] = 1 / 3;
+        }
+    }
+    const label = new Int32Array(n).fill(-1), stamp = new Int32Array(n).fill(-1), queue = new Int32Array(n);
+    let regionCount = 0, iterations = 0, converged = true;
+    for (let start = 0; start < n; start++) {
+        if (unknown[start] === 0 || label[start] >= 0)
+            continue;
+        // The region (8-connected), breadth first.
+        const region = regionCount++;
+        let head = 0, tail = 0;
+        queue[tail++] = start;
+        label[start] = region;
+        while (head < tail) {
+            const p = queue[head++], px = p % width, py = (p - px) / width;
+            for (let dy = -1; dy <= 1; dy++)
+                for (let dx = -1; dx <= 1; dx++) {
+                    const x = px + dx, y = py + dy;
+                    if (x < 0 || y < 0 || x >= width || y >= height)
+                        continue;
+                    const q = y * width + x;
+                    if (unknown[q] === 1 && label[q] < 0) {
+                        label[q] = region;
+                        queue[tail++] = q;
+                    }
+                }
+        }
+        // Fast: the mean ratios over the ring (each unclipped pixel once).
+        let r0 = 0, r1 = 0, r2 = 0, count = 0;
+        const R = NEIGHBOR_HIGHLIGHT_RING;
+        for (let k = 0; k < tail; k++) {
+            const p = queue[k], px = p % width, py = (p - px) / width;
+            for (let y = Math.max(0, py - R); y <= Math.min(height - 1, py + R); y++)
+                for (let x = Math.max(0, px - R); x <= Math.min(width - 1, px + R); x++) {
+                    const q = y * width + x;
+                    if (unknown[q] === 1 || stamp[q] === region)
+                        continue;
+                    stamp[q] = region;
+                    r0 += ratios[3 * q];
+                    r1 += ratios[3 * q + 1];
+                    r2 += ratios[3 * q + 2];
+                    count++;
+                }
+        }
+        const m0 = count > 0 ? r0 / count : 1 / 3, m1 = count > 0 ? r1 / count : 1 / 3, m2 = count > 0 ? r2 / count : 1 / 3;
+        for (let k = 0; k < tail; k++) {
+            const p = queue[k];
+            ratios[3 * p] = m0;
+            ratios[3 * p + 1] = m1;
+            ratios[3 * p + 2] = m2;
+        }
+        if (fast || count === 0)
+            continue;
+        // Full: successive over-relaxation of the discrete Laplace equation over the region, in the region's breadth-first order.
+        const { tolerance, maximumIterations } = NEIGHBOR_HIGHLIGHT_SOR;
+        let x0 = width, x1 = -1, y0 = height, y1 = -1;
+        for (let k = 0; k < tail; k++) {
+            const p = queue[k], px = p % width, py = (p - px) / width;
+            x0 = Math.min(x0, px);
+            x1 = Math.max(x1, px);
+            y0 = Math.min(y0, py);
+            y1 = Math.max(y1, py);
+        }
+        const omega = neighborHighlightOmega(x1 - x0 + 1, y1 - y0 + 1);
+        let sweep = 0, change = Infinity;
+        while (sweep < maximumIterations && change > tolerance) {
+            change = 0;
+            sweep++;
+            for (let k = 0; k < tail; k++) {
+                const p = queue[k], px = p % width, py = (p - px) / width;
+                let a0 = 0, a1 = 0, a2 = 0, m = 0;
+                for (let d = 0; d < 4; d++) {
+                    const x = px + (d === 0 ? -1 : d === 1 ? 1 : 0), y = py + (d === 2 ? -1 : d === 3 ? 1 : 0);
+                    if (x < 0 || y < 0 || x >= width || y >= height)
+                        continue;
+                    const q = y * width + x;
+                    a0 += ratios[3 * q];
+                    a1 += ratios[3 * q + 1];
+                    a2 += ratios[3 * q + 2];
+                    m++;
+                }
+                if (m === 0)
+                    continue;
+                for (let c = 0; c < 3; c++) {
+                    const old = ratios[3 * p + c], target = (c === 0 ? a0 : c === 1 ? a1 : a2) / m, next = old + omega * (target - old);
+                    ratios[3 * p + c] = next;
+                    const delta = Math.abs(next - old);
+                    if (delta > change)
+                        change = delta;
+                }
+            }
+        }
+        iterations = Math.max(iterations, sweep);
+        if (change > tolerance)
+            converged = false;
+    }
+    return { clipped, ratios, regionCount, iterations, converged };
+}
+/**
+ * The rebuild of one clipped pixel (step 4) into b, with its slopes ∂b/∂c (row-major 3 × 3, a.e.) when `slopes` is given.
+ * `c` are the normalized camera channels, `g` the white-balance gains, `rho` the pixel's ratios.
+ */
+export function rebuildClippedPixel(c0, c1, c2, g, rho, b, slopes) {
+    const H = NEIGHBOR_HIGHLIGHT_THRESHOLD, u0 = c0 * g[0], u1 = c1 * g[1], u2 = c2 * g[2];
+    const k0 = c0 < H, k1 = c1 < H, k2 = c2 < H;
+    if (slopes)
+        slopes.fill(0);
+    if (k0 && k1 && k2) {
+        b[0] = u0;
+        b[1] = u1;
+        b[2] = u2;
+        if (slopes) {
+            slopes[0] = g[0];
+            slopes[4] = g[1];
+            slopes[8] = g[2];
+        }
+        return;
+    }
+    const r0 = rho[0], r1 = rho[1], r2 = rho[2];
+    if (!k0 && !k1 && !k2) {
+        const q0 = r0 > 0 ? u0 / r0 : 0, q1 = r1 > 0 ? u1 / r1 : 0, q2 = r2 > 0 ? u2 / r2 : 0, m = q0 > q1 ? (q0 > q2 ? 0 : 2) : (q1 > q2 ? 1 : 2), s = m === 0 ? q0 : m === 1 ? q1 : q2;
+        b[0] = s * r0;
+        b[1] = s * r1;
+        b[2] = s * r2;
+        if (slopes) {
+            const rm = m === 0 ? r0 : m === 1 ? r1 : r2;
+            if (rm > 0)
+                for (let k = 0; k < 3; k++)
+                    slopes[k * 3 + m] = (k === 0 ? r0 : k === 1 ? r1 : r2) * g[m] / rm;
+        }
+        return;
+    }
+    const su = (k0 ? u0 : 0) + (k1 ? u1 : 0) + (k2 ? u2 : 0), sr = (k0 ? r0 : 0) + (k1 ? r1 : 0) + (k2 ? r2 : 0), s = sr > 0 ? su / sr : 0;
+    for (let k = 0; k < 3; k++) {
+        const kept = k === 0 ? k0 : k === 1 ? k1 : k2, uk = k === 0 ? u0 : k === 1 ? u1 : u2, rk = k === 0 ? r0 : k === 1 ? r1 : r2;
+        if (kept) {
+            b[k] = uk;
+            if (slopes)
+                slopes[k * 4] = g[k];
+            continue;
+        }
+        const estimate = s * rk;
+        if (estimate > uk) {
+            b[k] = estimate;
+            if (slopes && sr > 0)
+                for (let j = 0; j < 3; j++)
+                    if (j === 0 ? k0 : j === 1 ? k1 : k2)
+                        slopes[k * 3 + j] = rk * g[j] / sr;
+        }
+        else {
+            b[k] = uk;
+            if (slopes)
+                slopes[k * 4] = g[k];
+        }
+    }
+}
+/**
+ * Rebuilds the clipped pixels of a white-balanced image u (3 per pixel, in place; becomes b), with the normalized channels
+ * c = u/g deciding which channels are clipped.
+ */
+export function reconstructHighlightsFromNeighbors(u, width, height, gains, fast) {
+    const H = NEIGHBOR_HIGHLIGHT_THRESHOLD, g0 = gains[0], g1 = gains[1], g2 = gains[2];
+    const field = neighborHighlightRatios(width, height, (i) => u[3 * i] / g0 >= H || u[3 * i + 1] / g1 >= H || u[3 * i + 2] / g2 >= H, (i, out) => { out[0] = u[3 * i]; out[1] = u[3 * i + 1]; out[2] = u[3 * i + 2]; }, fast);
+    const b = new Float64Array(3);
+    let rebuilt = 0;
+    for (let i = 0; i < width * height; i++) {
+        if (field.clipped[i] === 0)
+            continue;
+        rebuildClippedPixel(u[3 * i] / g0, u[3 * i + 1] / g1, u[3 * i + 2] / g2, gains, field.ratios.subarray(3 * i, 3 * i + 3), b);
+        u[3 * i] = b[0];
+        u[3 * i + 1] = b[1];
+        u[3 * i + 2] = b[2];
+        rebuilt++;
+    }
+    return { method: fast ? "ring-average" : "harmonic-interpolation", reconstructedPixelCount: rebuilt, regionCount: field.regionCount, iterations: field.iterations, converged: field.converged };
+}
+//# sourceMappingURL=neighbor-highlight-reconstruction.js.map

package/dist/color/neighbor-highlight-reconstruction.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"neighbor-highlight-reconstruction.js","sourceRoot":"","sources":["../../src/color/neighbor-highlight-reconstruction.ts"],"names":[],"mappings":"AAAA,sCAAsC;AAEtC;;;;;;;;;;;;;;;;;;;;;;GAsBG;AAEH,qEAAqE;AACrE,MAAM,CAAC,MAAM,4BAA4B,GAAG,IAAI,CAAC;AACjD,+GAA+G;AAC/G,MAAM,CAAC,MAAM,yBAAyB,GAAG,CAAC,CAAC;AAC3C,gFAAgF;AAChF,MAAM,CAAC,MAAM,uBAAuB,GAAG,CAAC,CAAC;AACzC,wIAAwI;AACxI,MAAM,CAAC,MAAM,sBAAsB,GAAG,eAAe,CAAC,MAAM,CAAC,MAAM,CAAC,EAAE,SAAS,EAAE,IAAI,EAAE,iBAAiB,EAAE,KAAK,EAAE,CAAC,CAAC;AAEnH;;;;;;;;;GASG;AACH,MAAM,UAAU,sBAAsB,CAAC,CAAS,EAAE,CAAS;IACzD,MAAM,GAAG,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IAC5E,OAAO,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;AACzD,CAAC;AAYD;;;;GAIG;AACH,MAAM,UAAU,uBAAuB,CAAC,KAAa,EAAE,MAAc,EAAE,SAAiC,EAAE,OAA+C,EAAE,IAAa;IAEtK,MAAM,CAAC,GAAG,KAAK,GAAG,MAAM,EAAE,OAAO,GAAG,IAAI,UAAU,CAAC,CAAC,CAAC,EAAE,OAAO,GAAG,IAAI,UAAU,CAAC,CAAC,CAAC,EAAE,MAAM,GAAG,IAAI,YAAY,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,yBAAyB,CAAC;IAC7K,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;QAAE,IAAI,SAAS,CAAC,CAAC,CAAC;YAAE,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IAC7D,mEAAmE;IACnE,MAAM,IAAI,GAAG,IAAI,UAAU,CAAC,CAAC,CAAC,CAAC;IAC/B,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,IAAI,OAAO,CAAC,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC,KAAK,CAAC;gBAAE,IAAI,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,GAAG,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACxL,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE;QAAE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,EAAE,CAAC,EAAE;YAAE,IAAI,IAAI,CAAC,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC,KAAK,CAAC;gBAAE,KAAK,IAAI,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,IAAI,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE;oBAAE,OAAO,CAAC,EAAE,GAAG,KAAK,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;IACvM,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,IAAI,OAAO,CAAC,CAAC,CAAC,KAAK,CAAC;YAAE,SAAS;QAC/B,OAAO,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACd,MAAM,GAAG,GAAG,CAAC,CAAC,CAAC,CAAE,GAAG,CAAC,CAAC,CAAC,CAAE,GAAG,CAAC,CAAC,CAAC,CAAE,CAAC;QAClC,IAAI,GAAG,GAAG,CAAC,EAAE,CAAC;YAAC,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAE,GAAG,GAAG,CAAC;YAAC,MAAM,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAE,GAAG,GAAG,CAAC;YAAC,MAAM,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAE,GAAG,GAAG,CAAC;QAAC,CAAC;aAC1G,CAAC;YAAC,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;YAAC,MAAM,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;YAAC,MAAM,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QAAC,CAAC;IACvF,CAAC;IACD,MAAM,KAAK,GAAG,IAAI,UAAU,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,KAAK,GAAG,IAAI,UAAU,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,KAAK,GAAG,IAAI,UAAU,CAAC,CAAC,CAAC,CAAC;IACxG,IAAI,WAAW,GAAG,CAAC,EAAE,UAAU,GAAG,CAAC,EAAE,SAAS,GAAG,IAAI,CAAC;IACtD,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,CAAC,EAAE,KAAK,EAAE,EAAE,CAAC;QACvC,IAAI,OAAO,CAAC,KAAK,CAAC,KAAK,CAAC,IAAI,KAAK,CAAC,KAAK,CAAE,IAAI,CAAC;YAAE,SAAS;QACzD,2CAA2C;QAC3C,MAAM,MAAM,GAAG,WAAW,EAAE,CAAC;QAC7B,IAAI,IAAI,GAAG,CAAC,EAAE,IAAI,GAAG,CAAC,CAAC;QACvB,KAAK,CAAC,IAAI,EAAE,CAAC,GAAG,KAAK,CAAC;QAAC,KAAK,CAAC,KAAK,CAAC,GAAG,MAAM,CAAC;QAC7C,OAAO,IAAI,GAAG,IAAI,EAAE,CAAC;YACnB,MAAM,CAAC,GAAG,KAAK,CAAC,IAAI,EAAE,CAAE,EAAE,EAAE,GAAG,CAAC,GAAG,KAAK,EAAE,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,KAAK,CAAC;YAChE,KAAK,IAAI,EAAE,GAAG,CAAC,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,EAAE,EAAE;gBAAE,KAAK,IAAI,EAAE,GAAG,CAAC,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;oBACjE,MAAM,CAAC,GAAG,EAAE,GAAG,EAAE,EAAE,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC;oBAC/B,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,KAAK,IAAI,CAAC,IAAI,MAAM;wBAAE,SAAS;oBAC1D,MAAM,CAAC,GAAG,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC;oBACxB,IAAI,OAAO,CAAC,CAAC,CAAC,KAAK,CAAC,IAAI,KAAK,CAAC,CAAC,CAAE,GAAG,CAAC,EAAE,CAAC;wBAAC,KAAK,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC;wBAAC,KAAK,CAAC,IAAI,EAAE,CAAC,GAAG,CAAC,CAAC;oBAAC,CAAC;gBAClF,CAAC;QACH,CAAC;QACD,mEAAmE;QACnE,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,KAAK,GAAG,CAAC,CAAC;QACtC,MAAM,CAAC,GAAG,uBAAuB,CAAC;QAClC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;YAC9B,MAAM,CAAC,GAAG,KAAK,CAAC,CAAC,CAAE,EAAE,EAAE,GAAG,CAAC,GAAG,KAAK,EAAE,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,KAAK,CAAC;YAC3D,KAAK,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE;gBAAE,KAAK,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;oBAClJ,MAAM,CAAC,GAAG,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC;oBACxB,IAAI,OAAO,CAAC,CAAC,CAAC,KAAK,CAAC,IAAI,KAAK,CAAC,CAAC,CAAC,KAAK,MAAM;wBAAE,SAAS;oBACtD,KAAK,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC;oBAAC,EAAE,IAAI,MAAM,CAAC,CAAC,GAAG,CAAC,CAAE,CAAC;oBAAC,EAAE,IAAI,MAAM,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAE,CAAC;oBAAC,EAAE,IAAI,MAAM,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAE,CAAC;oBAAC,KAAK,EAAE,CAAC;gBACvG,CAAC;QACH,CAAC;QACD,MAAM,EAAE,GAAG,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,EAAE,GAAG,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,EAAE,GAAG,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACpH,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;YAAC,MAAM,CAAC,GAAG,KAAK,CAAC,CAAC,CAAE,CAAC;YAAC,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC;YAAC,MAAM,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC;YAAC,MAAM,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC;QAAC,CAAC;QAC3H,IAAI,IAAI,IAAI,KAAK,KAAK,CAAC;YAAE,SAAS;QAClC,0HAA0H;QAC1H,MAAM,EAAE,SAAS,EAAE,iBAAiB,EAAE,GAAG,sBAAsB,CAAC;QAChE,IAAI,EAAE,GAAG,KAAK,EAAE,EAAE,GAAG,CAAC,CAAC,EAAE,EAAE,GAAG,MAAM,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC;QAC9C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;YAAC,MAAM,CAAC,GAAG,KAAK,CAAC,CAAC,CAAE,EAAE,EAAE,GAAG,CAAC,GAAG,KAAK,EAAE,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,KAAK,CAAC;YAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;YAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;YAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;YAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;QAAC,CAAC;QAC1L,MAAM,KAAK,GAAG,sBAAsB,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC;QAC/D,IAAI,KAAK,GAAG,CAAC,EAAE,MAAM,GAAG,QAAQ,CAAC;QACjC,OAAO,KAAK,GAAG,iBAAiB,IAAI,MAAM,GAAG,SAAS,EAAE,CAAC;YACvD,MAAM,GAAG,CAAC,CAAC;YAAC,KAAK,EAAE,CAAC;YACpB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC9B,MAAM,CAAC,GAAG,KAAK,CAAC,CAAC,CAAE,EAAE,EAAE,GAAG,CAAC,GAAG,KAAK,EAAE,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,KAAK,CAAC;gBAC3D,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC;gBAClC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC3B,MAAM,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBAC3F,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,KAAK,IAAI,CAAC,IAAI,MAAM;wBAAE,SAAS;oBAC1D,MAAM,CAAC,GAAG,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC;oBACxB,EAAE,IAAI,MAAM,CAAC,CAAC,GAAG,CAAC,CAAE,CAAC;oBAAC,EAAE,IAAI,MAAM,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAE,CAAC;oBAAC,EAAE,IAAI,MAAM,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAE,CAAC;oBAAC,CAAC,EAAE,CAAC;gBAChF,CAAC;gBACD,IAAI,CAAC,KAAK,CAAC;oBAAE,SAAS;gBACtB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC3B,MAAM,GAAG,GAAG,MAAM,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAE,EAAE,MAAM,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,IAAI,GAAG,GAAG,GAAG,KAAK,GAAG,CAAC,MAAM,GAAG,GAAG,CAAC,CAAC;oBACrH,MAAM,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC;oBACzB,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,GAAG,GAAG,CAAC,CAAC;oBACnC,IAAI,KAAK,GAAG,MAAM;wBAAE,MAAM,GAAG,KAAK,CAAC;gBACrC,CAAC;YACH,CAAC;QACH,CAAC;QACD,UAAU,GAAG,IAAI,CAAC,GAAG,CAAC,UAAU,EAAE,KAAK,CAAC,CAAC;QACzC,IAAI,MAAM,GAAG,SAAS;YAAE,SAAS,GAAG,KAAK,CAAC;IAC5C,CAAC;IACD,OAAO,EAAE,OAAO,EAAE,MAAM,EAAE,WAAW,EAAE,UAAU,EAAE,SAAS,EAAE,CAAC;AACjE,CAAC;AAED;;;GAGG;AACH,MAAM,UAAU,mBAAmB,CAAC,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,CAAoB,EAAE,GAAsB,EAAE,CAA0B,EAAE,MAAqB;IACrK,MAAM,CAAC,GAAG,4BAA4B,EAAE,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,CAAE,EAAE,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,CAAE,EAAE,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,CAAE,CAAC;IAC1F,MAAM,EAAE,GAAG,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;IAC5C,IAAI,MAAM;QAAE,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAC3B,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,EAAE,CAAC;QACnB,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QAChC,IAAI,MAAM,EAAE,CAAC;YAAC,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAE,CAAC;YAAC,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAE,CAAC;YAAC,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAE,CAAC;QAAC,CAAC;QACxE,OAAO;IACT,CAAC;IACD,MAAM,EAAE,GAAG,GAAG,CAAC,CAAC,CAAE,EAAE,EAAE,GAAG,GAAG,CAAC,CAAC,CAAE,EAAE,EAAE,GAAG,GAAG,CAAC,CAAC,CAAE,CAAC;IAC/C,IAAI,CAAC,EAAE,IAAI,CAAC,EAAE,IAAI,CAAC,EAAE,EAAE,CAAC;QACtB,MAAM,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC;QACjL,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC;QAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC;QAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC;QAC5C,IAAI,MAAM,EAAE,CAAC;YAAC,MAAM,EAAE,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC;YAAC,IAAI,EAAE,GAAG,CAAC;gBAAE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;oBAAE,MAAM,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAE,GAAG,EAAE,CAAC;QAAC,CAAC;QACzK,OAAO;IACT,CAAC;IACD,MAAM,EAAE,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IACvI,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,MAAM,IAAI,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,EAAE,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,EAAE,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC;QAC5H,IAAI,IAAI,EAAE,CAAC;YAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;YAAC,IAAI,MAAM;gBAAE,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAE,CAAC;YAAC,SAAS;QAAC,CAAC;QACrE,MAAM,QAAQ,GAAG,CAAC,GAAG,EAAE,CAAC;QACxB,IAAI,QAAQ,GAAG,EAAE,EAAE,CAAC;YAClB,CAAC,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC;YAChB,IAAI,MAAM,IAAI,EAAE,GAAG,CAAC;gBAAE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;oBAAE,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE;wBAAE,MAAM,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,CAAE,GAAG,EAAE,CAAC;QAC9H,CAAC;aAAM,CAAC;YAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;YAAC,IAAI,MAAM;gBAAE,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAE,CAAC;QAAC,CAAC;IAC1D,CAAC;AACH,CAAC;AAED;;;GAGG;AACH,MAAM,UAAU,kCAAkC,CAAC,CAAe,EAAE,KAAa,EAAE,MAAc,EAAE,KAAwB,EAAE,IAAa;IACxI,MAAM,CAAC,GAAG,4BAA4B,EAAE,EAAE,GAAG,KAAK,CAAC,CAAC,CAAE,EAAE,EAAE,GAAG,KAAK,CAAC,CAAC,CAAE,EAAE,EAAE,GAAG,KAAK,CAAC,CAAC,CAAE,CAAC;IACvF,MAAM,KAAK,GAAG,uBAAuB,CAAC,KAAK,EAAE,MAAM,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAE,GAAG,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAE,GAAG,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAE,GAAG,EAAE,IAAI,CAAC,EACnI,CAAC,CAAC,EAAE,GAAG,EAAE,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAE,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC;IAC7F,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;IAC9B,IAAI,OAAO,GAAG,CAAC,CAAC;IAChB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,GAAG,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACxC,IAAI,KAAK,CAAC,OAAO,CAAC,CAAC,CAAC,KAAK,CAAC;YAAE,SAAS;QACrC,mBAAmB,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAE,GAAG,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAE,GAAG,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAE,GAAG,EAAE,EAAE,KAAK,EAAE,KAAK,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAC/H,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAE,CAAC;QAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAE,CAAC;QAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAE,CAAC;QAC7D,OAAO,EAAE,CAAC;IACZ,CAAC;IACD,OAAO,EAAE,MAAM,EAAE,IAAI,CAAC,CAAC,CAAC,cAAc,CAAC,CAAC,CAAC,wBAAwB,EAAE,uBAAuB,EAAE,OAAO,EAAE,WAAW,EAAE,KAAK,CAAC,WAAW,EAAE,UAAU,EAAE,KAAK,CAAC,UAAU,EAAE,SAAS,EAAE,KAAK,CAAC,SAAS,EAAE,CAAC;AAClM,CAAC","sourcesContent":["// SPDX-License-Identifier: Apache-2.0\n\n/**\n * #602 (2.22.0): highlight reconstruction from neighboring pixels, `highlightMode: \"reconstruct-neighbors\"` (opt-in; an estimate,\n * not a recovery of the clipped scene signal). The engine's own method, acting on the developed, white-balanced camera values\n * u = g·c of the whole image at use (no tables):\n *\n * 1. A pixel is clipped when one of its normalized camera channels c_k reaches `NEIGHBOR_HIGHLIGHT_THRESHOLD` of white; clipped\n *    pixels form regions (8-connected).\n * 2. Each unclipped pixel has the color ratios ρ = u/(u₀ + u₁ + u₂) (its chromaticity in white-balanced camera space). Pixels\n *    within `NEIGHBOR_HIGHLIGHT_MARGIN` pixels of a clipped one give no ratios: their demosaic reads clipped sites, which lowers\n *    the clipped channel (and noise selects them by the threshold); the regions below are the clipped pixels widened by it.\n * 3. The clipped region takes its ratios from its unclipped boundary. **Fast:** the mean of ρ over the unclipped pixels within\n *    `NEIGHBOR_HIGHLIGHT_RING` pixels (Chebyshev distance) of the region. **Full:** the discrete harmonic (Laplace) interpolation\n *    of ρ inward: ρ at each clipped pixel is the mean of its 4-neighbors', with the unclipped neighbors' ratios as Dirichlet\n *    boundary values (pixels outside the image are left out), solved by successive over-relaxation from the fast values. A region\n *    with no unclipped neighbor keeps the neutral ratios.\n * 4. A clipped pixel keeps its unclipped channels; its scale s = Σ_unclipped u_k / Σ_unclipped ρ_k, and each clipped channel\n *    becomes max(u_k, s·ρ_k) (never below what was recorded). Where every channel is clipped, s = max_k u_k/ρ_k.\n * The result is not clipped at the white level: the tone curve and gamut handling then render it.\n *\n * Mathematics only, general and free: the discrete Laplace equation and its Dirichlet problem, and successive over-relaxation\n * (e.g. G. Strang, MIT OpenCourseWare 18.086, \"Mathematical Methods for Engineers II\", finite differences for Laplace's equation\n * and iterative solvers, free at https://ocw.mit.edu/courses/18-086-mathematical-methods-for-engineers-ii-spring-2006/).\n */\n\n/** Normalized channel level at which a channel counts as clipped. */\nexport const NEIGHBOR_HIGHLIGHT_THRESHOLD = 0.98;\n/** Pixels this close (Chebyshev distance) to a clipped pixel give no ratios and are filled with the region. */\nexport const NEIGHBOR_HIGHLIGHT_MARGIN = 2;\n/** The fast method's ring width in pixels around a (widened) clipped region. */\nexport const NEIGHBOR_HIGHLIGHT_RING = 2;\n/** The full method's tolerance (the largest change of a sweep) and iteration cap; the relaxation factor is `neighborHighlightOmega`. */\nexport const NEIGHBOR_HIGHLIGHT_SOR = /* @__PURE__ */ Object.freeze({ tolerance: 1e-9, maximumIterations: 20000 });\n\n/**\n * #605 (2.23.0): the relaxation factor for a region whose members span a × b pixels: the optimal ω = 2/(1 + √(1 − ρ²)) of SOR\n * for the discrete Laplace equation on the a × b rectangle, whose Jacobi iteration has spectral radius\n * ρ = (cos(π/(a + 1)) + cos(π/(b + 1)))/2 (Young's theory of successive over-relaxation; e.g. G. Strang, MIT OpenCourseWare\n * 18.086, and R. J. LeVeque, \"Finite Difference Methods for Ordinary and Partial Differential Equations\", SIAM 2007, §4.2).\n * A region inside the rectangle has a Jacobi radius at most that of the rectangle, so this ω is at or above its optimum (where\n * SOR converges at the rate ω − 1). The system (each member's neighbor count times its ratio minus its neighbors') is\n * symmetric positive definite, so SOR converges for every ω in (0, 2) (Ostrowski–Reich). 2.22.0 used ω = 1.8 for every region:\n * 5,728 sweeps on a 260 × 260 region against about 400 (docs/DISPLAY_VIEW.md).\n */\nexport function neighborHighlightOmega(a: number, b: number): number {\n  const rho = (Math.cos(Math.PI / (a + 1)) + Math.cos(Math.PI / (b + 1))) / 2;\n  return 2 / (1 + Math.sqrt(Math.max(0, 1 - rho * rho)));\n}\n\nexport interface NeighborHighlightReport {\n  method: \"ring-average\" | \"harmonic-interpolation\";\n  /** Clipped pixels rebuilt, and the clipped regions. */\n  reconstructedPixelCount: number;\n  regionCount: number;\n  /** Full method: the largest number of relaxation sweeps a region needed, and whether every region met the tolerance. */\n  iterations: number;\n  converged: boolean;\n}\n\n/**\n * The region ratios of an image: per pixel the ratios ρ its rebuild uses (`ratios`, 3 per pixel; for an unclipped pixel its own),\n * and which pixels are clipped. `clippedAt(i)` tells whether pixel i is clipped; `valueAt(i, out)` writes pixel i's u.\n * Pixels are i = y·width + x.\n */\nexport function neighborHighlightRatios(width: number, height: number, clippedAt: (i: number) => boolean, valueAt: (i: number, out: Float64Array) => void, fast: boolean):\n  { clipped: Uint8Array; ratios: Float64Array; regionCount: number; iterations: number; converged: boolean } {\n  const n = width * height, clipped = new Uint8Array(n), unknown = new Uint8Array(n), ratios = new Float64Array(3 * n), u = new Float64Array(3), D = NEIGHBOR_HIGHLIGHT_MARGIN;\n  for (let i = 0; i < n; i++) if (clippedAt(i)) clipped[i] = 1;\n  // The clipped pixels widened by the margin (a separable dilation).\n  const rows = new Uint8Array(n);\n  for (let y = 0; y < height; y++) for (let x = 0; x < width; x++) if (clipped[y * width + x] === 1) rows.fill(1, y * width + Math.max(0, x - D), y * width + Math.min(width, x + D + 1));\n  for (let x = 0; x < width; x++) for (let y = 0; y < height; y++) if (rows[y * width + x] === 1) for (let yy = Math.max(0, y - D); yy <= Math.min(height - 1, y + D); yy++) unknown[yy * width + x] = 1;\n  for (let i = 0; i < n; i++) {\n    if (unknown[i] === 1) continue;\n    valueAt(i, u);\n    const sum = u[0]! + u[1]! + u[2]!;\n    if (sum > 0) { ratios[3 * i] = u[0]! / sum; ratios[3 * i + 1] = u[1]! / sum; ratios[3 * i + 2] = u[2]! / sum; }\n    else { ratios[3 * i] = 1 / 3; ratios[3 * i + 1] = 1 / 3; ratios[3 * i + 2] = 1 / 3; }\n  }\n  const label = new Int32Array(n).fill(-1), stamp = new Int32Array(n).fill(-1), queue = new Int32Array(n);\n  let regionCount = 0, iterations = 0, converged = true;\n  for (let start = 0; start < n; start++) {\n    if (unknown[start] === 0 || label[start]! >= 0) continue;\n    // The region (8-connected), breadth first.\n    const region = regionCount++;\n    let head = 0, tail = 0;\n    queue[tail++] = start; label[start] = region;\n    while (head < tail) {\n      const p = queue[head++]!, px = p % width, py = (p - px) / width;\n      for (let dy = -1; dy <= 1; dy++) for (let dx = -1; dx <= 1; dx++) {\n        const x = px + dx, y = py + dy;\n        if (x < 0 || y < 0 || x >= width || y >= height) continue;\n        const q = y * width + x;\n        if (unknown[q] === 1 && label[q]! < 0) { label[q] = region; queue[tail++] = q; }\n      }\n    }\n    // Fast: the mean ratios over the ring (each unclipped pixel once).\n    let r0 = 0, r1 = 0, r2 = 0, count = 0;\n    const R = NEIGHBOR_HIGHLIGHT_RING;\n    for (let k = 0; k < tail; k++) {\n      const p = queue[k]!, px = p % width, py = (p - px) / width;\n      for (let y = Math.max(0, py - R); y <= Math.min(height - 1, py + R); y++) for (let x = Math.max(0, px - R); x <= Math.min(width - 1, px + R); x++) {\n        const q = y * width + x;\n        if (unknown[q] === 1 || stamp[q] === region) continue;\n        stamp[q] = region; r0 += ratios[3 * q]!; r1 += ratios[3 * q + 1]!; r2 += ratios[3 * q + 2]!; count++;\n      }\n    }\n    const m0 = count > 0 ? r0 / count : 1 / 3, m1 = count > 0 ? r1 / count : 1 / 3, m2 = count > 0 ? r2 / count : 1 / 3;\n    for (let k = 0; k < tail; k++) { const p = queue[k]!; ratios[3 * p] = m0; ratios[3 * p + 1] = m1; ratios[3 * p + 2] = m2; }\n    if (fast || count === 0) continue;\n    // Full: successive over-relaxation of the discrete Laplace equation over the region, in the region's breadth-first order.\n    const { tolerance, maximumIterations } = NEIGHBOR_HIGHLIGHT_SOR;\n    let x0 = width, x1 = -1, y0 = height, y1 = -1;\n    for (let k = 0; k < tail; k++) { const p = queue[k]!, px = p % width, py = (p - px) / width; x0 = Math.min(x0, px); x1 = Math.max(x1, px); y0 = Math.min(y0, py); y1 = Math.max(y1, py); }\n    const omega = neighborHighlightOmega(x1 - x0 + 1, y1 - y0 + 1);\n    let sweep = 0, change = Infinity;\n    while (sweep < maximumIterations && change > tolerance) {\n      change = 0; sweep++;\n      for (let k = 0; k < tail; k++) {\n        const p = queue[k]!, px = p % width, py = (p - px) / width;\n        let a0 = 0, a1 = 0, a2 = 0, m = 0;\n        for (let d = 0; d < 4; d++) {\n          const x = px + (d === 0 ? -1 : d === 1 ? 1 : 0), y = py + (d === 2 ? -1 : d === 3 ? 1 : 0);\n          if (x < 0 || y < 0 || x >= width || y >= height) continue;\n          const q = y * width + x;\n          a0 += ratios[3 * q]!; a1 += ratios[3 * q + 1]!; a2 += ratios[3 * q + 2]!; m++;\n        }\n        if (m === 0) continue;\n        for (let c = 0; c < 3; c++) {\n          const old = ratios[3 * p + c]!, target = (c === 0 ? a0 : c === 1 ? a1 : a2) / m, next = old + omega * (target - old);\n          ratios[3 * p + c] = next;\n          const delta = Math.abs(next - old);\n          if (delta > change) change = delta;\n        }\n      }\n    }\n    iterations = Math.max(iterations, sweep);\n    if (change > tolerance) converged = false;\n  }\n  return { clipped, ratios, regionCount, iterations, converged };\n}\n\n/**\n * The rebuild of one clipped pixel (step 4) into b, with its slopes ∂b/∂c (row-major 3 × 3, a.e.) when `slopes` is given.\n * `c` are the normalized camera channels, `g` the white-balance gains, `rho` the pixel's ratios.\n */\nexport function rebuildClippedPixel(c0: number, c1: number, c2: number, g: ArrayLike<number>, rho: ArrayLike<number>, b: Float64Array | number[], slopes?: Float64Array): void {\n  const H = NEIGHBOR_HIGHLIGHT_THRESHOLD, u0 = c0 * g[0]!, u1 = c1 * g[1]!, u2 = c2 * g[2]!;\n  const k0 = c0 < H, k1 = c1 < H, k2 = c2 < H;\n  if (slopes) slopes.fill(0);\n  if (k0 && k1 && k2) {\n    b[0] = u0; b[1] = u1; b[2] = u2;\n    if (slopes) { slopes[0] = g[0]!; slopes[4] = g[1]!; slopes[8] = g[2]!; }\n    return;\n  }\n  const r0 = rho[0]!, r1 = rho[1]!, r2 = rho[2]!;\n  if (!k0 && !k1 && !k2) {\n    const q0 = r0 > 0 ? u0 / r0 : 0, q1 = r1 > 0 ? u1 / r1 : 0, q2 = r2 > 0 ? u2 / r2 : 0, m = q0 > q1 ? (q0 > q2 ? 0 : 2) : (q1 > q2 ? 1 : 2), s = m === 0 ? q0 : m === 1 ? q1 : q2;\n    b[0] = s * r0; b[1] = s * r1; b[2] = s * r2;\n    if (slopes) { const rm = m === 0 ? r0 : m === 1 ? r1 : r2; if (rm > 0) for (let k = 0; k < 3; k++) slopes[k * 3 + m] = (k === 0 ? r0 : k === 1 ? r1 : r2) * g[m]! / rm; }\n    return;\n  }\n  const su = (k0 ? u0 : 0) + (k1 ? u1 : 0) + (k2 ? u2 : 0), sr = (k0 ? r0 : 0) + (k1 ? r1 : 0) + (k2 ? r2 : 0), s = sr > 0 ? su / sr : 0;\n  for (let k = 0; k < 3; k++) {\n    const kept = k === 0 ? k0 : k === 1 ? k1 : k2, uk = k === 0 ? u0 : k === 1 ? u1 : u2, rk = k === 0 ? r0 : k === 1 ? r1 : r2;\n    if (kept) { b[k] = uk; if (slopes) slopes[k * 4] = g[k]!; continue; }\n    const estimate = s * rk;\n    if (estimate > uk) {\n      b[k] = estimate;\n      if (slopes && sr > 0) for (let j = 0; j < 3; j++) if (j === 0 ? k0 : j === 1 ? k1 : k2) slopes[k * 3 + j] = rk * g[j]! / sr;\n    } else { b[k] = uk; if (slopes) slopes[k * 4] = g[k]!; }\n  }\n}\n\n/**\n * Rebuilds the clipped pixels of a white-balanced image u (3 per pixel, in place; becomes b), with the normalized channels\n * c = u/g deciding which channels are clipped.\n */\nexport function reconstructHighlightsFromNeighbors(u: Float64Array, width: number, height: number, gains: ArrayLike<number>, fast: boolean): NeighborHighlightReport {\n  const H = NEIGHBOR_HIGHLIGHT_THRESHOLD, g0 = gains[0]!, g1 = gains[1]!, g2 = gains[2]!;\n  const field = neighborHighlightRatios(width, height, (i) => u[3 * i]! / g0 >= H || u[3 * i + 1]! / g1 >= H || u[3 * i + 2]! / g2 >= H,\n    (i, out) => { out[0] = u[3 * i]!; out[1] = u[3 * i + 1]!; out[2] = u[3 * i + 2]!; }, fast);\n  const b = new Float64Array(3);\n  let rebuilt = 0;\n  for (let i = 0; i < width * height; i++) {\n    if (field.clipped[i] === 0) continue;\n    rebuildClippedPixel(u[3 * i]! / g0, u[3 * i + 1]! / g1, u[3 * i + 2]! / g2, gains, field.ratios.subarray(3 * i, 3 * i + 3), b);\n    u[3 * i] = b[0]!; u[3 * i + 1] = b[1]!; u[3 * i + 2] = b[2]!;\n    rebuilt++;\n  }\n  return { method: fast ? \"ring-average\" : \"harmonic-interpolation\", reconstructedPixelCount: rebuilt, regionCount: field.regionCount, iterations: field.iterations, converged: field.converged };\n}\n"]}

package/dist/color/sensor-color-development.d.ts

@@ -33,6 +33,11 @@ export interface SensorColorDevelopment {
     highlightMode: RawHighlightMode;
     /** Develops one normalized camera pixel into out[offset..offset+2] (linear sRGB/D65); updates counts. */
     developInto(c0: number, c1: number, c2: number, out: Float64Array | number[], offset: number, counts: SensorColorDevelopmentCounts): void;
+    /**
+     * #602: the colour step alone, X·M·diag(w)·b, for white-balanced values b that a whole-image highlight step already formed
+     * (`highlightMode: "reconstruct-neighbors"`); the same products and order as the end of developInto.
+     */
+    developWhiteBalancedInto(b0: number, b1: number, b2: number, out: Float64Array | number[], offset: number): void;
     /** Convenience array form of developInto (counts discarded). */
     develop(values: readonly number[]): number[];
 }

package/dist/color/sensor-color-development.js

@@ -117,8 +117,14 @@ export function prepareSensorColorDevelopment(cameraChannelsToXyz, gains, highli
         out[offset + 1] = ((0 + x1[0] * X) + x1[1] * Y) + x1[2] * Z;
         out[offset + 2] = ((0 + x2[0] * X) + x2[1] * Y) + x2[2] * Z;
     }
+    function developWhiteBalancedInto(b0, b1, b2, out, offset) {
+        const X = ((0 + m0[0] * b0) + m0[1] * b1) + m0[2] * b2, Y = ((0 + m1[0] * b0) + m1[1] * b1) + m1[2] * b2, Z = ((0 + m2[0] * b0) + m2[1] * b1) + m2[2] * b2;
+        out[offset] = ((0 + x0[0] * X) + x0[1] * Y) + x0[2] * Z;
+        out[offset + 1] = ((0 + x1[0] * X) + x1[1] * Y) + x1[2] * Z;
+        out[offset + 2] = ((0 + x2[0] * X) + x2[1] * Y) + x2[2] * Z;
+    }
     return {
-        referenceWhiteXyz: model.referenceWhiteXyz, cameraNeutralD65: Object.freeze(neutral), appliedGains: Object.freeze(applied),
+        developWhiteBalancedInto, referenceWhiteXyz: model.referenceWhiteXyz, cameraNeutralD65: Object.freeze(neutral), appliedGains: Object.freeze(applied),
         asShotNeutral: Object.freeze(asShot), cameraToXyz: m, highlightMode, developInto,
         develop(values) {
             const out = [0, 0, 0];

package/dist/color/sensor-color-development.js.map

@@ -1 +1 @@
-{"version":3,"file":"sensor-color-development.js","sourceRoot":"","sources":["../../src/color/sensor-color-development.ts"],"names":[],"mappings":"AAAA,sCAAsC;AAEtC;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAsCG;AAEH,OAAO,EAAE,yBAAyB,EAAE,MAAM,gCAAgC,CAAC;AAC3E,OAAO,EAAE,wBAAwB,EAAE,MAAM,oBAAoB,CAAC;AAG9D,MAAM,CAAC,MAAM,8BAA8B,GAAG,eAAe,CAAC,MAAM,CAAC,MAAM,CAAC,EAAE,EAAE,EAAE,0BAA0B,EAAE,OAAO,EAAE,OAAO,EAAW,CAAC,CAAC;AAG3I,mFAAmF;AACnF,MAAM,CAAC,MAAM,kCAAkC,GAAG,GAAG,CAAC;AACtD,sGAAsG;AACtG,MAAM,0BAA0B,GAAG,KAAK,CAAC;AAEzC,qEAAqE;AACrE,SAAS,QAAQ,CAAC,CAAS,EAAE,MAAyB;IACpD,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,GAAC,CAAC,GAAC,MAAM,CAAC,CAAC,CAAE,EAAE,CAAC,CAAC,CAAC,CAAC;AACxE,CAAC;AACD,oGAAoG;AACpG,MAAM,UAAU,uBAAuB,CAAC,CAAS;IAC/C,MAAM,CAAC,CAAC,EAAC,CAAC,EAAC,CAAC,EAAC,CAAC,EAAC,CAAC,EAAC,CAAC,EAAC,CAAC,EAAC,CAAC,EAAC,CAAC,CAAC,GAAC,CAAC,CAAC,IAAI,EAAc,EAAE,GAAG,GAAC,CAAE,GAAC,CAAC,CAAE,GAAC,CAAE,GAAC,CAAE,GAAC,CAAE,CAAC,GAAC,CAAE,GAAC,CAAC,CAAE,GAAC,CAAE,GAAC,CAAE,GAAC,CAAE,CAAC,GAAC,CAAE,GAAC,CAAC,CAAE,GAAC,CAAE,GAAC,CAAE,GAAC,CAAE,CAAC,CAAC;IACvG,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAC,IAAI;QAAE,MAAM,IAAI,yBAAyB,CAAC,oEAAoE,CAAC,CAAC;IAC3J,MAAM,MAAM,GAAC,CAAC,CAAC,CAAE,GAAC,CAAE,GAAC,CAAE,GAAC,CAAE,EAAC,CAAE,GAAC,CAAE,GAAC,CAAE,GAAC,CAAE,EAAC,CAAE,GAAC,CAAE,GAAC,CAAE,GAAC,CAAE,CAAC,EAAC,CAAC,CAAE,GAAC,CAAE,GAAC,CAAE,GAAC,CAAE,EAAC,CAAE,GAAC,CAAE,GAAC,CAAE,GAAC,CAAE,EAAC,CAAE,GAAC,CAAE,GAAC,CAAE,GAAC,CAAE,CAAC,EAAC,CAAC,CAAE,GAAC,CAAE,GAAC,CAAE,GAAC,CAAE,EAAC,CAAE,GAAC,CAAE,GAAC,CAAE,GAAC,CAAE,EAAC,CAAE,GAAC,CAAE,GAAC,CAAE,GAAC,CAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,GAAC,GAAG,CAAC,CAAC,CAAC;IACjK,IAAI,MAAM,CAAC,IAAI,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAC,GAAG,CAAC;QAAE,MAAM,IAAI,yBAAyB,CAAC,sDAAsD,CAAC,CAAC;IACnK,OAAO,MAAM,CAAC;AAChB,CAAC;AAED,0GAA0G;AAC1G,MAAM,UAAU,yBAAyB,CAAC,WAAmB;IAC3D,MAAM,KAAK,GAAG,wBAAwB,EAAE,CAAC,iBAAiB,CAAC;IAC3D,OAAO,QAAQ,CAAC,WAAW,EAAE,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;AACtD,CAAC;AAsBD;;;;GAIG;AACH,MAAM,UAAU,6BAA6B,CAAC,mBAA2B,EAAE,KAA+B,EACxG,gBAAkC,MAAM;IACxC,IAAI,aAAa,KAAK,MAAM,IAAI,aAAa,KAAK,aAAa;QAAE,MAAM,IAAI,yBAAyB,CAAC,iCAAiC,CAAC,CAAC;IACxI,MAAM,KAAK,GAAG,wBAAwB,EAAE,CAAC;IACzC,MAAM,GAAG,GAAG,yBAAyB,CAAC,uBAAuB,CAAC,mBAAmB,CAAC,CAAC,CAAC;IACpF,IAAI,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;QAAE,MAAM,IAAI,yBAAyB,CAAC,uDAAuD,CAAC,CAAC;IACnJ,IAAI,OAAO,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,GAAC,GAAG,CAAC,CAAC,CAAE,CAAC,CAAC;IACxC,IAAI,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAC,CAAC,CAAC,IAAI,0BAA0B,CAAC;QAAE,OAAO,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAC3F,IAAI,KAAK,IAAI,CAAC,KAAK,CAAC,MAAM,KAAK,CAAC,IAAI,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC;QAAE,MAAM,IAAI,yBAAyB,CAAC,kDAAkD,CAAC,CAAC;IACjL,MAAM,OAAO,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,GAAC,CAAC,CAAC,CAAC;IAC7D,MAAM,MAAM,GAAG,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,OAAO,CAAC,CAAC,CAAE,GAAC,CAAC,CAAC,EAAE,SAAS,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,OAAO,CAAC,CAAC;IACnF,MAAM,CAAC,GAAG,mBAAmB,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAC,OAAO,CAAC,CAAC,CAAE,CAAC,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,cAAc,CAAC;IACvG,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,CAAE,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC,CAAE,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC,CAAE,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC,CAAE,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC,CAAE,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC,CAAE,CAAC;IAC7E,MAAM,EAAE,GAAG,OAAO,CAAC,CAAC,CAAE,EAAE,EAAE,GAAG,OAAO,CAAC,CAAC,CAAE,EAAE,EAAE,GAAG,OAAO,CAAC,CAAC,CAAE,EAAE,WAAW,GAAG,aAAa,KAAK,aAAa,CAAC;IAC1G,MAAM,KAAK,GAAG,kCAAkC,CAAC;IACjD,SAAS,WAAW,CAAC,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,GAA4B,EAAE,MAAc,EAAE,MAAoC;QACzI,MAAM,EAAE,GAAG,EAAE,GAAC,EAAE,EAAE,EAAE,GAAG,EAAE,GAAC,EAAE,EAAE,EAAE,GAAG,EAAE,GAAC,EAAE,CAAC;QACzC,IAAI,EAAU,EAAE,EAAU,EAAE,EAAU,CAAC;QACvC,MAAM,EAAE,GAAG,EAAE,GAAG,SAAS,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,EAAE,EAAE,EAAE,GAAG,EAAE,GAAG,SAAS,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,EAAE,EAAE,EAAE,GAAG,EAAE,GAAG,SAAS,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC;QACvH,IAAI,CAAC,WAAW,EAAE,CAAC;YACjB,EAAE,GAAG,EAAE,CAAC;YAAC,EAAE,GAAG,EAAE,CAAC;YAAC,EAAE,GAAG,EAAE,CAAC;YAC1B,MAAM,CAAC,8BAA8B,IAAI,MAAM,CAAC,EAAE,GAAG,SAAS,CAAC,GAAC,MAAM,CAAC,EAAE,GAAG,SAAS,CAAC,GAAC,MAAM,CAAC,EAAE,GAAG,SAAS,CAAC,CAAC;QAChH,CAAC;aAAM,CAAC;YACN,MAAM,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;YAC9D,IAAI,CAAC,GAAG,KAAK,EAAE,CAAC;gBACd,MAAM,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAC,KAAK,CAAC,GAAC,CAAC,CAAC,GAAC,KAAK,CAAC,EAAE,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;gBACpG,EAAE,GAAG,CAAC,CAAC,GAAC,CAAC,CAAC,GAAC,EAAE,GAAC,CAAC,GAAC,CAAC,CAAC;gBAAC,EAAE,GAAG,CAAC,CAAC,GAAC,CAAC,CAAC,GAAC,EAAE,GAAC,CAAC,GAAC,CAAC,CAAC;gBAAC,EAAE,GAAG,CAAC,CAAC,GAAC,CAAC,CAAC,GAAC,EAAE,GAAC,CAAC,GAAC,CAAC,CAAC;gBAAC,MAAM,CAAC,uBAAuB,EAAE,CAAC;YAC5F,CAAC;iBAAM,CAAC;gBAAC,EAAE,GAAG,EAAE,CAAC;gBAAC,EAAE,GAAG,EAAE,CAAC;gBAAC,EAAE,GAAG,EAAE,CAAC;YAAC,CAAC;QACvC,CAAC;QACD,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,EAAE,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,EAAE,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,EAAE,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,EAAE,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,EAAE,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,EAAE,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,EAAE,CAAC;QAChI,GAAG,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,CAAC,CAAC;QAC/C,GAAG,CAAC,MAAM,GAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,CAAC,CAAC;QACjD,GAAG,CAAC,MAAM,GAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,CAAC,CAAC;IACnD,CAAC;IACD,OAAO;QACL,iBAAiB,EAAE,KAAK,CAAC,iBAAiB,EAAE,gBAAgB,EAAE,MAAM,CAAC,MAAM,CAAC,OAAO,CAAC,EAAE,YAAY,EAAE,MAAM,CAAC,MAAM,CAAC,OAAO,CAAC;QAC1H,aAAa,EAAE,MAAM,CAAC,MAAM,CAAC,MAAM,CAAC,EAAE,WAAW,EAAE,CAAC,EAAE,aAAa,EAAE,WAAW;QAChF,OAAO,CAAC,MAAyB;YAC/B,MAAM,GAAG,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;YAAC,WAAW,CAAC,MAAM,CAAC,CAAC,CAAE,EAAE,MAAM,CAAC,CAAC,CAAE,EAAE,MAAM,CAAC,CAAC,CAAE,EAAE,GAAG,EAAE,CAAC,EAAE,EAAE,8BAA8B,EAAE,CAAC,EAAE,uBAAuB,EAAE,CAAC,EAAE,CAAC,CAAC;YAAC,OAAO,GAAG,CAAC;QAChK,CAAC;KACF,CAAC;AACJ,CAAC","sourcesContent":["// SPDX-License-Identifier: Apache-2.0\n\n/**\n * RAW colour development: turns normalized sensor RGB into linear sRGB, using a declared D65 camera\n * profile. Internal.\n *\n * The caller is responsible for the profile's evidence, for binding it to the CFA (colour filter\n * array) and colour basis, and for applying white balance (WB) only once. Only validated 3×3 matrices\n * are accepted, and their declared determinant and inverse bounds are kept.\n *\n * Development model (`sensor-color-development` 2.0.0). The input c is the normalized camera\n * channels: black level subtracted, divided by the effective white level and clipped at 1 per site,\n * all earlier in the pipeline.\n *\n *   b = highlight(diag(g)·c),   linear sRGB = X · M · diag(w) · b\n *\n * Where:\n * - M is the profile's camera→XYZ matrix and X is the XYZ→linear-sRGB matrix.\n * - w = M⁻¹·XYZ(D65), normalized to green. These are the camera-neutral coordinates of D65\n *   (DNG 1.7 \"camera neutral\").\n * - The WB gains g map the scene white to equal channel values. diag(w) then maps equal channel values\n *   to D65, for any invertible M. So a grey card lit by the gains' illuminant comes out neutral,\n *   whether or not M is white-normalized.\n * - Without a resolved WB intent, g = 1/w, and development reduces to X·M·c (a D65 scene white).\n * - The camera neutral of the scene white, written to DNG `AsShotNeutral`, is g_green/g.\n * This module does not chromatically adapt the scene, infer calibration, change geometry or write\n * container metadata.\n *\n * Highlight handling (after each site is clipped at the sensor white level):\n * - `clip` (default): every WB-scaled channel is clipped at the smallest gain. That is the level at\n *   which the least-amplified channel saturates, so fully clipped pixels stay neutral. This is the\n *   documented default of dcraw (`-H 0`, \"clip all highlights to solid white\") and of most raw\n *   converters.\n * - `reconstruct`: a desaturation blend defined by this engine. Let s be the largest normalized\n *   channel. The weight t = clamp((s − 0.8)/0.2, 0, 1) blends the unclipped WB result u toward the\n *   neutral n = max_k min(u_k, g_min): b = (1 − t)·u + t·n. Below 0.8 of white, colour is unchanged\n *   and unclipped; at white it is neutral. This approximation follows the idea of highlight\n *   desaturation (Guo, Cheng & Hu, \"Correcting over-exposure in photographs\", CVPR 2010). It does not\n *   recover the clipped scene signal.\n * @see docs/PHOTOGRAPHIC_EXPORT.md and docs/NATIVE_RAW_AND_PRINT.md for equations and limits.\n */\n\nimport { InvalidConfigurationError } from \"../core/configuration-error.js\";\nimport { resolveCaptureColorModel } from \"./capture-color.js\";\n\ntype Matrix = readonly (readonly number[])[];\nexport const SENSOR_COLOR_DEVELOPMENT_MODEL = /* @__PURE__ */ Object.freeze({ id: \"sensor-color-development\", version: \"2.0.0\" } as const);\n/** RAW development highlight policy; see the module notes. */\nexport type RawHighlightMode = \"clip\" | \"reconstruct\";\n/** Normalized channel level where `reconstruct` starts blending toward neutral. */\nexport const RAW_HIGHLIGHT_RECONSTRUCTION_START = 0.8;\n/** Profiles whose camera neutral is within this relative distance of [1,1,1] are white-normalized. */\nconst WHITE_NORMALIZED_TOLERANCE = 1e-12;\n\n/** Multiply in declared row order: ((0 + m0·v0) + m1·v1) + m2·v2. */\nfunction multiply(m: Matrix, values: readonly number[]): number[] {\n  return m.map((row) => row.reduce((sum, v, i) => sum+v*values[i]!, 0));\n}\n/** Accepts validated 3×3 matrices only, and keeps their declared determinant and inverse bounds. */\nexport function invertSensorColorMatrix(m: Matrix): number[][] {\n  const [a,b,c,d,e,f,g,h,i]=m.flat() as number[], det=a!*(e!*i!-f!*h!)-b!*(d!*i!-f!*g!)+c!*(d!*h!-e!*g!);\n  if (!Number.isFinite(det) || Math.abs(det)<1e-9) throw new InvalidConfigurationError(\"Export color matrix is singular or outside supported conditioning.\");\n  const result=[[e!*i!-f!*h!,c!*h!-b!*i!,b!*f!-c!*e!],[f!*g!-d!*i!,a!*i!-c!*g!,c!*d!-a!*f!],[d!*h!-e!*g!,b!*g!-a!*h!,a!*e!-b!*d!]].map((r) => r.map((v) => v/det));\n  if (result.flat().some((v) => !Number.isFinite(v) || Math.abs(v)>100)) throw new InvalidConfigurationError(\"Export inverse color matrix exceeds supported range.\");\n  return result;\n}\n\n/** Camera coordinates of the declared D65 white (not normalized). The caller decides container policy. */\nexport function sensorColorReferenceWhite(xyzToCamera: Matrix): number[] {\n  const white = resolveCaptureColorModel().referenceWhiteXyz;\n  return multiply(xyzToCamera, [white.x, 1, white.z]);\n}\n\n/** Per-pixel highlight/WB counters accumulated by {@link SensorColorDevelopment.developInto}. */\nexport interface SensorColorDevelopmentCounts { whiteBalanceClippedSampleCount: number; reconstructedPixelCount: number }\n\nexport interface SensorColorDevelopment {\n  referenceWhiteXyz: { x: number; y: 1; z: number };\n  /** Green-normalized camera-neutral coordinates of D65 (exactly [1,1,1] for white-normalized profiles). */\n  cameraNeutralD65: readonly number[];\n  /** Applied channel gains: the resolved WB gains, or 1/cameraNeutralD65 when WB is not required. */\n  appliedGains: readonly number[];\n  /** Green-normalized camera neutral of the assumed scene white: g_green/g (DNG AsShotNeutral). */\n  asShotNeutral: readonly number[];\n  /** X·M·diag(w) rows, applied as two products in the documented order. */\n  cameraToXyz: readonly (readonly number[])[];\n  highlightMode: RawHighlightMode;\n  /** Develops one normalized camera pixel into out[offset..offset+2] (linear sRGB/D65); updates counts. */\n  developInto(c0: number, c1: number, c2: number, out: Float64Array | number[], offset: number, counts: SensorColorDevelopmentCounts): void;\n  /** Convenience array form of developInto (counts discarded). */\n  develop(values: readonly number[]): number[];\n}\n\n/**\n * Internal development of validated, normalized sensor RGB with a declared D65 profile.\n * `gains` is the resolved WB intent, or null when WB is not required. Both the packed native path and\n * the bounded reference path use the returned operator, so there is one implementation.\n */\nexport function prepareSensorColorDevelopment(cameraChannelsToXyz: Matrix, gains: readonly number[] | null,\n  highlightMode: RawHighlightMode = \"clip\"): SensorColorDevelopment {\n  if (highlightMode !== \"clip\" && highlightMode !== \"reconstruct\") throw new InvalidConfigurationError(\"Unsupported RAW highlight mode.\");\n  const model = resolveCaptureColorModel();\n  const raw = sensorColorReferenceWhite(invertSensorColorMatrix(cameraChannelsToXyz));\n  if (raw.some((x) => !(x > 0) || !Number.isFinite(x))) throw new InvalidConfigurationError(\"Reference white requires positive camera coordinates.\");\n  let neutral = raw.map((x) => x/raw[1]!);\n  if (neutral.every((x) => Math.abs(x-1) <= WHITE_NORMALIZED_TOLERANCE)) neutral = [1, 1, 1];\n  if (gains && (gains.length !== 3 || gains.some((x) => !(x > 0) || !Number.isFinite(x)))) throw new InvalidConfigurationError(\"White-balance gains must be positive and finite.\");\n  const applied = gains ? [...gains] : neutral.map((x) => 1/x);\n  const asShot = applied.map((x) => applied[1]!/x), clipLevel = Math.min(...applied);\n  const m = cameraChannelsToXyz.map((row) => row.map((v, i) => v*neutral[i]!)), x = model.xyzToCameraRgb;\n  const m0 = m[0]!, m1 = m[1]!, m2 = m[2]!, x0 = x[0]!, x1 = x[1]!, x2 = x[2]!;\n  const g0 = applied[0]!, g1 = applied[1]!, g2 = applied[2]!, reconstruct = highlightMode === \"reconstruct\";\n  const start = RAW_HIGHLIGHT_RECONSTRUCTION_START;\n  function developInto(c0: number, c1: number, c2: number, out: Float64Array | number[], offset: number, counts: SensorColorDevelopmentCounts): void {\n    const u0 = c0*g0, u1 = c1*g1, u2 = c2*g2;\n    let b0: number, b1: number, b2: number;\n    const k0 = u0 > clipLevel ? clipLevel : u0, k1 = u1 > clipLevel ? clipLevel : u1, k2 = u2 > clipLevel ? clipLevel : u2;\n    if (!reconstruct) {\n      b0 = k0; b1 = k1; b2 = k2;\n      counts.whiteBalanceClippedSampleCount += Number(u0 > clipLevel)+Number(u1 > clipLevel)+Number(u2 > clipLevel);\n    } else {\n      const s = c0 > c1 ? (c0 > c2 ? c0 : c2) : (c1 > c2 ? c1 : c2);\n      if (s > start) {\n        const t = s >= 1 ? 1 : (s-start)/(1-start), n = k0 > k1 ? (k0 > k2 ? k0 : k2) : (k1 > k2 ? k1 : k2);\n        b0 = (1-t)*u0+t*n; b1 = (1-t)*u1+t*n; b2 = (1-t)*u2+t*n; counts.reconstructedPixelCount++;\n      } else { b0 = u0; b1 = u1; b2 = u2; }\n    }\n    const X = ((0+m0[0]!*b0)+m0[1]!*b1)+m0[2]!*b2, Y = ((0+m1[0]!*b0)+m1[1]!*b1)+m1[2]!*b2, Z = ((0+m2[0]!*b0)+m2[1]!*b1)+m2[2]!*b2;\n    out[offset] = ((0+x0[0]!*X)+x0[1]!*Y)+x0[2]!*Z;\n    out[offset+1] = ((0+x1[0]!*X)+x1[1]!*Y)+x1[2]!*Z;\n    out[offset+2] = ((0+x2[0]!*X)+x2[1]!*Y)+x2[2]!*Z;\n  }\n  return {\n    referenceWhiteXyz: model.referenceWhiteXyz, cameraNeutralD65: Object.freeze(neutral), appliedGains: Object.freeze(applied),\n    asShotNeutral: Object.freeze(asShot), cameraToXyz: m, highlightMode, developInto,\n    develop(values: readonly number[]): number[] {\n      const out = [0, 0, 0]; developInto(values[0]!, values[1]!, values[2]!, out, 0, { whiteBalanceClippedSampleCount: 0, reconstructedPixelCount: 0 }); return out;\n    }\n  };\n}\n"]}
+{"version":3,"file":"sensor-color-development.js","sourceRoot":"","sources":["../../src/color/sensor-color-development.ts"],"names":[],"mappings":"AAAA,sCAAsC;AAEtC;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAsCG;AAEH,OAAO,EAAE,yBAAyB,EAAE,MAAM,gCAAgC,CAAC;AAC3E,OAAO,EAAE,wBAAwB,EAAE,MAAM,oBAAoB,CAAC;AAG9D,MAAM,CAAC,MAAM,8BAA8B,GAAG,eAAe,CAAC,MAAM,CAAC,MAAM,CAAC,EAAE,EAAE,EAAE,0BAA0B,EAAE,OAAO,EAAE,OAAO,EAAW,CAAC,CAAC;AAG3I,mFAAmF;AACnF,MAAM,CAAC,MAAM,kCAAkC,GAAG,GAAG,CAAC;AACtD,sGAAsG;AACtG,MAAM,0BAA0B,GAAG,KAAK,CAAC;AAEzC,qEAAqE;AACrE,SAAS,QAAQ,CAAC,CAAS,EAAE,MAAyB;IACpD,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,GAAC,CAAC,GAAC,MAAM,CAAC,CAAC,CAAE,EAAE,CAAC,CAAC,CAAC,CAAC;AACxE,CAAC;AACD,oGAAoG;AACpG,MAAM,UAAU,uBAAuB,CAAC,CAAS;IAC/C,MAAM,CAAC,CAAC,EAAC,CAAC,EAAC,CAAC,EAAC,CAAC,EAAC,CAAC,EAAC,CAAC,EAAC,CAAC,EAAC,CAAC,EAAC,CAAC,CAAC,GAAC,CAAC,CAAC,IAAI,EAAc,EAAE,GAAG,GAAC,CAAE,GAAC,CAAC,CAAE,GAAC,CAAE,GAAC,CAAE,GAAC,CAAE,CAAC,GAAC,CAAE,GAAC,CAAC,CAAE,GAAC,CAAE,GAAC,CAAE,GAAC,CAAE,CAAC,GAAC,CAAE,GAAC,CAAC,CAAE,GAAC,CAAE,GAAC,CAAE,GAAC,CAAE,CAAC,CAAC;IACvG,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAC,IAAI;QAAE,MAAM,IAAI,yBAAyB,CAAC,oEAAoE,CAAC,CAAC;IAC3J,MAAM,MAAM,GAAC,CAAC,CAAC,CAAE,GAAC,CAAE,GAAC,CAAE,GAAC,CAAE,EAAC,CAAE,GAAC,CAAE,GAAC,CAAE,GAAC,CAAE,EAAC,CAAE,GAAC,CAAE,GAAC,CAAE,GAAC,CAAE,CAAC,EAAC,CAAC,CAAE,GAAC,CAAE,GAAC,CAAE,GAAC,CAAE,EAAC,CAAE,GAAC,CAAE,GAAC,CAAE,GAAC,CAAE,EAAC,CAAE,GAAC,CAAE,GAAC,CAAE,GAAC,CAAE,CAAC,EAAC,CAAC,CAAE,GAAC,CAAE,GAAC,CAAE,GAAC,CAAE,EAAC,CAAE,GAAC,CAAE,GAAC,CAAE,GAAC,CAAE,EAAC,CAAE,GAAC,CAAE,GAAC,CAAE,GAAC,CAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,GAAC,GAAG,CAAC,CAAC,CAAC;IACjK,IAAI,MAAM,CAAC,IAAI,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAC,GAAG,CAAC;QAAE,MAAM,IAAI,yBAAyB,CAAC,sDAAsD,CAAC,CAAC;IACnK,OAAO,MAAM,CAAC;AAChB,CAAC;AAED,0GAA0G;AAC1G,MAAM,UAAU,yBAAyB,CAAC,WAAmB;IAC3D,MAAM,KAAK,GAAG,wBAAwB,EAAE,CAAC,iBAAiB,CAAC;IAC3D,OAAO,QAAQ,CAAC,WAAW,EAAE,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;AACtD,CAAC;AA2BD;;;;GAIG;AACH,MAAM,UAAU,6BAA6B,CAAC,mBAA2B,EAAE,KAA+B,EACxG,gBAAkC,MAAM;IACxC,IAAI,aAAa,KAAK,MAAM,IAAI,aAAa,KAAK,aAAa;QAAE,MAAM,IAAI,yBAAyB,CAAC,iCAAiC,CAAC,CAAC;IACxI,MAAM,KAAK,GAAG,wBAAwB,EAAE,CAAC;IACzC,MAAM,GAAG,GAAG,yBAAyB,CAAC,uBAAuB,CAAC,mBAAmB,CAAC,CAAC,CAAC;IACpF,IAAI,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;QAAE,MAAM,IAAI,yBAAyB,CAAC,uDAAuD,CAAC,CAAC;IACnJ,IAAI,OAAO,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,GAAC,GAAG,CAAC,CAAC,CAAE,CAAC,CAAC;IACxC,IAAI,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAC,CAAC,CAAC,IAAI,0BAA0B,CAAC;QAAE,OAAO,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAC3F,IAAI,KAAK,IAAI,CAAC,KAAK,CAAC,MAAM,KAAK,CAAC,IAAI,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC;QAAE,MAAM,IAAI,yBAAyB,CAAC,kDAAkD,CAAC,CAAC;IACjL,MAAM,OAAO,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,GAAC,CAAC,CAAC,CAAC;IAC7D,MAAM,MAAM,GAAG,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,OAAO,CAAC,CAAC,CAAE,GAAC,CAAC,CAAC,EAAE,SAAS,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,OAAO,CAAC,CAAC;IACnF,MAAM,CAAC,GAAG,mBAAmB,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAC,OAAO,CAAC,CAAC,CAAE,CAAC,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,cAAc,CAAC;IACvG,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,CAAE,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC,CAAE,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC,CAAE,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC,CAAE,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC,CAAE,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC,CAAE,CAAC;IAC7E,MAAM,EAAE,GAAG,OAAO,CAAC,CAAC,CAAE,EAAE,EAAE,GAAG,OAAO,CAAC,CAAC,CAAE,EAAE,EAAE,GAAG,OAAO,CAAC,CAAC,CAAE,EAAE,WAAW,GAAG,aAAa,KAAK,aAAa,CAAC;IAC1G,MAAM,KAAK,GAAG,kCAAkC,CAAC;IACjD,SAAS,WAAW,CAAC,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,GAA4B,EAAE,MAAc,EAAE,MAAoC;QACzI,MAAM,EAAE,GAAG,EAAE,GAAC,EAAE,EAAE,EAAE,GAAG,EAAE,GAAC,EAAE,EAAE,EAAE,GAAG,EAAE,GAAC,EAAE,CAAC;QACzC,IAAI,EAAU,EAAE,EAAU,EAAE,EAAU,CAAC;QACvC,MAAM,EAAE,GAAG,EAAE,GAAG,SAAS,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,EAAE,EAAE,EAAE,GAAG,EAAE,GAAG,SAAS,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,EAAE,EAAE,EAAE,GAAG,EAAE,GAAG,SAAS,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC;QACvH,IAAI,CAAC,WAAW,EAAE,CAAC;YACjB,EAAE,GAAG,EAAE,CAAC;YAAC,EAAE,GAAG,EAAE,CAAC;YAAC,EAAE,GAAG,EAAE,CAAC;YAC1B,MAAM,CAAC,8BAA8B,IAAI,MAAM,CAAC,EAAE,GAAG,SAAS,CAAC,GAAC,MAAM,CAAC,EAAE,GAAG,SAAS,CAAC,GAAC,MAAM,CAAC,EAAE,GAAG,SAAS,CAAC,CAAC;QAChH,CAAC;aAAM,CAAC;YACN,MAAM,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;YAC9D,IAAI,CAAC,GAAG,KAAK,EAAE,CAAC;gBACd,MAAM,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAC,KAAK,CAAC,GAAC,CAAC,CAAC,GAAC,KAAK,CAAC,EAAE,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;gBACpG,EAAE,GAAG,CAAC,CAAC,GAAC,CAAC,CAAC,GAAC,EAAE,GAAC,CAAC,GAAC,CAAC,CAAC;gBAAC,EAAE,GAAG,CAAC,CAAC,GAAC,CAAC,CAAC,GAAC,EAAE,GAAC,CAAC,GAAC,CAAC,CAAC;gBAAC,EAAE,GAAG,CAAC,CAAC,GAAC,CAAC,CAAC,GAAC,EAAE,GAAC,CAAC,GAAC,CAAC,CAAC;gBAAC,MAAM,CAAC,uBAAuB,EAAE,CAAC;YAC5F,CAAC;iBAAM,CAAC;gBAAC,EAAE,GAAG,EAAE,CAAC;gBAAC,EAAE,GAAG,EAAE,CAAC;gBAAC,EAAE,GAAG,EAAE,CAAC;YAAC,CAAC;QACvC,CAAC;QACD,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,EAAE,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,EAAE,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,EAAE,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,EAAE,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,EAAE,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,EAAE,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,EAAE,CAAC;QAChI,GAAG,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,CAAC,CAAC;QAC/C,GAAG,CAAC,MAAM,GAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,CAAC,CAAC;QACjD,GAAG,CAAC,MAAM,GAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,CAAC,CAAC;IACnD,CAAC;IACD,SAAS,wBAAwB,CAAC,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,GAA4B,EAAE,MAAc;QAChH,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,EAAE,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,EAAE,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,EAAE,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,EAAE,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,EAAE,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,EAAE,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,EAAE,CAAC;QAChI,GAAG,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,CAAC,CAAC;QAC/C,GAAG,CAAC,MAAM,GAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,CAAC,CAAC;QACjD,GAAG,CAAC,MAAM,GAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC,CAAE,GAAC,CAAC,CAAC;IACnD,CAAC;IACD,OAAO;QACL,wBAAwB,EAAE,iBAAiB,EAAE,KAAK,CAAC,iBAAiB,EAAE,gBAAgB,EAAE,MAAM,CAAC,MAAM,CAAC,OAAO,CAAC,EAAE,YAAY,EAAE,MAAM,CAAC,MAAM,CAAC,OAAO,CAAC;QACpJ,aAAa,EAAE,MAAM,CAAC,MAAM,CAAC,MAAM,CAAC,EAAE,WAAW,EAAE,CAAC,EAAE,aAAa,EAAE,WAAW;QAChF,OAAO,CAAC,MAAyB;YAC/B,MAAM,GAAG,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;YAAC,WAAW,CAAC,MAAM,CAAC,CAAC,CAAE,EAAE,MAAM,CAAC,CAAC,CAAE,EAAE,MAAM,CAAC,CAAC,CAAE,EAAE,GAAG,EAAE,CAAC,EAAE,EAAE,8BAA8B,EAAE,CAAC,EAAE,uBAAuB,EAAE,CAAC,EAAE,CAAC,CAAC;YAAC,OAAO,GAAG,CAAC;QAChK,CAAC;KACF,CAAC;AACJ,CAAC","sourcesContent":["// SPDX-License-Identifier: Apache-2.0\n\n/**\n * RAW colour development: turns normalized sensor RGB into linear sRGB, using a declared D65 camera\n * profile. Internal.\n *\n * The caller is responsible for the profile's evidence, for binding it to the CFA (colour filter\n * array) and colour basis, and for applying white balance (WB) only once. Only validated 3×3 matrices\n * are accepted, and their declared determinant and inverse bounds are kept.\n *\n * Development model (`sensor-color-development` 2.0.0). The input c is the normalized camera\n * channels: black level subtracted, divided by the effective white level and clipped at 1 per site,\n * all earlier in the pipeline.\n *\n *   b = highlight(diag(g)·c),   linear sRGB = X · M · diag(w) · b\n *\n * Where:\n * - M is the profile's camera→XYZ matrix and X is the XYZ→linear-sRGB matrix.\n * - w = M⁻¹·XYZ(D65), normalized to green. These are the camera-neutral coordinates of D65\n *   (DNG 1.7 \"camera neutral\").\n * - The WB gains g map the scene white to equal channel values. diag(w) then maps equal channel values\n *   to D65, for any invertible M. So a grey card lit by the gains' illuminant comes out neutral,\n *   whether or not M is white-normalized.\n * - Without a resolved WB intent, g = 1/w, and development reduces to X·M·c (a D65 scene white).\n * - The camera neutral of the scene white, written to DNG `AsShotNeutral`, is g_green/g.\n * This module does not chromatically adapt the scene, infer calibration, change geometry or write\n * container metadata.\n *\n * Highlight handling (after each site is clipped at the sensor white level):\n * - `clip` (default): every WB-scaled channel is clipped at the smallest gain. That is the level at\n *   which the least-amplified channel saturates, so fully clipped pixels stay neutral. This is the\n *   documented default of dcraw (`-H 0`, \"clip all highlights to solid white\") and of most raw\n *   converters.\n * - `reconstruct`: a desaturation blend defined by this engine. Let s be the largest normalized\n *   channel. The weight t = clamp((s − 0.8)/0.2, 0, 1) blends the unclipped WB result u toward the\n *   neutral n = max_k min(u_k, g_min): b = (1 − t)·u + t·n. Below 0.8 of white, colour is unchanged\n *   and unclipped; at white it is neutral. This approximation follows the idea of highlight\n *   desaturation (Guo, Cheng & Hu, \"Correcting over-exposure in photographs\", CVPR 2010). It does not\n *   recover the clipped scene signal.\n * @see docs/PHOTOGRAPHIC_EXPORT.md and docs/NATIVE_RAW_AND_PRINT.md for equations and limits.\n */\n\nimport { InvalidConfigurationError } from \"../core/configuration-error.js\";\nimport { resolveCaptureColorModel } from \"./capture-color.js\";\n\ntype Matrix = readonly (readonly number[])[];\nexport const SENSOR_COLOR_DEVELOPMENT_MODEL = /* @__PURE__ */ Object.freeze({ id: \"sensor-color-development\", version: \"2.0.0\" } as const);\n/** RAW development highlight policy; see the module notes. */\nexport type RawHighlightMode = \"clip\" | \"reconstruct\";\n/** Normalized channel level where `reconstruct` starts blending toward neutral. */\nexport const RAW_HIGHLIGHT_RECONSTRUCTION_START = 0.8;\n/** Profiles whose camera neutral is within this relative distance of [1,1,1] are white-normalized. */\nconst WHITE_NORMALIZED_TOLERANCE = 1e-12;\n\n/** Multiply in declared row order: ((0 + m0·v0) + m1·v1) + m2·v2. */\nfunction multiply(m: Matrix, values: readonly number[]): number[] {\n  return m.map((row) => row.reduce((sum, v, i) => sum+v*values[i]!, 0));\n}\n/** Accepts validated 3×3 matrices only, and keeps their declared determinant and inverse bounds. */\nexport function invertSensorColorMatrix(m: Matrix): number[][] {\n  const [a,b,c,d,e,f,g,h,i]=m.flat() as number[], det=a!*(e!*i!-f!*h!)-b!*(d!*i!-f!*g!)+c!*(d!*h!-e!*g!);\n  if (!Number.isFinite(det) || Math.abs(det)<1e-9) throw new InvalidConfigurationError(\"Export color matrix is singular or outside supported conditioning.\");\n  const result=[[e!*i!-f!*h!,c!*h!-b!*i!,b!*f!-c!*e!],[f!*g!-d!*i!,a!*i!-c!*g!,c!*d!-a!*f!],[d!*h!-e!*g!,b!*g!-a!*h!,a!*e!-b!*d!]].map((r) => r.map((v) => v/det));\n  if (result.flat().some((v) => !Number.isFinite(v) || Math.abs(v)>100)) throw new InvalidConfigurationError(\"Export inverse color matrix exceeds supported range.\");\n  return result;\n}\n\n/** Camera coordinates of the declared D65 white (not normalized). The caller decides container policy. */\nexport function sensorColorReferenceWhite(xyzToCamera: Matrix): number[] {\n  const white = resolveCaptureColorModel().referenceWhiteXyz;\n  return multiply(xyzToCamera, [white.x, 1, white.z]);\n}\n\n/** Per-pixel highlight/WB counters accumulated by {@link SensorColorDevelopment.developInto}. */\nexport interface SensorColorDevelopmentCounts { whiteBalanceClippedSampleCount: number; reconstructedPixelCount: number }\n\nexport interface SensorColorDevelopment {\n  referenceWhiteXyz: { x: number; y: 1; z: number };\n  /** Green-normalized camera-neutral coordinates of D65 (exactly [1,1,1] for white-normalized profiles). */\n  cameraNeutralD65: readonly number[];\n  /** Applied channel gains: the resolved WB gains, or 1/cameraNeutralD65 when WB is not required. */\n  appliedGains: readonly number[];\n  /** Green-normalized camera neutral of the assumed scene white: g_green/g (DNG AsShotNeutral). */\n  asShotNeutral: readonly number[];\n  /** X·M·diag(w) rows, applied as two products in the documented order. */\n  cameraToXyz: readonly (readonly number[])[];\n  highlightMode: RawHighlightMode;\n  /** Develops one normalized camera pixel into out[offset..offset+2] (linear sRGB/D65); updates counts. */\n  developInto(c0: number, c1: number, c2: number, out: Float64Array | number[], offset: number, counts: SensorColorDevelopmentCounts): void;\n  /**\n   * #602: the colour step alone, X·M·diag(w)·b, for white-balanced values b that a whole-image highlight step already formed\n   * (`highlightMode: \"reconstruct-neighbors\"`); the same products and order as the end of developInto.\n   */\n  developWhiteBalancedInto(b0: number, b1: number, b2: number, out: Float64Array | number[], offset: number): void;\n  /** Convenience array form of developInto (counts discarded). */\n  develop(values: readonly number[]): number[];\n}\n\n/**\n * Internal development of validated, normalized sensor RGB with a declared D65 profile.\n * `gains` is the resolved WB intent, or null when WB is not required. Both the packed native path and\n * the bounded reference path use the returned operator, so there is one implementation.\n */\nexport function prepareSensorColorDevelopment(cameraChannelsToXyz: Matrix, gains: readonly number[] | null,\n  highlightMode: RawHighlightMode = \"clip\"): SensorColorDevelopment {\n  if (highlightMode !== \"clip\" && highlightMode !== \"reconstruct\") throw new InvalidConfigurationError(\"Unsupported RAW highlight mode.\");\n  const model = resolveCaptureColorModel();\n  const raw = sensorColorReferenceWhite(invertSensorColorMatrix(cameraChannelsToXyz));\n  if (raw.some((x) => !(x > 0) || !Number.isFinite(x))) throw new InvalidConfigurationError(\"Reference white requires positive camera coordinates.\");\n  let neutral = raw.map((x) => x/raw[1]!);\n  if (neutral.every((x) => Math.abs(x-1) <= WHITE_NORMALIZED_TOLERANCE)) neutral = [1, 1, 1];\n  if (gains && (gains.length !== 3 || gains.some((x) => !(x > 0) || !Number.isFinite(x)))) throw new InvalidConfigurationError(\"White-balance gains must be positive and finite.\");\n  const applied = gains ? [...gains] : neutral.map((x) => 1/x);\n  const asShot = applied.map((x) => applied[1]!/x), clipLevel = Math.min(...applied);\n  const m = cameraChannelsToXyz.map((row) => row.map((v, i) => v*neutral[i]!)), x = model.xyzToCameraRgb;\n  const m0 = m[0]!, m1 = m[1]!, m2 = m[2]!, x0 = x[0]!, x1 = x[1]!, x2 = x[2]!;\n  const g0 = applied[0]!, g1 = applied[1]!, g2 = applied[2]!, reconstruct = highlightMode === \"reconstruct\";\n  const start = RAW_HIGHLIGHT_RECONSTRUCTION_START;\n  function developInto(c0: number, c1: number, c2: number, out: Float64Array | number[], offset: number, counts: SensorColorDevelopmentCounts): void {\n    const u0 = c0*g0, u1 = c1*g1, u2 = c2*g2;\n    let b0: number, b1: number, b2: number;\n    const k0 = u0 > clipLevel ? clipLevel : u0, k1 = u1 > clipLevel ? clipLevel : u1, k2 = u2 > clipLevel ? clipLevel : u2;\n    if (!reconstruct) {\n      b0 = k0; b1 = k1; b2 = k2;\n      counts.whiteBalanceClippedSampleCount += Number(u0 > clipLevel)+Number(u1 > clipLevel)+Number(u2 > clipLevel);\n    } else {\n      const s = c0 > c1 ? (c0 > c2 ? c0 : c2) : (c1 > c2 ? c1 : c2);\n      if (s > start) {\n        const t = s >= 1 ? 1 : (s-start)/(1-start), n = k0 > k1 ? (k0 > k2 ? k0 : k2) : (k1 > k2 ? k1 : k2);\n        b0 = (1-t)*u0+t*n; b1 = (1-t)*u1+t*n; b2 = (1-t)*u2+t*n; counts.reconstructedPixelCount++;\n      } else { b0 = u0; b1 = u1; b2 = u2; }\n    }\n    const X = ((0+m0[0]!*b0)+m0[1]!*b1)+m0[2]!*b2, Y = ((0+m1[0]!*b0)+m1[1]!*b1)+m1[2]!*b2, Z = ((0+m2[0]!*b0)+m2[1]!*b1)+m2[2]!*b2;\n    out[offset] = ((0+x0[0]!*X)+x0[1]!*Y)+x0[2]!*Z;\n    out[offset+1] = ((0+x1[0]!*X)+x1[1]!*Y)+x1[2]!*Z;\n    out[offset+2] = ((0+x2[0]!*X)+x2[1]!*Y)+x2[2]!*Z;\n  }\n  function developWhiteBalancedInto(b0: number, b1: number, b2: number, out: Float64Array | number[], offset: number): void {\n    const X = ((0+m0[0]!*b0)+m0[1]!*b1)+m0[2]!*b2, Y = ((0+m1[0]!*b0)+m1[1]!*b1)+m1[2]!*b2, Z = ((0+m2[0]!*b0)+m2[1]!*b1)+m2[2]!*b2;\n    out[offset] = ((0+x0[0]!*X)+x0[1]!*Y)+x0[2]!*Z;\n    out[offset+1] = ((0+x1[0]!*X)+x1[1]!*Y)+x1[2]!*Z;\n    out[offset+2] = ((0+x2[0]!*X)+x2[1]!*Y)+x2[2]!*Z;\n  }\n  return {\n    developWhiteBalancedInto, referenceWhiteXyz: model.referenceWhiteXyz, cameraNeutralD65: Object.freeze(neutral), appliedGains: Object.freeze(applied),\n    asShotNeutral: Object.freeze(asShot), cameraToXyz: m, highlightMode, developInto,\n    develop(values: readonly number[]): number[] {\n      const out = [0, 0, 0]; developInto(values[0]!, values[1]!, values[2]!, out, 0, { whiteBalanceClippedSampleCount: 0, reconstructedPixelCount: 0 }); return out;\n    }\n  };\n}\n"]}

package/dist/color/spectral-spikiness.d.ts

@@ -0,0 +1,11 @@
+/**
+ * #590 (2.22.0): how much a light's spectrum can defeat the facade's spectral midpoint quadrature. The composite midpoint rule
+ * over [a, b] with step h has the error (b − a)·h²·f″(ξ)/24 (DLMF §3.5(i)); normalized by the integral it gives the relative
+ * error indicator R(h) = h²·κ/24 with κ = Σ_λ |S(λ + 1) − 2S(λ) + S(λ − 1)| / Σ_λ S(λ) on a 1 nm grid over 380–780 nm (nm⁻²),
+ * from the illuminant's spectral power distribution alone (no sensor channel). A light is "spiky at h" when R(h) exceeds
+ * {@link SPECTRAL_SPIKINESS_THRESHOLD} (a design choice: the bound allows more than 1 % error).
+ */
+/** R(h) above this marks a spiky light (owner decision 2026-11-10: the fluorescent and LED representative sources count). */
+export declare const SPECTRAL_SPIKINESS_THRESHOLD = 0.01;
+/** κ, nm⁻², of a spectral power distribution evaluated at 380, 381, …, 780 nm (0 for a zero spectrum). */
+export declare function spectralSpikiness(evaluate: (wavelengthNanometers: number) => number): number;

package/dist/color/spectral-spikiness.js

@@ -0,0 +1,23 @@
+// SPDX-License-Identifier: Apache-2.0
+/**
+ * #590 (2.22.0): how much a light's spectrum can defeat the facade's spectral midpoint quadrature. The composite midpoint rule
+ * over [a, b] with step h has the error (b − a)·h²·f″(ξ)/24 (DLMF §3.5(i)); normalized by the integral it gives the relative
+ * error indicator R(h) = h²·κ/24 with κ = Σ_λ |S(λ + 1) − 2S(λ) + S(λ − 1)| / Σ_λ S(λ) on a 1 nm grid over 380–780 nm (nm⁻²),
+ * from the illuminant's spectral power distribution alone (no sensor channel). A light is "spiky at h" when R(h) exceeds
+ * {@link SPECTRAL_SPIKINESS_THRESHOLD} (a design choice: the bound allows more than 1 % error).
+ */
+/** R(h) above this marks a spiky light (owner decision 2026-11-10: the fluorescent and LED representative sources count). */
+export const SPECTRAL_SPIKINESS_THRESHOLD = 0.01;
+/** κ, nm⁻², of a spectral power distribution evaluated at 380, 381, …, 780 nm (0 for a zero spectrum). */
+export function spectralSpikiness(evaluate) {
+    let a = evaluate(380), b = evaluate(381), sum = a + b, second = 0;
+    for (let nm = 382; nm <= 780; nm++) {
+        const c = evaluate(nm);
+        second += Math.abs(c - 2 * b + a);
+        sum += c;
+        a = b;
+        b = c;
+    }
+    return sum > 0 ? second / sum : 0;
+}
+//# sourceMappingURL=spectral-spikiness.js.map

package/dist/color/spectral-spikiness.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"spectral-spikiness.js","sourceRoot":"","sources":["../../src/color/spectral-spikiness.ts"],"names":[],"mappings":"AAAA,sCAAsC;AAEtC;;;;;;GAMG;AAEH,6HAA6H;AAC7H,MAAM,CAAC,MAAM,4BAA4B,GAAG,IAAI,CAAC;AAEjD,0GAA0G;AAC1G,MAAM,UAAU,iBAAiB,CAAC,QAAkD;IAClF,IAAI,CAAC,GAAG,QAAQ,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,CAAC,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,EAAE,MAAM,GAAG,CAAC,CAAC;IAClE,KAAK,IAAI,EAAE,GAAG,GAAG,EAAE,EAAE,IAAI,GAAG,EAAE,EAAE,EAAE,EAAE,CAAC;QAAC,MAAM,CAAC,GAAG,QAAQ,CAAC,EAAE,CAAC,CAAC;QAAC,MAAM,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QAAC,GAAG,IAAI,CAAC,CAAC;QAAC,CAAC,GAAG,CAAC,CAAC;QAAC,CAAC,GAAG,CAAC,CAAC;IAAC,CAAC;IAC1H,OAAO,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,MAAM,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;AACpC,CAAC","sourcesContent":["// SPDX-License-Identifier: Apache-2.0\n\n/**\n * #590 (2.22.0): how much a light's spectrum can defeat the facade's spectral midpoint quadrature. The composite midpoint rule\n * over [a, b] with step h has the error (b − a)·h²·f″(ξ)/24 (DLMF §3.5(i)); normalized by the integral it gives the relative\n * error indicator R(h) = h²·κ/24 with κ = Σ_λ |S(λ + 1) − 2S(λ) + S(λ − 1)| / Σ_λ S(λ) on a 1 nm grid over 380–780 nm (nm⁻²),\n * from the illuminant's spectral power distribution alone (no sensor channel). A light is \"spiky at h\" when R(h) exceeds\n * {@link SPECTRAL_SPIKINESS_THRESHOLD} (a design choice: the bound allows more than 1 % error).\n */\n\n/** R(h) above this marks a spiky light (owner decision 2026-11-10: the fluorescent and LED representative sources count). */\nexport const SPECTRAL_SPIKINESS_THRESHOLD = 0.01;\n\n/** κ, nm⁻², of a spectral power distribution evaluated at 380, 381, …, 780 nm (0 for a zero spectrum). */\nexport function spectralSpikiness(evaluate: (wavelengthNanometers: number) => number): number {\n  let a = evaluate(380), b = evaluate(381), sum = a + b, second = 0;\n  for (let nm = 382; nm <= 780; nm++) { const c = evaluate(nm); second += Math.abs(c - 2 * b + a); sum += c; a = b; b = c; }\n  return sum > 0 ? second / sum : 0;\n}\n"]}

package/dist/core/execution-error.d.ts

@@ -0,0 +1,21 @@
+/**
+ * The image-kernel backends' error class and its reasons (2.22.0: kept outside src/gpu/ so that the CPU's on-demand code loader,
+ * src/facade/on-demand-kits.ts, can report `quality-code-unavailable` without bringing GPU code into a bundle's first load;
+ * re-exported unchanged from src/gpu/image-kernel-execution.ts and the package root).
+ */
+/** Why a requested `webgpu-f32` backend ran the CPU reference instead. `webgpu-disabled`: `setWebGpuEnabled(false)` in this realm. */
+export type ImageKernelFallbackReason = "webgpu-unavailable" | "webgpu-disabled" | "adapter-unavailable" | "device-request-failed" | "insufficient-limits" | "unsupported-platform-endianness";
+/**
+ * Failure classes reported by {@link ImageKernelExecutionError}. `conformance-failed` (#481): a device failed the display view's
+ * run-time self-check with `gpu: "require"`. `quality-code-unavailable` (2.22.0, #594): on-demand full-quality code could not be
+ * loaded (for example a failed chunk download); `prepareQualityPaths()` rejects with it, a display-view render with `gpu: "require"`
+ * fails with it (with `gpu: "auto"` it runs the CPU reference), and a CPU request whose full-quality or opt-in code does not load
+ * fails with it (nothing is remembered: the next request tries again).
+ */
+export type ImageKernelFailureReason = ImageKernelFallbackReason | "device-lost" | "out-of-memory" | "gpu-validation" | "aborted" | "disposed" | "busy" | "unsupported-on-backend" | "unsupported-float32-domain" | "indeterminate-float32-classification" | "conformance-failed" | "quality-code-unavailable";
+/** Domain error for getting and running a backend. Input validation still throws the engine's existing errors. */
+export declare class ImageKernelExecutionError extends Error {
+    readonly code = "IMAGE_KERNEL_EXECUTION";
+    readonly reason: ImageKernelFailureReason;
+    constructor(reason: ImageKernelFailureReason, message: string);
+}

package/dist/core/execution-error.js

@@ -0,0 +1,12 @@
+// SPDX-License-Identifier: Apache-2.0
+/** Domain error for getting and running a backend. Input validation still throws the engine's existing errors. */
+export class ImageKernelExecutionError extends Error {
+    code = "IMAGE_KERNEL_EXECUTION";
+    reason;
+    constructor(reason, message) {
+        super(message);
+        this.name = "ImageKernelExecutionError";
+        this.reason = reason;
+    }
+}
+//# sourceMappingURL=execution-error.js.map

package/dist/core/execution-error.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"execution-error.js","sourceRoot":"","sources":["../../src/core/execution-error.ts"],"names":[],"mappings":"AAAA,sCAAsC;AAuBtC,kHAAkH;AAClH,MAAM,OAAO,yBAA0B,SAAQ,KAAK;IACzC,IAAI,GAAG,wBAAwB,CAAC;IAChC,MAAM,CAA2B;IAE1C,YAAY,MAAgC,EAAE,OAAe;QAC3D,KAAK,CAAC,OAAO,CAAC,CAAC;QACf,IAAI,CAAC,IAAI,GAAG,2BAA2B,CAAC;QACxC,IAAI,CAAC,MAAM,GAAG,MAAM,CAAC;IACvB,CAAC;CACF","sourcesContent":["// SPDX-License-Identifier: Apache-2.0\n\n/**\n * The image-kernel backends' error class and its reasons (2.22.0: kept outside src/gpu/ so that the CPU's on-demand code loader,\n * src/facade/on-demand-kits.ts, can report `quality-code-unavailable` without bringing GPU code into a bundle's first load;\n * re-exported unchanged from src/gpu/image-kernel-execution.ts and the package root).\n */\n\n/** Why a requested `webgpu-f32` backend ran the CPU reference instead. `webgpu-disabled`: `setWebGpuEnabled(false)` in this realm. */\nexport type ImageKernelFallbackReason = \"webgpu-unavailable\" | \"webgpu-disabled\" | \"adapter-unavailable\" | \"device-request-failed\" |\n  \"insufficient-limits\" | \"unsupported-platform-endianness\";\n\n/**\n * Failure classes reported by {@link ImageKernelExecutionError}. `conformance-failed` (#481): a device failed the display view's\n * run-time self-check with `gpu: \"require\"`. `quality-code-unavailable` (2.22.0, #594): on-demand full-quality code could not be\n * loaded (for example a failed chunk download); `prepareQualityPaths()` rejects with it, a display-view render with `gpu: \"require\"`\n * fails with it (with `gpu: \"auto\"` it runs the CPU reference), and a CPU request whose full-quality or opt-in code does not load\n * fails with it (nothing is remembered: the next request tries again).\n */\nexport type ImageKernelFailureReason = ImageKernelFallbackReason | \"device-lost\" | \"out-of-memory\" | \"gpu-validation\" |\n  \"aborted\" | \"disposed\" | \"busy\" | \"unsupported-on-backend\" | \"unsupported-float32-domain\" | \"indeterminate-float32-classification\" | \"conformance-failed\" |\n  \"quality-code-unavailable\";\n\n/** Domain error for getting and running a backend. Input validation still throws the engine's existing errors. */\nexport class ImageKernelExecutionError extends Error {\n  readonly code = \"IMAGE_KERNEL_EXECUTION\";\n  readonly reason: ImageKernelFailureReason;\n\n  constructor(reason: ImageKernelFailureReason, message: string) {\n    super(message);\n    this.name = \"ImageKernelExecutionError\";\n    this.reason = reason;\n  }\n}\n\n"]}

package/dist/core/version.d.ts

@@ -6,4 +6,4 @@
  * own versions. Release checks confirm it matches package.json and the packed package.
  * @see docs/RELEASE_1_0_1.md
  */
-export declare const ENGINE_API_VERSION: "2.21.0";
+export declare const ENGINE_API_VERSION: "2.23.0";

package/dist/core/version.js

@@ -7,5 +7,5 @@
  * own versions. Release checks confirm it matches package.json and the packed package.
  * @see docs/RELEASE_1_0_1.md
  */
-export const ENGINE_API_VERSION = "2.21.0";
+export const ENGINE_API_VERSION = "2.23.0";
 //# sourceMappingURL=version.js.map

package/dist/core/version.js.map

@@ -1 +1 @@
-{"version":3,"file":"version.js","sourceRoot":"","sources":["../../src/core/version.ts"],"names":[],"mappings":"AAAA,sCAAsC;AAEtC;;;;;;;GAOG;AACH,MAAM,CAAC,MAAM,kBAAkB,GAAG,QAAiB,CAAC","sourcesContent":["// SPDX-License-Identifier: Apache-2.0\n\n/**\n * Version of the engine release. Since 1.0.1 it equals the npm package version.\n *\n * New captures and plans record this version as their creator. Archived records keep the\n * version they were made with. Proof-of-concept (POC), schema and model contracts have their\n * own versions. Release checks confirm it matches package.json and the packed package.\n * @see docs/RELEASE_1_0_1.md\n */\nexport const ENGINE_API_VERSION = \"2.21.0\" as const;\n"]}
+{"version":3,"file":"version.js","sourceRoot":"","sources":["../../src/core/version.ts"],"names":[],"mappings":"AAAA,sCAAsC;AAEtC;;;;;;;GAOG;AACH,MAAM,CAAC,MAAM,kBAAkB,GAAG,QAAiB,CAAC","sourcesContent":["// SPDX-License-Identifier: Apache-2.0\n\n/**\n * Version of the engine release. Since 1.0.1 it equals the npm package version.\n *\n * New captures and plans record this version as their creator. Archived records keep the\n * version they were made with. Proof-of-concept (POC), schema and model contracts have their\n * own versions. Release checks confirm it matches package.json and the packed package.\n * @see docs/RELEASE_1_0_1.md\n */\nexport const ENGINE_API_VERSION = \"2.23.0\" as const;\n"]}

package/dist/facade/catalogue-lens-signal.d.ts

@@ -3,6 +3,28 @@ import { type ResolveLensSampledPsfInput } from "../optics/lens-psf-profile.js";
 import type { DisplayViewScene } from "./expected-signal.js";
 /** Fewest grid cells between neighboring profile field nodes (10 since #589: with 4, 0.4 % of quick frames failed the run-time check). */
 export declare const CATALOGUE_FAST_GRID_MINIMUM_SUBDIVISIONS = 10;
+/**
+ * #603 (2.22.0): the fewest cells on a fast request (on screen: 2.20.0's grid, its run-time fallbacks to the exact method included);
+ * the full grid, {@link CATALOGUE_FAST_GRID_MINIMUM_SUBDIVISIONS}, on an exact one (the exports and the calculator).
+ */
+export declare const CATALOGUE_FAST_GRID_SPEED_SUBDIVISIONS = 4;
+/**
+ * #603 (owner decision 2026-11-11): the lenses that keep the full grid on fast requests too, because the 4-cell grid's run-time
+ * fallbacks to the exact method made their fast frames slower on average. Measured November 11, 2026 (arm64, Node.js 24; every
+ * catalogue lens but the phone lenses, wide open and at max(f/8, widest), Siemens star and color patches, whole frames 192, 384 and
+ * 768 pixels wide, two runs each; docs/LENS_CATALOGUE.md, "Fast method"): only this lens fell back, 4 of its 8 frames at 192 pixels
+ * (f/2.8), 11.6 s each, against 0.17 s for the 10-cell grid.
+ */
+export declare const CATALOGUE_FAST_GRID_FULL_ON_SCREEN: readonly {
+    lensId: string;
+    fallbacksAtFourCells: {
+        frames: number;
+        fellBack: number;
+        widthPixels: number;
+    };
+}[];
+/** #603: the fewest grid cells for a lens on a fast or an exact request. */
+export declare function catalogueFastGridSubdivisions(lensId: string, fast: boolean): number;
 /** Largest move of a node kernel's intensity centroid across one grid cell, in sites (sets finer grids where lateral color moves the PSF). */
 export declare const CATALOGUE_FAST_GRID_CENTROID_STEP_SITES = 0.05;
 /**
@@ -44,6 +66,8 @@ export interface CatalogueLensSignalInput {
     /** The planar provider's plane distance and frame-unit scale. */
     distanceM: number;
     metresPerFrameUnit: number;
+    /** #603: the fewest grid cells between neighbouring field nodes (default {@link CATALOGUE_FAST_GRID_MINIMUM_SUBDIVISIONS}). */
+    gridSubdivisions?: number;
     /** Internal test hook: resolve every site's kernels exactly (the exact method's sum, slowly) instead of blending the grid. */
     exactKernels?: boolean;
 }
@@ -66,7 +90,7 @@ export interface CatalogueLensSignal {
  * ⌈Δc / GRID_CENTROID_STEP_SITES⌉) equal cells, Δc the largest centroid move between the interval's two nodes over the other
  * axis's nodes and the wavelength nodes.
  */
-export declare function catalogueFastGridAxes(psf: CatalogueLensSignalInput["psf"]): {
+export declare function catalogueFastGridAxes(psf: CatalogueLensSignalInput["psf"], minimum?: number): {
     x: number[];
     y: number[];
 };
@@ -87,6 +111,7 @@ export declare function catalogueFastAssumptions(run: {
     rows: number;
     resolved: number;
     verificationRelativeDifference: number;
+    subdivisions: number;
 }): {
     id: string;
     value: Record<string, unknown>;