@photivra/engine 2.18.0 → 2.20.0

package/dist/processing/development-detail.d.ts

@@ -1,3 +1,4 @@
+import { type SdrPixelCounts } from "../output/native-capture-sdr.js";
 import { type SdrRenderingProfile } from "../output/sdr-rendering.js";
 /** Version of the `development-detail` model (diffraction correction and unsharp mask). */
 export declare const DEVELOPMENT_DETAIL_MODEL_VERSION: "0.1.0";
@@ -86,6 +87,15 @@ export declare function applyDiffractionCorrection(rgb: Float64Array, width: num
     samples: Float64Array;
     diagnostics: DiffractionCorrectionDiagnostics;
 };
+/** Runs a row-by-row detail stage to its end and returns its result. */
+export declare function drainDetailRows<T>(rows: Generator<number, T, void>): T;
+/**
+ * `applyDiffractionCorrection` in place, one output row per step (#571); each step yields the number of rows done. Output row y
+ * reads input rows y − R … y + R, so it waits in a ring of R + 1 rows and replaces input row y once row y + R is done. Every sample
+ * is the same sum in the same order as in a separate output image, so the result is bit-identical, with one image in memory
+ * instead of two. `rgb` must hold width × height pixels.
+ */
+export declare function correctDiffractionInPlace(rgb: Float64Array, width: number, height: number, filter: Pick<DiffractionCorrectionFilter, "weights" | "settings">): Generator<number, DiffractionCorrectionDiagnostics, void>;
 /** Normalized Gaussian weights exp(−i²/(2σ²)), i = −r..r with r = ⌈3σ⌉. */
 export declare function unsharpMaskKernel(radiusPixels: number): Float64Array;
 /** What the unsharp mask did. Counts are samples (one per channel). */
@@ -112,6 +122,15 @@ export declare function applyUnsharpMask(encoded: Float64Array, width: number, h
     samples: Float64Array;
     diagnostics: UnsharpMaskDiagnostics;
 };
