@photivra/engine 1.4.0 → 2.0.0

package/dist/optics/complex-pupil-psf.js

@@ -1,31 +1,32 @@
 // SPDX-License-Identifier: Apache-2.0
 /**
- * Module boundary and integration notes.
- * Validate an untrusted declaration and return the normalized typed contract. Unknown enum values,
- * missing required fields and incompatible scientific data fail at this boundary. Complex pupil
- * evaluation sums declared amplitude and phase on a bounded regular pupil grid. PSF shape is
- * normalized separately from pupil throughput; coherent phase, pupil energy and spectral/field
- * applicability cannot be replaced by a scalar blur radius. Numeric profiles require owned or
- * reusable-data provenance.
- * Evaluate a declared complex pupil with explicit amplitude/phase and sampling support; normalize PSF
- * shape and retain pupil throughput as a separate quantity. Complex pupil evaluation sums declared
- * amplitude and phase on a bounded regular pupil grid. PSF shape is normalized separately from pupil
- * throughput; coherent phase, pupil energy and spectral/field applicability cannot be replaced by a
- * scalar blur radius. Numeric profiles require owned or reusable-data provenance.
- * @see docs/PSF_FOUNDATION.md for equations, coordinate/unit conventions, blockers and support limits.
+ * PSF (point-spread function) from a declared complex pupil: amplitude and phase sampled on a grid.
+ *
+ * The profile parser validates untrusted input and returns the normalized typed contract. Unknown
+ * enum values, missing required fields and incompatible scientific data are rejected.
+ *
+ * The complex pupil is evaluated by summing the declared amplitude and phase over a bounded, regular
+ * pupil grid. The PSF shape is normalized separately from the pupil throughput. A single blur radius
+ * cannot stand in for coherent phase, pupil energy, or the wavelength and field range where the
+ * profile applies. Numeric profiles need provenance showing the data is owned or reusable.
+ * @see docs/PSF_FOUNDATION.md for the equations, units, coordinate conventions, known blockers and
+ * support limits.
  */
 import { approximationResult } from "../core/calculation-result.js";
 import { InvalidConfigurationError } from "../core/configuration-error.js";
 import { parseEvidenceList } from "../core/evidence-provenance.js";
 import { InvalidScientificInputError } from "../core/validation.js";
 import { parseFocusPlane } from "./focus-state.js";
+import { requireInputKeys, rejectUnknownFields, rejectUnknownVariantFields, UNCERTAINTY_VARIANT_FIELDS } from "../core/record-validation.js";
 export const LENS_COMPLEX_PUPIL_PROFILE_SCHEMA_VERSION = "0.1.0";