+/**
+ * The unsharp mask one output row per step (#571); each step yields the number of rows done. The horizontal Gaussian pass of
+ * each input row is kept in a ring of 2r + 1 rows, and `prepareRow(y)`, when given, runs once on input row y before it is
+ * first read (`renderSdrWithUnsharpMaskInPlace` encodes the row there). Output row y goes to `out`: the clipped value, or
+ * Math.round(value · scale) with a `scale`. Every sum is the same, in the same order, as in whole-image passes, so the
+ * result is bit-identical, with 2r + 1 rows of working memory instead of three images. `settings` must be parsed and the
+ * raster sizes must match.
+ */
+export declare function unsharpMaskRows(encoded: Float64Array, width: number, height: number, settings: UnsharpMaskSettings, out: Float64Array | Uint8Array | Uint16Array, scale?: number, prepareRow?: (y: number) => void): Generator<number, UnsharpMaskDiagnostics, void>;
 /** SDR counters in the shape of the SDR rendering diagnostics. */
 export interface DetailSdrCounts {
     negative: number;
@@ -130,3 +149,13 @@ export declare function renderSdrWithUnsharpMask(rgb: Float64Array, width: numbe
     counts: DetailSdrCounts;
     sharpening: UnsharpMaskDiagnostics;
 };
+/**
+ * `renderSdrWithUnsharpMask` one output row per step (#571), encoding each row of `rgb` in place just before the mask first reads
+ * it (`rgb` then holds the encoded values). Bit-identical to the whole-image function, with no second image in memory. `settings`
+ * must be parsed and `rgb` must hold width × height pixels.
+ */
+export declare function renderSdrWithUnsharpMaskInPlace(rgb: Float64Array, width: number, height: number, profile: SdrRenderingProfile, settings: UnsharpMaskSettings): Generator<number, {
+    codes: Uint8Array | Uint16Array;
+    counts: SdrPixelCounts;
+    sharpening: UnsharpMaskDiagnostics;
+}, void>;

package/dist/processing/development-detail.js

@@ -12,12 +12,13 @@
  * Development order with these stages: (exposure gain) → demosaic → colour and white balance → diffraction correction
  * (linear) → rendering exposure and tone curve → gamut handling → sRGB transfer → unsharp mask → quantization (encode).
  *
- * Image edges are handled by clamping (edge pixels repeated). Both stages work on a whole in-memory image.
+ * Image edges are handled by clamping (edge pixels repeated). Both stages need the whole image; since #571 they also run row by
+ * row in place (`correctDiffractionInPlace`, `renderSdrWithUnsharpMaskInPlace`), with the same sums in the same order.
  * @see docs/DIGITAL_OPTICS_FOUNDATION.md for the equations, sources, trade-offs and limits.
  */
-import { encodeSrgbComponent } from "../color/srgb-transfer.js";
 import { InvalidConfigurationError } from "../core/configuration-error.js";
 import { InvalidScientificResultError } from "../core/validation.js";
+import { renderSdrRunInto } from "../output/native-capture-sdr.js";
 import { prepareSdrPixelTone } from "../output/sdr-rendering.js";
 /** Version of the `development-detail` model (diffraction correction and unsharp mask). */
 export const DEVELOPMENT_DETAIL_MODEL_VERSION = "0.1.0";
@@ -158,38 +159,62 @@ export function evaluateDiffractionCorrectionResponse(weights, radius, fx, fy) {
 export function applyDiffractionCorrection(rgb, width, height, filter) {
     if (rgb.length !== width * height * 3)
         throw new InvalidConfigurationError("Diffraction correction raster size mismatch.");
-    const R = filter.settings.kernelRadiusPixels, size = 2 * R + 1, w = filter.weights, out = new Float64Array(rgb.length);
+    const samples = Float64Array.from(rgb);
+    return { samples, diagnostics: drainDetailRows(correctDiffractionInPlace(samples, width, height, filter)) };
+}
+/** Runs a row-by-row detail stage to its end and returns its result. */
+export function drainDetailRows(rows) {
+    let step = rows.next();
+    while (!step.done)
+        step = rows.next();
+    return step.value;
+}
+/**
+ * `applyDiffractionCorrection` in place, one output row per step (#571); each step yields the number of rows done. Output row y
+ * reads input rows y − R … y + R, so it waits in a ring of R + 1 rows and replaces input row y once row y + R is done. Every sample
+ * is the same sum in the same order as in a separate output image, so the result is bit-identical, with one image in memory
+ * instead of two. `rgb` must hold width × height pixels.
+ */
+export function* correctDiffractionInPlace(rgb, width, height, filter) {
+    const R = filter.settings.kernelRadiusPixels, size = 2 * R + 1, w = filter.weights, rowLength = width * 3, ring = new Float64Array((R + 1) * rowLength);
     const xs = new Int32Array(width + 2 * R);
     for (let i = 0; i < xs.length; i++)
         xs[i] = Math.min(width - 1, Math.max(0, i - R)) * 3;
     let negative = 0, above = 0;
-    for (let y = 0; y < height; y++) {
-        const o = y * width * 3;
+    for (let y = 0; y <= height + R; y++) {
+        // Output row y − R − 1 holds the slot of row y; input row y − R − 1 is read by no later output row, so it takes its place.
+        const o = (y % (R + 1)) * rowLength;
+        if (y > R)
+            rgb.set(ring.subarray(o, o + rowLength), (y - R - 1) * rowLength);
+        if (y >= height)
+            continue;
+        ring.fill(0, o, o + rowLength);
         for (let dy = -R; dy <= R; dy++) {
-            const row = Math.min(height - 1, Math.max(0, y + dy)) * width * 3, wr = (dy + R) * size;
+            const row = Math.min(height - 1, Math.max(0, y + dy)) * rowLength, wr = (dy + R) * size;
             for (let dx = 0; dx < size; dx++) {
                 const weight = w[wr + dx];
                 if (weight === 0)
                     continue;
                 for (let x = 0; x < width; x++) {
                     const s = row + xs[x + dx], d = o + 3 * x;
-                    out[d] = out[d] + weight * rgb[s];
-                    out[d + 1] = out[d + 1] + weight * rgb[s + 1];
-                    out[d + 2] = out[d + 2] + weight * rgb[s + 2];
+                    ring[d] = ring[d] + weight * rgb[s];
+                    ring[d + 1] = ring[d + 1] + weight * rgb[s + 1];
+                    ring[d + 2] = ring[d + 2] + weight * rgb[s + 2];
                 }
             }
         }
+        for (let i = o; i < o + rowLength; i++) {
+            const v = ring[i];
+            if (!Number.isFinite(v))
+                throw new InvalidScientificResultError("Diffraction correction produced a non-finite sample.");
+            if (v < 0)
+                negative++;
+            else if (v > 1)
+                above++;
+        }
+        yield y + 1;
     }
-    for (let i = 0; i < out.length; i++) {
-        const v = out[i];
-        if (!Number.isFinite(v))
-            throw new InvalidScientificResultError("Diffraction correction produced a non-finite sample.");
-        if (v < 0)
-            negative++;
-        else if (v > 1)
-            above++;
-    }
-    return { samples: out, diagnostics: { negativeSampleCount: negative, aboveWhiteSampleCount: above } };
+    return { negativeSampleCount: negative, aboveWhiteSampleCount: above };
 }
 /** Normalized Gaussian weights exp(−i²/(2σ²)), i = −r..r with r = ⌈3σ⌉. */
 export function unsharpMaskKernel(radiusPixels) {
@@ -213,64 +238,79 @@ export function applyUnsharpMask(encoded, width, height, settings) {
     const s = parseUnsharpMaskSettings(settings);
     if (encoded.length !== width * height * 3)
         throw new InvalidConfigurationError("Unsharp mask raster size mismatch.");
-    const k = unsharpMaskKernel(s.radiusPixels), r = (k.length - 1) / 2, horizontal = new Float64Array(encoded.length), blurred = new Float64Array(encoded.length);
-    for (let y = 0; y < height; y++)
-        for (let x = 0; x < width; x++) {
-            const d = (y * width + x) * 3;
-            let a = 0, b = 0, c = 0;
-            for (let i = -r; i <= r; i++) {
-                const sx = (y * width + Math.min(width - 1, Math.max(0, x + i))) * 3, w = k[i + r];
-                a += w * encoded[sx];
-                b += w * encoded[sx + 1];
-                c += w * encoded[sx + 2];
+    const samples = new Float64Array(encoded.length);
+    return { samples, diagnostics: drainDetailRows(unsharpMaskRows(encoded, width, height, s, samples)) };
+}
+/**
+ * The unsharp mask one output row per step (#571); each step yields the number of rows done. The horizontal Gaussian pass of
+ * each input row is kept in a ring of 2r + 1 rows, and `prepareRow(y)`, when given, runs once on input row y before it is
+ * first read (`renderSdrWithUnsharpMaskInPlace` encodes the row there). Output row y goes to `out`: the clipped value, or
+ * Math.round(value · scale) with a `scale`. Every sum is the same, in the same order, as in whole-image passes, so the
+ * result is bit-identical, with 2r + 1 rows of working memory instead of three images. `settings` must be parsed and the
+ * raster sizes must match.
+ */
+export function* unsharpMaskRows(encoded, width, height, settings, out, scale, prepareRow) {
+    const k = unsharpMaskKernel(settings.radiusPixels), r = (k.length - 1) / 2, span = 2 * r + 1, rowLength = width * 3, ring = new Float64Array(span * rowLength);
+    const amount = settings.amount, threshold = settings.threshold;
+    let sharpened = 0, held = 0, low = 0, high = 0, increase = 0, decrease = 0, ready = 0;
+    for (let y = 0; y < height; y++) {
+        // Rows up to y + r, prepared and horizontally blurred.
+        for (const last = Math.min(height - 1, y + r); ready <= last; ready++) {
+            prepareRow?.(ready);
+            const source = ready * rowLength, slot = (ready % span) * rowLength;
+            for (let x = 0; x < width; x++) {
+                let a = 0, b = 0, c = 0;
+                for (let i = -r; i <= r; i++) {
+                    const sx = source + Math.min(width - 1, Math.max(0, x + i)) * 3, w = k[i + r];
+                    a += w * encoded[sx];
+                    b += w * encoded[sx + 1];
+                    c += w * encoded[sx + 2];
+                }
+                ring[slot + 3 * x] = a;
+                ring[slot + 3 * x + 1] = b;
+                ring[slot + 3 * x + 2] = c;
             }
-            horizontal[d] = a;
-            horizontal[d + 1] = b;
-            horizontal[d + 2] = c;
         }
-    for (let y = 0; y < height; y++)
+        const o = y * rowLength;
         for (let x = 0; x < width; x++) {
-            const d = (y * width + x) * 3;
             let a = 0, b = 0, c = 0;
             for (let i = -r; i <= r; i++) {
-                const sy = (Math.min(height - 1, Math.max(0, y + i)) * width + x) * 3, w = k[i + r];
-                a += w * horizontal[sy];
-                b += w * horizontal[sy + 1];
-                c += w * horizontal[sy + 2];
+                const sy = (Math.min(height - 1, Math.max(0, y + i)) % span) * rowLength + 3 * x, w = k[i + r];
+                a += w * ring[sy];
+                b += w * ring[sy + 1];
+                c += w * ring[sy + 2];
+            }
+            for (let channel = 0; channel < 3; channel++) {
+                const index = o + 3 * x + channel, value = encoded[index], detail = value - (channel === 0 ? a : channel === 1 ? b : c);
+                let v = value;
+                if (detail === 0 || Math.abs(detail) < threshold) {
+                    if (detail !== 0)
+                        held++;
+                }
+                else {
+                    const change = amount * detail;
+                    v = value + change;
+                    sharpened++;
+                    if (change > increase)
+                        increase = change;
+                    if (-change > decrease)
+                        decrease = -change;
+                    if (v < 0) {
+                        v = 0;
+                        low++;
+                    }
+                    else if (v > 1) {
+                        v = 1;
+                        high++;
+                    }
+                }
+                out[index] = scale === undefined ? v : Math.round(v * scale);
             }
-            blurred[d] = a;
-            blurred[d + 1] = b;
-            blurred[d + 2] = c;
-        }
-    const out = new Float64Array(encoded.length);
-    let sharpened = 0, held = 0, low = 0, high = 0, increase = 0, decrease = 0;
-    for (let i = 0; i < encoded.length; i++) {
-        const x = encoded[i], detail = x - blurred[i];
-        if (detail === 0 || Math.abs(detail) < s.threshold) {
-            if (detail !== 0)
-                held++;
-            out[i] = x;
-            continue;
-        }
-        const change = s.amount * detail;
-        let v = x + change;
-        sharpened++;
-        if (change > increase)
-            increase = change;
-        if (-change > decrease)
-            decrease = -change;
-        if (v < 0) {
-            v = 0;
-            low++;
-        }
-        else if (v > 1) {
-            v = 1;
-            high++;
         }
-        out[i] = v;
+        yield y + 1;
     }
-    return { samples: out, diagnostics: { kernelRadiusPixels: r, sharpenedSampleCount: sharpened, thresholdHeldSampleCount: held, clippedLowSampleCount: low,
-            clippedHighSampleCount: high, maximumIncrease: increase, maximumDecrease: decrease } };
+    return { kernelRadiusPixels: r, sharpenedSampleCount: sharpened, thresholdHeldSampleCount: held, clippedLowSampleCount: low,
+        clippedHighSampleCount: high, maximumIncrease: increase, maximumDecrease: decrease };
 }
 /**
  * SDR rendering with an unsharp mask between the sRGB transfer and quantization. Per pixel it runs the reference renderer's
@@ -280,46 +320,27 @@ export function applyUnsharpMask(encoded, width, height, settings) {
 export function renderSdrWithUnsharpMask(rgb, width, height, profile, settings) {
     if (rgb.length !== width * height * 3)
         throw new InvalidConfigurationError("SDR raster size mismatch.");
-    const scale = 2 ** profile.renderingExposureEv, tone = prepareSdrPixelTone(profile), scratch = new Float64Array(3), encoded = new Float64Array(rgb.length);
+    const { codes, counts: { negative, above, changed, low, high }, sharpening } = drainDetailRows(renderSdrWithUnsharpMaskInPlace(Float64Array.from(rgb), width, height, profile, parseUnsharpMaskSettings(settings)));
+    return { codes, counts: { negative, above, changed, low, high }, sharpening };
+}
+/**
+ * `renderSdrWithUnsharpMask` one output row per step (#571), encoding each row of `rgb` in place just before the mask first reads
+ * it (`rgb` then holds the encoded values). Bit-identical to the whole-image function, with no second image in memory. `settings`
+ * must be parsed and `rgb` must hold width × height pixels.
+ */
+export function* renderSdrWithUnsharpMaskInPlace(rgb, width, height, profile, settings) {
+    const scale = 2 ** profile.renderingExposureEv, tone = prepareSdrPixelTone(profile), scratch = new Float64Array(3);
     if (!Number.isFinite(scale) || scale <= 0)
         throw new InvalidConfigurationError("Rendering scale is unrepresentable.");
-    const reject = profile.gamutHandling === "reject-out-of-range", counts = { negative: 0, above: 0, changed: 0, low: 0, high: 0 };
-    for (let i = 0; i < rgb.length; i += 3) {
-        const r = rgb[i], g = rgb[i + 1], b = rgb[i + 2];
-        const er = (r === 0 ? 0 : r) * scale, eg = (g === 0 ? 0 : g) * scale, eb = (b === 0 ? 0 : b) * scale;
-        if (!Number.isFinite(er) || !Number.isFinite(eg) || !Number.isFinite(eb))
-            throw new InvalidScientificResultError("Rendering exposure produced a non-finite sample.");
-        tone(er, eg, eb, scratch);
-        for (let c = 0; c < 3; c++) {
-            const exposed = c === 0 ? er : c === 1 ? eg : eb;
-            if (exposed < 0)
-                counts.negative++;
-            if (exposed > 1)
-                counts.above++;
-            let mapped = scratch[c];
-            if (mapped !== exposed)
-                counts.changed++;
-            if (mapped === 0)
-                mapped = 0;
-            if (mapped < 0 || mapped > 1) {
-                if (reject)
-                    throw new InvalidConfigurationError("Tone output exceeds the declared SDR gamut range.");
-                if (mapped < 0) {
-                    counts.low++;
-                    mapped = 0;
-                }
-                else {
-                    counts.high++;
-                    mapped = 1;
-                }
-            }
-            encoded[i + c] = encodeSrgbComponent(mapped);
-        }
-    }
-    const sharpened = applyUnsharpMask(encoded, width, height, settings), maximum = 2 ** profile.bitDepth - 1;
+    const counts = { negative: 0, above: 0, changed: 0, low: 0, high: 0, outOfGamut: false };
+    // The native SDR kernel per row, unquantized (maximum 0), in place: the reference renderer's operations in its order.
+    const encodeRow = (y) => {
+        renderSdrRunInto(rgb, y * width * 3, width, profile, tone, scale, 0, rgb, y * width * 3, scratch, counts);
+        if (counts.outOfGamut)
+            throw new InvalidConfigurationError("Tone output exceeds the declared SDR gamut range.");
+    };
     const codes = profile.bitDepth === 8 ? new Uint8Array(rgb.length) : new Uint16Array(rgb.length);
-    for (let i = 0; i < codes.length; i++)
-        codes[i] = Math.round(sharpened.samples[i] * maximum);
-    return { codes, counts, sharpening: sharpened.diagnostics };
+    const sharpening = yield* unsharpMaskRows(rgb, width, height, settings, codes, 2 ** profile.bitDepth - 1, encodeRow);
+    return { codes, counts, sharpening };
 }
 //# sourceMappingURL=development-detail.js.map

package/dist/processing/development-detail.js.map

@@ -1 +1 @@
-{"version":3,"file":"development-detail.js","sourceRoot":"","sources":["../../src/processing/development-detail.ts"],"names":[],"mappings":"AAAA,sCAAsC;AAEtC;;;;;;;;;;;;;;;GAeG;AAEH,OAAO,EAAE,mBAAmB,EAAE,MAAM,2BAA2B,CAAC;AAChE,OAAO,EAAE,yBAAyB,EAAE,MAAM,gCAAgC,CAAC;AAC3E,OAAO,EAAE,4BAA4B,EAAE,MAAM,uBAAuB,CAAC;AACrE,OAAO,EAAE,mBAAmB,EAA4B,MAAM,4BAA4B,CAAC;AAE3F,2FAA2F;AAC3F,MAAM,CAAC,MAAM,gCAAgC,GAAG,OAAgB,CAAC;AAOjE,6GAA6G;AAC7G,MAAM,CAAC,MAAM,yBAAyB,GAAG,eAAe,CAAC,MAAM,CAAC,MAAM,CAAC;IACrE,0BAA0B,EAAE,CAAC,EAAE,oBAAoB,EAAE,CAAC;IACtD,mBAAmB,EAAE,GAAG,EAAE,mBAAmB,EAAE,GAAG;IAClD,qBAAqB,EAAE,IAAI,EAAE,qBAAqB,EAAE,CAAC;IACrD,WAAW,EAAE,EAAE,EAAE,yBAAyB,EAAE,EAAE;CAC/C,CAAC,CAAC;AAEH,yFAAyF;AACzF,MAAM,4BAA4B,GAAG,kBAAkB,GAAG,IAAI,CAAC,EAAE,CAAC;AAElE,SAAS,QAAQ,CAAC,KAAc,EAAE,IAAY,EAAE,OAAe,EAAE,OAAe,EAAE,gBAAgB,GAAG,KAAK;IACxG,IAAI,OAAO,KAAK,KAAK,QAAQ,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,IAAI,CAAC,gBAAgB,CAAC,CAAC,CAAC,KAAK,IAAI,OAAO,CAAC,CAAC,CAAC,KAAK,GAAG,OAAO,CAAC,IAAI,KAAK,GAAG,OAAO,EAAE,CAAC;QACvI,MAAM,IAAI,yBAAyB,CAAC,IAAI,GAAG,8BAA8B,GAAG,CAAC,gBAAgB,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,OAAO,GAAG,IAAI,GAAG,OAAO,GAAG,IAAI,CAAC,CAAC;IAChJ,CAAC;IACD,OAAO,KAAK,CAAC;AACf,CAAC;AAED,uCAAuC;AACvC,MAAM,UAAU,wBAAwB,CAAC,KAA0B,EAAE,IAAI,GAAG,YAAY;IACtF,OAAO,EAAE,YAAY,EAAE,QAAQ,CAAC,KAAK,CAAC,YAAY,EAAE,IAAI,GAAG,eAAe,EAAE,CAAC,EAAE,yBAAyB,CAAC,0BAA0B,EAAE,IAAI,CAAC;QACxI,MAAM,EAAE,QAAQ,CAAC,KAAK,CAAC,MAAM,EAAE,IAAI,GAAG,SAAS,EAAE,CAAC,EAAE,yBAAyB,CAAC,oBAAoB,EAAE,IAAI,CAAC;QACzG,SAAS,EAAE,QAAQ,CAAC,KAAK,CAAC,SAAS,EAAE,IAAI,GAAG,YAAY,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC;AACtE,CAAC;AACD,iDAAiD;AACjD,MAAM,UAAU,kCAAkC,CAAC,KAAoC,EAAE,IAAI,GAAG,uBAAuB;IACrH,MAAM,MAAM,GAAG,KAAK,CAAC,kBAAkB,CAAC;IACxC,IAAI,OAAO,MAAM,KAAK,QAAQ,IAAI,CAAC,MAAM,CAAC,aAAa,CAAC,MAAM,CAAC,IAAI,MAAM,GAAG,CAAC,IAAI,MAAM,GAAG,yBAAyB,CAAC,yBAAyB,EAAE,CAAC;QAC9I,MAAM,IAAI,yBAAyB,CAAC,IAAI,GAAG,mDAAmD,GAAG,yBAAyB,CAAC,yBAAyB,GAAG,GAAG,CAAC,CAAC;IAC9J,CAAC;IACD,OAAO,EAAE,YAAY,EAAE,QAAQ,CAAC,KAAK,CAAC,YAAY,EAAE,IAAI,GAAG,eAAe,EAAE,yBAAyB,CAAC,mBAAmB,EAAE,yBAAyB,CAAC,mBAAmB,CAAC;QACvK,cAAc,EAAE,QAAQ,CAAC,KAAK,CAAC,cAAc,EAAE,IAAI,GAAG,iBAAiB,EAAE,yBAAyB,CAAC,qBAAqB,EAAE,yBAAyB,CAAC,qBAAqB,CAAC;QAC1K,WAAW,EAAE,QAAQ,CAAC,KAAK,CAAC,WAAW,EAAE,IAAI,GAAG,cAAc,EAAE,CAAC,EAAE,yBAAyB,CAAC,WAAW,CAAC,EAAE,kBAAkB,EAAE,MAAM,EAAE,CAAC;AAC5I,CAAC;AACD,iHAAiH;AACjH,MAAM,UAAU,wCAAwC,CAAC,cAAsB,EAAE,YAAoB,EAAE,qBAA6B;IAClI,MAAM,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,4BAA4B,GAAG,YAAY,GAAG,IAAI,GAAG,cAAc,GAAG,qBAAqB,CAAC,CAAC;IACrH,OAAO,IAAI,CAAC,GAAG,CAAC,yBAAyB,CAAC,yBAAyB,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC;AAC3F,CAAC;AAED;;;GAGG;AACH,MAAM,UAAU,2BAA2B,CAAC,mBAA2B;IACrE,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,mBAAmB,CAAC,CAAC;IACxC,IAAI,CAAC,IAAI,CAAC;QAAE,OAAO,CAAC,CAAC;IACrB,OAAO,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;AACjE,CAAC;AAyBD;;;;;;GAMG;AACH,MAAM,UAAU,iCAAiC,CAAC,KAAyG;IACzJ,MAAM,QAAQ,GAAG,kCAAkC,CAAC,KAAK,CAAC,QAAQ,CAAC,CAAC;IACpE,MAAM,cAAc,GAAG,QAAQ,CAAC,KAAK,CAAC,cAAc,EAAE,gBAAgB,EAAE,CAAC,EAAE,GAAG,EAAE,IAAI,CAAC,EAAE,KAAK,GAAG,QAAQ,CAAC,KAAK,CAAC,qBAAqB,EAAE,uBAAuB,EAAE,CAAC,EAAE,GAAG,EAAE,IAAI,CAAC,CAAC;IAC5K,MAAM,MAAM,GAAG,KAAK,GAAG,CAAC,QAAQ,CAAC,YAAY,GAAG,IAAI,GAAG,cAAc,CAAC,EAAE,KAAK,GAAG,QAAQ,CAAC,cAAc,EAAE,KAAK,GAAG,QAAQ,CAAC,WAAW,CAAC;IACtI,MAAM,MAAM,GAAG,CAAC,GAAW,EAAU,EAAE;QACrC,MAAM,CAAC,GAAG,2BAA2B,CAAC,GAAG,GAAG,MAAM,CAAC,CAAC;QACpD,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;IAC1E,CAAC,CAAC;IACF,MAAM,CAAC,GAAG,QAAQ,CAAC,kBAAkB,EAAE,IAAI,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IACxD,IAAI,CAAC,GAAG,EAAE,CAAC;IACX,OAAO,CAAC,GAAG,CAAC,GAAG,IAAI;QAAE,CAAC,IAAI,CAAC,CAAC;IAC5B,kIAAkI;IAClI,MAAM,QAAQ,GAAG,IAAI,YAAY,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IAC/C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;QAAE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE;YAAE,QAAQ,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;IACvH,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACxC,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;QAC9B,MAAM,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;QAC1C,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;YAC9B,MAAM,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;YAC1E,IAAI,CAAC,KAAK,CAAC;gBAAE,SAAS;YACtB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE;gBAAE,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAE,IAAI,CAAC,GAAG,QAAQ,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAE,CAAC;QACvF,CAAC;IACH,CAAC;IACD,MAAM,OAAO,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACpD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE;QAAE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YACzD,IAAI,GAAG,GAAG,CAAC,CAAC;YACZ,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE;gBAAE,GAAG,IAAI,QAAQ,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAE,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAE,CAAC;YACzF,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,KAAK,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAChG,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC;QACnD,CAAC;IACD,MAAM,OAAO,GAAG,IAAI,YAAY,CAAC,IAAI,GAAG,IAAI,CAAC,CAAC;IAC9C,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,KAAK,IAAI,EAAE,GAAG,CAAC,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,EAAE,EAAE;QAAE,KAAK,IAAI,EAAE,GAAG,CAAC,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;YACjE,MAAM,CAAC,GAAG,OAAO,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAE,CAAC;YAC1D,OAAO,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,IAAI,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;YAAC,KAAK,IAAI,CAAC,CAAC;QACpD,CAAC;IACD,IAAI,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC;QAAE,MAAM,IAAI,4BAA4B,CAAC,0DAA0D,CAAC,CAAC;IAChJ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC,EAAE;QAAE,OAAO,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,CAAE,GAAG,KAAK,CAAC;IAC1E,IAAI,cAAc,GAAG,CAAC,CAAC;IACvB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,IAAI,EAAE,CAAC,EAAE;QAAE,cAAc,GAAG,IAAI,CAAC,GAAG,CAAC,cAAc,EAAE,MAAM,CAAC,CAAC,GAAG,IAAI,GAAG,MAAM,CAAC,CAAC,CAAC;IACrG,IAAI,gBAAgB,GAAG,CAAC,EAAE,MAAM,GAAkB,IAAI,EAAE,WAAW,GAAG,CAAC,CAAC;IACxE,MAAM,IAAI,GAAG,EAAE,CAAC;IAChB,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,IAAI,IAAI,GAAG,CAAC,EAAE,EAAE,EAAE;QAAE,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,IAAI,IAAI,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;YAC7E,MAAM,EAAE,GAAG,EAAE,GAAG,IAAI,EAAE,EAAE,GAAG,EAAE,GAAG,IAAI,EAAE,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,qCAAqC,CAAC,OAAO,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;YACrH,gBAAgB,GAAG,IAAI,CAAC,GAAG,CAAC,gBAAgB,EAAE,QAAQ,CAAC,CAAC;YACxD,IAAI,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,IAAI,MAAM;gBAAE,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,IAAI,CAAC,EAAE,QAAQ,CAAC,CAAC;QAC7E,CAAC;IACD,MAAM,IAAI,GAAG,GAAG,CAAC;IACjB,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,IAAI,EAAE,EAAE,EAAE;QAAE,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,IAAI,EAAE,EAAE,EAAE;YAAE,IAAI,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,IAAI,MAAM;gBAAE,WAAW,EAAE,CAAC;IAC9J,OAAO,EAAE,QAAQ,EAAE,cAAc,EAAE,qBAAqB,EAAE,KAAK,EAAE,oBAAoB,EAAE,MAAM,EAAE,+BAA+B,EAAE,WAAW,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC;QACzJ,cAAc,EAAE,gBAAgB,EAAE,2BAA2B,EAAE,MAAM,EAAE,OAAO,EAAE,CAAC;AACrF,CAAC;AACD,iHAAiH;AACjH,MAAM,UAAU,qCAAqC,CAAC,OAAqB,EAAE,MAAc,EAAE,EAAU,EAAE,EAAU;IACjH,MAAM,IAAI,GAAG,CAAC,GAAG,MAAM,GAAG,CAAC,CAAC;IAC5B,IAAI,GAAG,GAAG,CAAC,CAAC;IACZ,KAAK,IAAI,EAAE,GAAG,CAAC,MAAM,EAAE,EAAE,IAAI,MAAM,EAAE,EAAE,EAAE;QAAE,KAAK,IAAI,EAAE,GAAG,CAAC,MAAM,EAAE,EAAE,IAAI,MAAM,EAAE,EAAE,EAAE;YAAE,GAAG,IAAI,OAAO,CAAC,CAAC,EAAE,GAAG,MAAM,CAAC,GAAG,IAAI,GAAG,EAAE,GAAG,MAAM,CAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC;IACxL,OAAO,GAAG,CAAC;AACb,CAAC;AASD;;GAEG;AACH,MAAM,UAAU,0BAA0B,CAAC,GAAiB,EAAE,KAAa,EAAE,MAAc,EAAE,MAAiE;IAE5J,IAAI,GAAG,CAAC,MAAM,KAAK,KAAK,GAAG,MAAM,GAAG,CAAC;QAAE,MAAM,IAAI,yBAAyB,CAAC,8CAA8C,CAAC,CAAC;IAC3H,MAAM,CAAC,GAAG,MAAM,CAAC,QAAQ,CAAC,kBAAkB,EAAE,IAAI,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,OAAO,EAAE,GAAG,GAAG,IAAI,YAAY,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC;IACvH,MAAM,EAAE,GAAG,IAAI,UAAU,CAAC,KAAK,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;IACzC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,MAAM,EAAE,CAAC,EAAE;QAAE,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IACxF,IAAI,QAAQ,GAAG,CAAC,EAAE,KAAK,GAAG,CAAC,CAAC;IAC5B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QAChC,MAAM,CAAC,GAAG,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC;QACxB,KAAK,IAAI,EAAE,GAAG,CAAC,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;YAChC,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,KAAK,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC;YACxF,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,IAAI,EAAE,EAAE,EAAE,EAAE,CAAC;gBACjC,MAAM,MAAM,GAAG,CAAC,CAAC,EAAE,GAAG,EAAE,CAAE,CAAC;gBAC3B,IAAI,MAAM,KAAK,CAAC;oBAAE,SAAS;gBAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC/B,MAAM,CAAC,GAAG,GAAG,GAAG,EAAE,CAAC,CAAC,GAAG,EAAE,CAAE,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;oBAC3C,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAE,GAAG,MAAM,GAAG,GAAG,CAAC,CAAC,CAAE,CAAC;oBAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAE,GAAG,MAAM,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAE,CAAC;oBAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAE,GAAG,MAAM,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAE,CAAC;gBACxI,CAAC;YACH,CAAC;QACH,CAAC;IACH,CAAC;IACD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACpC,MAAM,CAAC,GAAG,GAAG,CAAC,CAAC,CAAE,CAAC;QAClB,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC;YAAE,MAAM,IAAI,4BAA4B,CAAC,sDAAsD,CAAC,CAAC;QACxH,IAAI,CAAC,GAAG,CAAC;YAAE,QAAQ,EAAE,CAAC;aAAM,IAAI,CAAC,GAAG,CAAC;YAAE,KAAK,EAAE,CAAC;IACjD,CAAC;IACD,OAAO,EAAE,OAAO,EAAE,GAAG,EAAE,WAAW,EAAE,EAAE,mBAAmB,EAAE,QAAQ,EAAE,qBAAqB,EAAE,KAAK,EAAE,EAAE,CAAC;AACxG,CAAC;AAED,2EAA2E;AAC3E,MAAM,UAAU,iBAAiB,CAAC,YAAoB;IACpD,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,YAAY,CAAC,EAAE,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;IACvE,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,KAAK,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAAC,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,YAAY,GAAG,YAAY,CAAC,CAAC,CAAC;QAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;QAAC,KAAK,IAAI,CAAC,CAAC;IAAC,CAAC;IAC3H,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,CAAC,EAAE;QAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAE,GAAG,KAAK,CAAC;IACxD,OAAO,CAAC,CAAC;AACX,CAAC;AAgBD;;;;GAIG;AACH,MAAM,UAAU,gBAAgB,CAAC,OAAqB,EAAE,KAAa,EAAE,MAAc,EAAE,QAA6B;IAElH,MAAM,CAAC,GAAG,wBAAwB,CAAC,QAAQ,CAAC,CAAC;IAC7C,IAAI,OAAO,CAAC,MAAM,KAAK,KAAK,GAAG,MAAM,GAAG,CAAC;QAAE,MAAM,IAAI,yBAAyB,CAAC,oCAAoC,CAAC,CAAC;IACrH,MAAM,CAAC,GAAG,iBAAiB,CAAC,CAAC,CAAC,YAAY,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,UAAU,GAAG,IAAI,YAAY,CAAC,OAAO,CAAC,MAAM,CAAC,EAAE,OAAO,GAAG,IAAI,YAAY,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC;IAC/J,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,EAAE,CAAC;YAChE,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;YAC9B,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC;YACxB,KAAK,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC7B,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAE,CAAC;gBACpF,CAAC,IAAI,CAAC,GAAG,OAAO,CAAC,EAAE,CAAE,CAAC;gBAAC,CAAC,IAAI,CAAC,GAAG,OAAO,CAAC,EAAE,GAAG,CAAC,CAAE,CAAC;gBAAC,CAAC,IAAI,CAAC,GAAG,OAAO,CAAC,EAAE,GAAG,CAAC,CAAE,CAAC;YAC9E,CAAC;YACD,UAAU,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YAAC,UAAU,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;YAAC,UAAU,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;QAClE,CAAC;IACD,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,EAAE,CAAC;YAChE,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;YAC9B,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC;YACxB,KAAK,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC7B,MAAM,EAAE,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAE,CAAC;gBACrF,CAAC,IAAI,CAAC,GAAG,UAAU,CAAC,EAAE,CAAE,CAAC;gBAAC,CAAC,IAAI,CAAC,GAAG,UAAU,CAAC,EAAE,GAAG,CAAC,CAAE,CAAC;gBAAC,CAAC,IAAI,CAAC,GAAG,UAAU,CAAC,EAAE,GAAG,CAAC,CAAE,CAAC;YACvF,CAAC;YACD,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YAAC,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;YAAC,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;QACzD,CAAC;IACD,MAAM,GAAG,GAAG,IAAI,YAAY,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC;IAC7C,IAAI,SAAS,GAAG,CAAC,EAAE,IAAI,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,IAAI,GAAG,CAAC,EAAE,QAAQ,GAAG,CAAC,EAAE,QAAQ,GAAG,CAAC,CAAC;IAC3E,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACxC,MAAM,CAAC,GAAG,OAAO,CAAC,CAAC,CAAE,EAAE,MAAM,GAAG,CAAC,GAAG,OAAO,CAAC,CAAC,CAAE,CAAC;QAChD,IAAI,MAAM,KAAK,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,SAAS,EAAE,CAAC;YAAC,IAAI,MAAM,KAAK,CAAC;gBAAE,IAAI,EAAE,CAAC;YAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YAAC,SAAS;QAAC,CAAC;QACvG,MAAM,MAAM,GAAG,CAAC,CAAC,MAAM,GAAG,MAAM,CAAC;QACjC,IAAI,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC;QACnB,SAAS,EAAE,CAAC;QACZ,IAAI,MAAM,GAAG,QAAQ;YAAE,QAAQ,GAAG,MAAM,CAAC;QACzC,IAAI,CAAC,MAAM,GAAG,QAAQ;YAAE,QAAQ,GAAG,CAAC,MAAM,CAAC;QAC3C,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;YAAC,CAAC,GAAG,CAAC,CAAC;YAAC,GAAG,EAAE,CAAC;QAAC,CAAC;aAAM,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;YAAC,CAAC,GAAG,CAAC,CAAC;YAAC,IAAI,EAAE,CAAC;QAAC,CAAC;QAC/D,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IACb,CAAC;IACD,OAAO,EAAE,OAAO,EAAE,GAAG,EAAE,WAAW,EAAE,EAAE,kBAAkB,EAAE,CAAC,EAAE,oBAAoB,EAAE,SAAS,EAAE,wBAAwB,EAAE,IAAI,EAAE,qBAAqB,EAAE,GAAG;YACtJ,sBAAsB,EAAE,IAAI,EAAE,eAAe,EAAE,QAAQ,EAAE,eAAe,EAAE,QAAQ,EAAE,EAAE,CAAC;AAC3F,CAAC;AAID;;;;GAIG;AACH,MAAM,UAAU,wBAAwB,CAAC,GAAiB,EAAE,KAAa,EAAE,MAAc,EAAE,OAA4B,EAAE,QAA6B;IAEpJ,IAAI,GAAG,CAAC,MAAM,KAAK,KAAK,GAAG,MAAM,GAAG,CAAC;QAAE,MAAM,IAAI,yBAAyB,CAAC,2BAA2B,CAAC,CAAC;IACxG,MAAM,KAAK,GAAG,CAAC,IAAI,OAAO,CAAC,mBAAmB,EAAE,IAAI,GAAG,mBAAmB,CAAC,OAAO,CAAC,EAAE,OAAO,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,EAAE,OAAO,GAAG,IAAI,YAAY,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC;IAC3J,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,IAAI,KAAK,IAAI,CAAC;QAAE,MAAM,IAAI,yBAAyB,CAAC,qCAAqC,CAAC,CAAC;IACtH,MAAM,MAAM,GAAG,OAAO,CAAC,aAAa,KAAK,qBAAqB,EAAE,MAAM,GAAoB,EAAE,QAAQ,EAAE,CAAC,EAAE,KAAK,EAAE,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,IAAI,EAAE,CAAC,EAAE,CAAC;IACjJ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,MAAM,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC;QACvC,MAAM,CAAC,GAAG,GAAG,CAAC,CAAC,CAAE,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAE,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAE,CAAC;QACpD,MAAM,EAAE,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,EAAE,EAAE,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,EAAE,EAAE,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;QACrG,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC;YAAE,MAAM,IAAI,4BAA4B,CAAC,kDAAkD,CAAC,CAAC;QACrK,IAAI,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,OAAO,CAAC,CAAC;QAC1B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,MAAM,OAAO,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC;YACjD,IAAI,OAAO,GAAG,CAAC;gBAAE,MAAM,CAAC,QAAQ,EAAE,CAAC;YACnC,IAAI,OAAO,GAAG,CAAC;gBAAE,MAAM,CAAC,KAAK,EAAE,CAAC;YAChC,IAAI,MAAM,GAAG,OAAO,CAAC,CAAC,CAAE,CAAC;YACzB,IAAI,MAAM,KAAK,OAAO;gBAAE,MAAM,CAAC,OAAO,EAAE,CAAC;YACzC,IAAI,MAAM,KAAK,CAAC;gBAAE,MAAM,GAAG,CAAC,CAAC;YAC7B,IAAI,MAAM,GAAG,CAAC,IAAI,MAAM,GAAG,CAAC,EAAE,CAAC;gBAC7B,IAAI,MAAM;oBAAE,MAAM,IAAI,yBAAyB,CAAC,mDAAmD,CAAC,CAAC;gBACrG,IAAI,MAAM,GAAG,CAAC,EAAE,CAAC;oBAAC,MAAM,CAAC,GAAG,EAAE,CAAC;oBAAC,MAAM,GAAG,CAAC,CAAC;gBAAC,CAAC;qBAAM,CAAC;oBAAC,MAAM,CAAC,IAAI,EAAE,CAAC;oBAAC,MAAM,GAAG,CAAC,CAAC;gBAAC,CAAC;YACnF,CAAC;YACD,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,mBAAmB,CAAC,MAAM,CAAC,CAAC;QAC/C,CAAC;IACH,CAAC;IACD,MAAM,SAAS,GAAG,gBAAgB,CAAC,OAAO,EAAE,KAAK,EAAE,MAAM,EAAE,QAAQ,CAAC,EAAE,OAAO,GAAG,CAAC,IAAI,OAAO,CAAC,QAAQ,GAAG,CAAC,CAAC;IAC1G,MAAM,KAAK,GAAG,OAAO,CAAC,QAAQ,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,UAAU,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,IAAI,WAAW,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC;IAChG,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,MAAM,EAAE,CAAC,EAAE;QAAE,KAAK,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,SAAS,CAAC,OAAO,CAAC,CAAC,CAAE,GAAG,OAAO,CAAC,CAAC;IAC9F,OAAO,EAAE,KAAK,EAAE,MAAM,EAAE,UAAU,EAAE,SAAS,CAAC,WAAW,EAAE,CAAC;AAC9D,CAAC","sourcesContent":["// SPDX-License-Identifier: Apache-2.0\n\n/**\n * Development detail stages (#552), both opt-in and off by default:\n *\n * - **Diffraction correction**: a bounded, regularized (Wiener–Tikhonov) deconvolution of the developed linear RGB with the\n *   diffraction OTF of an ideal circular pupil at the capture's working f-number. The OTF is zero beyond the cut-off\n *   f_c = 1/(λ·N_w), so detail beyond it never reached the sensor; the filter's design gain there is 0 and it cannot, and\n *   does not claim to, restore it.\n * - **Sharpening**: an unsharp mask (radius, amount, threshold) on the display-encoded values, after the tone curve and\n *   gamut handling and before quantization. It raises edge contrast and creates halos; it adds no optical detail.\n *\n * Development order with these stages: (exposure gain) → demosaic → colour and white balance → diffraction correction\n * (linear) → rendering exposure and tone curve → gamut handling → sRGB transfer → unsharp mask → quantization (encode).\n *\n * Image edges are handled by clamping (edge pixels repeated). Both stages work on a whole in-memory image.\n * @see docs/DIGITAL_OPTICS_FOUNDATION.md for the equations, sources, trade-offs and limits.\n */\n\nimport { encodeSrgbComponent } from \"../color/srgb-transfer.js\";\nimport { InvalidConfigurationError } from \"../core/configuration-error.js\";\nimport { InvalidScientificResultError } from \"../core/validation.js\";\nimport { prepareSdrPixelTone, type SdrRenderingProfile } from \"../output/sdr-rendering.js\";\n\n/** Version of the `development-detail` model (diffraction correction and unsharp mask). */\nexport const DEVELOPMENT_DETAIL_MODEL_VERSION = \"0.1.0\" as const;\n\n/** Unsharp mask settings: Gaussian standard deviation in output pixels, gain on the detail signal, and threshold in encoded units [0, 1]. */\nexport interface UnsharpMaskSettings { radiusPixels: number; amount: number; threshold: number }\n/** Diffraction correction settings: design wavelength, Wiener noise-to-signal regularization, gain bound and FIR radius in pixels. */\nexport interface DiffractionCorrectionSettings { wavelengthNm: number; regularization: number; maximumGain: number; kernelRadiusPixels: number }\n\n/** Accepted ranges. Bounds are engineering limits on cost and on noise gain, not quality recommendations. */\nexport const DEVELOPMENT_DETAIL_LIMITS = /* @__PURE__ */ Object.freeze({\n  maximumUnsharpRadiusPixels: 8, maximumUnsharpAmount: 5,\n  minimumWavelengthNm: 380, maximumWavelengthNm: 780,\n  minimumRegularization: 1e-4, maximumRegularization: 1,\n  maximumGain: 10, maximumKernelRadiusPixels: 16\n});\n\n/** Third dark ring of the Airy pattern, in units of λ·N (j₁,₃/π, j₁,₃ = 10.1734681…). */\nconst AIRY_THIRD_ZERO_PER_LAMBDA_N = 10.173468135062722 / Math.PI;\n\nfunction finiteIn(value: unknown, name: string, minimum: number, maximum: number, minimumExclusive = false): number {\n  if (typeof value !== \"number\" || !Number.isFinite(value) || (minimumExclusive ? value <= minimum : value < minimum) || value > maximum) {\n    throw new InvalidConfigurationError(name + \" must be a finite number in \" + (minimumExclusive ? \"(\" : \"[\") + minimum + \", \" + maximum + \"].\");\n  }\n  return value;\n}\n\n/** Validates unsharp mask settings. */\nexport function parseUnsharpMaskSettings(value: UnsharpMaskSettings, path = \"sharpening\"): UnsharpMaskSettings {\n  return { radiusPixels: finiteIn(value.radiusPixels, path + \".radiusPixels\", 0, DEVELOPMENT_DETAIL_LIMITS.maximumUnsharpRadiusPixels, true),\n    amount: finiteIn(value.amount, path + \".amount\", 0, DEVELOPMENT_DETAIL_LIMITS.maximumUnsharpAmount, true),\n    threshold: finiteIn(value.threshold, path + \".threshold\", 0, 1) };\n}\n/** Validates diffraction correction settings. */\nexport function parseDiffractionCorrectionSettings(value: DiffractionCorrectionSettings, path = \"diffractionCorrection\"): DiffractionCorrectionSettings {\n  const radius = value.kernelRadiusPixels;\n  if (typeof radius !== \"number\" || !Number.isSafeInteger(radius) || radius < 1 || radius > DEVELOPMENT_DETAIL_LIMITS.maximumKernelRadiusPixels) {\n    throw new InvalidConfigurationError(path + \".kernelRadiusPixels must be an integer from 1 to \" + DEVELOPMENT_DETAIL_LIMITS.maximumKernelRadiusPixels + \".\");\n  }\n  return { wavelengthNm: finiteIn(value.wavelengthNm, path + \".wavelengthNm\", DEVELOPMENT_DETAIL_LIMITS.minimumWavelengthNm, DEVELOPMENT_DETAIL_LIMITS.maximumWavelengthNm),\n    regularization: finiteIn(value.regularization, path + \".regularization\", DEVELOPMENT_DETAIL_LIMITS.minimumRegularization, DEVELOPMENT_DETAIL_LIMITS.maximumRegularization),\n    maximumGain: finiteIn(value.maximumGain, path + \".maximumGain\", 1, DEVELOPMENT_DETAIL_LIMITS.maximumGain), kernelRadiusPixels: radius };\n}\n/** Default FIR radius: the Airy pattern's third dark ring (3.2383·λ·N_w) in pixels, rounded up, within 1..16. */\nexport function defaultDiffractionCorrectionRadiusPixels(workingFNumber: number, wavelengthNm: number, pixelPitchMicrometers: number): number {\n  const sites = Math.ceil(AIRY_THIRD_ZERO_PER_LAMBDA_N * wavelengthNm / 1000 * workingFNumber / pixelPitchMicrometers);\n  return Math.min(DEVELOPMENT_DETAIL_LIMITS.maximumKernelRadiusPixels, Math.max(1, sites));\n}\n\n/**\n * Incoherent OTF of an ideal, aberration-free circular pupil at normalized frequency s = f/f_c:\n * H(s) = (2/π)(acos s − s·√(1 − s²)) for s < 1, and 0 for s ≥ 1 (MIT OCW 2.71, lectures 21–22; Goodman, §7.3).\n */\nexport function circularPupilDiffractionMtf(normalizedFrequency: number): number {\n  const s = Math.abs(normalizedFrequency);\n  if (s >= 1) return 0;\n  return 2 / Math.PI * (Math.acos(s) - s * Math.sqrt(1 - s * s));\n}\n\n/** A designed diffraction correction filter and what its response does. Frequencies are in cycles per pixel. */\nexport interface DiffractionCorrectionFilter {\n  settings: DiffractionCorrectionSettings;\n  workingFNumber: number;\n  pixelPitchMicrometers: number;\n  /** f_c · pitch = pitch/(λ·N_w): the diffraction cut-off. The sampling Nyquist frequency is 0.5. */\n  cutoffCyclesPerPixel: number;\n  /** Fraction of the Nyquist square |fx|, |fy| ≤ 0.5 lying beyond the cut-off, where the lens transmitted nothing. */\n  nyquistAreaBeyondCutoffFraction: number;\n  /** Largest gain of the designed (untruncated) response: min(maximumGain, max over s of H(1+α)/(H² + α)). */\n  designPeakGain: number;\n  /** Largest |response| of the windowed, truncated FIR over a 1/64-cycle grid of the Nyquist square. */\n  responsePeakGain: number;\n  /**\n   * Largest |response| of the FIR at grid frequencies beyond the cut-off (null when the whole Nyquist square is inside it).\n   * The design is 0 there; a non-zero value is truncation leakage that passes noise and aliased signal, never lost detail,\n   * because the optics delivered no signal at those frequencies.\n   */\n  responseMaximumBeyondCutoff: number | null;\n  /** Row-major (2R+1)² weights summing to 1; weights[(dy + R)(2R + 1) + dx + R]. */\n  weights: Float64Array;\n}\n\n/**\n * Designs the bounded Wiener–Tikhonov correction. In frequency (normalized s = ρ/f_c, ρ the radial frequency):\n *   G(s) = min(maximumGain, H(s)·(1 + α)/(H(s)² + α)),   G(0) = 1, and G(s) = 0 for s ≥ 1,\n * where α is the regularization (an assumed noise-to-signal power ratio, Wiener's K). The response is sampled on an M × M\n * grid (M a power of two ≥ max(64, 4(2R + 1))), inverse transformed, truncated to the (2R + 1)² window with a radial raised\n * cosine taper and renormalized to unit sum. The achieved response is then measured, not assumed.\n */\nexport function designDiffractionCorrectionFilter(input: { workingFNumber: number; pixelPitchMicrometers: number; settings: DiffractionCorrectionSettings }): DiffractionCorrectionFilter {\n  const settings = parseDiffractionCorrectionSettings(input.settings);\n  const workingFNumber = finiteIn(input.workingFNumber, \"workingFNumber\", 0, 1e4, true), pitch = finiteIn(input.pixelPitchMicrometers, \"pixelPitchMicrometers\", 0, 1e4, true);\n  const cutoff = pitch / (settings.wavelengthNm / 1000 * workingFNumber), alpha = settings.regularization, bound = settings.maximumGain;\n  const gainAt = (rho: number): number => {\n    const h = circularPupilDiffractionMtf(rho / cutoff);\n    return h === 0 ? 0 : Math.min(bound, h * (1 + alpha) / (h * h + alpha));\n  };\n  const R = settings.kernelRadiusPixels, size = 2 * R + 1;\n  let M = 64;\n  while (M < 4 * size) M *= 2;\n  // A(ky, x) = Σ_kx G cos(2π kx x/M); h(x, y) = (1/M²) Σ_ky cos(2π ky y/M) A(ky, x). G is even in each axis, so h is real and even.\n  const cosTable = new Float64Array(M * (R + 1));\n  for (let k = 0; k < M; k++) for (let x = 0; x <= R; x++) cosTable[k * (R + 1) + x] = Math.cos(2 * Math.PI * k * x / M);\n  const A = new Float64Array(M * (R + 1));\n  for (let ky = 0; ky < M; ky++) {\n    const fy = (ky < M / 2 ? ky : ky - M) / M;\n    for (let kx = 0; kx < M; kx++) {\n      const fx = (kx < M / 2 ? kx : kx - M) / M, g = gainAt(Math.hypot(fx, fy));\n      if (g === 0) continue;\n      for (let x = 0; x <= R; x++) A[ky * (R + 1) + x]! += g * cosTable[kx * (R + 1) + x]!;\n    }\n  }\n  const quarter = new Float64Array((R + 1) * (R + 1));\n  for (let y = 0; y <= R; y++) for (let x = 0; x <= R; x++) {\n    let sum = 0;\n    for (let ky = 0; ky < M; ky++) sum += cosTable[ky * (R + 1) + y]! * A[ky * (R + 1) + x]!;\n    const r = Math.hypot(x, y), taper = r < R + 1 ? 0.5 * (1 + Math.cos(Math.PI * r / (R + 1))) : 0;\n    quarter[y * (R + 1) + x] = sum / (M * M) * taper;\n  }\n  const weights = new Float64Array(size * size);\n  let total = 0;\n  for (let dy = -R; dy <= R; dy++) for (let dx = -R; dx <= R; dx++) {\n    const w = quarter[Math.abs(dy) * (R + 1) + Math.abs(dx)]!;\n    weights[(dy + R) * size + dx + R] = w; total += w;\n  }\n  if (!(total > 0) || !Number.isFinite(total)) throw new InvalidScientificResultError(\"Diffraction correction filter has no finite DC response.\");\n  for (let i = 0; i < weights.length; i++) weights[i] = weights[i]! / total;\n  let designPeakGain = 0;\n  for (let i = 0; i <= 4096; i++) designPeakGain = Math.max(designPeakGain, gainAt(i / 4096 * cutoff));\n  let responsePeakGain = 0, beyond: number | null = null, beyondCount = 0;\n  const GRID = 64;\n  for (let iy = 0; iy <= GRID / 2; iy++) for (let ix = 0; ix <= GRID / 2; ix++) {\n    const fx = ix / GRID, fy = iy / GRID, response = Math.abs(evaluateDiffractionCorrectionResponse(weights, R, fx, fy));\n    responsePeakGain = Math.max(responsePeakGain, response);\n    if (Math.hypot(fx, fy) >= cutoff) beyond = Math.max(beyond ?? 0, response);\n  }\n  const AREA = 256;\n  for (let iy = 0; iy < AREA; iy++) for (let ix = 0; ix < AREA; ix++) if (Math.hypot((ix + 0.5) / (2 * AREA), (iy + 0.5) / (2 * AREA)) >= cutoff) beyondCount++;\n  return { settings, workingFNumber, pixelPitchMicrometers: pitch, cutoffCyclesPerPixel: cutoff, nyquistAreaBeyondCutoffFraction: beyondCount / (AREA * AREA),\n    designPeakGain, responsePeakGain, responseMaximumBeyondCutoff: beyond, weights };\n}\n/** Frequency response Σ w(x, y)·cos(2π(fx·x + fy·y)) of a centred (2R + 1)² FIR at (fx, fy) cycles per pixel. */\nexport function evaluateDiffractionCorrectionResponse(weights: Float64Array, radius: number, fx: number, fy: number): number {\n  const size = 2 * radius + 1;\n  let sum = 0;\n  for (let dy = -radius; dy <= radius; dy++) for (let dx = -radius; dx <= radius; dx++) sum += weights[(dy + radius) * size + dx + radius]! * Math.cos(2 * Math.PI * (fx * dx + fy * dy));\n  return sum;\n}\n\n/** What the correction did to the linear samples. */\nexport interface DiffractionCorrectionDiagnostics {\n  /** Corrected samples below 0 (ringing below black). They are clipped by the later gamut handling. */\n  negativeSampleCount: number;\n  /** Corrected samples above 1 (above the RAW white level). */\n  aboveWhiteSampleCount: number;\n}\n/**\n * Applies the FIR to interleaved linear RGB (width × height × 3), each channel separately, edges clamped. Returns new samples.\n */\nexport function applyDiffractionCorrection(rgb: Float64Array, width: number, height: number, filter: Pick<DiffractionCorrectionFilter, \"weights\" | \"settings\">):\n  { samples: Float64Array; diagnostics: DiffractionCorrectionDiagnostics } {\n  if (rgb.length !== width * height * 3) throw new InvalidConfigurationError(\"Diffraction correction raster size mismatch.\");\n  const R = filter.settings.kernelRadiusPixels, size = 2 * R + 1, w = filter.weights, out = new Float64Array(rgb.length);\n  const xs = new Int32Array(width + 2 * R);\n  for (let i = 0; i < xs.length; i++) xs[i] = Math.min(width - 1, Math.max(0, i - R)) * 3;\n  let negative = 0, above = 0;\n  for (let y = 0; y < height; y++) {\n    const o = y * width * 3;\n    for (let dy = -R; dy <= R; dy++) {\n      const row = Math.min(height - 1, Math.max(0, y + dy)) * width * 3, wr = (dy + R) * size;\n      for (let dx = 0; dx < size; dx++) {\n        const weight = w[wr + dx]!;\n        if (weight === 0) continue;\n        for (let x = 0; x < width; x++) {\n          const s = row + xs[x + dx]!, d = o + 3 * x;\n          out[d] = out[d]! + weight * rgb[s]!; out[d + 1] = out[d + 1]! + weight * rgb[s + 1]!; out[d + 2] = out[d + 2]! + weight * rgb[s + 2]!;\n        }\n      }\n    }\n  }\n  for (let i = 0; i < out.length; i++) {\n    const v = out[i]!;\n    if (!Number.isFinite(v)) throw new InvalidScientificResultError(\"Diffraction correction produced a non-finite sample.\");\n    if (v < 0) negative++; else if (v > 1) above++;\n  }\n  return { samples: out, diagnostics: { negativeSampleCount: negative, aboveWhiteSampleCount: above } };\n}\n\n/** Normalized Gaussian weights exp(−i²/(2σ²)), i = −r..r with r = ⌈3σ⌉. */\nexport function unsharpMaskKernel(radiusPixels: number): Float64Array {\n  const r = Math.ceil(3 * radiusPixels), k = new Float64Array(2 * r + 1);\n  let total = 0;\n  for (let i = -r; i <= r; i++) { const v = Math.exp(-i * i / (2 * radiusPixels * radiusPixels)); k[i + r] = v; total += v; }\n  for (let i = 0; i < k.length; i++) k[i] = k[i]! / total;\n  return k;\n}\n/** What the unsharp mask did. Counts are samples (one per channel). */\nexport interface UnsharpMaskDiagnostics {\n  /** Half-width of the separable Gaussian, ⌈3σ⌉ pixels. */\n  kernelRadiusPixels: number;\n  /** Samples changed by the mask (|detail| ≥ threshold and detail ≠ 0). */\n  sharpenedSampleCount: number;\n  /** Samples left unchanged because 0 < |detail| < threshold. */\n  thresholdHeldSampleCount: number;\n  /** Sharpened samples clipped at 0 or 1 (the halo did not fit in the code range). */\n  clippedLowSampleCount: number;\n  clippedHighSampleCount: number;\n  /** Largest increase and decrease amount·detail before clipping, in encoded units: the halo strength. */\n  maximumIncrease: number;\n  maximumDecrease: number;\n}\n/**\n * Unsharp mask on interleaved encoded RGB in [0, 1] (each channel separately, edges clamped):\n *   d = x − (G_σ ∗ x);   y = x + amount·d where |d| ≥ threshold, else y = x;   then y is clipped to [0, 1].\n * Returns new samples.\n */\nexport function applyUnsharpMask(encoded: Float64Array, width: number, height: number, settings: UnsharpMaskSettings):\n  { samples: Float64Array; diagnostics: UnsharpMaskDiagnostics } {\n  const s = parseUnsharpMaskSettings(settings);\n  if (encoded.length !== width * height * 3) throw new InvalidConfigurationError(\"Unsharp mask raster size mismatch.\");\n  const k = unsharpMaskKernel(s.radiusPixels), r = (k.length - 1) / 2, horizontal = new Float64Array(encoded.length), blurred = new Float64Array(encoded.length);\n  for (let y = 0; y < height; y++) for (let x = 0; x < width; x++) {\n    const d = (y * width + x) * 3;\n    let a = 0, b = 0, c = 0;\n    for (let i = -r; i <= r; i++) {\n      const sx = (y * width + Math.min(width - 1, Math.max(0, x + i))) * 3, w = k[i + r]!;\n      a += w * encoded[sx]!; b += w * encoded[sx + 1]!; c += w * encoded[sx + 2]!;\n    }\n    horizontal[d] = a; horizontal[d + 1] = b; horizontal[d + 2] = c;\n  }\n  for (let y = 0; y < height; y++) for (let x = 0; x < width; x++) {\n    const d = (y * width + x) * 3;\n    let a = 0, b = 0, c = 0;\n    for (let i = -r; i <= r; i++) {\n      const sy = (Math.min(height - 1, Math.max(0, y + i)) * width + x) * 3, w = k[i + r]!;\n      a += w * horizontal[sy]!; b += w * horizontal[sy + 1]!; c += w * horizontal[sy + 2]!;\n    }\n    blurred[d] = a; blurred[d + 1] = b; blurred[d + 2] = c;\n  }\n  const out = new Float64Array(encoded.length);\n  let sharpened = 0, held = 0, low = 0, high = 0, increase = 0, decrease = 0;\n  for (let i = 0; i < encoded.length; i++) {\n    const x = encoded[i]!, detail = x - blurred[i]!;\n    if (detail === 0 || Math.abs(detail) < s.threshold) { if (detail !== 0) held++; out[i] = x; continue; }\n    const change = s.amount * detail;\n    let v = x + change;\n    sharpened++;\n    if (change > increase) increase = change;\n    if (-change > decrease) decrease = -change;\n    if (v < 0) { v = 0; low++; } else if (v > 1) { v = 1; high++; }\n    out[i] = v;\n  }\n  return { samples: out, diagnostics: { kernelRadiusPixels: r, sharpenedSampleCount: sharpened, thresholdHeldSampleCount: held, clippedLowSampleCount: low,\n    clippedHighSampleCount: high, maximumIncrease: increase, maximumDecrease: decrease } };\n}\n\n/** SDR counters in the shape of the SDR rendering diagnostics. */\nexport interface DetailSdrCounts { negative: number; above: number; changed: number; low: number; high: number }\n/**\n * SDR rendering with an unsharp mask between the sRGB transfer and quantization. Per pixel it runs the reference renderer's\n * operations in its order (−0 → +0, rendering exposure 2^EV with reference white 1, the shared tone operator, gamut\n * handling, sRGB transfer); the whole encoded image is then sharpened and rounded to nearest (ties up).\n */\nexport function renderSdrWithUnsharpMask(rgb: Float64Array, width: number, height: number, profile: SdrRenderingProfile, settings: UnsharpMaskSettings):\n  { codes: Uint8Array | Uint16Array; counts: DetailSdrCounts; sharpening: UnsharpMaskDiagnostics } {\n  if (rgb.length !== width * height * 3) throw new InvalidConfigurationError(\"SDR raster size mismatch.\");\n  const scale = 2 ** profile.renderingExposureEv, tone = prepareSdrPixelTone(profile), scratch = new Float64Array(3), encoded = new Float64Array(rgb.length);\n  if (!Number.isFinite(scale) || scale <= 0) throw new InvalidConfigurationError(\"Rendering scale is unrepresentable.\");\n  const reject = profile.gamutHandling === \"reject-out-of-range\", counts: DetailSdrCounts = { negative: 0, above: 0, changed: 0, low: 0, high: 0 };\n  for (let i = 0; i < rgb.length; i += 3) {\n    const r = rgb[i]!, g = rgb[i + 1]!, b = rgb[i + 2]!;\n    const er = (r === 0 ? 0 : r) * scale, eg = (g === 0 ? 0 : g) * scale, eb = (b === 0 ? 0 : b) * scale;\n    if (!Number.isFinite(er) || !Number.isFinite(eg) || !Number.isFinite(eb)) throw new InvalidScientificResultError(\"Rendering exposure produced a non-finite sample.\");\n    tone(er, eg, eb, scratch);\n    for (let c = 0; c < 3; c++) {\n      const exposed = c === 0 ? er : c === 1 ? eg : eb;\n      if (exposed < 0) counts.negative++;\n      if (exposed > 1) counts.above++;\n      let mapped = scratch[c]!;\n      if (mapped !== exposed) counts.changed++;\n      if (mapped === 0) mapped = 0;\n      if (mapped < 0 || mapped > 1) {\n        if (reject) throw new InvalidConfigurationError(\"Tone output exceeds the declared SDR gamut range.\");\n        if (mapped < 0) { counts.low++; mapped = 0; } else { counts.high++; mapped = 1; }\n      }\n      encoded[i + c] = encodeSrgbComponent(mapped);\n    }\n  }\n  const sharpened = applyUnsharpMask(encoded, width, height, settings), maximum = 2 ** profile.bitDepth - 1;\n  const codes = profile.bitDepth === 8 ? new Uint8Array(rgb.length) : new Uint16Array(rgb.length);\n  for (let i = 0; i < codes.length; i++) codes[i] = Math.round(sharpened.samples[i]! * maximum);\n  return { codes, counts, sharpening: sharpened.diagnostics };\n}\n"]}
+{"version":3,"file":"development-detail.js","sourceRoot":"","sources":["../../src/processing/development-detail.ts"],"names":[],"mappings":"AAAA,sCAAsC;AAEtC;;;;;;;;;;;;;;;;GAgBG;AAEH,OAAO,EAAE,yBAAyB,EAAE,MAAM,gCAAgC,CAAC;AAC3E,OAAO,EAAE,4BAA4B,EAAE,MAAM,uBAAuB,CAAC;AACrE,OAAO,EAAE,gBAAgB,EAAuB,MAAM,iCAAiC,CAAC;AACxF,OAAO,EAAE,mBAAmB,EAA4B,MAAM,4BAA4B,CAAC;AAE3F,2FAA2F;AAC3F,MAAM,CAAC,MAAM,gCAAgC,GAAG,OAAgB,CAAC;AAOjE,6GAA6G;AAC7G,MAAM,CAAC,MAAM,yBAAyB,GAAG,eAAe,CAAC,MAAM,CAAC,MAAM,CAAC;IACrE,0BAA0B,EAAE,CAAC,EAAE,oBAAoB,EAAE,CAAC;IACtD,mBAAmB,EAAE,GAAG,EAAE,mBAAmB,EAAE,GAAG;IAClD,qBAAqB,EAAE,IAAI,EAAE,qBAAqB,EAAE,CAAC;IACrD,WAAW,EAAE,EAAE,EAAE,yBAAyB,EAAE,EAAE;CAC/C,CAAC,CAAC;AAEH,yFAAyF;AACzF,MAAM,4BAA4B,GAAG,kBAAkB,GAAG,IAAI,CAAC,EAAE,CAAC;AAElE,SAAS,QAAQ,CAAC,KAAc,EAAE,IAAY,EAAE,OAAe,EAAE,OAAe,EAAE,gBAAgB,GAAG,KAAK;IACxG,IAAI,OAAO,KAAK,KAAK,QAAQ,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,IAAI,CAAC,gBAAgB,CAAC,CAAC,CAAC,KAAK,IAAI,OAAO,CAAC,CAAC,CAAC,KAAK,GAAG,OAAO,CAAC,IAAI,KAAK,GAAG,OAAO,EAAE,CAAC;QACvI,MAAM,IAAI,yBAAyB,CAAC,IAAI,GAAG,8BAA8B,GAAG,CAAC,gBAAgB,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,OAAO,GAAG,IAAI,GAAG,OAAO,GAAG,IAAI,CAAC,CAAC;IAChJ,CAAC;IACD,OAAO,KAAK,CAAC;AACf,CAAC;AAED,uCAAuC;AACvC,MAAM,UAAU,wBAAwB,CAAC,KAA0B,EAAE,IAAI,GAAG,YAAY;IACtF,OAAO,EAAE,YAAY,EAAE,QAAQ,CAAC,KAAK,CAAC,YAAY,EAAE,IAAI,GAAG,eAAe,EAAE,CAAC,EAAE,yBAAyB,CAAC,0BAA0B,EAAE,IAAI,CAAC;QACxI,MAAM,EAAE,QAAQ,CAAC,KAAK,CAAC,MAAM,EAAE,IAAI,GAAG,SAAS,EAAE,CAAC,EAAE,yBAAyB,CAAC,oBAAoB,EAAE,IAAI,CAAC;QACzG,SAAS,EAAE,QAAQ,CAAC,KAAK,CAAC,SAAS,EAAE,IAAI,GAAG,YAAY,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC;AACtE,CAAC;AACD,iDAAiD;AACjD,MAAM,UAAU,kCAAkC,CAAC,KAAoC,EAAE,IAAI,GAAG,uBAAuB;IACrH,MAAM,MAAM,GAAG,KAAK,CAAC,kBAAkB,CAAC;IACxC,IAAI,OAAO,MAAM,KAAK,QAAQ,IAAI,CAAC,MAAM,CAAC,aAAa,CAAC,MAAM,CAAC,IAAI,MAAM,GAAG,CAAC,IAAI,MAAM,GAAG,yBAAyB,CAAC,yBAAyB,EAAE,CAAC;QAC9I,MAAM,IAAI,yBAAyB,CAAC,IAAI,GAAG,mDAAmD,GAAG,yBAAyB,CAAC,yBAAyB,GAAG,GAAG,CAAC,CAAC;IAC9J,CAAC;IACD,OAAO,EAAE,YAAY,EAAE,QAAQ,CAAC,KAAK,CAAC,YAAY,EAAE,IAAI,GAAG,eAAe,EAAE,yBAAyB,CAAC,mBAAmB,EAAE,yBAAyB,CAAC,mBAAmB,CAAC;QACvK,cAAc,EAAE,QAAQ,CAAC,KAAK,CAAC,cAAc,EAAE,IAAI,GAAG,iBAAiB,EAAE,yBAAyB,CAAC,qBAAqB,EAAE,yBAAyB,CAAC,qBAAqB,CAAC;QAC1K,WAAW,EAAE,QAAQ,CAAC,KAAK,CAAC,WAAW,EAAE,IAAI,GAAG,cAAc,EAAE,CAAC,EAAE,yBAAyB,CAAC,WAAW,CAAC,EAAE,kBAAkB,EAAE,MAAM,EAAE,CAAC;AAC5I,CAAC;AACD,iHAAiH;AACjH,MAAM,UAAU,wCAAwC,CAAC,cAAsB,EAAE,YAAoB,EAAE,qBAA6B;IAClI,MAAM,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,4BAA4B,GAAG,YAAY,GAAG,IAAI,GAAG,cAAc,GAAG,qBAAqB,CAAC,CAAC;IACrH,OAAO,IAAI,CAAC,GAAG,CAAC,yBAAyB,CAAC,yBAAyB,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC;AAC3F,CAAC;AAED;;;GAGG;AACH,MAAM,UAAU,2BAA2B,CAAC,mBAA2B;IACrE,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,mBAAmB,CAAC,CAAC;IACxC,IAAI,CAAC,IAAI,CAAC;QAAE,OAAO,CAAC,CAAC;IACrB,OAAO,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;AACjE,CAAC;AAyBD;;;;;;GAMG;AACH,MAAM,UAAU,iCAAiC,CAAC,KAAyG;IACzJ,MAAM,QAAQ,GAAG,kCAAkC,CAAC,KAAK,CAAC,QAAQ,CAAC,CAAC;IACpE,MAAM,cAAc,GAAG,QAAQ,CAAC,KAAK,CAAC,cAAc,EAAE,gBAAgB,EAAE,CAAC,EAAE,GAAG,EAAE,IAAI,CAAC,EAAE,KAAK,GAAG,QAAQ,CAAC,KAAK,CAAC,qBAAqB,EAAE,uBAAuB,EAAE,CAAC,EAAE,GAAG,EAAE,IAAI,CAAC,CAAC;IAC5K,MAAM,MAAM,GAAG,KAAK,GAAG,CAAC,QAAQ,CAAC,YAAY,GAAG,IAAI,GAAG,cAAc,CAAC,EAAE,KAAK,GAAG,QAAQ,CAAC,cAAc,EAAE,KAAK,GAAG,QAAQ,CAAC,WAAW,CAAC;IACtI,MAAM,MAAM,GAAG,CAAC,GAAW,EAAU,EAAE;QACrC,MAAM,CAAC,GAAG,2BAA2B,CAAC,GAAG,GAAG,MAAM,CAAC,CAAC;QACpD,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;IAC1E,CAAC,CAAC;IACF,MAAM,CAAC,GAAG,QAAQ,CAAC,kBAAkB,EAAE,IAAI,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IACxD,IAAI,CAAC,GAAG,EAAE,CAAC;IACX,OAAO,CAAC,GAAG,CAAC,GAAG,IAAI;QAAE,CAAC,IAAI,CAAC,CAAC;IAC5B,kIAAkI;IAClI,MAAM,QAAQ,GAAG,IAAI,YAAY,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IAC/C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;QAAE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE;YAAE,QAAQ,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;IACvH,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACxC,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;QAC9B,MAAM,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;QAC1C,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;YAC9B,MAAM,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;YAC1E,IAAI,CAAC,KAAK,CAAC;gBAAE,SAAS;YACtB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE;gBAAE,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAE,IAAI,CAAC,GAAG,QAAQ,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAE,CAAC;QACvF,CAAC;IACH,CAAC;IACD,MAAM,OAAO,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACpD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE;QAAE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YACzD,IAAI,GAAG,GAAG,CAAC,CAAC;YACZ,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE;gBAAE,GAAG,IAAI,QAAQ,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAE,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAE,CAAC;YACzF,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,KAAK,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAChG,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC;QACnD,CAAC;IACD,MAAM,OAAO,GAAG,IAAI,YAAY,CAAC,IAAI,GAAG,IAAI,CAAC,CAAC;IAC9C,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,KAAK,IAAI,EAAE,GAAG,CAAC,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,EAAE,EAAE;QAAE,KAAK,IAAI,EAAE,GAAG,CAAC,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;YACjE,MAAM,CAAC,GAAG,OAAO,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAE,CAAC;YAC1D,OAAO,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,IAAI,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;YAAC,KAAK,IAAI,CAAC,CAAC;QACpD,CAAC;IACD,IAAI,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC;QAAE,MAAM,IAAI,4BAA4B,CAAC,0DAA0D,CAAC,CAAC;IAChJ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC,EAAE;QAAE,OAAO,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,CAAE,GAAG,KAAK,CAAC;IAC1E,IAAI,cAAc,GAAG,CAAC,CAAC;IACvB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,IAAI,EAAE,CAAC,EAAE;QAAE,cAAc,GAAG,IAAI,CAAC,GAAG,CAAC,cAAc,EAAE,MAAM,CAAC,CAAC,GAAG,IAAI,GAAG,MAAM,CAAC,CAAC,CAAC;IACrG,IAAI,gBAAgB,GAAG,CAAC,EAAE,MAAM,GAAkB,IAAI,EAAE,WAAW,GAAG,CAAC,CAAC;IACxE,MAAM,IAAI,GAAG,EAAE,CAAC;IAChB,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,IAAI,IAAI,GAAG,CAAC,EAAE,EAAE,EAAE;QAAE,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,IAAI,IAAI,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;YAC7E,MAAM,EAAE,GAAG,EAAE,GAAG,IAAI,EAAE,EAAE,GAAG,EAAE,GAAG,IAAI,EAAE,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,qCAAqC,CAAC,OAAO,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;YACrH,gBAAgB,GAAG,IAAI,CAAC,GAAG,CAAC,gBAAgB,EAAE,QAAQ,CAAC,CAAC;YACxD,IAAI,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,IAAI,MAAM;gBAAE,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,IAAI,CAAC,EAAE,QAAQ,CAAC,CAAC;QAC7E,CAAC;IACD,MAAM,IAAI,GAAG,GAAG,CAAC;IACjB,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,IAAI,EAAE,EAAE,EAAE;QAAE,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,IAAI,EAAE,EAAE,EAAE;YAAE,IAAI,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,IAAI,MAAM;gBAAE,WAAW,EAAE,CAAC;IAC9J,OAAO,EAAE,QAAQ,EAAE,cAAc,EAAE,qBAAqB,EAAE,KAAK,EAAE,oBAAoB,EAAE,MAAM,EAAE,+BAA+B,EAAE,WAAW,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC;QACzJ,cAAc,EAAE,gBAAgB,EAAE,2BAA2B,EAAE,MAAM,EAAE,OAAO,EAAE,CAAC;AACrF,CAAC;AACD,iHAAiH;AACjH,MAAM,UAAU,qCAAqC,CAAC,OAAqB,EAAE,MAAc,EAAE,EAAU,EAAE,EAAU;IACjH,MAAM,IAAI,GAAG,CAAC,GAAG,MAAM,GAAG,CAAC,CAAC;IAC5B,IAAI,GAAG,GAAG,CAAC,CAAC;IACZ,KAAK,IAAI,EAAE,GAAG,CAAC,MAAM,EAAE,EAAE,IAAI,MAAM,EAAE,EAAE,EAAE;QAAE,KAAK,IAAI,EAAE,GAAG,CAAC,MAAM,EAAE,EAAE,IAAI,MAAM,EAAE,EAAE,EAAE;YAAE,GAAG,IAAI,OAAO,CAAC,CAAC,EAAE,GAAG,MAAM,CAAC,GAAG,IAAI,GAAG,EAAE,GAAG,MAAM,CAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC;IACxL,OAAO,GAAG,CAAC;AACb,CAAC;AASD;;GAEG;AACH,MAAM,UAAU,0BAA0B,CAAC,GAAiB,EAAE,KAAa,EAAE,MAAc,EAAE,MAAiE;IAE5J,IAAI,GAAG,CAAC,MAAM,KAAK,KAAK,GAAG,MAAM,GAAG,CAAC;QAAE,MAAM,IAAI,yBAAyB,CAAC,8CAA8C,CAAC,CAAC;IAC3H,MAAM,OAAO,GAAG,YAAY,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;IACvC,OAAO,EAAE,OAAO,EAAE,WAAW,EAAE,eAAe,CAAC,yBAAyB,CAAC,OAAO,EAAE,KAAK,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC,EAAE,CAAC;AAC9G,CAAC;AACD,wEAAwE;AACxE,MAAM,UAAU,eAAe,CAAI,IAAgC;IACjE,IAAI,IAAI,GAAG,IAAI,CAAC,IAAI,EAAE,CAAC;IACvB,OAAO,CAAC,IAAI,CAAC,IAAI;QAAE,IAAI,GAAG,IAAI,CAAC,IAAI,EAAE,CAAC;IACtC,OAAO,IAAI,CAAC,KAAK,CAAC;AACpB,CAAC;AACD;;;;;GAKG;AACH,MAAM,SAAS,CAAC,CAAC,yBAAyB,CAAC,GAAiB,EAAE,KAAa,EAAE,MAAc,EAAE,MAAiE;IAE5J,MAAM,CAAC,GAAG,MAAM,CAAC,QAAQ,CAAC,kBAAkB,EAAE,IAAI,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,OAAO,EAAE,SAAS,GAAG,KAAK,GAAG,CAAC,EAAE,IAAI,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC;IACxJ,MAAM,EAAE,GAAG,IAAI,UAAU,CAAC,KAAK,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;IACzC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,MAAM,EAAE,CAAC,EAAE;QAAE,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IACxF,IAAI,QAAQ,GAAG,CAAC,EAAE,KAAK,GAAG,CAAC,CAAC;IAC5B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,MAAM,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QACrC,2HAA2H;QAC3H,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,SAAS,CAAC;QACpC,IAAI,CAAC,GAAG,CAAC;YAAE,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC,EAAE,CAAC,GAAG,SAAS,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC;QAC7E,IAAI,CAAC,IAAI,MAAM;YAAE,SAAS;QAC1B,IAAI,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,GAAG,SAAS,CAAC,CAAC;QAC/B,KAAK,IAAI,EAAE,GAAG,CAAC,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;YAChC,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,SAAS,EAAE,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC;YACxF,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,IAAI,EAAE,EAAE,EAAE,EAAE,CAAC;gBACjC,MAAM,MAAM,GAAG,CAAC,CAAC,EAAE,GAAG,EAAE,CAAE,CAAC;gBAC3B,IAAI,MAAM,KAAK,CAAC;oBAAE,SAAS;gBAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC/B,MAAM,CAAC,GAAG,GAAG,GAAG,EAAE,CAAC,CAAC,GAAG,EAAE,CAAE,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;oBAC3C,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAE,GAAG,MAAM,GAAG,GAAG,CAAC,CAAC,CAAE,CAAC;oBAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAE,GAAG,MAAM,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAE,CAAC;oBAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAE,GAAG,MAAM,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAE,CAAC;gBAC9I,CAAC;YACH,CAAC;QACH,CAAC;QACD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,SAAS,EAAE,CAAC,EAAE,EAAE,CAAC;YACvC,MAAM,CAAC,GAAG,IAAI,CAAC,CAAC,CAAE,CAAC;YACnB,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC;gBAAE,MAAM,IAAI,4BAA4B,CAAC,sDAAsD,CAAC,CAAC;YACxH,IAAI,CAAC,GAAG,CAAC;gBAAE,QAAQ,EAAE,CAAC;iBAAM,IAAI,CAAC,GAAG,CAAC;gBAAE,KAAK,EAAE,CAAC;QACjD,CAAC;QACD,MAAM,CAAC,GAAG,CAAC,CAAC;IACd,CAAC;IACD,OAAO,EAAE,mBAAmB,EAAE,QAAQ,EAAE,qBAAqB,EAAE,KAAK,EAAE,CAAC;AACzE,CAAC;AAED,2EAA2E;AAC3E,MAAM,UAAU,iBAAiB,CAAC,YAAoB;IACpD,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,YAAY,CAAC,EAAE,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;IACvE,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,KAAK,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAAC,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,YAAY,GAAG,YAAY,CAAC,CAAC,CAAC;QAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;QAAC,KAAK,IAAI,CAAC,CAAC;IAAC,CAAC;IAC3H,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,CAAC,EAAE;QAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAE,GAAG,KAAK,CAAC;IACxD,OAAO,CAAC,CAAC;AACX,CAAC;AAgBD;;;;GAIG;AACH,MAAM,UAAU,gBAAgB,CAAC,OAAqB,EAAE,KAAa,EAAE,MAAc,EAAE,QAA6B;IAElH,MAAM,CAAC,GAAG,wBAAwB,CAAC,QAAQ,CAAC,CAAC;IAC7C,IAAI,OAAO,CAAC,MAAM,KAAK,KAAK,GAAG,MAAM,GAAG,CAAC;QAAE,MAAM,IAAI,yBAAyB,CAAC,oCAAoC,CAAC,CAAC;IACrH,MAAM,OAAO,GAAG,IAAI,YAAY,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC;IACjD,OAAO,EAAE,OAAO,EAAE,WAAW,EAAE,eAAe,CAAC,eAAe,CAAC,OAAO,EAAE,KAAK,EAAE,MAAM,EAAE,CAAC,EAAE,OAAO,CAAC,CAAC,EAAE,CAAC;AACxG,CAAC;AACD;;;;;;;GAOG;AACH,MAAM,SAAS,CAAC,CAAC,eAAe,CAAC,OAAqB,EAAE,KAAa,EAAE,MAAc,EAAE,QAA6B,EAAE,GAA4C,EAChK,KAAc,EAAE,UAAgC;IAChD,MAAM,CAAC,GAAG,iBAAiB,CAAC,QAAQ,CAAC,YAAY,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,IAAI,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,SAAS,GAAG,KAAK,GAAG,CAAC,EAAE,IAAI,GAAG,IAAI,YAAY,CAAC,IAAI,GAAG,SAAS,CAAC,CAAC;IAC/J,MAAM,MAAM,GAAG,QAAQ,CAAC,MAAM,EAAE,SAAS,GAAG,QAAQ,CAAC,SAAS,CAAC;IAC/D,IAAI,SAAS,GAAG,CAAC,EAAE,IAAI,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,IAAI,GAAG,CAAC,EAAE,QAAQ,GAAG,CAAC,EAAE,QAAQ,GAAG,CAAC,EAAE,KAAK,GAAG,CAAC,CAAC;IACtF,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QAChC,uDAAuD;QACvD,KAAK,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,EAAE,KAAK,IAAI,IAAI,EAAE,KAAK,EAAE,EAAE,CAAC;YACtE,UAAU,EAAE,CAAC,KAAK,CAAC,CAAC;YACpB,MAAM,MAAM,GAAG,KAAK,GAAG,SAAS,EAAE,IAAI,GAAG,CAAC,KAAK,GAAG,IAAI,CAAC,GAAG,SAAS,CAAC;YACpE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC/B,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC;gBACxB,KAAK,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC7B,MAAM,EAAE,GAAG,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAE,CAAC;oBAC/E,CAAC,IAAI,CAAC,GAAG,OAAO,CAAC,EAAE,CAAE,CAAC;oBAAC,CAAC,IAAI,CAAC,GAAG,OAAO,CAAC,EAAE,GAAG,CAAC,CAAE,CAAC;oBAAC,CAAC,IAAI,CAAC,GAAG,OAAO,CAAC,EAAE,GAAG,CAAC,CAAE,CAAC;gBAC9E,CAAC;gBACD,IAAI,CAAC,IAAI,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;gBAAC,IAAI,CAAC,IAAI,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;gBAAC,IAAI,CAAC,IAAI,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;YACjF,CAAC;QACH,CAAC;QACD,MAAM,CAAC,GAAG,CAAC,GAAG,SAAS,CAAC;QACxB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;YAC/B,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC;YACxB,KAAK,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC7B,MAAM,EAAE,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,SAAS,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAE,CAAC;gBAChG,CAAC,IAAI,CAAC,GAAG,IAAI,CAAC,EAAE,CAAE,CAAC;gBAAC,CAAC,IAAI,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAE,CAAC;gBAAC,CAAC,IAAI,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAE,CAAC;YACrE,CAAC;YACD,KAAK,IAAI,OAAO,GAAG,CAAC,EAAE,OAAO,GAAG,CAAC,EAAE,OAAO,EAAE,EAAE,CAAC;gBAC7C,MAAM,KAAK,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,OAAO,EAAE,KAAK,GAAG,OAAO,CAAC,KAAK,CAAE,EAAE,MAAM,GAAG,KAAK,GAAG,CAAC,OAAO,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,OAAO,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACzH,IAAI,CAAC,GAAG,KAAK,CAAC;gBACd,IAAI,MAAM,KAAK,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,GAAG,SAAS,EAAE,CAAC;oBAAC,IAAI,MAAM,KAAK,CAAC;wBAAE,IAAI,EAAE,CAAC;gBAAC,CAAC;qBAC1E,CAAC;oBACJ,MAAM,MAAM,GAAG,MAAM,GAAG,MAAM,CAAC;oBAC/B,CAAC,GAAG,KAAK,GAAG,MAAM,CAAC;oBACnB,SAAS,EAAE,CAAC;oBACZ,IAAI,MAAM,GAAG,QAAQ;wBAAE,QAAQ,GAAG,MAAM,CAAC;oBACzC,IAAI,CAAC,MAAM,GAAG,QAAQ;wBAAE,QAAQ,GAAG,CAAC,MAAM,CAAC;oBAC3C,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;wBAAC,CAAC,GAAG,CAAC,CAAC;wBAAC,GAAG,EAAE,CAAC;oBAAC,CAAC;yBAAM,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;wBAAC,CAAC,GAAG,CAAC,CAAC;wBAAC,IAAI,EAAE,CAAC;oBAAC,CAAC;gBACjE,CAAC;gBACD,GAAG,CAAC,KAAK,CAAC,GAAG,KAAK,KAAK,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC;YAC/D,CAAC;QACH,CAAC;QACD,MAAM,CAAC,GAAG,CAAC,CAAC;IACd,CAAC;IACD,OAAO,EAAE,kBAAkB,EAAE,CAAC,EAAE,oBAAoB,EAAE,SAAS,EAAE,wBAAwB,EAAE,IAAI,EAAE,qBAAqB,EAAE,GAAG;QACzH,sBAAsB,EAAE,IAAI,EAAE,eAAe,EAAE,QAAQ,EAAE,eAAe,EAAE,QAAQ,EAAE,CAAC;AACzF,CAAC;AAID;;;;GAIG;AACH,MAAM,UAAU,wBAAwB,CAAC,GAAiB,EAAE,KAAa,EAAE,MAAc,EAAE,OAA4B,EAAE,QAA6B;IAEpJ,IAAI,GAAG,CAAC,MAAM,KAAK,KAAK,GAAG,MAAM,GAAG,CAAC;QAAE,MAAM,IAAI,yBAAyB,CAAC,2BAA2B,CAAC,CAAC;IACxG,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,EAAE,QAAQ,EAAE,KAAK,EAAE,OAAO,EAAE,GAAG,EAAE,IAAI,EAAE,EAAE,UAAU,EAAE,GAC1E,eAAe,CAAC,+BAA+B,CAAC,YAAY,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,KAAK,EAAE,MAAM,EAAE,OAAO,EAAE,wBAAwB,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;IACvI,OAAO,EAAE,KAAK,EAAE,MAAM,EAAE,EAAE,QAAQ,EAAE,KAAK,EAAE,OAAO,EAAE,GAAG,EAAE,IAAI,EAAE,EAAE,UAAU,EAAE,CAAC;AAChF,CAAC;AACD;;;;GAIG;AACH,MAAM,SAAS,CAAC,CAAC,+BAA+B,CAAC,GAAiB,EAAE,KAAa,EAAE,MAAc,EAAE,OAA4B,EAAE,QAA6B;IAE5J,MAAM,KAAK,GAAG,CAAC,IAAI,OAAO,CAAC,mBAAmB,EAAE,IAAI,GAAG,mBAAmB,CAAC,OAAO,CAAC,EAAE,OAAO,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;IACnH,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,IAAI,KAAK,IAAI,CAAC;QAAE,MAAM,IAAI,yBAAyB,CAAC,qCAAqC,CAAC,CAAC;IACtH,MAAM,MAAM,GAAmB,EAAE,QAAQ,EAAE,CAAC,EAAE,KAAK,EAAE,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,IAAI,EAAE,CAAC,EAAE,UAAU,EAAE,KAAK,EAAE,CAAC;IACzG,sHAAsH;IACtH,MAAM,SAAS,GAAG,CAAC,CAAS,EAAQ,EAAE;QACpC,gBAAgB,CAAC,GAAG,EAAE,CAAC,GAAG,KAAK,GAAG,CAAC,EAAE,KAAK,EAAE,OAAO,EAAE,IAAI,EAAE,KAAK,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,KAAK,GAAG,CAAC,EAAE,OAAO,EAAE,MAAM,CAAC,CAAC;QAC1G,IAAI,MAAM,CAAC,UAAU;YAAE,MAAM,IAAI,yBAAyB,CAAC,mDAAmD,CAAC,CAAC;IAClH,CAAC,CAAC;IACF,MAAM,KAAK,GAAG,OAAO,CAAC,QAAQ,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,UAAU,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,IAAI,WAAW,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC;IAChG,MAAM,UAAU,GAAG,KAAK,CAAC,CAAC,eAAe,CAAC,GAAG,EAAE,KAAK,EAAE,MAAM,EAAE,QAAQ,EAAE,KAAK,EAAE,CAAC,IAAI,OAAO,CAAC,QAAQ,GAAG,CAAC,EAAE,SAAS,CAAC,CAAC;IACrH,OAAO,EAAE,KAAK,EAAE,MAAM,EAAE,UAAU,EAAE,CAAC;AACvC,CAAC","sourcesContent":["// SPDX-License-Identifier: Apache-2.0\n\n/**\n * Development detail stages (#552), both opt-in and off by default:\n *\n * - **Diffraction correction**: a bounded, regularized (Wiener–Tikhonov) deconvolution of the developed linear RGB with the\n *   diffraction OTF of an ideal circular pupil at the capture's working f-number. The OTF is zero beyond the cut-off\n *   f_c = 1/(λ·N_w), so detail beyond it never reached the sensor; the filter's design gain there is 0 and it cannot, and\n *   does not claim to, restore it.\n * - **Sharpening**: an unsharp mask (radius, amount, threshold) on the display-encoded values, after the tone curve and\n *   gamut handling and before quantization. It raises edge contrast and creates halos; it adds no optical detail.\n *\n * Development order with these stages: (exposure gain) → demosaic → colour and white balance → diffraction correction\n * (linear) → rendering exposure and tone curve → gamut handling → sRGB transfer → unsharp mask → quantization (encode).\n *\n * Image edges are handled by clamping (edge pixels repeated). Both stages need the whole image; since #571 they also run row by\n * row in place (`correctDiffractionInPlace`, `renderSdrWithUnsharpMaskInPlace`), with the same sums in the same order.\n * @see docs/DIGITAL_OPTICS_FOUNDATION.md for the equations, sources, trade-offs and limits.\n */\n\nimport { InvalidConfigurationError } from \"../core/configuration-error.js\";\nimport { InvalidScientificResultError } from \"../core/validation.js\";\nimport { renderSdrRunInto, type SdrPixelCounts } from \"../output/native-capture-sdr.js\";\nimport { prepareSdrPixelTone, type SdrRenderingProfile } from \"../output/sdr-rendering.js\";\n\n/** Version of the `development-detail` model (diffraction correction and unsharp mask). */\nexport const DEVELOPMENT_DETAIL_MODEL_VERSION = \"0.1.0\" as const;\n\n/** Unsharp mask settings: Gaussian standard deviation in output pixels, gain on the detail signal, and threshold in encoded units [0, 1]. */\nexport interface UnsharpMaskSettings { radiusPixels: number; amount: number; threshold: number }\n/** Diffraction correction settings: design wavelength, Wiener noise-to-signal regularization, gain bound and FIR radius in pixels. */\nexport interface DiffractionCorrectionSettings { wavelengthNm: number; regularization: number; maximumGain: number; kernelRadiusPixels: number }\n\n/** Accepted ranges. Bounds are engineering limits on cost and on noise gain, not quality recommendations. */\nexport const DEVELOPMENT_DETAIL_LIMITS = /* @__PURE__ */ Object.freeze({\n  maximumUnsharpRadiusPixels: 8, maximumUnsharpAmount: 5,\n  minimumWavelengthNm: 380, maximumWavelengthNm: 780,\n  minimumRegularization: 1e-4, maximumRegularization: 1,\n  maximumGain: 10, maximumKernelRadiusPixels: 16\n});\n\n/** Third dark ring of the Airy pattern, in units of λ·N (j₁,₃/π, j₁,₃ = 10.1734681…). */\nconst AIRY_THIRD_ZERO_PER_LAMBDA_N = 10.173468135062722 / Math.PI;\n\nfunction finiteIn(value: unknown, name: string, minimum: number, maximum: number, minimumExclusive = false): number {\n  if (typeof value !== \"number\" || !Number.isFinite(value) || (minimumExclusive ? value <= minimum : value < minimum) || value > maximum) {\n    throw new InvalidConfigurationError(name + \" must be a finite number in \" + (minimumExclusive ? \"(\" : \"[\") + minimum + \", \" + maximum + \"].\");\n  }\n  return value;\n}\n\n/** Validates unsharp mask settings. */\nexport function parseUnsharpMaskSettings(value: UnsharpMaskSettings, path = \"sharpening\"): UnsharpMaskSettings {\n  return { radiusPixels: finiteIn(value.radiusPixels, path + \".radiusPixels\", 0, DEVELOPMENT_DETAIL_LIMITS.maximumUnsharpRadiusPixels, true),\n    amount: finiteIn(value.amount, path + \".amount\", 0, DEVELOPMENT_DETAIL_LIMITS.maximumUnsharpAmount, true),\n    threshold: finiteIn(value.threshold, path + \".threshold\", 0, 1) };\n}\n/** Validates diffraction correction settings. */\nexport function parseDiffractionCorrectionSettings(value: DiffractionCorrectionSettings, path = \"diffractionCorrection\"): DiffractionCorrectionSettings {\n  const radius = value.kernelRadiusPixels;\n  if (typeof radius !== \"number\" || !Number.isSafeInteger(radius) || radius < 1 || radius > DEVELOPMENT_DETAIL_LIMITS.maximumKernelRadiusPixels) {\n    throw new InvalidConfigurationError(path + \".kernelRadiusPixels must be an integer from 1 to \" + DEVELOPMENT_DETAIL_LIMITS.maximumKernelRadiusPixels + \".\");\n  }\n  return { wavelengthNm: finiteIn(value.wavelengthNm, path + \".wavelengthNm\", DEVELOPMENT_DETAIL_LIMITS.minimumWavelengthNm, DEVELOPMENT_DETAIL_LIMITS.maximumWavelengthNm),\n    regularization: finiteIn(value.regularization, path + \".regularization\", DEVELOPMENT_DETAIL_LIMITS.minimumRegularization, DEVELOPMENT_DETAIL_LIMITS.maximumRegularization),\n    maximumGain: finiteIn(value.maximumGain, path + \".maximumGain\", 1, DEVELOPMENT_DETAIL_LIMITS.maximumGain), kernelRadiusPixels: radius };\n}\n/** Default FIR radius: the Airy pattern's third dark ring (3.2383·λ·N_w) in pixels, rounded up, within 1..16. */\nexport function defaultDiffractionCorrectionRadiusPixels(workingFNumber: number, wavelengthNm: number, pixelPitchMicrometers: number): number {\n  const sites = Math.ceil(AIRY_THIRD_ZERO_PER_LAMBDA_N * wavelengthNm / 1000 * workingFNumber / pixelPitchMicrometers);\n  return Math.min(DEVELOPMENT_DETAIL_LIMITS.maximumKernelRadiusPixels, Math.max(1, sites));\n}\n\n/**\n * Incoherent OTF of an ideal, aberration-free circular pupil at normalized frequency s = f/f_c:\n * H(s) = (2/π)(acos s − s·√(1 − s²)) for s < 1, and 0 for s ≥ 1 (MIT OCW 2.71, lectures 21–22; Goodman, §7.3).\n */\nexport function circularPupilDiffractionMtf(normalizedFrequency: number): number {\n  const s = Math.abs(normalizedFrequency);\n  if (s >= 1) return 0;\n  return 2 / Math.PI * (Math.acos(s) - s * Math.sqrt(1 - s * s));\n}\n\n/** A designed diffraction correction filter and what its response does. Frequencies are in cycles per pixel. */\nexport interface DiffractionCorrectionFilter {\n  settings: DiffractionCorrectionSettings;\n  workingFNumber: number;\n  pixelPitchMicrometers: number;\n  /** f_c · pitch = pitch/(λ·N_w): the diffraction cut-off. The sampling Nyquist frequency is 0.5. */\n  cutoffCyclesPerPixel: number;\n  /** Fraction of the Nyquist square |fx|, |fy| ≤ 0.5 lying beyond the cut-off, where the lens transmitted nothing. */\n  nyquistAreaBeyondCutoffFraction: number;\n  /** Largest gain of the designed (untruncated) response: min(maximumGain, max over s of H(1+α)/(H² + α)). */\n  designPeakGain: number;\n  /** Largest |response| of the windowed, truncated FIR over a 1/64-cycle grid of the Nyquist square. */\n  responsePeakGain: number;\n  /**\n   * Largest |response| of the FIR at grid frequencies beyond the cut-off (null when the whole Nyquist square is inside it).\n   * The design is 0 there; a non-zero value is truncation leakage that passes noise and aliased signal, never lost detail,\n   * because the optics delivered no signal at those frequencies.\n   */\n  responseMaximumBeyondCutoff: number | null;\n  /** Row-major (2R+1)² weights summing to 1; weights[(dy + R)(2R + 1) + dx + R]. */\n  weights: Float64Array;\n}\n\n/**\n * Designs the bounded Wiener–Tikhonov correction. In frequency (normalized s = ρ/f_c, ρ the radial frequency):\n *   G(s) = min(maximumGain, H(s)·(1 + α)/(H(s)² + α)),   G(0) = 1, and G(s) = 0 for s ≥ 1,\n * where α is the regularization (an assumed noise-to-signal power ratio, Wiener's K). The response is sampled on an M × M\n * grid (M a power of two ≥ max(64, 4(2R + 1))), inverse transformed, truncated to the (2R + 1)² window with a radial raised\n * cosine taper and renormalized to unit sum. The achieved response is then measured, not assumed.\n */\nexport function designDiffractionCorrectionFilter(input: { workingFNumber: number; pixelPitchMicrometers: number; settings: DiffractionCorrectionSettings }): DiffractionCorrectionFilter {\n  const settings = parseDiffractionCorrectionSettings(input.settings);\n  const workingFNumber = finiteIn(input.workingFNumber, \"workingFNumber\", 0, 1e4, true), pitch = finiteIn(input.pixelPitchMicrometers, \"pixelPitchMicrometers\", 0, 1e4, true);\n  const cutoff = pitch / (settings.wavelengthNm / 1000 * workingFNumber), alpha = settings.regularization, bound = settings.maximumGain;\n  const gainAt = (rho: number): number => {\n    const h = circularPupilDiffractionMtf(rho / cutoff);\n    return h === 0 ? 0 : Math.min(bound, h * (1 + alpha) / (h * h + alpha));\n  };\n  const R = settings.kernelRadiusPixels, size = 2 * R + 1;\n  let M = 64;\n  while (M < 4 * size) M *= 2;\n  // A(ky, x) = Σ_kx G cos(2π kx x/M); h(x, y) = (1/M²) Σ_ky cos(2π ky y/M) A(ky, x). G is even in each axis, so h is real and even.\n  const cosTable = new Float64Array(M * (R + 1));\n  for (let k = 0; k < M; k++) for (let x = 0; x <= R; x++) cosTable[k * (R + 1) + x] = Math.cos(2 * Math.PI * k * x / M);\n  const A = new Float64Array(M * (R + 1));\n  for (let ky = 0; ky < M; ky++) {\n    const fy = (ky < M / 2 ? ky : ky - M) / M;\n    for (let kx = 0; kx < M; kx++) {\n      const fx = (kx < M / 2 ? kx : kx - M) / M, g = gainAt(Math.hypot(fx, fy));\n      if (g === 0) continue;\n      for (let x = 0; x <= R; x++) A[ky * (R + 1) + x]! += g * cosTable[kx * (R + 1) + x]!;\n    }\n  }\n  const quarter = new Float64Array((R + 1) * (R + 1));\n  for (let y = 0; y <= R; y++) for (let x = 0; x <= R; x++) {\n    let sum = 0;\n    for (let ky = 0; ky < M; ky++) sum += cosTable[ky * (R + 1) + y]! * A[ky * (R + 1) + x]!;\n    const r = Math.hypot(x, y), taper = r < R + 1 ? 0.5 * (1 + Math.cos(Math.PI * r / (R + 1))) : 0;\n    quarter[y * (R + 1) + x] = sum / (M * M) * taper;\n  }\n  const weights = new Float64Array(size * size);\n  let total = 0;\n  for (let dy = -R; dy <= R; dy++) for (let dx = -R; dx <= R; dx++) {\n    const w = quarter[Math.abs(dy) * (R + 1) + Math.abs(dx)]!;\n    weights[(dy + R) * size + dx + R] = w; total += w;\n  }\n  if (!(total > 0) || !Number.isFinite(total)) throw new InvalidScientificResultError(\"Diffraction correction filter has no finite DC response.\");\n  for (let i = 0; i < weights.length; i++) weights[i] = weights[i]! / total;\n  let designPeakGain = 0;\n  for (let i = 0; i <= 4096; i++) designPeakGain = Math.max(designPeakGain, gainAt(i / 4096 * cutoff));\n  let responsePeakGain = 0, beyond: number | null = null, beyondCount = 0;\n  const GRID = 64;\n  for (let iy = 0; iy <= GRID / 2; iy++) for (let ix = 0; ix <= GRID / 2; ix++) {\n    const fx = ix / GRID, fy = iy / GRID, response = Math.abs(evaluateDiffractionCorrectionResponse(weights, R, fx, fy));\n    responsePeakGain = Math.max(responsePeakGain, response);\n    if (Math.hypot(fx, fy) >= cutoff) beyond = Math.max(beyond ?? 0, response);\n  }\n  const AREA = 256;\n  for (let iy = 0; iy < AREA; iy++) for (let ix = 0; ix < AREA; ix++) if (Math.hypot((ix + 0.5) / (2 * AREA), (iy + 0.5) / (2 * AREA)) >= cutoff) beyondCount++;\n  return { settings, workingFNumber, pixelPitchMicrometers: pitch, cutoffCyclesPerPixel: cutoff, nyquistAreaBeyondCutoffFraction: beyondCount / (AREA * AREA),\n    designPeakGain, responsePeakGain, responseMaximumBeyondCutoff: beyond, weights };\n}\n/** Frequency response Σ w(x, y)·cos(2π(fx·x + fy·y)) of a centred (2R + 1)² FIR at (fx, fy) cycles per pixel. */\nexport function evaluateDiffractionCorrectionResponse(weights: Float64Array, radius: number, fx: number, fy: number): number {\n  const size = 2 * radius + 1;\n  let sum = 0;\n  for (let dy = -radius; dy <= radius; dy++) for (let dx = -radius; dx <= radius; dx++) sum += weights[(dy + radius) * size + dx + radius]! * Math.cos(2 * Math.PI * (fx * dx + fy * dy));\n  return sum;\n}\n\n/** What the correction did to the linear samples. */\nexport interface DiffractionCorrectionDiagnostics {\n  /** Corrected samples below 0 (ringing below black). They are clipped by the later gamut handling. */\n  negativeSampleCount: number;\n  /** Corrected samples above 1 (above the RAW white level). */\n  aboveWhiteSampleCount: number;\n}\n/**\n * Applies the FIR to interleaved linear RGB (width × height × 3), each channel separately, edges clamped. Returns new samples.\n */\nexport function applyDiffractionCorrection(rgb: Float64Array, width: number, height: number, filter: Pick<DiffractionCorrectionFilter, \"weights\" | \"settings\">):\n  { samples: Float64Array; diagnostics: DiffractionCorrectionDiagnostics } {\n  if (rgb.length !== width * height * 3) throw new InvalidConfigurationError(\"Diffraction correction raster size mismatch.\");\n  const samples = Float64Array.from(rgb);\n  return { samples, diagnostics: drainDetailRows(correctDiffractionInPlace(samples, width, height, filter)) };\n}\n/** Runs a row-by-row detail stage to its end and returns its result. */\nexport function drainDetailRows<T>(rows: Generator<number, T, void>): T {\n  let step = rows.next();\n  while (!step.done) step = rows.next();\n  return step.value;\n}\n/**\n * `applyDiffractionCorrection` in place, one output row per step (#571); each step yields the number of rows done. Output row y\n * reads input rows y − R … y + R, so it waits in a ring of R + 1 rows and replaces input row y once row y + R is done. Every sample\n * is the same sum in the same order as in a separate output image, so the result is bit-identical, with one image in memory\n * instead of two. `rgb` must hold width × height pixels.\n */\nexport function* correctDiffractionInPlace(rgb: Float64Array, width: number, height: number, filter: Pick<DiffractionCorrectionFilter, \"weights\" | \"settings\">):\n  Generator<number, DiffractionCorrectionDiagnostics, void> {\n  const R = filter.settings.kernelRadiusPixels, size = 2 * R + 1, w = filter.weights, rowLength = width * 3, ring = new Float64Array((R + 1) * rowLength);\n  const xs = new Int32Array(width + 2 * R);\n  for (let i = 0; i < xs.length; i++) xs[i] = Math.min(width - 1, Math.max(0, i - R)) * 3;\n  let negative = 0, above = 0;\n  for (let y = 0; y <= height + R; y++) {\n    // Output row y − R − 1 holds the slot of row y; input row y − R − 1 is read by no later output row, so it takes its place.\n    const o = (y % (R + 1)) * rowLength;\n    if (y > R) rgb.set(ring.subarray(o, o + rowLength), (y - R - 1) * rowLength);\n    if (y >= height) continue;\n    ring.fill(0, o, o + rowLength);\n    for (let dy = -R; dy <= R; dy++) {\n      const row = Math.min(height - 1, Math.max(0, y + dy)) * rowLength, wr = (dy + R) * size;\n      for (let dx = 0; dx < size; dx++) {\n        const weight = w[wr + dx]!;\n        if (weight === 0) continue;\n        for (let x = 0; x < width; x++) {\n          const s = row + xs[x + dx]!, d = o + 3 * x;\n          ring[d] = ring[d]! + weight * rgb[s]!; ring[d + 1] = ring[d + 1]! + weight * rgb[s + 1]!; ring[d + 2] = ring[d + 2]! + weight * rgb[s + 2]!;\n        }\n      }\n    }\n    for (let i = o; i < o + rowLength; i++) {\n      const v = ring[i]!;\n      if (!Number.isFinite(v)) throw new InvalidScientificResultError(\"Diffraction correction produced a non-finite sample.\");\n      if (v < 0) negative++; else if (v > 1) above++;\n    }\n    yield y + 1;\n  }\n  return { negativeSampleCount: negative, aboveWhiteSampleCount: above };\n}\n\n/** Normalized Gaussian weights exp(−i²/(2σ²)), i = −r..r with r = ⌈3σ⌉. */\nexport function unsharpMaskKernel(radiusPixels: number): Float64Array {\n  const r = Math.ceil(3 * radiusPixels), k = new Float64Array(2 * r + 1);\n  let total = 0;\n  for (let i = -r; i <= r; i++) { const v = Math.exp(-i * i / (2 * radiusPixels * radiusPixels)); k[i + r] = v; total += v; }\n  for (let i = 0; i < k.length; i++) k[i] = k[i]! / total;\n  return k;\n}\n/** What the unsharp mask did. Counts are samples (one per channel). */\nexport interface UnsharpMaskDiagnostics {\n  /** Half-width of the separable Gaussian, ⌈3σ⌉ pixels. */\n  kernelRadiusPixels: number;\n  /** Samples changed by the mask (|detail| ≥ threshold and detail ≠ 0). */\n  sharpenedSampleCount: number;\n  /** Samples left unchanged because 0 < |detail| < threshold. */\n  thresholdHeldSampleCount: number;\n  /** Sharpened samples clipped at 0 or 1 (the halo did not fit in the code range). */\n  clippedLowSampleCount: number;\n  clippedHighSampleCount: number;\n  /** Largest increase and decrease amount·detail before clipping, in encoded units: the halo strength. */\n  maximumIncrease: number;\n  maximumDecrease: number;\n}\n/**\n * Unsharp mask on interleaved encoded RGB in [0, 1] (each channel separately, edges clamped):\n *   d = x − (G_σ ∗ x);   y = x + amount·d where |d| ≥ threshold, else y = x;   then y is clipped to [0, 1].\n * Returns new samples.\n */\nexport function applyUnsharpMask(encoded: Float64Array, width: number, height: number, settings: UnsharpMaskSettings):\n  { samples: Float64Array; diagnostics: UnsharpMaskDiagnostics } {\n  const s = parseUnsharpMaskSettings(settings);\n  if (encoded.length !== width * height * 3) throw new InvalidConfigurationError(\"Unsharp mask raster size mismatch.\");\n  const samples = new Float64Array(encoded.length);\n  return { samples, diagnostics: drainDetailRows(unsharpMaskRows(encoded, width, height, s, samples)) };\n}\n/**\n * The unsharp mask one output row per step (#571); each step yields the number of rows done. The horizontal Gaussian pass of\n * each input row is kept in a ring of 2r + 1 rows, and `prepareRow(y)`, when given, runs once on input row y before it is\n * first read (`renderSdrWithUnsharpMaskInPlace` encodes the row there). Output row y goes to `out`: the clipped value, or\n * Math.round(value · scale) with a `scale`. Every sum is the same, in the same order, as in whole-image passes, so the\n * result is bit-identical, with 2r + 1 rows of working memory instead of three images. `settings` must be parsed and the\n * raster sizes must match.\n */\nexport function* unsharpMaskRows(encoded: Float64Array, width: number, height: number, settings: UnsharpMaskSettings, out: Float64Array | Uint8Array | Uint16Array,\n  scale?: number, prepareRow?: (y: number) => void): Generator<number, UnsharpMaskDiagnostics, void> {\n  const k = unsharpMaskKernel(settings.radiusPixels), r = (k.length - 1) / 2, span = 2 * r + 1, rowLength = width * 3, ring = new Float64Array(span * rowLength);\n  const amount = settings.amount, threshold = settings.threshold;\n  let sharpened = 0, held = 0, low = 0, high = 0, increase = 0, decrease = 0, ready = 0;\n  for (let y = 0; y < height; y++) {\n    // Rows up to y + r, prepared and horizontally blurred.\n    for (const last = Math.min(height - 1, y + r); ready <= last; ready++) {\n      prepareRow?.(ready);\n      const source = ready * rowLength, slot = (ready % span) * rowLength;\n      for (let x = 0; x < width; x++) {\n        let a = 0, b = 0, c = 0;\n        for (let i = -r; i <= r; i++) {\n          const sx = source + Math.min(width - 1, Math.max(0, x + i)) * 3, w = k[i + r]!;\n          a += w * encoded[sx]!; b += w * encoded[sx + 1]!; c += w * encoded[sx + 2]!;\n        }\n        ring[slot + 3 * x] = a; ring[slot + 3 * x + 1] = b; ring[slot + 3 * x + 2] = c;\n      }\n    }\n    const o = y * rowLength;\n    for (let x = 0; x < width; x++) {\n      let a = 0, b = 0, c = 0;\n      for (let i = -r; i <= r; i++) {\n        const sy = (Math.min(height - 1, Math.max(0, y + i)) % span) * rowLength + 3 * x, w = k[i + r]!;\n        a += w * ring[sy]!; b += w * ring[sy + 1]!; c += w * ring[sy + 2]!;\n      }\n      for (let channel = 0; channel < 3; channel++) {\n        const index = o + 3 * x + channel, value = encoded[index]!, detail = value - (channel === 0 ? a : channel === 1 ? b : c);\n        let v = value;\n        if (detail === 0 || Math.abs(detail) < threshold) { if (detail !== 0) held++; }\n        else {\n          const change = amount * detail;\n          v = value + change;\n          sharpened++;\n          if (change > increase) increase = change;\n          if (-change > decrease) decrease = -change;\n          if (v < 0) { v = 0; low++; } else if (v > 1) { v = 1; high++; }\n        }\n        out[index] = scale === undefined ? v : Math.round(v * scale);\n      }\n    }\n    yield y + 1;\n  }\n  return { kernelRadiusPixels: r, sharpenedSampleCount: sharpened, thresholdHeldSampleCount: held, clippedLowSampleCount: low,\n    clippedHighSampleCount: high, maximumIncrease: increase, maximumDecrease: decrease };\n}\n\n/** SDR counters in the shape of the SDR rendering diagnostics. */\nexport interface DetailSdrCounts { negative: number; above: number; changed: number; low: number; high: number }\n/**\n * SDR rendering with an unsharp mask between the sRGB transfer and quantization. Per pixel it runs the reference renderer's\n * operations in its order (−0 → +0, rendering exposure 2^EV with reference white 1, the shared tone operator, gamut\n * handling, sRGB transfer); the whole encoded image is then sharpened and rounded to nearest (ties up).\n */\nexport function renderSdrWithUnsharpMask(rgb: Float64Array, width: number, height: number, profile: SdrRenderingProfile, settings: UnsharpMaskSettings):\n  { codes: Uint8Array | Uint16Array; counts: DetailSdrCounts; sharpening: UnsharpMaskDiagnostics } {\n  if (rgb.length !== width * height * 3) throw new InvalidConfigurationError(\"SDR raster size mismatch.\");\n  const { codes, counts: { negative, above, changed, low, high }, sharpening } =\n    drainDetailRows(renderSdrWithUnsharpMaskInPlace(Float64Array.from(rgb), width, height, profile, parseUnsharpMaskSettings(settings)));\n  return { codes, counts: { negative, above, changed, low, high }, sharpening };\n}\n/**\n * `renderSdrWithUnsharpMask` one output row per step (#571), encoding each row of `rgb` in place just before the mask first reads\n * it (`rgb` then holds the encoded values). Bit-identical to the whole-image function, with no second image in memory. `settings`\n * must be parsed and `rgb` must hold width × height pixels.\n */\nexport function* renderSdrWithUnsharpMaskInPlace(rgb: Float64Array, width: number, height: number, profile: SdrRenderingProfile, settings: UnsharpMaskSettings):\n  Generator<number, { codes: Uint8Array | Uint16Array; counts: SdrPixelCounts; sharpening: UnsharpMaskDiagnostics }, void> {\n  const scale = 2 ** profile.renderingExposureEv, tone = prepareSdrPixelTone(profile), scratch = new Float64Array(3);\n  if (!Number.isFinite(scale) || scale <= 0) throw new InvalidConfigurationError(\"Rendering scale is unrepresentable.\");\n  const counts: SdrPixelCounts = { negative: 0, above: 0, changed: 0, low: 0, high: 0, outOfGamut: false };\n  // The native SDR kernel per row, unquantized (maximum 0), in place: the reference renderer's operations in its order.\n  const encodeRow = (y: number): void => {\n    renderSdrRunInto(rgb, y * width * 3, width, profile, tone, scale, 0, rgb, y * width * 3, scratch, counts);\n    if (counts.outOfGamut) throw new InvalidConfigurationError(\"Tone output exceeds the declared SDR gamut range.\");\n  };\n  const codes = profile.bitDepth === 8 ? new Uint8Array(rgb.length) : new Uint16Array(rgb.length);\n  const sharpening = yield* unsharpMaskRows(rgb, width, height, settings, codes, 2 ** profile.bitDepth - 1, encodeRow);\n  return { codes, counts, sharpening };\n}\n"]}

package/dist/schema/illumination-temporal.js

@@ -12,48 +12,10 @@
 import { InvalidConfigurationError } from "../core/configuration-error.js";
 import { parseEvidenceList } from "../core/evidence-provenance.js";
 import { InvalidScientificInputError, requirePositiveInteger } from "../core/validation.js";
-import { requireInputKeys, rejectUnknownFields, rejectUnknownVariantFields, UNCERTAINTY_VARIANT_FIELDS, requireFixedFieldsWhenPresent } from "../core/record-validation.js";
+import { requireInputKeys, rejectUnknownVariantFields, requireFixedFieldsWhenPresent, requireConfigRecord as requireRecord, requireTrimmedConfigString as requireNonEmptyString, requireConfigFinite as requireFinite, parseConfigScientificStatus as parseScientificStatus, requireConfigFiniteAboveZero as requirePositiveFinite, requireConfigFiniteAtLeastZero as requireNonNegativeFinite } from "../core/record-validation.js";
+import { parseUncertainty } from "../optics/lens-psf-profile.js";
+import { positiveModulo } from "../optics/complex-pupil-psf.js";
 export const SCENE_ILLUMINATION_TEMPORAL_PROFILE_SCHEMA_VERSION = "0.1.0";
-function requireRecord(value, path, allowedKeys) {
-    if (typeof value !== "object" ||
-        value === null ||
-        Array.isArray(value)) {
-        throw new InvalidConfigurationError(path + " must be an object.");
-    }
-    if (allowedKeys !== null)
-        rejectUnknownFields(value, allowedKeys, path);
-    return value;
-}
-function requireNonEmptyString(value, path) {
-    if (typeof value !== "string" ||
-        value.trim().length === 0) {
-        throw new InvalidConfigurationError(path + " must be a non-empty string.");
-    }
-    return value.trim();
-}
-function requireFinite(value, path) {
-    if (typeof value !== "number" ||
-        !Number.isFinite(value)) {
-        throw new InvalidConfigurationError(path + " must be finite.");
-    }
-    return value;
-}
-function requirePositiveFinite(value, path) {
-    const parsed = requireFinite(value, path);
-    if (parsed <= 0) {
-        throw new InvalidConfigurationError(path +
-            " must be greater than zero.");
-    }
-    return parsed;
-}
-function requireNonNegativeFinite(value, path) {
-    const parsed = requireFinite(value, path);
-    if (parsed < 0) {
-        throw new InvalidConfigurationError(path +
-            " must be greater than or equal to zero.");
-    }
-    return parsed;
-}
 function requireReusableNumericEvidence(value, path) {
     const evidence = parseEvidenceList(value, path);
     if (!evidence.some((entry) => entry.reuseStatus ===
@@ -65,31 +27,6 @@ function requireReusableNumericEvidence(value, path) {
     }
     return evidence;
 }
-function parseScientificStatus(value, path) {
-    if (value !== "calibrated" &&
-        value !== "approximation") {
-        throw new InvalidConfigurationError(path + " is invalid.");
-    }
-    return value;
-}
-function parseUncertainty(value, path) {
-    const record = requireRecord(value, path, null);
-    rejectUnknownVariantFields(record, "kind", UNCERTAINTY_VARIANT_FIELDS, path);
-    if (record.kind === "relative") {
-        return {
-            kind: "relative",
-            fraction: requireNonNegativeFinite(record.fraction, path + ".fraction"),
-            basis: requireNonEmptyString(record.basis, path + ".basis")
-        };
-    }
-    if (record.kind === "not-quantified") {
-        return {
-            kind: "not-quantified",
-            limitation: requireNonEmptyString(record.limitation, path + ".limitation")
-        };
-    }
-    throw new InvalidConfigurationError(path + ".kind is invalid.");
-}
 function parseRegistrationUncertainty(value, path) {
     const record = requireRecord(value, path, null);
     rejectUnknownVariantFields(record, "kind", { "absolute-seconds": ["kind", "plusMinusSeconds", "basis"], "not-quantified": ["kind", "limitation"] }, path);
@@ -394,12 +331,6 @@ function interpolateWaveform(samples, timeSeconds) {
     return last
         .relativeMagnitudeMultiplier;
 }
-function positiveModulo(value, modulus) {
-    const remainder = value % modulus;
-    return remainder < 0
-        ? remainder + modulus
-        : remainder;
-}
 function evaluateWaveform(waveform, waveformTimeSeconds) {
     if (waveform.kind ===
         "aperiodic-relative-multiplier") {