-function requireRecord(value, path) {
+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) {
@@ -86,7 +87,8 @@ function parseStatus(value, path) {
     return value;
 }
 function parseUncertainty(value, path) {
-    const record = requireRecord(value, path);
+    const record = requireRecord(value, path, null);
+    rejectUnknownVariantFields(record, "kind", UNCERTAINTY_VARIANT_FIELDS, path);
     if (record.kind === "relative") {
         const fraction = requireFinite(record.fraction, path + ".fraction");
         if (fraction < 0) {
@@ -125,7 +127,7 @@ function parseLimitations(value, path) {
     return parsed;
 }
 function parseGrid(value) {
-    const record = requireRecord(value, "lensComplexPupil.grid");
+    const record = requireRecord(value, "lensComplexPupil.grid", ["widthSamples", "heightSamples", "centerSampleX", "centerSampleY", "relativeAmplitude", "opticalPathDifferenceMicrometers", "pupilSamplePitchMmX", "pupilSamplePitchMmY"]);
     const widthSamples = requirePositiveSafeInteger(record.widthSamples, "lensComplexPupil.grid.widthSamples");
     const heightSamples = requirePositiveSafeInteger(record.heightSamples, "lensComplexPupil.grid.heightSamples");
     const total = widthSamples *
@@ -179,20 +181,20 @@ function parseGrid(value) {
     };
 }
 /**
- * Validate an untrusted declaration and return the normalized typed contract. Unknown enum values,
- * missing required fields and incompatible scientific data fail at this boundary.
+ * Validates an untrusted declaration and returns the normalized typed contract. Unknown enum values,
+ * missing required fields and incompatible scientific data are rejected here.
  *
- * Complex pupil evaluation sums declared amplitude and phase on a bounded regular pupil grid. PSF
- * shape is normalized separately from pupil throughput; coherent phase, pupil energy and
- * spectral/field applicability cannot be replaced by a scalar blur radius. Numeric profiles require
- * owned or reusable-data provenance.
- * @param value - unknown. Treated as untrusted data; static typing alone is not validation.
- * @returns LensComplexPupilProfile. Return shape and scientific status are explicit; no calibration is inferred from successful execution.
+ * The complex pupil is evaluated by summing the declared amplitude and phase over a bounded, regular
+ * pupil grid. The PSF shape is normalized separately from the pupil throughput. A single blur radius
+ * cannot stand in for coherent phase, pupil energy, or the wavelength and field range where the
+ * profile applies. Numeric profiles need provenance showing the data is owned or reusable.
+ * @param value - unknown. Untrusted input; TypeScript types alone are not validation.
+ * @returns LensComplexPupilProfile. The result states its shape and scientific status; a successful call does not imply calibration.
  *
  * @see docs/API_REFERENCE.md for the root export and exact type graph.
  */
 export function parseLensComplexPupilProfile(value) {
-    const record = requireRecord(value, "lensComplexPupil");
+    const record = requireRecord(value, "lensComplexPupil", ["schemaVersion", "opticalDomain", "representation", "pupilCoordinateSystem", "imageFieldAxes", "amplitudeMeaning", "wavefrontMeaning", "throughputOwnership", "kernelEnergyNormalization", "propagationModel", "sensorOpticalStackIncluded", "sensorSamplingIncluded", "reconstructionIncluded", "strayLightIncluded", "scientificStatus", "uncertainty", "context", "responseIncludes", "profileId", "profileVersion", "relativePupilThroughputFactor", "grid", "evidence", "limitations"]);
     if (record.schemaVersion !==
         LENS_COMPLEX_PUPIL_PROFILE_SCHEMA_VERSION) {
         throw new InvalidConfigurationError('lensComplexPupil.schemaVersion must be "' +
@@ -236,9 +238,9 @@ export function parseLensComplexPupilProfile(value) {
             "not-quantified") {
         throw new InvalidConfigurationError("Calibrated complex-pupil profiles must declare quantified relative uncertainty.");
     }
-    const context = requireRecord(record.context, "lensComplexPupil.context");
-    const field = requireRecord(context.fieldPointMm, "lensComplexPupil.context.fieldPointMm");
-    const response = requireRecord(record.responseIncludes, "lensComplexPupil.responseIncludes");
+    const context = requireRecord(record.context, "lensComplexPupil.context", ["fieldPointMm", "focalLengthMm", "focus", "apertureFNumber", "wavelengthNm", "signedDefocusImagePlaneMicrometers"]);
+    const field = requireRecord(context.fieldPointMm, "lensComplexPupil.context.fieldPointMm", ["x", "y"]);
+    const response = requireRecord(record.responseIncludes, "lensComplexPupil.responseIncludes", ["diffraction", "aberration", "defocus", "pupilClippingShape"]);
     if (response.diffraction !==
         true ||
         typeof response.aberration !==
@@ -306,59 +308,151 @@ function pupilComplexValue(amplitude, opdMicrometers, wavelengthMicrometers) {
             Math.sin(phase)
     };
 }
-function calculateDiscreteFraunhoferKernel(profile, propagationDistanceMm) {
+function supportSpan(grid, axis) {
+    let minimum = Number.POSITIVE_INFINITY;
+    let maximum = Number.NEGATIVE_INFINITY;
+    grid.relativeAmplitude.forEach((amplitude, index) => {
+        if (amplitude > 0) {
+            const coordinate = axis === "x"
+                ? index % grid.widthSamples
+                : Math.floor(index / grid.widthSamples);
+            minimum = Math.min(minimum, coordinate);
+            maximum = Math.max(maximum, coordinate);
+        }
+    });
+    return maximum - minimum + 1;
+}
+/**
+ * Sampling guard for the discrete Fraunhofer reference (Ban et al.,
+ * arXiv:2511.06936, eqs. 6-9, free: N*dp = lambda*f/dc, Q = lambda*F/dc and
+ * Q >= 2. Also see: Goodman, Introduction to Fourier Optics, sampling of
+ * DFT/FFT Fraunhofer propagation).
+ *
+ * A DFT (discrete Fourier transform) of size M over pupil pitch dp gives an
+ * image-plane pitch of lambda*z/(M*dp).
+ *
+ * The intensity PSF is the transform of the pupil autocorrelation. That
+ * autocorrelation is twice as wide as the pupil. So Nyquist sampling of the
+ * intensity needs Q = lambda*(z/D)/pitch = M*dp/D >= 2.
+ *
+ * For a pupil spanning S samples, the discrete autocorrelation has 2S-1 lags.
+ * So M >= 2S-1 is the exact condition for no circular wrap-around, and it
+ * implies Q = M/(S-1) > 2. The declared grid is zero padded to that size.
+ *
+ * The non-zero support must also fit inside the declared f/N cone. Its
+ * diameter D may not exceed z/N for the resolved propagation distance z; that
+ * is, the image-side cone D/z may not be wider than 1/N. With the nominal
+ * focal-plane reference (z = f) this is the entrance-aperture diameter f/N. An
+ * exit pupil declared at its own distance z uses that distance
+ * (docs/PSF_FOUNDATION.md).
+ */
+function planPupilSampling(profile, propagationDistanceMm) {
+    const grid = profile.grid;
+    const supportX = supportSpan(grid, "x");
+    const supportY = supportSpan(grid, "y");
+    const diameterX = (supportX - 1) * grid.pupilSamplePitchMmX;
+    const diameterY = (supportY - 1) * grid.pupilSamplePitchMmY;
+    const declaredDiameter = propagationDistanceMm / profile.context.apertureFNumber;
+    // Relative allowance for floating-point rounding only.
+    const limit = declaredDiameter * (1 + 1e-9);
+    if (diameterX > limit || diameterY > limit) {
+        throw new InvalidScientificInputError("Complex-pupil amplitude support spans " +
+            diameterX.toPrecision(6) +
+            " x " +
+            diameterY.toPrecision(6) +
+            " mm between outermost non-zero samples, wider than the declared f/N cone diameter " +
+            declaredDiameter.toPrecision(6) +
+            " mm (propagationDistanceMm / apertureFNumber).");
+    }
+    const transformWidth = Math.max(grid.widthSamples, 2 * supportX - 1);
+    const transformHeight = Math.max(grid.heightSamples, 2 * supportY - 1);
+    return {
+        sampling: {
+            criterion: "discrete-intensity-alias-free-goodman-q-above-2",
+            pupilSupportSamplesX: supportX,
+            pupilSupportSamplesY: supportY,
+            pupilSupportDiameterMmX: diameterX,
+            pupilSupportDiameterMmY: diameterY,
+            declaredPupilDiameterMm: declaredDiameter,
+            pupilGridWidthSamples: grid.widthSamples,
+            pupilGridHeightSamples: grid.heightSamples,
+            transformWidthSamples: transformWidth,
+            transformHeightSamples: transformHeight,
+            zeroPaddingApplied: transformWidth !== grid.widthSamples ||
+                transformHeight !== grid.heightSamples,
+            samplingFactorQX: supportX > 1 ? transformWidth / (supportX - 1) : null,
+            samplingFactorQY: supportY > 1 ? transformHeight / (supportY - 1) : null
+        },
+        // Symmetric padding keeps zero frequency at the declared centre, shifted
+        // by the leading padding.
+        outputCenterX: grid.centerSampleX +
+            Math.floor((transformWidth - grid.widthSamples) / 2),
+        outputCenterY: grid.centerSampleY +
+            Math.floor((transformHeight - grid.heightSamples) / 2)
+    };
+}
+function twiddleTable(size) {
+    const cos = [];
+    const sin = [];
+    for (let index = 0; index < size; index += 1) {
+        const angle = (-2 * Math.PI * index) / size;
+        cos.push(Math.cos(angle));
+        sin.push(Math.sin(angle));
+    }
+    return { cos, sin };
+}
+function positiveModulo(value, size) {
+    const remainder = value % size;
+    return remainder < 0 ? remainder + size : remainder;
+}
+function calculateDiscreteFraunhoferKernel(profile, propagationDistanceMm, plan) {
     const grid = profile.grid;
     const width = grid.widthSamples;
     const height = grid.heightSamples;
-    const wavelengthMicrometers = profile.context
-        .wavelengthNm /
-        1000;
-    const wavelengthMm = profile.context
-        .wavelengthNm /
-        1_000_000;
-    const pupil = grid.relativeAmplitude.map((amplitude, index) => pupilComplexValue(amplitude, grid
-        .opticalPathDifferenceMicrometers[index], wavelengthMicrometers));
-    const intensity = new Array(width * height).fill(0);
-    for (let outY = 0; outY < height; outY += 1) {
-        const ky = outY -
-            grid.centerSampleY;
-        for (let outX = 0; outX < width; outX += 1) {
-            const kx = outX -
-                grid.centerSampleX;
+    const outWidth = plan.sampling.transformWidthSamples;
+    const outHeight = plan.sampling.transformHeightSamples;
+    const wavelengthMicrometers = profile.context.wavelengthNm / 1000;
+    const wavelengthMm = profile.context.wavelengthNm / 1_000_000;
+    const pupil = grid.relativeAmplitude.map((amplitude, index) => pupilComplexValue(amplitude, grid.opticalPathDifferenceMicrometers[index], wavelengthMicrometers));
+    // Separable zero-padded DFT: pupil rows to image kx first, then columns.
+    // Padding samples are zero, so only declared pupil samples are summed.
+    const twiddleX = twiddleTable(outWidth);
+    const twiddleY = twiddleTable(outHeight);
+    const rowRe = new Float64Array(height * outWidth);
+    const rowIm = new Float64Array(height * outWidth);
+    for (let y = 0; y < height; y += 1) {
+        for (let outX = 0; outX < outWidth; outX += 1) {
+            const kx = outX - plan.outputCenterX;
+            let re = 0;
+            let im = 0;
+            for (let x = 0; x < width; x += 1) {
+                const source = pupil[y * width + x];
+                const t = positiveModulo(kx * (x - grid.centerSampleX), outWidth);
+                const cos = twiddleX.cos[t];
+                const sin = twiddleX.sin[t];
+                re += source.re * cos - source.im * sin;
+                im += source.re * sin + source.im * cos;
+            }
+            rowRe[y * outWidth + outX] = re;
+            rowIm[y * outWidth + outX] = im;
+        }
+    }
+    const intensity = new Array(outWidth * outHeight).fill(0);
+    for (let outY = 0; outY < outHeight; outY += 1) {
+        const ky = outY - plan.outputCenterY;
+        for (let outX = 0; outX < outWidth; outX += 1) {
             let re = 0;
             let im = 0;
             for (let y = 0; y < height; y += 1) {
-                const py = y -
-                    grid.centerSampleY;
-                for (let x = 0; x < width; x += 1) {
-                    const px = x -
-                        grid.centerSampleX;
-                    const source = pupil[y * width +
-                        x];
-                    const angle = -2 *
-                        Math.PI *
-                        (kx * px /
-                            width +
-                            ky * py /
-                                height);
-                    const cos = Math.cos(angle);
-                    const sin = Math.sin(angle);
-                    re +=
-                        source.re *
-                            cos -
-                            source.im *
-                                sin;
-                    im +=
-                        source.re *
-                            sin +
-                            source.im *
-                                cos;
-                }
+                const sourceRe = rowRe[y * outWidth + outX];
+                const sourceIm = rowIm[y * outWidth + outX];
+                const t = positiveModulo(ky * (y - grid.centerSampleY), outHeight);
+                const cos = twiddleY.cos[t];
+                const sin = twiddleY.sin[t];
+                re += sourceRe * cos - sourceIm * sin;
+                im += sourceRe * sin + sourceIm * cos;
             }
-            intensity[outY * width +
-                outX] =
-                re * re +
-                    im * im;
+            intensity[outY * outWidth + outX] = re * re + im * im;
         }
     }
     const sum = intensity.reduce((acc, value) => acc + value, 0);
@@ -367,16 +461,10 @@ function calculateDiscreteFraunhoferKernel(profile, propagationDistanceMm) {
         throw new InvalidScientificInputError("Complex-pupil propagation must produce finite positive PSF energy.");
     }
     const normalized = intensity.map((value) => value / sum);
-    const samplePitchMmX = wavelengthMm *
-        propagationDistanceMm /
-        (width *
-            grid
-                .pupilSamplePitchMmX);
-    const samplePitchMmY = wavelengthMm *
-        propagationDistanceMm /
-        (height *
-            grid
-                .pupilSamplePitchMmY);
+    const samplePitchMmX = wavelengthMm * propagationDistanceMm /
+        (outWidth * grid.pupilSamplePitchMmX);
+    const samplePitchMmY = wavelengthMm * propagationDistanceMm /
+        (outHeight * grid.pupilSamplePitchMmY);
     if (!Number.isFinite(samplePitchMmX) ||
         !Number.isFinite(samplePitchMmY) ||
         samplePitchMmX <= 0 ||
@@ -384,35 +472,37 @@ function calculateDiscreteFraunhoferKernel(profile, propagationDistanceMm) {
         throw new InvalidScientificInputError("Complex-pupil PSF sample pitch must remain finite and positive.");
     }
     return {
-        widthSamples: width,
-        heightSamples: height,
+        widthSamples: outWidth,
+        heightSamples: outHeight,
         samplePitchMicrometersX: samplePitchMmX *
             1000,
         samplePitchMicrometersY: samplePitchMmY *
             1000,
-        centerSampleX: grid.centerSampleX,
-        centerSampleY: grid.centerSampleY,
+        centerSampleX: plan.outputCenterX,
+        centerSampleY: plan.outputCenterY,
         normalizedIntensity: normalized
     };
 }
 /**
- * Evaluate a declared complex pupil with explicit amplitude/phase and sampling support; normalize PSF
- * shape and retain pupil throughput as a separate quantity.
+ * Evaluates a declared complex pupil with explicit amplitude, phase and sampling support. Returns the
+ * normalized PSF shape, and the pupil throughput as a separate quantity.
  *
- * Complex pupil evaluation sums declared amplitude and phase on a bounded regular pupil grid. PSF
- * shape is normalized separately from pupil throughput; coherent phase, pupil energy and
- * spectral/field applicability cannot be replaced by a scalar blur radius. Numeric profiles require
- * owned or reusable-data provenance.
- * @param input - CalculateLensComplexPupilPsfInput. See the linked contract for coordinate, unit and profile binding semantics.
- * @returns CalculationResult<LensComplexPupilPsf>. Return shape and scientific status are explicit; no calibration is inferred from successful execution.
+ * The complex pupil is evaluated by summing the declared amplitude and phase over a bounded, regular
+ * pupil grid. The PSF shape is normalized separately from the pupil throughput. A single blur radius
+ * cannot stand in for coherent phase, pupil energy, or the wavelength and field range where the
+ * profile applies. Numeric profiles need provenance showing the data is owned or reusable.
+ * @param input - CalculateLensComplexPupilPsfInput. See that type for coordinates, units and profile binding.
+ * @returns CalculationResult<LensComplexPupilPsf>. The result states its shape and scientific status; a successful call does not imply calibration.
  *
  * @see docs/API_REFERENCE.md for the root export and exact type graph.
  */
 export function calculateLensComplexPupilPsf(input) {
+    requireInputKeys(input, "calculateLensComplexPupilPsf", ["profile", "propagationDistanceMm"]);
     const profile = parseLensComplexPupilProfile(input.profile);
     const propagationDistanceMm = input.propagationDistanceMm === undefined ? profile.context.focalLengthMm : input.propagationDistanceMm;
     requirePositiveFinite(propagationDistanceMm, "propagationDistanceMm");
-    const kernel = calculateDiscreteFraunhoferKernel(profile, propagationDistanceMm);
+    const plan = planPupilSampling(profile, propagationDistanceMm);
+    const kernel = calculateDiscreteFraunhoferKernel(profile, propagationDistanceMm, plan);
     return approximationResult({
         propagationDistanceMm,
         propagationDistanceSource: input.propagationDistanceMm === undefined ? "nominal-focal-plane-reference" : "explicit",
@@ -430,6 +520,7 @@ export function calculateLensComplexPupilPsf(input) {
             }
         },
         kernel,
+        sampling: plan.sampling,
         relativePupilThroughputFactor: profile
             .relativePupilThroughputFactor,
         throughputOwnership: "separate-relative-pupil-throughput-factor",
@@ -451,8 +542,9 @@ export function calculateLensComplexPupilPsf(input) {
             pupilProfile: profile.evidence
         },
         uncertainty: profile.uncertainty
-    }, "lens-complex-pupil-fraunhofer-psf", "1.1.0", [
+    }, "lens-complex-pupil-fraunhofer-psf", "2.0.0", [
         "The reference evaluator performs deterministic scalar Fraunhofer propagation of one explicit complex pupil amplitude plus optical-path-difference grid.",
+        "The pupil is zero padded so the DFT spans at least 2S-1 samples for a non-zero support of S samples per axis: the discrete intensity PSF then has no wrap-around aliasing and Goodman's sampling factor Q = lambda*N/pitch exceeds 2. Pupil support wider than the declared f/N aperture is rejected.",
         "Physical sample pitch uses the explicit propagation distance when supplied; omission preserves the nominal focal-plane reference. Finite focus metadata does not resolve the pupil-to-image geometry automatically.",
         "Diffraction and aberration phase are evaluated together in the same pupil propagation; callers must not stack an independent diffraction blur over this result.",
         "The output PSF is unit-energy normalized. The explicit relativePupilThroughputFactor is not multiplied into the kernel and remains owned by the separate optical-throughput path.",

package/dist/optics/complex-pupil-psf.js.map

@@ -1 +1 @@
-{"version":3,"file":"complex-pupil-psf.js","sourceRoot":"","sources":["../../src/optics/complex-pupil-psf.ts"],"names":[],"mappings":"AAAA,sCAAsC;AAEtC;;;;;;;;;;;;;;GAcG;AAEH,OAAO,EACL,mBAAmB,EAEpB,MAAM,+BAA+B,CAAC;AACvC,OAAO,EAAE,yBAAyB,EAAE,MAAM,gCAAgC,CAAC;AAC3E,OAAO,EACL,iBAAiB,EAElB,MAAM,gCAAgC,CAAC;AACxC,OAAO,EAAE,2BAA2B,EAAE,MAAM,uBAAuB,CAAC;AACpE,OAAO,EACL,eAAe,EAEhB,MAAM,kBAAkB,CAAC;AAS1B,MAAM,CAAC,MAAM,yCAAyC,GACpD,OAAgB,CAAC;AAgInB,SAAS,aAAa,CACpB,KAAc,EACd,IAAY;IAEZ,IACE,OAAO,KAAK,KAAK,QAAQ;QACzB,KAAK,KAAK,IAAI;QACd,KAAK,CAAC,OAAO,CAAC,KAAK,CAAC,EACpB,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,IAAI,GAAG,qBAAqB,CAC7B,CAAC;IACJ,CAAC;IACD,OAAO,KAAsB,CAAC;AAChC,CAAC;AAED,SAAS,qBAAqB,CAC5B,KAAc,EACd,IAAY;IAEZ,IACE,OAAO,KAAK,KAAK,QAAQ;QACzB,KAAK,CAAC,IAAI,EAAE,CAAC,MAAM,KAAK,CAAC,EACzB,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,IAAI;YACF,8BAA8B,CACjC,CAAC;IACJ,CAAC;IACD,OAAO,KAAK,CAAC,IAAI,EAAE,CAAC;AACtB,CAAC;AAED,SAAS,aAAa,CACpB,KAAc,EACd,IAAY;IAEZ,IACE,OAAO,KAAK,KAAK,QAAQ;QACzB,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,EACvB,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,IAAI,GAAG,kBAAkB,CAC1B,CAAC;IACJ,CAAC;IACD,OAAO,KAAK,CAAC;AACf,CAAC;AAED,SAAS,qBAAqB,CAC5B,KAAc,EACd,IAAY;IAEZ,MAAM,MAAM,GACV,aAAa,CAAC,KAAK,EAAE,IAAI,CAAC,CAAC;IAC7B,IAAI,MAAM,IAAI,CAAC,EAAE,CAAC;QAChB,MAAM,IAAI,yBAAyB,CACjC,IAAI;YACF,6BAA6B,CAChC,CAAC;IACJ,CAAC;IACD,OAAO,MAAM,CAAC;AAChB,CAAC;AAED,SAAS,mBAAmB,CAC1B,KAAc,EACd,IAAY;IAEZ,MAAM,MAAM,GACV,aAAa,CAAC,KAAK,EAAE,IAAI,CAAC,CAAC;IAC7B,IACE,MAAM,GAAG,CAAC;QACV,MAAM,GAAG,CAAC,EACV,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,IAAI;YACF,kCAAkC,CACrC,CAAC;IACJ,CAAC;IACD,OAAO,MAAM,CAAC;AAChB,CAAC;AAED,SAAS,0BAA0B,CACjC,KAAc,EACd,IAAY;IAEZ,IACE,OAAO,KAAK,KAAK,QAAQ;QACzB,CAAC,MAAM,CAAC,aAAa,CAAC,KAAK,CAAC;QAC5B,KAAK,IAAI,CAAC,EACV,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,IAAI;YACF,mCAAmC,CACtC,CAAC;IACJ,CAAC;IACD,OAAO,KAAK,CAAC;AACf,CAAC;AAED,SAAS,6BAA6B,CACpC,KAAc,EACd,IAAY;IAEZ,IACE,OAAO,KAAK,KAAK,QAAQ;QACzB,CAAC,MAAM,CAAC,aAAa,CAAC,KAAK,CAAC;QAC5B,KAAK,GAAG,CAAC,EACT,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,IAAI;YACF,uCAAuC,CAC1C,CAAC;IACJ,CAAC;IACD,OAAO,KAAK,CAAC;AACf,CAAC;AAED,SAAS,WAAW,CAClB,KAAc,EACd,IAAY;IAEZ,IACE,KAAK,KAAK,YAAY;QACtB,KAAK,KAAK,eAAe,EACzB,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,IAAI,GAAG,cAAc,CACtB,CAAC;IACJ,CAAC;IACD,OAAO,KAAK,CAAC;AACf,CAAC;AAED,SAAS,gBAAgB,CACvB,KAAc,EACd,IAAY;IAEZ,MAAM,MAAM,GACV,aAAa,CAAC,KAAK,EAAE,IAAI,CAAC,CAAC;IAC7B,IAAI,MAAM,CAAC,IAAI,KAAK,UAAU,EAAE,CAAC;QAC/B,MAAM,QAAQ,GACZ,aAAa,CACX,MAAM,CAAC,QAAQ,EACf,IAAI,GAAG,WAAW,CACnB,CAAC;QACJ,IAAI,QAAQ,GAAG,CAAC,EAAE,CAAC;YACjB,MAAM,IAAI,yBAAyB,CACjC,IAAI;gBACF,gCAAgC,CACnC,CAAC;QACJ,CAAC;QACD,OAAO;YACL,IAAI,EAAE,UAAU;YAChB,QAAQ;YACR,KAAK,EACH,qBAAqB,CACnB,MAAM,CAAC,KAAK,EACZ,IAAI,GAAG,QAAQ,CAChB;SACJ,CAAC;IACJ,CAAC;IACD,IACE,MAAM,CAAC,IAAI;QACX,gBAAgB,EAChB,CAAC;QACD,OAAO;YACL,IAAI,EAAE,gBAAgB;YACtB,UAAU,EACR,qBAAqB,CACnB,MAAM,CAAC,UAAU,EACjB,IAAI;gBACF,aAAa,CAChB;SACJ,CAAC;IACJ,CAAC;IACD,MAAM,IAAI,yBAAyB,CACjC,IAAI,GAAG,mBAAmB,CAC3B,CAAC;AACJ,CAAC;AAED,SAAS,gBAAgB,CACvB,KAAc,EACd,IAAY;IAEZ,IAAI,CAAC,KAAK,CAAC,OAAO,CAAC,KAAK,CAAC,EAAE,CAAC;QAC1B,MAAM,IAAI,yBAAyB,CACjC,IAAI,GAAG,oBAAoB,CAC5B,CAAC;IACJ,CAAC;IACD,MAAM,MAAM,GACV,KAAK,CAAC,GAAG,CACP,CAAC,KAAK,EAAE,KAAK,EAAE,EAAE,CACf,qBAAqB,CACnB,KAAK,EACL,IAAI;QACF,GAAG;QACH,KAAK;QACL,GAAG,CACN,CACJ,CAAC;IACJ,IACE,IAAI,GAAG,CAAC,MAAM,CAAC,CAAC,IAAI;QACpB,MAAM,CAAC,MAAM,EACb,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,IAAI;YACF,+BAA+B,CAClC,CAAC;IACJ,CAAC;IACD,OAAO,MAAM,CAAC;AAChB,CAAC;AAED,SAAS,SAAS,CAChB,KAAc;IAEd,MAAM,MAAM,GACV,aAAa,CACX,KAAK,EACL,uBAAuB,CACxB,CAAC;IACJ,MAAM,YAAY,GAChB,0BAA0B,CACxB,MAAM,CAAC,YAAY,EACnB,oCAAoC,CACrC,CAAC;IACJ,MAAM,aAAa,GACjB,0BAA0B,CACxB,MAAM,CAAC,aAAa,EACpB,qCAAqC,CACtC,CAAC;IACJ,MAAM,KAAK,GACT,YAAY;QACZ,aAAa,CAAC;IAChB,IACE,CAAC,MAAM,CAAC,aAAa,CAAC,KAAK,CAAC;QAC5B,KAAK,GAAG,IAAI,EACZ,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,0DAA0D,CAC3D,CAAC;IACJ,CAAC;IACD,MAAM,aAAa,GACjB,6BAA6B,CAC3B,MAAM,CAAC,aAAa,EACpB,qCAAqC,CACtC,CAAC;IACJ,MAAM,aAAa,GACjB,6BAA6B,CAC3B,MAAM,CAAC,aAAa,EACpB,qCAAqC,CACtC,CAAC;IACJ,IACE,aAAa;QACX,YAAY;QACd,aAAa;YACX,aAAa,EACf,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,+DAA+D,CAChE,CAAC;IACJ,CAAC;IACD,IACE,CAAC,KAAK,CAAC,OAAO,CACZ,MAAM,CAAC,iBAAiB,CACzB;QACD,MAAM,CAAC,iBAAiB;aACrB,MAAM,KAAK,KAAK,EACnB,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,yFAAyF,CAC1F,CAAC;IACJ,CAAC;IACD,IACE,CAAC,KAAK,CAAC,OAAO,CACZ,MAAM;SACH,gCAAgC,CACpC;QACD,MAAM;aACH,gCAAgC;aAChC,MAAM,KAAK,KAAK,EACnB,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,wGAAwG,CACzG,CAAC;IACJ,CAAC;IAED,MAAM,iBAAiB,GACrB,MAAM,CAAC,iBAAiB,CAAC,GAAG,CAC1B,CAAC,KAAK,EAAE,KAAK,EAAE,EAAE,CACf,mBAAmB,CACjB,KAAK,EACL,0CAA0C;QACxC,KAAK;QACL,GAAG,CACN,CACJ,CAAC;IACJ,IACE,CAAC,iBAAiB,CAAC,IAAI,CACrB,CAAC,KAAK,EAAE,EAAE,CACR,KAAK,GAAG,CAAC,CACZ,EACD,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,6EAA6E,CAC9E,CAAC;IACJ,CAAC;IACD,MAAM,GAAG,GACP,MAAM;SACH,gCAAgC;SAChC,GAAG,CACF,CAAC,KAAK,EAAE,KAAK,EAAE,EAAE,CACf,aAAa,CACX,KAAK,EACL,yDAAyD;QACvD,KAAK;QACL,GAAG,CACN,CACJ,CAAC;IAEN,OAAO;QACL,YAAY;QACZ,aAAa;QACb,mBAAmB,EACjB,qBAAqB,CACnB,MAAM;aACH,mBAAmB,EACtB,2CAA2C,CAC5C;QACH,mBAAmB,EACjB,qBAAqB,CACnB,MAAM;aACH,mBAAmB,EACtB,2CAA2C,CAC5C;QACH,aAAa;QACb,aAAa;QACb,iBAAiB;QACjB,gCAAgC,EAC9B,GAAG;KACN,CAAC;AACJ,CAAC;AAED;;;;;;;;;;;;GAYG;AACH,MAAM,UAAU,4BAA4B,CAC1C,KAAc;IAEd,MAAM,MAAM,GACV,aAAa,CACX,KAAK,EACL,kBAAkB,CACnB,CAAC;IACJ,IACE,MAAM,CAAC,aAAa;QACpB,yCAAyC,EACzC,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,0CAA0C;YACxC,yCAAyC;YACzC,IAAI,CACP,CAAC;IACJ,CAAC;IACD,IACE,MAAM,CAAC,aAAa;QAClB,gCAAgC;QAClC,MAAM,CAAC,cAAc;YACnB,kCAAkC;QACpC,MAAM,CAAC,qBAAqB;YAC1B,2CAA2C;QAC7C,MAAM,CAAC,cAAc;YACnB,iBAAiB;QACnB,MAAM,CAAC,gBAAgB;YACrB,wCAAwC;QAC1C,MAAM,CAAC,gBAAgB;YACrB,qCAAqC;QACvC,MAAM,CAAC,mBAAmB;YACxB,2CAA2C;QAC7C,MAAM,CAAC,yBAAyB;YAC9B,mBAAmB;QACrB,MAAM,CAAC,gBAAgB;YACrB,sCAAsC,EACxC,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,0EAA0E,CAC3E,CAAC;IACJ,CAAC;IACD,IACE,MAAM,CAAC,0BAA0B;QAC/B,KAAK;QACP,MAAM,CAAC,sBAAsB;YAC3B,KAAK;QACP,MAAM,CAAC,sBAAsB;YAC3B,KAAK;QACP,MAAM,CAAC,kBAAkB;YACvB,KAAK,EACP,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,mFAAmF,CACpF,CAAC;IACJ,CAAC;IACD,MAAM,MAAM,GACV,WAAW,CACT,MAAM,CAAC,gBAAgB,EACvB,mCAAmC,CACpC,CAAC;IACJ,MAAM,WAAW,GACf,gBAAgB,CACd,MAAM,CAAC,WAAW,EAClB,8BAA8B,CAC/B,CAAC;IACJ,IACE,MAAM,KAAK,YAAY;QACvB,WAAW,CAAC,IAAI;YACd,gBAAgB,EAClB,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,iFAAiF,CAClF,CAAC;IACJ,CAAC;IACD,MAAM,OAAO,GACX,aAAa,CACX,MAAM,CAAC,OAAO,EACd,0BAA0B,CAC3B,CAAC;IACJ,MAAM,KAAK,GACT,aAAa,CACX,OAAO,CAAC,YAAY,EACpB,uCAAuC,CACxC,CAAC;IACJ,MAAM,QAAQ,GACZ,aAAa,CACX,MAAM,CAAC,gBAAgB,EACvB,mCAAmC,CACpC,CAAC;IACJ,IACE,QAAQ,CAAC,WAAW;QACpB,IAAI;QACJ,OAAO,QAAQ,CAAC,UAAU;YACxB,SAAS;QACX,OAAO,QAAQ,CAAC,OAAO;YACrB,SAAS;QACX,OAAO,QAAQ,CAAC,kBAAkB;YAChC,SAAS,EACX,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,4GAA4G,CAC7G,CAAC;IACJ,CAAC;IAED,OAAO;QACL,aAAa,EACX,yCAAyC;QAC3C,SAAS,EACP,qBAAqB,CACnB,MAAM,CAAC,SAAS,EAChB,4BAA4B,CAC7B;QACH,cAAc,EACZ,qBAAqB,CACnB,MAAM,CAAC,cAAc,EACrB,iCAAiC,CAClC;QACH,gBAAgB,EACd,MAAM;QACR,aAAa,EACX,gCAAgC;QAClC,cAAc,EACZ,kCAAkC;QACpC,qBAAqB,EACnB,2CAA2C;QAC7C,cAAc,EACZ,iBAAiB;QACnB,gBAAgB,EACd,wCAAwC;QAC1C,gBAAgB,EACd,qCAAqC;QACvC,mBAAmB,EACjB,2CAA2C;QAC7C,yBAAyB,EACvB,mBAAmB;QACrB,gBAAgB,EACd,sCAAsC;QACxC,OAAO,EAAE;YACP,aAAa,EACX,qBAAqB,CACnB,OAAO,CAAC,aAAa,EACrB,wCAAwC,CACzC;YACH,KAAK,EACH,eAAe,CACb,OAAO,CAAC,KAAK,CACd;YACH,eAAe,EACb,qBAAqB,CACnB,OAAO,CAAC,eAAe,EACvB,0CAA0C,CAC3C;YACH,YAAY,EAAE;gBACZ,CAAC,EACC,aAAa,CACX,KAAK,CAAC,CAAC,EACP,yCAAyC,CAC1C;gBACH,CAAC,EACC,aAAa,CACX,KAAK,CAAC,CAAC,EACP,yCAAyC,CAC1C;aACJ;YACD,YAAY,EACV,qBAAqB,CACnB,OAAO,CAAC,YAAY,EACpB,uCAAuC,CACxC;YACH,kCAAkC,EAChC,aAAa,CACX,OAAO;iBACJ,kCAAkC,EACrC,6DAA6D,CAC9D;SACJ;QACD,gBAAgB,EAAE;YAChB,WAAW,EAAE,IAAI;YACjB,UAAU,EACR,QAAQ,CAAC,UAAqB;YAChC,OAAO,EACL,QAAQ,CAAC,OAAkB;YAC7B,kBAAkB,EAChB,QAAQ;iBACL,kBAA6B;SACnC;QACD,6BAA6B,EAC3B,mBAAmB,CACjB,MAAM;aACH,6BAA6B,EAChC,gDAAgD,CACjD;QACH,IAAI,EACF,SAAS,CAAC,MAAM,CAAC,IAAI,CAAC;QACxB,0BAA0B,EACxB,KAAK;QACP,sBAAsB,EACpB,KAAK;QACP,sBAAsB,EACpB,KAAK;QACP,kBAAkB,EAChB,KAAK;QACP,QAAQ,EACN,iBAAiB,CACf,MAAM,CAAC,QAAQ,EACf,2BAA2B,CAC5B;QACH,WAAW;QACX,WAAW,EACT,gBAAgB,CACd,MAAM,CAAC,WAAW,EAClB,8BAA8B,CAC/B;KACJ,CAAC;AACJ,CAAC;AAOD,SAAS,iBAAiB,CACxB,SAAiB,EACjB,cAAsB,EACtB,qBACQ;IAER,MAAM,KAAK,GACT,CAAC;QACD,IAAI,CAAC,EAAE;QACP,cAAc;QACd,qBAAqB,CAAC;IACxB,OAAO;QACL,EAAE,EACA,SAAS;YACT,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC;QACjB,EAAE,EACA,SAAS;YACT,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC;KAClB,CAAC;AACJ,CAAC;AAED,SAAS,iCAAiC,CACxC,OACyB,EACzB,qBAA6B;IAE7B,MAAM,IAAI,GACR,OAAO,CAAC,IAAI,CAAC;IACf,MAAM,KAAK,GACT,IAAI,CAAC,YAAY,CAAC;IACpB,MAAM,MAAM,GACV,IAAI,CAAC,aAAa,CAAC;IACrB,MAAM,qBAAqB,GACzB,OAAO,CAAC,OAAO;SACZ,YAAY;QACf,IAAI,CAAC;IACP,MAAM,YAAY,GAChB,OAAO,CAAC,OAAO;SACZ,YAAY;QACf,SAAS,CAAC;IACZ,MAAM,KAAK,GACT,IAAI,CAAC,iBAAiB,CAAC,GAAG,CACxB,CAAC,SAAS,EAAE,KAAK,EAAE,EAAE,CACnB,iBAAiB,CACf,SAAS,EACT,IAAI;SACD,gCAAgC,CAC/B,KAAK,CACL,EACJ,qBAAqB,CACtB,CACJ,CAAC;IAEJ,MAAM,SAAS,GACb,IAAI,KAAK,CACP,KAAK,GAAG,MAAM,CACf,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAEZ,KACE,IAAI,IAAI,GAAG,CAAC,EACZ,IAAI,GAAG,MAAM,EACb,IAAI,IAAI,CAAC,EACT,CAAC;QACD,MAAM,EAAE,GACN,IAAI;YACJ,IAAI,CAAC,aAAa,CAAC;QACrB,KACE,IAAI,IAAI,GAAG,CAAC,EACZ,IAAI,GAAG,KAAK,EACZ,IAAI,IAAI,CAAC,EACT,CAAC;YACD,MAAM,EAAE,GACN,IAAI;gBACJ,IAAI,CAAC,aAAa,CAAC;YACrB,IAAI,EAAE,GAAG,CAAC,CAAC;YACX,IAAI,EAAE,GAAG,CAAC,CAAC;YAEX,KACE,IAAI,CAAC,GAAG,CAAC,EACT,CAAC,GAAG,MAAM,EACV,CAAC,IAAI,CAAC,EACN,CAAC;gBACD,MAAM,EAAE,GACN,CAAC;oBACD,IAAI,CAAC,aAAa,CAAC;gBACrB,KACE,IAAI,CAAC,GAAG,CAAC,EACT,CAAC,GAAG,KAAK,EACT,CAAC,IAAI,CAAC,EACN,CAAC;oBACD,MAAM,EAAE,GACN,CAAC;wBACD,IAAI,CAAC,aAAa,CAAC;oBACrB,MAAM,MAAM,GACV,KAAK,CACH,CAAC,GAAG,KAAK;wBACT,CAAC,CACD,CAAC;oBACL,MAAM,KAAK,GACT,CAAC,CAAC;wBACF,IAAI,CAAC,EAAE;wBACP,CACE,EAAE,GAAG,EAAE;4BACL,KAAK;4BACP,EAAE,GAAG,EAAE;gCACL,MAAM,CACT,CAAC;oBACJ,MAAM,GAAG,GACP,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;oBAClB,MAAM,GAAG,GACP,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;oBAClB,EAAE;wBACA,MAAM,CAAC,EAAE;4BACP,GAAG;4BACL,MAAM,CAAC,EAAE;gCACP,GAAG,CAAC;oBACR,EAAE;wBACA,MAAM,CAAC,EAAE;4BACP,GAAG;4BACL,MAAM,CAAC,EAAE;gCACP,GAAG,CAAC;gBACV,CAAC;YACH,CAAC;YAED,SAAS,CACP,IAAI,GAAG,KAAK;gBACZ,IAAI,CACL;gBACC,EAAE,GAAG,EAAE;oBACP,EAAE,GAAG,EAAE,CAAC;QACZ,CAAC;IACH,CAAC;IAED,MAAM,GAAG,GACP,SAAS,CAAC,MAAM,CACd,CAAC,GAAG,EAAE,KAAK,EAAE,EAAE,CACb,GAAG,GAAG,KAAK,EACb,CAAC,CACF,CAAC;IACJ,IACE,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC;QACrB,GAAG,IAAI,CAAC,EACR,CAAC;QACD,MAAM,IAAI,2BAA2B,CACnC,oEAAoE,CACrE,CAAC;IACJ,CAAC;IACD,MAAM,UAAU,GACd,SAAS,CAAC,GAAG,CACX,CAAC,KAAK,EAAE,EAAE,CACR,KAAK,GAAG,GAAG,CACd,CAAC;IAEJ,MAAM,cAAc,GAClB,YAAY;QACZ,qBAAqB;QACrB,CACE,KAAK;YACL,IAAI;iBACD,mBAAmB,CACvB,CAAC;IACJ,MAAM,cAAc,GAClB,YAAY;QACZ,qBAAqB;QACrB,CACE,MAAM;YACN,IAAI;iBACD,mBAAmB,CACvB,CAAC;IAEJ,IACE,CAAC,MAAM,CAAC,QAAQ,CACd,cAAc,CACf;QACD,CAAC,MAAM,CAAC,QAAQ,CACd,cAAc,CACf;QACD,cAAc,IAAI,CAAC;QACnB,cAAc,IAAI,CAAC,EACnB,CAAC;QACD,MAAM,IAAI,2BAA2B,CACnC,iEAAiE,CAClE,CAAC;IACJ,CAAC;IAED,OAAO;QACL,YAAY,EACV,KAAK;QACP,aAAa,EACX,MAAM;QACR,uBAAuB,EACrB,cAAc;YACd,IAAI;QACN,uBAAuB,EACrB,cAAc;YACd,IAAI;QACN,aAAa,EACX,IAAI,CAAC,aAAa;QACpB,aAAa,EACX,IAAI,CAAC,aAAa;QACpB,mBAAmB,EACjB,UAAU;KACb,CAAC;AACJ,CAAC;AAED;;;;;;;;;;;;GAYG;AACH,MAAM,UAAU,4BAA4B,CAC1C,KACmC;IAEnC,MAAM,OAAO,GACX,4BAA4B,CAC1B,KAAK,CAAC,OAAO,CACd,CAAC;IACJ,MAAM,qBAAqB,GAAG,KAAK,CAAC,qBAAqB,KAAK,SAAS,CAAC,CAAC,CAAC,OAAO,CAAC,OAAO,CAAC,aAAa,CAAC,CAAC,CAAC,KAAK,CAAC,qBAAqB,CAAC;IACtI,qBAAqB,CAAC,qBAAqB,EAAE,uBAAuB,CAAC,CAAC;IACtE,MAAM,MAAM,GACV,iCAAiC,CAC/B,OAAO,EACP,qBAAqB,CACtB,CAAC;IAEJ,OAAO,mBAAmB,CACxB;QACE,qBAAqB;QACrB,yBAAyB,EAAE,KAAK,CAAC,qBAAqB,KAAK,SAAS,CAAC,CAAC,CAAC,+BAA+B,CAAC,CAAC,CAAC,UAAU;QACnH,SAAS,EACP,OAAO,CAAC,SAAS;QACnB,cAAc,EACZ,OAAO,CAAC,cAAc;QACxB,gBAAgB,EACd,OAAO,CAAC,gBAAgB;QAC1B,OAAO,EAAE;YACP,GAAG,OAAO,CAAC,OAAO;YAClB,KAAK,EAAE;gBACL,GAAG,OAAO,CAAC,OAAO,CAAC,KAAK;aACzB;YACD,YAAY,EAAE;gBACZ,GAAG,OAAO,CAAC,OAAO;qBACf,YAAY;aAChB;SACF;QACD,MAAM;QACN,6BAA6B,EAC3B,OAAO;aACJ,6BAA6B;QAClC,mBAAmB,EACjB,2CAA2C;QAC7C,wCAAwC,EACtC,IAAI;QACN,0BAA0B,EACxB,OAAO;aACJ,gBAAgB;aAChB,kBAAkB;QACvB,sCAAsC,EACpC,KAAK;QACP,yBAAyB,EACvB,IAAI;QACN,0BAA0B,EACxB,IAAI;QACN,2BAA2B,EACzB,KAAK;QACP,oCAAoC,EAClC,KAAK;QACP,0BAA0B,EACxB,KAAK;QACP,sBAAsB,EACpB,KAAK;QACP,sBAAsB,EACpB,KAAK;QACP,kBAAkB,EAChB,KAAK;QACP,4BAA4B,EAC1B,KAAK;QACP,iBAAiB,EAAE;YACjB,YAAY,EACV,OAAO,CAAC,QAAQ;SACnB;QACD,WAAW,EACT,OAAO,CAAC,WAAW;KACtB,EACD,mCAAmC,EACnC,OAAO,EACP;QACE,yJAAyJ;QACzJ,qNAAqN;QACrN,iKAAiK;QACjK,mLAAmL;QACnL,+HAA+H;QAC/H,qHAAqH;QACrH,wJAAwJ;QACxJ,uFAAuF;KACxF,CACF,CAAC;AACJ,CAAC"}
+{"version":3,"file":"complex-pupil-psf.js","sourceRoot":"","sources":["../../src/optics/complex-pupil-psf.ts"],"names":[],"mappings":"AAAA,sCAAsC;AAEtC;;;;;;;;;;;;GAYG;AAEH,OAAO,EACL,mBAAmB,EAEpB,MAAM,+BAA+B,CAAC;AACvC,OAAO,EAAE,yBAAyB,EAAE,MAAM,gCAAgC,CAAC;AAC3E,OAAO,EACL,iBAAiB,EAElB,MAAM,gCAAgC,CAAC;AACxC,OAAO,EAAE,2BAA2B,EAAE,MAAM,uBAAuB,CAAC;AACpE,OAAO,EACL,eAAe,EAEhB,MAAM,kBAAkB,CAAC;AAM1B,OAAO,EAAE,gBAAgB,EAAE,mBAAmB,EAAE,0BAA0B,EAAE,0BAA0B,EAAE,MAAM,8BAA8B,CAAC;AAI7I,MAAM,CAAC,MAAM,yCAAyC,GACpD,OAAgB,CAAC;AAkInB,SAAS,aAAa,CACpB,KAAc,EACd,IAAY,EACZ,WAAqC;IAErC,IACE,OAAO,KAAK,KAAK,QAAQ;QACzB,KAAK,KAAK,IAAI;QACd,KAAK,CAAC,OAAO,CAAC,KAAK,CAAC,EACpB,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,IAAI,GAAG,qBAAqB,CAC7B,CAAC;IACJ,CAAC;IACD,IAAI,WAAW,KAAK,IAAI;QAAE,mBAAmB,CAAC,KAAK,EAAE,WAAW,EAAE,IAAI,CAAC,CAAC;IACxE,OAAO,KAAsB,CAAC;AAChC,CAAC;AAED,SAAS,qBAAqB,CAC5B,KAAc,EACd,IAAY;IAEZ,IACE,OAAO,KAAK,KAAK,QAAQ;QACzB,KAAK,CAAC,IAAI,EAAE,CAAC,MAAM,KAAK,CAAC,EACzB,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,IAAI;YACF,8BAA8B,CACjC,CAAC;IACJ,CAAC;IACD,OAAO,KAAK,CAAC,IAAI,EAAE,CAAC;AACtB,CAAC;AAED,SAAS,aAAa,CACpB,KAAc,EACd,IAAY;IAEZ,IACE,OAAO,KAAK,KAAK,QAAQ;QACzB,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,EACvB,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,IAAI,GAAG,kBAAkB,CAC1B,CAAC;IACJ,CAAC;IACD,OAAO,KAAK,CAAC;AACf,CAAC;AAED,SAAS,qBAAqB,CAC5B,KAAc,EACd,IAAY;IAEZ,MAAM,MAAM,GACV,aAAa,CAAC,KAAK,EAAE,IAAI,CAAC,CAAC;IAC7B,IAAI,MAAM,IAAI,CAAC,EAAE,CAAC;QAChB,MAAM,IAAI,yBAAyB,CACjC,IAAI;YACF,6BAA6B,CAChC,CAAC;IACJ,CAAC;IACD,OAAO,MAAM,CAAC;AAChB,CAAC;AAED,SAAS,mBAAmB,CAC1B,KAAc,EACd,IAAY;IAEZ,MAAM,MAAM,GACV,aAAa,CAAC,KAAK,EAAE,IAAI,CAAC,CAAC;IAC7B,IACE,MAAM,GAAG,CAAC;QACV,MAAM,GAAG,CAAC,EACV,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,IAAI;YACF,kCAAkC,CACrC,CAAC;IACJ,CAAC;IACD,OAAO,MAAM,CAAC;AAChB,CAAC;AAED,SAAS,0BAA0B,CACjC,KAAc,EACd,IAAY;IAEZ,IACE,OAAO,KAAK,KAAK,QAAQ;QACzB,CAAC,MAAM,CAAC,aAAa,CAAC,KAAK,CAAC;QAC5B,KAAK,IAAI,CAAC,EACV,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,IAAI;YACF,mCAAmC,CACtC,CAAC;IACJ,CAAC;IACD,OAAO,KAAK,CAAC;AACf,CAAC;AAED,SAAS,6BAA6B,CACpC,KAAc,EACd,IAAY;IAEZ,IACE,OAAO,KAAK,KAAK,QAAQ;QACzB,CAAC,MAAM,CAAC,aAAa,CAAC,KAAK,CAAC;QAC5B,KAAK,GAAG,CAAC,EACT,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,IAAI;YACF,uCAAuC,CAC1C,CAAC;IACJ,CAAC;IACD,OAAO,KAAK,CAAC;AACf,CAAC;AAED,SAAS,WAAW,CAClB,KAAc,EACd,IAAY;IAEZ,IACE,KAAK,KAAK,YAAY;QACtB,KAAK,KAAK,eAAe,EACzB,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,IAAI,GAAG,cAAc,CACtB,CAAC;IACJ,CAAC;IACD,OAAO,KAAK,CAAC;AACf,CAAC;AAED,SAAS,gBAAgB,CACvB,KAAc,EACd,IAAY;IAEZ,MAAM,MAAM,GACV,aAAa,CAAC,KAAK,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;IACnC,0BAA0B,CAAC,MAAM,EAAE,MAAM,EAAE,0BAA0B,EAAE,IAAI,CAAC,CAAC;IAC7E,IAAI,MAAM,CAAC,IAAI,KAAK,UAAU,EAAE,CAAC;QAC/B,MAAM,QAAQ,GACZ,aAAa,CACX,MAAM,CAAC,QAAQ,EACf,IAAI,GAAG,WAAW,CACnB,CAAC;QACJ,IAAI,QAAQ,GAAG,CAAC,EAAE,CAAC;YACjB,MAAM,IAAI,yBAAyB,CACjC,IAAI;gBACF,gCAAgC,CACnC,CAAC;QACJ,CAAC;QACD,OAAO;YACL,IAAI,EAAE,UAAU;YAChB,QAAQ;YACR,KAAK,EACH,qBAAqB,CACnB,MAAM,CAAC,KAAK,EACZ,IAAI,GAAG,QAAQ,CAChB;SACJ,CAAC;IACJ,CAAC;IACD,IACE,MAAM,CAAC,IAAI;QACX,gBAAgB,EAChB,CAAC;QACD,OAAO;YACL,IAAI,EAAE,gBAAgB;YACtB,UAAU,EACR,qBAAqB,CACnB,MAAM,CAAC,UAAU,EACjB,IAAI;gBACF,aAAa,CAChB;SACJ,CAAC;IACJ,CAAC;IACD,MAAM,IAAI,yBAAyB,CACjC,IAAI,GAAG,mBAAmB,CAC3B,CAAC;AACJ,CAAC;AAED,SAAS,gBAAgB,CACvB,KAAc,EACd,IAAY;IAEZ,IAAI,CAAC,KAAK,CAAC,OAAO,CAAC,KAAK,CAAC,EAAE,CAAC;QAC1B,MAAM,IAAI,yBAAyB,CACjC,IAAI,GAAG,oBAAoB,CAC5B,CAAC;IACJ,CAAC;IACD,MAAM,MAAM,GACV,KAAK,CAAC,GAAG,CACP,CAAC,KAAK,EAAE,KAAK,EAAE,EAAE,CACf,qBAAqB,CACnB,KAAK,EACL,IAAI;QACF,GAAG;QACH,KAAK;QACL,GAAG,CACN,CACJ,CAAC;IACJ,IACE,IAAI,GAAG,CAAC,MAAM,CAAC,CAAC,IAAI;QACpB,MAAM,CAAC,MAAM,EACb,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,IAAI;YACF,+BAA+B,CAClC,CAAC;IACJ,CAAC;IACD,OAAO,MAAM,CAAC;AAChB,CAAC;AAED,SAAS,SAAS,CAChB,KAAc;IAEd,MAAM,MAAM,GACV,aAAa,CACX,KAAK,EACL,uBAAuB,EAAE,CAAC,cAAc,EAAE,eAAe,EAAE,eAAe,EAAE,eAAe,EAAE,mBAAmB,EAAE,kCAAkC,EAAE,qBAAqB,EAAE,qBAAqB,CAAC,CACpM,CAAC;IACJ,MAAM,YAAY,GAChB,0BAA0B,CACxB,MAAM,CAAC,YAAY,EACnB,oCAAoC,CACrC,CAAC;IACJ,MAAM,aAAa,GACjB,0BAA0B,CACxB,MAAM,CAAC,aAAa,EACpB,qCAAqC,CACtC,CAAC;IACJ,MAAM,KAAK,GACT,YAAY;QACZ,aAAa,CAAC;IAChB,IACE,CAAC,MAAM,CAAC,aAAa,CAAC,KAAK,CAAC;QAC5B,KAAK,GAAG,IAAI,EACZ,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,0DAA0D,CAC3D,CAAC;IACJ,CAAC;IACD,MAAM,aAAa,GACjB,6BAA6B,CAC3B,MAAM,CAAC,aAAa,EACpB,qCAAqC,CACtC,CAAC;IACJ,MAAM,aAAa,GACjB,6BAA6B,CAC3B,MAAM,CAAC,aAAa,EACpB,qCAAqC,CACtC,CAAC;IACJ,IACE,aAAa;QACX,YAAY;QACd,aAAa;YACX,aAAa,EACf,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,+DAA+D,CAChE,CAAC;IACJ,CAAC;IACD,IACE,CAAC,KAAK,CAAC,OAAO,CACZ,MAAM,CAAC,iBAAiB,CACzB;QACD,MAAM,CAAC,iBAAiB;aACrB,MAAM,KAAK,KAAK,EACnB,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,yFAAyF,CAC1F,CAAC;IACJ,CAAC;IACD,IACE,CAAC,KAAK,CAAC,OAAO,CACZ,MAAM;SACH,gCAAgC,CACpC;QACD,MAAM;aACH,gCAAgC;aAChC,MAAM,KAAK,KAAK,EACnB,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,wGAAwG,CACzG,CAAC;IACJ,CAAC;IAED,MAAM,iBAAiB,GACrB,MAAM,CAAC,iBAAiB,CAAC,GAAG,CAC1B,CAAC,KAAK,EAAE,KAAK,EAAE,EAAE,CACf,mBAAmB,CACjB,KAAK,EACL,0CAA0C;QACxC,KAAK;QACL,GAAG,CACN,CACJ,CAAC;IACJ,IACE,CAAC,iBAAiB,CAAC,IAAI,CACrB,CAAC,KAAK,EAAE,EAAE,CACR,KAAK,GAAG,CAAC,CACZ,EACD,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,6EAA6E,CAC9E,CAAC;IACJ,CAAC;IACD,MAAM,GAAG,GACP,MAAM;SACH,gCAAgC;SAChC,GAAG,CACF,CAAC,KAAK,EAAE,KAAK,EAAE,EAAE,CACf,aAAa,CACX,KAAK,EACL,yDAAyD;QACvD,KAAK;QACL,GAAG,CACN,CACJ,CAAC;IAEN,OAAO;QACL,YAAY;QACZ,aAAa;QACb,mBAAmB,EACjB,qBAAqB,CACnB,MAAM;aACH,mBAAmB,EACtB,2CAA2C,CAC5C;QACH,mBAAmB,EACjB,qBAAqB,CACnB,MAAM;aACH,mBAAmB,EACtB,2CAA2C,CAC5C;QACH,aAAa;QACb,aAAa;QACb,iBAAiB;QACjB,gCAAgC,EAC9B,GAAG;KACN,CAAC;AACJ,CAAC;AAED;;;;;;;;;;;;GAYG;AACH,MAAM,UAAU,4BAA4B,CAC1C,KAAc;IAEd,MAAM,MAAM,GACV,aAAa,CACX,KAAK,EACL,kBAAkB,EAAE,CAAC,eAAe,EAAE,eAAe,EAAE,gBAAgB,EAAE,uBAAuB,EAAE,gBAAgB,EAAE,kBAAkB,EAAE,kBAAkB,EAAE,qBAAqB,EAAE,2BAA2B,EAAE,kBAAkB,EAAE,4BAA4B,EAAE,wBAAwB,EAAE,wBAAwB,EAAE,oBAAoB,EAAE,kBAAkB,EAAE,aAAa,EAAE,SAAS,EAAE,kBAAkB,EAAE,WAAW,EAAE,gBAAgB,EAAE,+BAA+B,EAAE,MAAM,EAAE,UAAU,EAAE,aAAa,CAAC,CACjf,CAAC;IACJ,IACE,MAAM,CAAC,aAAa;QACpB,yCAAyC,EACzC,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,0CAA0C;YACxC,yCAAyC;YACzC,IAAI,CACP,CAAC;IACJ,CAAC;IACD,IACE,MAAM,CAAC,aAAa;QAClB,gCAAgC;QAClC,MAAM,CAAC,cAAc;YACnB,kCAAkC;QACpC,MAAM,CAAC,qBAAqB;YAC1B,2CAA2C;QAC7C,MAAM,CAAC,cAAc;YACnB,iBAAiB;QACnB,MAAM,CAAC,gBAAgB;YACrB,wCAAwC;QAC1C,MAAM,CAAC,gBAAgB;YACrB,qCAAqC;QACvC,MAAM,CAAC,mBAAmB;YACxB,2CAA2C;QAC7C,MAAM,CAAC,yBAAyB;YAC9B,mBAAmB;QACrB,MAAM,CAAC,gBAAgB;YACrB,sCAAsC,EACxC,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,0EAA0E,CAC3E,CAAC;IACJ,CAAC;IACD,IACE,MAAM,CAAC,0BAA0B;QAC/B,KAAK;QACP,MAAM,CAAC,sBAAsB;YAC3B,KAAK;QACP,MAAM,CAAC,sBAAsB;YAC3B,KAAK;QACP,MAAM,CAAC,kBAAkB;YACvB,KAAK,EACP,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,mFAAmF,CACpF,CAAC;IACJ,CAAC;IACD,MAAM,MAAM,GACV,WAAW,CACT,MAAM,CAAC,gBAAgB,EACvB,mCAAmC,CACpC,CAAC;IACJ,MAAM,WAAW,GACf,gBAAgB,CACd,MAAM,CAAC,WAAW,EAClB,8BAA8B,CAC/B,CAAC;IACJ,IACE,MAAM,KAAK,YAAY;QACvB,WAAW,CAAC,IAAI;YACd,gBAAgB,EAClB,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,iFAAiF,CAClF,CAAC;IACJ,CAAC;IACD,MAAM,OAAO,GACX,aAAa,CACX,MAAM,CAAC,OAAO,EACd,0BAA0B,EAAE,CAAC,cAAc,EAAE,eAAe,EAAE,OAAO,EAAE,iBAAiB,EAAE,cAAc,EAAE,oCAAoC,CAAC,CAChJ,CAAC;IACJ,MAAM,KAAK,GACT,aAAa,CACX,OAAO,CAAC,YAAY,EACpB,uCAAuC,EAAE,CAAC,GAAG,EAAE,GAAG,CAAC,CACpD,CAAC;IACJ,MAAM,QAAQ,GACZ,aAAa,CACX,MAAM,CAAC,gBAAgB,EACvB,mCAAmC,EAAE,CAAC,aAAa,EAAE,YAAY,EAAE,SAAS,EAAE,oBAAoB,CAAC,CACpG,CAAC;IACJ,IACE,QAAQ,CAAC,WAAW;QACpB,IAAI;QACJ,OAAO,QAAQ,CAAC,UAAU;YACxB,SAAS;QACX,OAAO,QAAQ,CAAC,OAAO;YACrB,SAAS;QACX,OAAO,QAAQ,CAAC,kBAAkB;YAChC,SAAS,EACX,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,4GAA4G,CAC7G,CAAC;IACJ,CAAC;IAED,OAAO;QACL,aAAa,EACX,yCAAyC;QAC3C,SAAS,EACP,qBAAqB,CACnB,MAAM,CAAC,SAAS,EAChB,4BAA4B,CAC7B;QACH,cAAc,EACZ,qBAAqB,CACnB,MAAM,CAAC,cAAc,EACrB,iCAAiC,CAClC;QACH,gBAAgB,EACd,MAAM;QACR,aAAa,EACX,gCAAgC;QAClC,cAAc,EACZ,kCAAkC;QACpC,qBAAqB,EACnB,2CAA2C;QAC7C,cAAc,EACZ,iBAAiB;QACnB,gBAAgB,EACd,wCAAwC;QAC1C,gBAAgB,EACd,qCAAqC;QACvC,mBAAmB,EACjB,2CAA2C;QAC7C,yBAAyB,EACvB,mBAAmB;QACrB,gBAAgB,EACd,sCAAsC;QACxC,OAAO,EAAE;YACP,aAAa,EACX,qBAAqB,CACnB,OAAO,CAAC,aAAa,EACrB,wCAAwC,CACzC;YACH,KAAK,EACH,eAAe,CACb,OAAO,CAAC,KAAK,CACd;YACH,eAAe,EACb,qBAAqB,CACnB,OAAO,CAAC,eAAe,EACvB,0CAA0C,CAC3C;YACH,YAAY,EAAE;gBACZ,CAAC,EACC,aAAa,CACX,KAAK,CAAC,CAAC,EACP,yCAAyC,CAC1C;gBACH,CAAC,EACC,aAAa,CACX,KAAK,CAAC,CAAC,EACP,yCAAyC,CAC1C;aACJ;YACD,YAAY,EACV,qBAAqB,CACnB,OAAO,CAAC,YAAY,EACpB,uCAAuC,CACxC;YACH,kCAAkC,EAChC,aAAa,CACX,OAAO;iBACJ,kCAAkC,EACrC,6DAA6D,CAC9D;SACJ;QACD,gBAAgB,EAAE;YAChB,WAAW,EAAE,IAAI;YACjB,UAAU,EACR,QAAQ,CAAC,UAAqB;YAChC,OAAO,EACL,QAAQ,CAAC,OAAkB;YAC7B,kBAAkB,EAChB,QAAQ;iBACL,kBAA6B;SACnC;QACD,6BAA6B,EAC3B,mBAAmB,CACjB,MAAM;aACH,6BAA6B,EAChC,gDAAgD,CACjD;QACH,IAAI,EACF,SAAS,CAAC,MAAM,CAAC,IAAI,CAAC;QACxB,0BAA0B,EACxB,KAAK;QACP,sBAAsB,EACpB,KAAK;QACP,sBAAsB,EACpB,KAAK;QACP,kBAAkB,EAChB,KAAK;QACP,QAAQ,EACN,iBAAiB,CACf,MAAM,CAAC,QAAQ,EACf,2BAA2B,CAC5B;QACH,WAAW;QACX,WAAW,EACT,gBAAgB,CACd,MAAM,CAAC,WAAW,EAClB,8BAA8B,CAC/B;KACJ,CAAC;AACJ,CAAC;AAOD,SAAS,iBAAiB,CACxB,SAAiB,EACjB,cAAsB,EACtB,qBACQ;IAER,MAAM,KAAK,GACT,CAAC;QACD,IAAI,CAAC,EAAE;QACP,cAAc;QACd,qBAAqB,CAAC;IACxB,OAAO;QACL,EAAE,EACA,SAAS;YACT,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC;QACjB,EAAE,EACA,SAAS;YACT,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC;KAClB,CAAC;AACJ,CAAC;AAyCD,SAAS,WAAW,CAClB,IAA0B,EAC1B,IAAe;IAEf,IAAI,OAAO,GAAG,MAAM,CAAC,iBAAiB,CAAC;IACvC,IAAI,OAAO,GAAG,MAAM,CAAC,iBAAiB,CAAC;IACvC,IAAI,CAAC,iBAAiB,CAAC,OAAO,CAAC,CAAC,SAAS,EAAE,KAAK,EAAE,EAAE;QAClD,IAAI,SAAS,GAAG,CAAC,EAAE,CAAC;YAClB,MAAM,UAAU,GACd,IAAI,KAAK,GAAG;gBACV,CAAC,CAAC,KAAK,GAAG,IAAI,CAAC,YAAY;gBAC3B,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,KAAK,GAAG,IAAI,CAAC,YAAY,CAAC,CAAC;YAC5C,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,EAAE,UAAU,CAAC,CAAC;YACxC,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,EAAE,UAAU,CAAC,CAAC;QAC1C,CAAC;IACH,CAAC,CAAC,CAAC;IACH,OAAO,OAAO,GAAG,OAAO,GAAG,CAAC,CAAC;AAC/B,CAAC;AAED;;;;;;;;;;;;;;;;;;;;;;;GAuBG;AACH,SAAS,iBAAiB,CACxB,OAAgC,EAChC,qBAA6B;IAE7B,MAAM,IAAI,GAAG,OAAO,CAAC,IAAI,CAAC;IAC1B,MAAM,QAAQ,GAAG,WAAW,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC;IACxC,MAAM,QAAQ,GAAG,WAAW,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC;IACxC,MAAM,SAAS,GAAG,CAAC,QAAQ,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,mBAAmB,CAAC;IAC5D,MAAM,SAAS,GAAG,CAAC,QAAQ,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,mBAAmB,CAAC;IAC5D,MAAM,gBAAgB,GACpB,qBAAqB,GAAG,OAAO,CAAC,OAAO,CAAC,eAAe,CAAC;IAC1D,uDAAuD;IACvD,MAAM,KAAK,GAAG,gBAAgB,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC;IAC5C,IAAI,SAAS,GAAG,KAAK,IAAI,SAAS,GAAG,KAAK,EAAE,CAAC;QAC3C,MAAM,IAAI,2BAA2B,CACnC,wCAAwC;YACtC,SAAS,CAAC,WAAW,CAAC,CAAC,CAAC;YACxB,KAAK;YACL,SAAS,CAAC,WAAW,CAAC,CAAC,CAAC;YACxB,oFAAoF;YACpF,gBAAgB,CAAC,WAAW,CAAC,CAAC,CAAC;YAC/B,gDAAgD,CACnD,CAAC;IACJ,CAAC;IACD,MAAM,cAAc,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,YAAY,EAAE,CAAC,GAAG,QAAQ,GAAG,CAAC,CAAC,CAAC;IACrE,MAAM,eAAe,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,aAAa,EAAE,CAAC,GAAG,QAAQ,GAAG,CAAC,CAAC,CAAC;IACvE,OAAO;QACL,QAAQ,EAAE;YACR,SAAS,EAAE,iDAAiD;YAC5D,oBAAoB,EAAE,QAAQ;YAC9B,oBAAoB,EAAE,QAAQ;YAC9B,uBAAuB,EAAE,SAAS;YAClC,uBAAuB,EAAE,SAAS;YAClC,uBAAuB,EAAE,gBAAgB;YACzC,qBAAqB,EAAE,IAAI,CAAC,YAAY;YACxC,sBAAsB,EAAE,IAAI,CAAC,aAAa;YAC1C,qBAAqB,EAAE,cAAc;YACrC,sBAAsB,EAAE,eAAe;YACvC,kBAAkB,EAChB,cAAc,KAAK,IAAI,CAAC,YAAY;gBACpC,eAAe,KAAK,IAAI,CAAC,aAAa;YACxC,gBAAgB,EACd,QAAQ,GAAG,CAAC,CAAC,CAAC,CAAC,cAAc,GAAG,CAAC,QAAQ,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI;YACvD,gBAAgB,EACd,QAAQ,GAAG,CAAC,CAAC,CAAC,CAAC,eAAe,GAAG,CAAC,QAAQ,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI;SACzD;QACD,yEAAyE;QACzE,0BAA0B;QAC1B,aAAa,EACX,IAAI,CAAC,aAAa;YAClB,IAAI,CAAC,KAAK,CAAC,CAAC,cAAc,GAAG,IAAI,CAAC,YAAY,CAAC,GAAG,CAAC,CAAC;QACtD,aAAa,EACX,IAAI,CAAC,aAAa;YAClB,IAAI,CAAC,KAAK,CAAC,CAAC,eAAe,GAAG,IAAI,CAAC,aAAa,CAAC,GAAG,CAAC,CAAC;KACzD,CAAC;AACJ,CAAC;AAED,SAAS,YAAY,CAAC,IAAY;IAChC,MAAM,GAAG,GAAa,EAAE,CAAC;IACzB,MAAM,GAAG,GAAa,EAAE,CAAC;IACzB,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,IAAI,EAAE,KAAK,IAAI,CAAC,EAAE,CAAC;QAC7C,MAAM,KAAK,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,KAAK,CAAC,GAAG,IAAI,CAAC;QAC5C,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC;QAC1B,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC;IAC5B,CAAC;IACD,OAAO,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;AACtB,CAAC;AAED,SAAS,cAAc,CAAC,KAAa,EAAE,IAAY;IACjD,MAAM,SAAS,GAAG,KAAK,GAAG,IAAI,CAAC;IAC/B,OAAO,SAAS,GAAG,CAAC,CAAC,CAAC,CAAC,SAAS,GAAG,IAAI,CAAC,CAAC,CAAC,SAAS,CAAC;AACtD,CAAC;AAED,SAAS,iCAAiC,CACxC,OACyB,EACzB,qBAA6B,EAC7B,IAAuB;IAEvB,MAAM,IAAI,GAAG,OAAO,CAAC,IAAI,CAAC;IAC1B,MAAM,KAAK,GAAG,IAAI,CAAC,YAAY,CAAC;IAChC,MAAM,MAAM,GAAG,IAAI,CAAC,aAAa,CAAC;IAClC,MAAM,QAAQ,GAAG,IAAI,CAAC,QAAQ,CAAC,qBAAqB,CAAC;IACrD,MAAM,SAAS,GAAG,IAAI,CAAC,QAAQ,CAAC,sBAAsB,CAAC;IACvD,MAAM,qBAAqB,GAAG,OAAO,CAAC,OAAO,CAAC,YAAY,GAAG,IAAI,CAAC;IAClE,MAAM,YAAY,GAAG,OAAO,CAAC,OAAO,CAAC,YAAY,GAAG,SAAS,CAAC;IAC9D,MAAM,KAAK,GACT,IAAI,CAAC,iBAAiB,CAAC,GAAG,CACxB,CAAC,SAAS,EAAE,KAAK,EAAE,EAAE,CACnB,iBAAiB,CACf,SAAS,EACT,IAAI,CAAC,gCAAgC,CAAC,KAAK,CAAE,EAC7C,qBAAqB,CACtB,CACJ,CAAC;IAEJ,yEAAyE;IACzE,uEAAuE;IACvE,MAAM,QAAQ,GAAG,YAAY,CAAC,QAAQ,CAAC,CAAC;IACxC,MAAM,QAAQ,GAAG,YAAY,CAAC,SAAS,CAAC,CAAC;IACzC,MAAM,KAAK,GAAG,IAAI,YAAY,CAAC,MAAM,GAAG,QAAQ,CAAC,CAAC;IAClD,MAAM,KAAK,GAAG,IAAI,YAAY,CAAC,MAAM,GAAG,QAAQ,CAAC,CAAC;IAClD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC;QACnC,KAAK,IAAI,IAAI,GAAG,CAAC,EAAE,IAAI,GAAG,QAAQ,EAAE,IAAI,IAAI,CAAC,EAAE,CAAC;YAC9C,MAAM,EAAE,GAAG,IAAI,GAAG,IAAI,CAAC,aAAa,CAAC;YACrC,IAAI,EAAE,GAAG,CAAC,CAAC;YACX,IAAI,EAAE,GAAG,CAAC,CAAC;YACX,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC;gBAClC,MAAM,MAAM,GAAG,KAAK,CAAC,CAAC,GAAG,KAAK,GAAG,CAAC,CAAE,CAAC;gBACrC,MAAM,CAAC,GAAG,cAAc,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,aAAa,CAAC,EAAE,QAAQ,CAAC,CAAC;gBAClE,MAAM,GAAG,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAE,CAAC;gBAC7B,MAAM,GAAG,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAE,CAAC;gBAC7B,EAAE,IAAI,MAAM,CAAC,EAAE,GAAG,GAAG,GAAG,MAAM,CAAC,EAAE,GAAG,GAAG,CAAC;gBACxC,EAAE,IAAI,MAAM,CAAC,EAAE,GAAG,GAAG,GAAG,MAAM,CAAC,EAAE,GAAG,GAAG,CAAC;YAC1C,CAAC;YACD,KAAK,CAAC,CAAC,GAAG,QAAQ,GAAG,IAAI,CAAC,GAAG,EAAE,CAAC;YAChC,KAAK,CAAC,CAAC,GAAG,QAAQ,GAAG,IAAI,CAAC,GAAG,EAAE,CAAC;QAClC,CAAC;IACH,CAAC;IAED,MAAM,SAAS,GAAG,IAAI,KAAK,CAAS,QAAQ,GAAG,SAAS,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAClE,KAAK,IAAI,IAAI,GAAG,CAAC,EAAE,IAAI,GAAG,SAAS,EAAE,IAAI,IAAI,CAAC,EAAE,CAAC;QAC/C,MAAM,EAAE,GAAG,IAAI,GAAG,IAAI,CAAC,aAAa,CAAC;QACrC,KAAK,IAAI,IAAI,GAAG,CAAC,EAAE,IAAI,GAAG,QAAQ,EAAE,IAAI,IAAI,CAAC,EAAE,CAAC;YAC9C,IAAI,EAAE,GAAG,CAAC,CAAC;YACX,IAAI,EAAE,GAAG,CAAC,CAAC;YACX,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC;gBACnC,MAAM,QAAQ,GAAG,KAAK,CAAC,CAAC,GAAG,QAAQ,GAAG,IAAI,CAAE,CAAC;gBAC7C,MAAM,QAAQ,GAAG,KAAK,CAAC,CAAC,GAAG,QAAQ,GAAG,IAAI,CAAE,CAAC;gBAC7C,MAAM,CAAC,GAAG,cAAc,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,aAAa,CAAC,EAAE,SAAS,CAAC,CAAC;gBACnE,MAAM,GAAG,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAE,CAAC;gBAC7B,MAAM,GAAG,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAE,CAAC;gBAC7B,EAAE,IAAI,QAAQ,GAAG,GAAG,GAAG,QAAQ,GAAG,GAAG,CAAC;gBACtC,EAAE,IAAI,QAAQ,GAAG,GAAG,GAAG,QAAQ,GAAG,GAAG,CAAC;YACxC,CAAC;YACD,SAAS,CAAC,IAAI,GAAG,QAAQ,GAAG,IAAI,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QACxD,CAAC;IACH,CAAC;IAED,MAAM,GAAG,GACP,SAAS,CAAC,MAAM,CACd,CAAC,GAAG,EAAE,KAAK,EAAE,EAAE,CACb,GAAG,GAAG,KAAK,EACb,CAAC,CACF,CAAC;IACJ,IACE,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC;QACrB,GAAG,IAAI,CAAC,EACR,CAAC;QACD,MAAM,IAAI,2BAA2B,CACnC,oEAAoE,CACrE,CAAC;IACJ,CAAC;IACD,MAAM,UAAU,GACd,SAAS,CAAC,GAAG,CACX,CAAC,KAAK,EAAE,EAAE,CACR,KAAK,GAAG,GAAG,CACd,CAAC;IAEJ,MAAM,cAAc,GAClB,YAAY,GAAG,qBAAqB;QACpC,CAAC,QAAQ,GAAG,IAAI,CAAC,mBAAmB,CAAC,CAAC;IACxC,MAAM,cAAc,GAClB,YAAY,GAAG,qBAAqB;QACpC,CAAC,SAAS,GAAG,IAAI,CAAC,mBAAmB,CAAC,CAAC;IAEzC,IACE,CAAC,MAAM,CAAC,QAAQ,CACd,cAAc,CACf;QACD,CAAC,MAAM,CAAC,QAAQ,CACd,cAAc,CACf;QACD,cAAc,IAAI,CAAC;QACnB,cAAc,IAAI,CAAC,EACnB,CAAC;QACD,MAAM,IAAI,2BAA2B,CACnC,iEAAiE,CAClE,CAAC;IACJ,CAAC;IAED,OAAO;QACL,YAAY,EAAE,QAAQ;QACtB,aAAa,EAAE,SAAS;QACxB,uBAAuB,EACrB,cAAc;YACd,IAAI;QACN,uBAAuB,EACrB,cAAc;YACd,IAAI;QACN,aAAa,EAAE,IAAI,CAAC,aAAa;QACjC,aAAa,EAAE,IAAI,CAAC,aAAa;QACjC,mBAAmB,EACjB,UAAU;KACb,CAAC;AACJ,CAAC;AAED;;;;;;;;;;;;GAYG;AACH,MAAM,UAAU,4BAA4B,CAC1C,KACmC;IAEnC,gBAAgB,CAAC,KAAK,EAAE,8BAA8B,EAAE,CAAC,SAAS,EAAE,uBAAuB,CAAC,CAAC,CAAC;IAC9F,MAAM,OAAO,GACX,4BAA4B,CAC1B,KAAK,CAAC,OAAO,CACd,CAAC;IACJ,MAAM,qBAAqB,GAAG,KAAK,CAAC,qBAAqB,KAAK,SAAS,CAAC,CAAC,CAAC,OAAO,CAAC,OAAO,CAAC,aAAa,CAAC,CAAC,CAAC,KAAK,CAAC,qBAAqB,CAAC;IACtI,qBAAqB,CAAC,qBAAqB,EAAE,uBAAuB,CAAC,CAAC;IACtE,MAAM,IAAI,GAAG,iBAAiB,CAAC,OAAO,EAAE,qBAAqB,CAAC,CAAC;IAC/D,MAAM,MAAM,GACV,iCAAiC,CAC/B,OAAO,EACP,qBAAqB,EACrB,IAAI,CACL,CAAC;IAEJ,OAAO,mBAAmB,CACxB;QACE,qBAAqB;QACrB,yBAAyB,EAAE,KAAK,CAAC,qBAAqB,KAAK,SAAS,CAAC,CAAC,CAAC,+BAA+B,CAAC,CAAC,CAAC,UAAU;QACnH,SAAS,EACP,OAAO,CAAC,SAAS;QACnB,cAAc,EACZ,OAAO,CAAC,cAAc;QACxB,gBAAgB,EACd,OAAO,CAAC,gBAAgB;QAC1B,OAAO,EAAE;YACP,GAAG,OAAO,CAAC,OAAO;YAClB,KAAK,EAAE;gBACL,GAAG,OAAO,CAAC,OAAO,CAAC,KAAK;aACzB;YACD,YAAY,EAAE;gBACZ,GAAG,OAAO,CAAC,OAAO;qBACf,YAAY;aAChB;SACF;QACD,MAAM;QACN,QAAQ,EAAE,IAAI,CAAC,QAAQ;QACvB,6BAA6B,EAC3B,OAAO;aACJ,6BAA6B;QAClC,mBAAmB,EACjB,2CAA2C;QAC7C,wCAAwC,EACtC,IAAI;QACN,0BAA0B,EACxB,OAAO;aACJ,gBAAgB;aAChB,kBAAkB;QACvB,sCAAsC,EACpC,KAAK;QACP,yBAAyB,EACvB,IAAI;QACN,0BAA0B,EACxB,IAAI;QACN,2BAA2B,EACzB,KAAK;QACP,oCAAoC,EAClC,KAAK;QACP,0BAA0B,EACxB,KAAK;QACP,sBAAsB,EACpB,KAAK;QACP,sBAAsB,EACpB,KAAK;QACP,kBAAkB,EAChB,KAAK;QACP,4BAA4B,EAC1B,KAAK;QACP,iBAAiB,EAAE;YACjB,YAAY,EACV,OAAO,CAAC,QAAQ;SACnB;QACD,WAAW,EACT,OAAO,CAAC,WAAW;KACtB,EACD,mCAAmC,EACnC,OAAO,EACP;QACE,yJAAyJ;QACzJ,uSAAuS;QACvS,qNAAqN;QACrN,iKAAiK;QACjK,mLAAmL;QACnL,+HAA+H;QAC/H,qHAAqH;QACrH,wJAAwJ;QACxJ,uFAAuF;KACxF,CACF,CAAC;AACJ,CAAC","sourcesContent":["// SPDX-License-Identifier: Apache-2.0\n\n/**\n * PSF (point-spread function) from a declared complex pupil: amplitude and phase sampled on a grid.\n *\n * The profile parser validates untrusted input and returns the normalized typed contract. Unknown\n * enum values, missing required fields and incompatible scientific data are rejected.\n *\n * The complex pupil is evaluated by summing the declared amplitude and phase over a bounded, regular\n * pupil grid. The PSF shape is normalized separately from the pupil throughput. A single blur radius\n * cannot stand in for coherent phase, pupil energy, or the wavelength and field range where the\n * profile applies. Numeric profiles need provenance showing the data is owned or reusable.\n * @see docs/PSF_FOUNDATION.md for the equations, units, coordinate conventions, known blockers and\n * support limits.\n */\n\nimport {\n  approximationResult,\n  type CalculationResult\n} from \"../core/calculation-result.js\";\nimport { InvalidConfigurationError } from \"../core/configuration-error.js\";\nimport {\n  parseEvidenceList,\n  type EvidenceProvenance\n} from \"../core/evidence-provenance.js\";\nimport { InvalidScientificInputError } from \"../core/validation.js\";\nimport {\n  parseFocusPlane,\n  type FocusPlane\n} from \"./focus-state.js\";\nimport type {\n  LensPsfKernel,\n  LensPsfScientificStatus,\n  LensPsfUncertainty\n} from \"./lens-psf-profile.js\";\nimport { requireInputKeys, rejectUnknownFields, rejectUnknownVariantFields, UNCERTAINTY_VARIANT_FIELDS } from \"../core/record-validation.js\";\n\ntype UnknownRecord = Record<string, unknown>;\n\nexport const LENS_COMPLEX_PUPIL_PROFILE_SCHEMA_VERSION =\n  \"0.1.0\" as const;\n\nexport interface LensComplexPupilGrid {\n  widthSamples: number;\n  heightSamples: number;\n  pupilSamplePitchMmX: number;\n  pupilSamplePitchMmY: number;\n  centerSampleX: number;\n  centerSampleY: number;\n  relativeAmplitude:\n    readonly number[];\n  opticalPathDifferenceMicrometers:\n    readonly number[];\n}\n\nexport interface LensComplexPupilProfile {\n  schemaVersion:\n    typeof LENS_COMPLEX_PUPIL_PROFILE_SCHEMA_VERSION;\n  profileId: string;\n  profileVersion: string;\n  scientificStatus:\n    LensPsfScientificStatus;\n  opticalDomain:\n    \"lens-primary-optical-path-only\";\n  representation:\n    \"complex-pupil-amplitude-plus-opd\";\n  pupilCoordinateSystem:\n    \"pupil-plane-metric-aligned-to-image-plane\";\n  imageFieldAxes:\n    \"+X right, +Y up\";\n  amplitudeMeaning:\n    \"relative-complex-pupil-amplitude-shape\";\n  wavefrontMeaning:\n    \"optical-path-difference-micrometers\";\n  throughputOwnership:\n    \"separate-relative-pupil-throughput-factor\";\n  kernelEnergyNormalization:\n    \"unit-energy-shape\";\n  propagationModel:\n    \"scalar-fraunhofer-discrete-reference\";\n  context: {\n    focalLengthMm: number;\n    focus: FocusPlane;\n    apertureFNumber: number;\n    fieldPointMm: {\n      x: number;\n      y: number;\n    };\n    wavelengthNm: number;\n    signedDefocusImagePlaneMicrometers:\n      number;\n  };\n  responseIncludes: {\n    diffraction: true;\n    aberration: boolean;\n    defocus: boolean;\n    pupilClippingShape:\n      boolean;\n  };\n  relativePupilThroughputFactor:\n    number;\n  grid: LensComplexPupilGrid;\n  sensorOpticalStackIncluded:\n    false;\n  sensorSamplingIncluded:\n    false;\n  reconstructionIncluded:\n    false;\n  strayLightIncluded: false;\n  evidence:\n    readonly EvidenceProvenance[];\n  uncertainty:\n    LensPsfUncertainty;\n  limitations: readonly string[];\n}\n\nexport interface CalculateLensComplexPupilPsfInput {\n  profile:\n    LensComplexPupilProfile;\n  /** Physical propagation distance from pupil to image, in mm. If omitted, the\n   * nominal focal-plane reference is used. At finite focus the caller must\n   * resolve the supported principal-plane geometry; profile metadata does not. */\n  propagationDistanceMm?: number;\n}\n\nexport interface LensComplexPupilPsf {\n  propagationDistanceMm: number;\n  propagationDistanceSource: \"explicit\" | \"nominal-focal-plane-reference\";\n  profileId: string;\n  profileVersion: string;\n  scientificStatus:\n    LensPsfScientificStatus;\n  context:\n    LensComplexPupilProfile[\"context\"];\n  kernel: LensPsfKernel;\n  /** Nyquist sampling diagnostics and the applied zero padding. */\n  sampling: LensComplexPupilPsfSampling;\n  relativePupilThroughputFactor:\n    number;\n  throughputOwnership:\n    \"separate-relative-pupil-throughput-factor\";\n  diffractionAndAberrationJointlyEvaluated:\n    true;\n  pupilClippingShapeIncluded:\n    boolean;\n  pupilClippingThroughputAppliedToKernel:\n    false;\n  scalarFraunhoferReference:\n    true;\n  full2dOrientationPreserved:\n    true;\n  geometricDistortionModified:\n    false;\n  lateralChromaticPositionShiftApplied:\n    false;\n  sensorOpticalStackIncluded:\n    false;\n  sensorSamplingIncluded: false;\n  reconstructionIncluded: false;\n  strayLightIncluded: false;\n  scalarSharpnessScoreProduced:\n    false;\n  componentEvidence: {\n    pupilProfile:\n      readonly EvidenceProvenance[];\n  };\n  uncertainty:\n    LensPsfUncertainty;\n}\n\nfunction requireRecord(\n  value: unknown,\n  path: string,\n  allowedKeys: readonly string[] | null\n): UnknownRecord {\n  if (\n    typeof value !== \"object\" ||\n    value === null ||\n    Array.isArray(value)\n  ) {\n    throw new InvalidConfigurationError(\n      path + \" must be an object.\"\n    );\n  }\n  if (allowedKeys !== null) rejectUnknownFields(value, allowedKeys, path);\n  return value as UnknownRecord;\n}\n\nfunction requireNonEmptyString(\n  value: unknown,\n  path: string\n): string {\n  if (\n    typeof value !== \"string\" ||\n    value.trim().length === 0\n  ) {\n    throw new InvalidConfigurationError(\n      path +\n        \" must be a non-empty string.\"\n    );\n  }\n  return value.trim();\n}\n\nfunction requireFinite(\n  value: unknown,\n  path: string\n): number {\n  if (\n    typeof value !== \"number\" ||\n    !Number.isFinite(value)\n  ) {\n    throw new InvalidConfigurationError(\n      path + \" must be finite.\"\n    );\n  }\n  return value;\n}\n\nfunction requirePositiveFinite(\n  value: unknown,\n  path: string\n): number {\n  const parsed =\n    requireFinite(value, path);\n  if (parsed <= 0) {\n    throw new InvalidConfigurationError(\n      path +\n        \" must be greater than zero.\"\n    );\n  }\n  return parsed;\n}\n\nfunction requireUnitInterval(\n  value: unknown,\n  path: string\n): number {\n  const parsed =\n    requireFinite(value, path);\n  if (\n    parsed < 0 ||\n    parsed > 1\n  ) {\n    throw new InvalidConfigurationError(\n      path +\n        \" must lie from zero through one.\"\n    );\n  }\n  return parsed;\n}\n\nfunction requirePositiveSafeInteger(\n  value: unknown,\n  path: string\n): number {\n  if (\n    typeof value !== \"number\" ||\n    !Number.isSafeInteger(value) ||\n    value <= 0\n  ) {\n    throw new InvalidConfigurationError(\n      path +\n        \" must be a positive safe integer.\"\n    );\n  }\n  return value;\n}\n\nfunction requireNonNegativeSafeInteger(\n  value: unknown,\n  path: string\n): number {\n  if (\n    typeof value !== \"number\" ||\n    !Number.isSafeInteger(value) ||\n    value < 0\n  ) {\n    throw new InvalidConfigurationError(\n      path +\n        \" must be a non-negative safe integer.\"\n    );\n  }\n  return value;\n}\n\nfunction parseStatus(\n  value: unknown,\n  path: string\n): LensPsfScientificStatus {\n  if (\n    value !== \"calibrated\" &&\n    value !== \"approximation\"\n  ) {\n    throw new InvalidConfigurationError(\n      path + \" is invalid.\"\n    );\n  }\n  return value;\n}\n\nfunction parseUncertainty(\n  value: unknown,\n  path: string\n): LensPsfUncertainty {\n  const record =\n    requireRecord(value, path, null);\n  rejectUnknownVariantFields(record, \"kind\", UNCERTAINTY_VARIANT_FIELDS, path);\n  if (record.kind === \"relative\") {\n    const fraction =\n      requireFinite(\n        record.fraction,\n        path + \".fraction\"\n      );\n    if (fraction < 0) {\n      throw new InvalidConfigurationError(\n        path +\n          \".fraction must be nonnegative.\"\n      );\n    }\n    return {\n      kind: \"relative\",\n      fraction,\n      basis:\n        requireNonEmptyString(\n          record.basis,\n          path + \".basis\"\n        )\n    };\n  }\n  if (\n    record.kind ===\n    \"not-quantified\"\n  ) {\n    return {\n      kind: \"not-quantified\",\n      limitation:\n        requireNonEmptyString(\n          record.limitation,\n          path +\n            \".limitation\"\n        )\n    };\n  }\n  throw new InvalidConfigurationError(\n    path + \".kind is invalid.\"\n  );\n}\n\nfunction parseLimitations(\n  value: unknown,\n  path: string\n): readonly string[] {\n  if (!Array.isArray(value)) {\n    throw new InvalidConfigurationError(\n      path + \" must be an array.\"\n    );\n  }\n  const parsed =\n    value.map(\n      (entry, index) =>\n        requireNonEmptyString(\n          entry,\n          path +\n            \"[\" +\n            index +\n            \"]\"\n        )\n    );\n  if (\n    new Set(parsed).size !==\n    parsed.length\n  ) {\n    throw new InvalidConfigurationError(\n      path +\n        \" must not contain duplicates.\"\n    );\n  }\n  return parsed;\n}\n\nfunction parseGrid(\n  value: unknown\n): LensComplexPupilGrid {\n  const record =\n    requireRecord(\n      value,\n      \"lensComplexPupil.grid\", [\"widthSamples\", \"heightSamples\", \"centerSampleX\", \"centerSampleY\", \"relativeAmplitude\", \"opticalPathDifferenceMicrometers\", \"pupilSamplePitchMmX\", \"pupilSamplePitchMmY\"]\n    );\n  const widthSamples =\n    requirePositiveSafeInteger(\n      record.widthSamples,\n      \"lensComplexPupil.grid.widthSamples\"\n    );\n  const heightSamples =\n    requirePositiveSafeInteger(\n      record.heightSamples,\n      \"lensComplexPupil.grid.heightSamples\"\n    );\n  const total =\n    widthSamples *\n    heightSamples;\n  if (\n    !Number.isSafeInteger(total) ||\n    total > 4096\n  ) {\n    throw new InvalidConfigurationError(\n      \"lensComplexPupil grid must contain at most 4096 samples.\"\n    );\n  }\n  const centerSampleX =\n    requireNonNegativeSafeInteger(\n      record.centerSampleX,\n      \"lensComplexPupil.grid.centerSampleX\"\n    );\n  const centerSampleY =\n    requireNonNegativeSafeInteger(\n      record.centerSampleY,\n      \"lensComplexPupil.grid.centerSampleY\"\n    );\n  if (\n    centerSampleX >=\n      widthSamples ||\n    centerSampleY >=\n      heightSamples\n  ) {\n    throw new InvalidConfigurationError(\n      \"lensComplexPupil grid center must lie within the pupil array.\"\n    );\n  }\n  if (\n    !Array.isArray(\n      record.relativeAmplitude\n    ) ||\n    record.relativeAmplitude\n      .length !== total\n  ) {\n    throw new InvalidConfigurationError(\n      \"lensComplexPupil.grid.relativeAmplitude length must equal widthSamples × heightSamples.\"\n    );\n  }\n  if (\n    !Array.isArray(\n      record\n        .opticalPathDifferenceMicrometers\n    ) ||\n    record\n      .opticalPathDifferenceMicrometers\n      .length !== total\n  ) {\n    throw new InvalidConfigurationError(\n      \"lensComplexPupil.grid.opticalPathDifferenceMicrometers length must equal widthSamples × heightSamples.\"\n    );\n  }\n\n  const relativeAmplitude =\n    record.relativeAmplitude.map(\n      (entry, index) =>\n        requireUnitInterval(\n          entry,\n          \"lensComplexPupil.grid.relativeAmplitude[\" +\n            index +\n            \"]\"\n        )\n    );\n  if (\n    !relativeAmplitude.some(\n      (value) =>\n        value > 0\n    )\n  ) {\n    throw new InvalidConfigurationError(\n      \"lensComplexPupil must contain at least one non-zero pupil-amplitude sample.\"\n    );\n  }\n  const opd =\n    record\n      .opticalPathDifferenceMicrometers\n      .map(\n        (entry, index) =>\n          requireFinite(\n            entry,\n            \"lensComplexPupil.grid.opticalPathDifferenceMicrometers[\" +\n              index +\n              \"]\"\n          )\n      );\n\n  return {\n    widthSamples,\n    heightSamples,\n    pupilSamplePitchMmX:\n      requirePositiveFinite(\n        record\n          .pupilSamplePitchMmX,\n        \"lensComplexPupil.grid.pupilSamplePitchMmX\"\n      ),\n    pupilSamplePitchMmY:\n      requirePositiveFinite(\n        record\n          .pupilSamplePitchMmY,\n        \"lensComplexPupil.grid.pupilSamplePitchMmY\"\n      ),\n    centerSampleX,\n    centerSampleY,\n    relativeAmplitude,\n    opticalPathDifferenceMicrometers:\n      opd\n  };\n}\n\n/**\n * Validates an untrusted declaration and returns the normalized typed contract. Unknown enum values,\n * missing required fields and incompatible scientific data are rejected here.\n *\n * The complex pupil is evaluated by summing the declared amplitude and phase over a bounded, regular\n * pupil grid. The PSF shape is normalized separately from the pupil throughput. A single blur radius\n * cannot stand in for coherent phase, pupil energy, or the wavelength and field range where the\n * profile applies. Numeric profiles need provenance showing the data is owned or reusable.\n * @param value - unknown. Untrusted input; TypeScript types alone are not validation.\n * @returns LensComplexPupilProfile. The result states its shape and scientific status; a successful call does not imply calibration.\n *\n * @see docs/API_REFERENCE.md for the root export and exact type graph.\n */\nexport function parseLensComplexPupilProfile(\n  value: unknown\n): LensComplexPupilProfile {\n  const record =\n    requireRecord(\n      value,\n      \"lensComplexPupil\", [\"schemaVersion\", \"opticalDomain\", \"representation\", \"pupilCoordinateSystem\", \"imageFieldAxes\", \"amplitudeMeaning\", \"wavefrontMeaning\", \"throughputOwnership\", \"kernelEnergyNormalization\", \"propagationModel\", \"sensorOpticalStackIncluded\", \"sensorSamplingIncluded\", \"reconstructionIncluded\", \"strayLightIncluded\", \"scientificStatus\", \"uncertainty\", \"context\", \"responseIncludes\", \"profileId\", \"profileVersion\", \"relativePupilThroughputFactor\", \"grid\", \"evidence\", \"limitations\"]\n    );\n  if (\n    record.schemaVersion !==\n    LENS_COMPLEX_PUPIL_PROFILE_SCHEMA_VERSION\n  ) {\n    throw new InvalidConfigurationError(\n      'lensComplexPupil.schemaVersion must be \"' +\n        LENS_COMPLEX_PUPIL_PROFILE_SCHEMA_VERSION +\n        '\".'\n    );\n  }\n  if (\n    record.opticalDomain !==\n      \"lens-primary-optical-path-only\" ||\n    record.representation !==\n      \"complex-pupil-amplitude-plus-opd\" ||\n    record.pupilCoordinateSystem !==\n      \"pupil-plane-metric-aligned-to-image-plane\" ||\n    record.imageFieldAxes !==\n      \"+X right, +Y up\" ||\n    record.amplitudeMeaning !==\n      \"relative-complex-pupil-amplitude-shape\" ||\n    record.wavefrontMeaning !==\n      \"optical-path-difference-micrometers\" ||\n    record.throughputOwnership !==\n      \"separate-relative-pupil-throughput-factor\" ||\n    record.kernelEnergyNormalization !==\n      \"unit-energy-shape\" ||\n    record.propagationModel !==\n      \"scalar-fraunhofer-discrete-reference\"\n  ) {\n    throw new InvalidConfigurationError(\n      \"lensComplexPupil representation/coordinate/energy semantics are invalid.\"\n    );\n  }\n  if (\n    record.sensorOpticalStackIncluded !==\n      false ||\n    record.sensorSamplingIncluded !==\n      false ||\n    record.reconstructionIncluded !==\n      false ||\n    record.strayLightIncluded !==\n      false\n  ) {\n    throw new InvalidConfigurationError(\n      \"V1 complex pupil profiles must exclude sensor/reconstruction/stray-light effects.\"\n    );\n  }\n  const status =\n    parseStatus(\n      record.scientificStatus,\n      \"lensComplexPupil.scientificStatus\"\n    );\n  const uncertainty =\n    parseUncertainty(\n      record.uncertainty,\n      \"lensComplexPupil.uncertainty\"\n    );\n  if (\n    status === \"calibrated\" &&\n    uncertainty.kind ===\n      \"not-quantified\"\n  ) {\n    throw new InvalidConfigurationError(\n      \"Calibrated complex-pupil profiles must declare quantified relative uncertainty.\"\n    );\n  }\n  const context =\n    requireRecord(\n      record.context,\n      \"lensComplexPupil.context\", [\"fieldPointMm\", \"focalLengthMm\", \"focus\", \"apertureFNumber\", \"wavelengthNm\", \"signedDefocusImagePlaneMicrometers\"]\n    );\n  const field =\n    requireRecord(\n      context.fieldPointMm,\n      \"lensComplexPupil.context.fieldPointMm\", [\"x\", \"y\"]\n    );\n  const response =\n    requireRecord(\n      record.responseIncludes,\n      \"lensComplexPupil.responseIncludes\", [\"diffraction\", \"aberration\", \"defocus\", \"pupilClippingShape\"]\n    );\n  if (\n    response.diffraction !==\n    true ||\n    typeof response.aberration !==\n      \"boolean\" ||\n    typeof response.defocus !==\n      \"boolean\" ||\n    typeof response.pupilClippingShape !==\n      \"boolean\"\n  ) {\n    throw new InvalidConfigurationError(\n      \"lensComplexPupil.responseIncludes is invalid; diffraction must be included in this pupil propagation path.\"\n    );\n  }\n\n  return {\n    schemaVersion:\n      LENS_COMPLEX_PUPIL_PROFILE_SCHEMA_VERSION,\n    profileId:\n      requireNonEmptyString(\n        record.profileId,\n        \"lensComplexPupil.profileId\"\n      ),\n    profileVersion:\n      requireNonEmptyString(\n        record.profileVersion,\n        \"lensComplexPupil.profileVersion\"\n      ),\n    scientificStatus:\n      status,\n    opticalDomain:\n      \"lens-primary-optical-path-only\",\n    representation:\n      \"complex-pupil-amplitude-plus-opd\",\n    pupilCoordinateSystem:\n      \"pupil-plane-metric-aligned-to-image-plane\",\n    imageFieldAxes:\n      \"+X right, +Y up\",\n    amplitudeMeaning:\n      \"relative-complex-pupil-amplitude-shape\",\n    wavefrontMeaning:\n      \"optical-path-difference-micrometers\",\n    throughputOwnership:\n      \"separate-relative-pupil-throughput-factor\",\n    kernelEnergyNormalization:\n      \"unit-energy-shape\",\n    propagationModel:\n      \"scalar-fraunhofer-discrete-reference\",\n    context: {\n      focalLengthMm:\n        requirePositiveFinite(\n          context.focalLengthMm,\n          \"lensComplexPupil.context.focalLengthMm\"\n        ),\n      focus:\n        parseFocusPlane(\n          context.focus\n        ),\n      apertureFNumber:\n        requirePositiveFinite(\n          context.apertureFNumber,\n          \"lensComplexPupil.context.apertureFNumber\"\n        ),\n      fieldPointMm: {\n        x:\n          requireFinite(\n            field.x,\n            \"lensComplexPupil.context.fieldPointMm.x\"\n          ),\n        y:\n          requireFinite(\n            field.y,\n            \"lensComplexPupil.context.fieldPointMm.y\"\n          )\n      },\n      wavelengthNm:\n        requirePositiveFinite(\n          context.wavelengthNm,\n          \"lensComplexPupil.context.wavelengthNm\"\n        ),\n      signedDefocusImagePlaneMicrometers:\n        requireFinite(\n          context\n            .signedDefocusImagePlaneMicrometers,\n          \"lensComplexPupil.context.signedDefocusImagePlaneMicrometers\"\n        )\n    },\n    responseIncludes: {\n      diffraction: true,\n      aberration:\n        response.aberration as boolean,\n      defocus:\n        response.defocus as boolean,\n      pupilClippingShape:\n        response\n          .pupilClippingShape as boolean\n    },\n    relativePupilThroughputFactor:\n      requireUnitInterval(\n        record\n          .relativePupilThroughputFactor,\n        \"lensComplexPupil.relativePupilThroughputFactor\"\n      ),\n    grid:\n      parseGrid(record.grid),\n    sensorOpticalStackIncluded:\n      false,\n    sensorSamplingIncluded:\n      false,\n    reconstructionIncluded:\n      false,\n    strayLightIncluded:\n      false,\n    evidence:\n      parseEvidenceList(\n        record.evidence,\n        \"lensComplexPupil.evidence\"\n      ),\n    uncertainty,\n    limitations:\n      parseLimitations(\n        record.limitations,\n        \"lensComplexPupil.limitations\"\n      )\n  };\n}\n\ninterface ComplexValue {\n  re: number;\n  im: number;\n}\n\nfunction pupilComplexValue(\n  amplitude: number,\n  opdMicrometers: number,\n  wavelengthMicrometers:\n    number\n): ComplexValue {\n  const phase =\n    2 *\n    Math.PI *\n    opdMicrometers /\n    wavelengthMicrometers;\n  return {\n    re:\n      amplitude *\n      Math.cos(phase),\n    im:\n      amplitude *\n      Math.sin(phase)\n  };\n}\n\nexport interface LensComplexPupilPsfSampling {\n  criterion: \"discrete-intensity-alias-free-goodman-q-above-2\";\n  /** Number of pupil samples spanned by non-zero amplitude, per axis. */\n  pupilSupportSamplesX: number;\n  pupilSupportSamplesY: number;\n  /**\n   * Sampled pupil diameter D = (support samples - 1) x pupil pitch: the\n   * centre-to-centre span of the outermost non-zero samples, in millimetres.\n   */\n  pupilSupportDiameterMmX: number;\n  pupilSupportDiameterMmY: number;\n  /**\n   * Widest support allowed by the declared f-number at the resolved propagation distance:\n   * propagationDistanceMm / apertureFNumber, in millimetres. For the nominal focal-plane reference\n   * this is focalLengthMm / apertureFNumber.\n   */\n  declaredPupilDiameterMm: number;\n  /** Declared pupil grid size before zero padding. */\n  pupilGridWidthSamples: number;\n  pupilGridHeightSamples: number;\n  /** DFT size after zero padding; equals the output kernel size. */\n  transformWidthSamples: number;\n  transformHeightSamples: number;\n  zeroPaddingApplied: boolean;\n  /**\n   * Achieved Goodman sampling factor Q = lambda*(z/D)/imagePitch =\n   * transform samples / (support samples - 1). Always greater than 2. Null\n   * when the support is a single sample on that axis (no finite diameter).\n   */\n  samplingFactorQX: number | null;\n  samplingFactorQY: number | null;\n}\n\ninterface PupilSamplingPlan {\n  sampling: LensComplexPupilPsfSampling;\n  outputCenterX: number;\n  outputCenterY: number;\n}\n\nfunction supportSpan(\n  grid: LensComplexPupilGrid,\n  axis: \"x\" | \"y\"\n): number {\n  let minimum = Number.POSITIVE_INFINITY;\n  let maximum = Number.NEGATIVE_INFINITY;\n  grid.relativeAmplitude.forEach((amplitude, index) => {\n    if (amplitude > 0) {\n      const coordinate =\n        axis === \"x\"\n          ? index % grid.widthSamples\n          : Math.floor(index / grid.widthSamples);\n      minimum = Math.min(minimum, coordinate);\n      maximum = Math.max(maximum, coordinate);\n    }\n  });\n  return maximum - minimum + 1;\n}\n\n/**\n * Sampling guard for the discrete Fraunhofer reference (Ban et al.,\n * arXiv:2511.06936, eqs. 6-9, free: N*dp = lambda*f/dc, Q = lambda*F/dc and\n * Q >= 2. Also see: Goodman, Introduction to Fourier Optics, sampling of\n * DFT/FFT Fraunhofer propagation).\n *\n * A DFT (discrete Fourier transform) of size M over pupil pitch dp gives an\n * image-plane pitch of lambda*z/(M*dp).\n *\n * The intensity PSF is the transform of the pupil autocorrelation. That\n * autocorrelation is twice as wide as the pupil. So Nyquist sampling of the\n * intensity needs Q = lambda*(z/D)/pitch = M*dp/D >= 2.\n *\n * For a pupil spanning S samples, the discrete autocorrelation has 2S-1 lags.\n * So M >= 2S-1 is the exact condition for no circular wrap-around, and it\n * implies Q = M/(S-1) > 2. The declared grid is zero padded to that size.\n *\n * The non-zero support must also fit inside the declared f/N cone. Its\n * diameter D may not exceed z/N for the resolved propagation distance z; that\n * is, the image-side cone D/z may not be wider than 1/N. With the nominal\n * focal-plane reference (z = f) this is the entrance-aperture diameter f/N. An\n * exit pupil declared at its own distance z uses that distance\n * (docs/PSF_FOUNDATION.md).\n */\nfunction planPupilSampling(\n  profile: LensComplexPupilProfile,\n  propagationDistanceMm: number\n): PupilSamplingPlan {\n  const grid = profile.grid;\n  const supportX = supportSpan(grid, \"x\");\n  const supportY = supportSpan(grid, \"y\");\n  const diameterX = (supportX - 1) * grid.pupilSamplePitchMmX;\n  const diameterY = (supportY - 1) * grid.pupilSamplePitchMmY;\n  const declaredDiameter =\n    propagationDistanceMm / profile.context.apertureFNumber;\n  // Relative allowance for floating-point rounding only.\n  const limit = declaredDiameter * (1 + 1e-9);\n  if (diameterX > limit || diameterY > limit) {\n    throw new InvalidScientificInputError(\n      \"Complex-pupil amplitude support spans \" +\n        diameterX.toPrecision(6) +\n        \" x \" +\n        diameterY.toPrecision(6) +\n        \" mm between outermost non-zero samples, wider than the declared f/N cone diameter \" +\n        declaredDiameter.toPrecision(6) +\n        \" mm (propagationDistanceMm / apertureFNumber).\"\n    );\n  }\n  const transformWidth = Math.max(grid.widthSamples, 2 * supportX - 1);\n  const transformHeight = Math.max(grid.heightSamples, 2 * supportY - 1);\n  return {\n    sampling: {\n      criterion: \"discrete-intensity-alias-free-goodman-q-above-2\",\n      pupilSupportSamplesX: supportX,\n      pupilSupportSamplesY: supportY,\n      pupilSupportDiameterMmX: diameterX,\n      pupilSupportDiameterMmY: diameterY,\n      declaredPupilDiameterMm: declaredDiameter,\n      pupilGridWidthSamples: grid.widthSamples,\n      pupilGridHeightSamples: grid.heightSamples,\n      transformWidthSamples: transformWidth,\n      transformHeightSamples: transformHeight,\n      zeroPaddingApplied:\n        transformWidth !== grid.widthSamples ||\n        transformHeight !== grid.heightSamples,\n      samplingFactorQX:\n        supportX > 1 ? transformWidth / (supportX - 1) : null,\n      samplingFactorQY:\n        supportY > 1 ? transformHeight / (supportY - 1) : null\n    },\n    // Symmetric padding keeps zero frequency at the declared centre, shifted\n    // by the leading padding.\n    outputCenterX:\n      grid.centerSampleX +\n      Math.floor((transformWidth - grid.widthSamples) / 2),\n    outputCenterY:\n      grid.centerSampleY +\n      Math.floor((transformHeight - grid.heightSamples) / 2)\n  };\n}\n\nfunction twiddleTable(size: number): { cos: number[]; sin: number[] } {\n  const cos: number[] = [];\n  const sin: number[] = [];\n  for (let index = 0; index < size; index += 1) {\n    const angle = (-2 * Math.PI * index) / size;\n    cos.push(Math.cos(angle));\n    sin.push(Math.sin(angle));\n  }\n  return { cos, sin };\n}\n\nfunction positiveModulo(value: number, size: number): number {\n  const remainder = value % size;\n  return remainder < 0 ? remainder + size : remainder;\n}\n\nfunction calculateDiscreteFraunhoferKernel(\n  profile:\n    LensComplexPupilProfile,\n  propagationDistanceMm: number,\n  plan: PupilSamplingPlan\n): LensPsfKernel {\n  const grid = profile.grid;\n  const width = grid.widthSamples;\n  const height = grid.heightSamples;\n  const outWidth = plan.sampling.transformWidthSamples;\n  const outHeight = plan.sampling.transformHeightSamples;\n  const wavelengthMicrometers = profile.context.wavelengthNm / 1000;\n  const wavelengthMm = profile.context.wavelengthNm / 1_000_000;\n  const pupil: ComplexValue[] =\n    grid.relativeAmplitude.map(\n      (amplitude, index) =>\n        pupilComplexValue(\n          amplitude,\n          grid.opticalPathDifferenceMicrometers[index]!,\n          wavelengthMicrometers\n        )\n    );\n\n  // Separable zero-padded DFT: pupil rows to image kx first, then columns.\n  // Padding samples are zero, so only declared pupil samples are summed.\n  const twiddleX = twiddleTable(outWidth);\n  const twiddleY = twiddleTable(outHeight);\n  const rowRe = new Float64Array(height * outWidth);\n  const rowIm = new Float64Array(height * outWidth);\n  for (let y = 0; y < height; y += 1) {\n    for (let outX = 0; outX < outWidth; outX += 1) {\n      const kx = outX - plan.outputCenterX;\n      let re = 0;\n      let im = 0;\n      for (let x = 0; x < width; x += 1) {\n        const source = pupil[y * width + x]!;\n        const t = positiveModulo(kx * (x - grid.centerSampleX), outWidth);\n        const cos = twiddleX.cos[t]!;\n        const sin = twiddleX.sin[t]!;\n        re += source.re * cos - source.im * sin;\n        im += source.re * sin + source.im * cos;\n      }\n      rowRe[y * outWidth + outX] = re;\n      rowIm[y * outWidth + outX] = im;\n    }\n  }\n\n  const intensity = new Array<number>(outWidth * outHeight).fill(0);\n  for (let outY = 0; outY < outHeight; outY += 1) {\n    const ky = outY - plan.outputCenterY;\n    for (let outX = 0; outX < outWidth; outX += 1) {\n      let re = 0;\n      let im = 0;\n      for (let y = 0; y < height; y += 1) {\n        const sourceRe = rowRe[y * outWidth + outX]!;\n        const sourceIm = rowIm[y * outWidth + outX]!;\n        const t = positiveModulo(ky * (y - grid.centerSampleY), outHeight);\n        const cos = twiddleY.cos[t]!;\n        const sin = twiddleY.sin[t]!;\n        re += sourceRe * cos - sourceIm * sin;\n        im += sourceRe * sin + sourceIm * cos;\n      }\n      intensity[outY * outWidth + outX] = re * re + im * im;\n    }\n  }\n\n  const sum =\n    intensity.reduce(\n      (acc, value) =>\n        acc + value,\n      0\n    );\n  if (\n    !Number.isFinite(sum) ||\n    sum <= 0\n  ) {\n    throw new InvalidScientificInputError(\n      \"Complex-pupil propagation must produce finite positive PSF energy.\"\n    );\n  }\n  const normalized =\n    intensity.map(\n      (value) =>\n        value / sum\n    );\n\n  const samplePitchMmX =\n    wavelengthMm * propagationDistanceMm /\n    (outWidth * grid.pupilSamplePitchMmX);\n  const samplePitchMmY =\n    wavelengthMm * propagationDistanceMm /\n    (outHeight * grid.pupilSamplePitchMmY);\n\n  if (\n    !Number.isFinite(\n      samplePitchMmX\n    ) ||\n    !Number.isFinite(\n      samplePitchMmY\n    ) ||\n    samplePitchMmX <= 0 ||\n    samplePitchMmY <= 0\n  ) {\n    throw new InvalidScientificInputError(\n      \"Complex-pupil PSF sample pitch must remain finite and positive.\"\n    );\n  }\n\n  return {\n    widthSamples: outWidth,\n    heightSamples: outHeight,\n    samplePitchMicrometersX:\n      samplePitchMmX *\n      1000,\n    samplePitchMicrometersY:\n      samplePitchMmY *\n      1000,\n    centerSampleX: plan.outputCenterX,\n    centerSampleY: plan.outputCenterY,\n    normalizedIntensity:\n      normalized\n  };\n}\n\n/**\n * Evaluates a declared complex pupil with explicit amplitude, phase and sampling support. Returns the\n * normalized PSF shape, and the pupil throughput as a separate quantity.\n *\n * The complex pupil is evaluated by summing the declared amplitude and phase over a bounded, regular\n * pupil grid. The PSF shape is normalized separately from the pupil throughput. A single blur radius\n * cannot stand in for coherent phase, pupil energy, or the wavelength and field range where the\n * profile applies. Numeric profiles need provenance showing the data is owned or reusable.\n * @param input - CalculateLensComplexPupilPsfInput. See that type for coordinates, units and profile binding.\n * @returns CalculationResult<LensComplexPupilPsf>. The result states its shape and scientific status; a successful call does not imply calibration.\n *\n * @see docs/API_REFERENCE.md for the root export and exact type graph.\n */\nexport function calculateLensComplexPupilPsf(\n  input:\n    CalculateLensComplexPupilPsfInput\n): CalculationResult<LensComplexPupilPsf> {\n  requireInputKeys(input, \"calculateLensComplexPupilPsf\", [\"profile\", \"propagationDistanceMm\"]);\n  const profile =\n    parseLensComplexPupilProfile(\n      input.profile\n    );\n  const propagationDistanceMm = input.propagationDistanceMm === undefined ? profile.context.focalLengthMm : input.propagationDistanceMm;\n  requirePositiveFinite(propagationDistanceMm, \"propagationDistanceMm\");\n  const plan = planPupilSampling(profile, propagationDistanceMm);\n  const kernel =\n    calculateDiscreteFraunhoferKernel(\n      profile,\n      propagationDistanceMm,\n      plan\n    );\n\n  return approximationResult(\n    {\n      propagationDistanceMm,\n      propagationDistanceSource: input.propagationDistanceMm === undefined ? \"nominal-focal-plane-reference\" : \"explicit\",\n      profileId:\n        profile.profileId,\n      profileVersion:\n        profile.profileVersion,\n      scientificStatus:\n        profile.scientificStatus,\n      context: {\n        ...profile.context,\n        focus: {\n          ...profile.context.focus\n        },\n        fieldPointMm: {\n          ...profile.context\n            .fieldPointMm\n        }\n      },\n      kernel,\n      sampling: plan.sampling,\n      relativePupilThroughputFactor:\n        profile\n          .relativePupilThroughputFactor,\n      throughputOwnership:\n        \"separate-relative-pupil-throughput-factor\",\n      diffractionAndAberrationJointlyEvaluated:\n        true,\n      pupilClippingShapeIncluded:\n        profile\n          .responseIncludes\n          .pupilClippingShape,\n      pupilClippingThroughputAppliedToKernel:\n        false,\n      scalarFraunhoferReference:\n        true,\n      full2dOrientationPreserved:\n        true,\n      geometricDistortionModified:\n        false,\n      lateralChromaticPositionShiftApplied:\n        false,\n      sensorOpticalStackIncluded:\n        false,\n      sensorSamplingIncluded:\n        false,\n      reconstructionIncluded:\n        false,\n      strayLightIncluded:\n        false,\n      scalarSharpnessScoreProduced:\n        false,\n      componentEvidence: {\n        pupilProfile:\n          profile.evidence\n      },\n      uncertainty:\n        profile.uncertainty\n    },\n    \"lens-complex-pupil-fraunhofer-psf\",\n    \"2.0.0\",\n    [\n      \"The reference evaluator performs deterministic scalar Fraunhofer propagation of one explicit complex pupil amplitude plus optical-path-difference grid.\",\n      \"The pupil is zero padded so the DFT spans at least 2S-1 samples for a non-zero support of S samples per axis: the discrete intensity PSF then has no wrap-around aliasing and Goodman's sampling factor Q = lambda*N/pitch exceeds 2. Pupil support wider than the declared f/N aperture is rejected.\",\n      \"Physical sample pitch uses the explicit propagation distance when supplied; omission preserves the nominal focal-plane reference. Finite focus metadata does not resolve the pupil-to-image geometry automatically.\",\n      \"Diffraction and aberration phase are evaluated together in the same pupil propagation; callers must not stack an independent diffraction blur over this result.\",\n      \"The output PSF is unit-energy normalized. The explicit relativePupilThroughputFactor is not multiplied into the kernel and remains owned by the separate optical-throughput path.\",\n      \"Pupil amplitude shape can represent mechanical clipping/cat-eye geometry without silently encoding unknown transmission loss.\",\n      \"The pupil and image-plane axes preserve full 2D orientation; no radial or sagittal/tangential averaging is imposed.\",\n      \"This scalar Fraunhofer grid is a bounded reference approximation, not a full vector electromagnetic, polarization, stray-light or sensor-system model.\",\n      \"Sensor OLPF, microlens, CFA, sampling and reconstruction are downstream and excluded.\"\n    ]\n  );\n}\n"]}

package/dist/optics/depth-of-field.d.ts

@@ -1,12 +1,17 @@
 /**
- * Module boundary and integration notes.
- * Calculates thin-lens depth-of-field limits using the conventional hyperfocal-distance formulation.
- * The model is geometric and does not include diffraction, aberrations, focus breathing, pupil
- * magnification, or macro/high-magnification effects.
- * Calculates the geometric defocus-circle diameter at the sensor plane for a subject plane away from
- * the selected focus plane.
- * @see docs/PHYSICS_FOUNDATION.md for equations, coordinate/unit conventions, blockers and support
- * limits.
+ * Thin-lens depth of field (DoF) and defocus blur.
+ *
+ * Depth-of-field limits are exact for a thin lens. They use the hyperfocal distance
+ * H = f^2/(Nc) + f:
+ *   near = s(H - f)/(H + s - 2f)
+ *   far  = s(H - f)/(H - s), infinite (null) when s >= H
+ * The model is geometric. It does not include diffraction, aberrations, focus breathing, pupil
+ * magnification or macro (high-magnification) effects.
+ *
+ * The defocus circle is the geometric blur-circle diameter at the sensor for a subject plane that is
+ * not at the focus distance.
+ * @see docs/PHYSICS_FOUNDATION.md for the equations, units, coordinate conventions, known blockers
+ * and support limits.
  */
 import { type CalculationResult } from "../core/calculation-result.js";
 export interface CalculateDepthOfFieldInput {
@@ -28,11 +33,16 @@ export interface DepthOfField {
     totalDepthOfFieldM: number | null;
 }
 /**
- * Calculates thin-lens depth-of-field limits using the conventional
- * hyperfocal-distance formulation.
+ * Calculates exact thin-lens depth-of-field limits. With hyperfocal
+ * distance H = f^2/(Nc) + f and focus distance s, the near limit is
+ * s(H - f)/(H + s - 2f) and the far limit is s(H - f)/(H - s), which is
+ * infinite (`null`) when s >= H. At both limits the geometric defocus circle
+ * from {@link calculateDefocusCircle} equals the supplied circle of confusion.
+ * The circle of confusion must be smaller than the aperture diameter f/N;
+ * otherwise the near limit would fall inside the focal length.
  *
- * The model is geometric and does not include diffraction, aberrations,
- * focus breathing, pupil magnification, or macro/high-magnification effects.
+ * The model is geometric. It does not include diffraction, aberrations,
+ * focus breathing, pupil magnification or macro (high-magnification) effects.
  *
  * @param input Lens, focus, aperture, and acceptable blur-circle inputs.
  * @returns Hyperfocal, near, and far depth-of-field distances.
@@ -54,13 +64,14 @@ export interface DefocusCircle {
 }
 /**
  * Calculates the geometric defocus-circle diameter at the sensor plane for a
- * subject plane away from the selected focus plane.
+ * subject plane away from the selected focus plane, using the stable closed
+ * form c = f^2 |s - d| / (N d (s - f)) (model 1.1.0).
  *
  * @param input Lens, aperture, focus-plane, and subject-plane distances.
  * @returns Geometric blur-circle diameter.
  */
 export declare function calculateDefocusCircle(input: CalculateDefocusCircleInput): CalculationResult<DefocusCircle>;
-/** @internal Immutable request-local identity; never a root public cache/API. */
+/** @internal Immutable identity for one request. Not a public cache or API. */
 export interface ResolvedDefocusFocusProjection {
     readonly focalLengthMm: number;
     readonly focusDistanceM: number;
@@ -68,4 +79,3 @@ export interface ResolvedDefocusFocusProjection {
 }
 /** @internal Reuses a projection already calculated for this exact POC request. */
 export declare function calculateDefocusCircleUsingFocusProjection(input: CalculateDefocusCircleInput, projection?: ResolvedDefocusFocusProjection): CalculationResult<DefocusCircle>;
-//# sourceMappingURL=depth-of-field.d.ts.map