@photivra/engine 2.9.0 → 2.10.0

package/dist/optics/irradiance-convolution.js.map

@@ -1 +1 @@
-{"version":3,"file":"irradiance-convolution.js","sourceRoot":"","sources":["../../src/optics/irradiance-convolution.ts"],"names":[],"mappings":"AAAA,sCAAsC;AAEtC;;;;;;;;;;;;;;;;;;;;;;;;;;GA0BG;AACH,OAAO,EAAE,yBAAyB,EAAE,MAAM,gCAAgC,CAAC;AAC3E,OAAO,EAAE,wBAAwB,EAAE,MAAM,8BAA8B,CAAC;AACxE,OAAO,EAAE,8BAA8B,EAAE,kBAAkB,EAAE,MAAM,yBAAyB,CAAC;AAE7F,MAAM,CAAC,MAAM,qCAAqC,GAAG,OAAgB,CAAC;AACtE,6HAA6H;AAC7H,MAAM,CAAC,MAAM,6BAA6B,GAAG,MAAM,CAAC,MAAM,CAAC;IACzD,aAAa,EAAE,UAAU;IACzB,aAAa,EAAE,8BAA8B;IAC7C,gBAAgB,EAAE,MAAM;IACxB,mBAAmB,EAAE,EAAE;IACvB,iBAAiB,EAAE,KAAK;IACxB,oBAAoB,EAAE,SAAS;IAC/B,wGAAwG;IACxG,gBAAgB,EAAE,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE;IAC9B,sFAAsF;IACtF,sBAAsB,EAAE,CAAC,IAAI,GAAG;CACjC,CAAC,CAAC;AA0CH,SAAS,eAAe,CAAC,KAAc,EAAE,OAAe;IACtD,IAAI,OAAO,KAAK,KAAK,QAAQ,IAAI,CAAC,MAAM,CAAC,aAAa,CAAC,KAAK,CAAC,IAAI,KAAK,IAAI,CAAC,IAAI,KAAK,GAAG,OAAO,EAAE,CAAC;QAC/F,MAAM,IAAI,yBAAyB,CAAC,0CAA0C,CAAC,CAAC;IAClF,CAAC;IACD,OAAO,KAAK,CAAC;AACf,CAAC;AAED,2FAA2F;AAC3F,MAAM,UAAU,kCAAkC,CAAC,KAAiC;IAClF,MAAM,CAAC,GAAG,wBAAwB,CAAC,KAAK,EAAE,CAAC,eAAe,EAAE,OAAO,EAAE,QAAQ,EAAE,SAAS,EAAE,cAAc,EAAE,YAAY;QACpH,aAAa,EAAE,SAAS,EAAE,UAAU,EAAE,eAAe,CAAC,EAAE,uCAAuC,CAAC,CAAC;IACnG,MAAM,CAAC,GAAG,6BAA6B,CAAC;IACxC,IAAI,CAAC,CAAC,aAAa,KAAK,qCAAqC,IAAI,CAAC,CAAC,QAAQ,KAAK,mBAAmB,EAAE,CAAC;QACpG,MAAM,IAAI,yBAAyB,CAAC,wDAAwD,CAAC,CAAC;IAChG,CAAC;IACD,MAAM,KAAK,GAAG,eAAe,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,gBAAgB,CAAC,EAAE,MAAM,GAAG,eAAe,CAAC,CAAC,CAAC,MAAM,EAAE,CAAC,CAAC,gBAAgB,CAAC,CAAC;IACnH,MAAM,aAAa,GAAG,kBAAkB,CAAC,CAAC,CAAC,aAAa,EAAE,CAAC,CAAC,aAAa,EAAE,CAAC,CAAC,aAAa,EACxF,oEAAoE,GAAG,CAAC,CAAC,aAAa,GAAG,GAAG,CAAC,CAAC;IAChG,IAAI,KAAK,GAAG,MAAM,GAAG,aAAa;QAAE,MAAM,IAAI,yBAAyB,CAAC,sDAAsD,CAAC,CAAC;IAChI,MAAM,MAAM,GAAG,CAAC,CAAC,YAAY,CAAC;IAC9B,IAAI,OAAO,MAAM,KAAK,QAAQ,IAAI,CAAC,MAAM,CAAC,aAAa,CAAC,MAAM,CAAC,IAAI,MAAM,GAAG,CAAC,IAAI,MAAM,GAAG,CAAC,CAAC,mBAAmB,EAAE,CAAC;QAChH,MAAM,IAAI,yBAAyB,CAAC,+CAA+C,CAAC,CAAC;IACvF,CAAC;IACD,MAAM,UAAU,GAAG,eAAe,CAAC,CAAC,CAAC,UAAU,EAAE,KAAK,CAAC,EAAE,WAAW,GAAG,eAAe,CAAC,CAAC,CAAC,WAAW,EAAE,MAAM,CAAC,CAAC;IAC9G,MAAM,YAAY,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,GAAG,UAAU,CAAC,EAAE,SAAS,GAAG,IAAI,CAAC,IAAI,CAAC,MAAM,GAAG,WAAW,CAAC,CAAC;IAChG,MAAM,UAAU,GAAG,CAAC,GAAG,MAAM,GAAG,CAAC,EAAE,QAAQ,GAAG,UAAU,GAAG,UAAU,EAAE,UAAU,GAAG,YAAY,GAAG,SAAS,CAAC;IAC7G,IAAI,UAAU,GAAG,CAAC,CAAC,iBAAiB,IAAI,UAAU,GAAG,QAAQ,GAAG,CAAC,CAAC,oBAAoB,EAAE,CAAC;QACvF,MAAM,IAAI,yBAAyB,CAAC,4DAA4D,CAAC,CAAC;IACpG,CAAC;IACD,MAAM,OAAO,GAAG,CAAC,CAAC,OAAO,EAAE,OAAO,GAAG,CAAC,CAAC,OAAO,CAAC;IAC/C,IAAI,CAAC,CAAC,OAAO,YAAY,YAAY,IAAI,OAAO,YAAY,YAAY,CAAC,IAAI,OAAO,CAAC,MAAM,KAAK,KAAK,GAAG,MAAM;QAC1G,CAAC,CAAC,OAAO,YAAY,YAAY,IAAI,OAAO,YAAY,YAAY,CAAC,IAAI,OAAO,CAAC,MAAM,KAAK,UAAU,GAAG,QAAQ,EAAE,CAAC;QACtH,MAAM,IAAI,yBAAyB,CAAC,mEAAmE,CAAC,CAAC;IAC3G,CAAC;IACD,IAAI,sBAAsB,GAAG,CAAC,CAAC;IAC/B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACxC,MAAM,CAAC,GAAG,OAAO,CAAC,CAAC,CAAE,CAAC;QACtB,uCAAuC;QACvC,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,sBAAsB,CAAC,EAAE,CAAC;YAC/C,MAAM,IAAI,yBAAyB,CAAC,oFAAoF,CAAC,CAAC;QAC5H,CAAC;QACD,IAAI,CAAC,GAAG,sBAAsB;YAAE,sBAAsB,GAAG,CAAC,CAAC;IAC7D,CAAC;IACD,IAAI,gBAAgB,GAAG,CAAC,CAAC;IACzB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,UAAU,EAAE,CAAC,EAAE,EAAE,CAAC;QACpC,IAAI,GAAG,GAAG,CAAC,CAAC;QACZ,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,QAAQ,EAAE,GAAG,GAAG,CAAC,GAAG,QAAQ,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;YAC5D,MAAM,CAAC,GAAG,OAAO,CAAC,CAAC,CAAE,CAAC;YACtB,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;gBAAE,MAAM,IAAI,yBAAyB,CAAC,sDAAsD,CAAC,CAAC;YACrH,GAAG,IAAI,CAAC,CAAC;QACX,CAAC;QACD,IAAI,GAAG,GAAG,CAAC,CAAC,gBAAgB;YAAE,MAAM,IAAI,yBAAyB,CAAC,qDAAqD,CAAC,CAAC;QACzH,IAAI,GAAG,GAAG,gBAAgB;YAAE,gBAAgB,GAAG,GAAG,CAAC;IACrD,CAAC;IACD,OAAO,MAAM,CAAC,MAAM,CAAC,EAAE,KAAK,EAAE,MAAM,EAAE,YAAY,EAAE,MAAM,EAAE,UAAU,EAAE,QAAQ,EAAE,UAAU,EAAE,WAAW,EAAE,YAAY,EAAE,SAAS;QAChI,gBAAgB,EAAE,sBAAsB,EAAE,CAAC,CAAC;AAChD,CAAC;AAED;;;;;GAKG;AACH,MAAM,UAAU,sBAAsB,CAAC,IAA+B,EAAE,OAAoC,EAC1G,OAAoC,EAAE,EAAU,EAAE,EAAU,EAAE,MAAoB;IAClF,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,YAAY,EAAE,CAAC,EAAE,UAAU,EAAE,QAAQ,EAAE,UAAU,EAAE,WAAW,EAAE,YAAY,EAAE,SAAS,EAAE,GAAG,IAAI,CAAC;IACxH,KAAK,IAAI,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;QAC7B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;YAC/B,IAAI,GAAG,GAAG,CAAC,CAAC;YACZ,KAAK,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC7B,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;gBACjB,IAAI,EAAE,GAAG,CAAC,IAAI,EAAE,IAAI,MAAM;oBAAE,SAAS;gBACrC,MAAM,YAAY,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,GAAG,WAAW,CAAC,EAAE,SAAS,GAAG,CAAC,CAAC,GAAG,YAAY,CAAC;gBAC1F,MAAM,SAAS,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,UAAU,GAAG,CAAC,EAAE,SAAS,GAAG,EAAE,GAAG,KAAK,CAAC;gBACnE,KAAK,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC7B,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;oBACjB,IAAI,EAAE,GAAG,CAAC,IAAI,EAAE,IAAI,KAAK;wBAAE,SAAS;oBACpC,MAAM,KAAK,GAAG,YAAY,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,GAAG,UAAU,CAAC,EAAE,YAAY,GAAG,CAAC,CAAC,CAAC;oBACrF,GAAG,IAAI,OAAO,CAAC,SAAS,GAAG,EAAE,CAAE,GAAG,OAAO,CAAC,KAAK,GAAG,QAAQ,GAAG,SAAS,GAAG,CAAC,CAAE,CAAC;gBAC/E,CAAC;YACH,CAAC;YACD,MAAM,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC;QACrC,CAAC;IACH,CAAC;AACH,CAAC;AAED,yGAAyG;AACzG,MAAM,UAAU,8BAA8B,CAAC,KAAiC;IAC9E,MAAM,IAAI,GAAG,kCAAkC,CAAC,KAAK,CAAC,CAAC;IACvD,MAAM,MAAM,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC;IAC1D,sBAAsB,CAAC,IAAI,EAAE,KAAK,CAAC,OAAO,EAAE,KAAK,CAAC,OAAO,EAAE,CAAC,EAAE,IAAI,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IACnF,OAAO,EAAE,IAAI,EAAE,MAAM,EAAE,CAAC;AAC1B,CAAC","sourcesContent":["// SPDX-License-Identifier: Apache-2.0\n\n/**\n * Blurs a sampled image-plane quantity with caller-supplied PSF (point-spread function) kernels, one\n * kernel per patch of the image. The quantity must be linear and additive (for example irradiance).\n * Each patch is treated as shift-invariant (\"isoplanatic\"): the same PSF applies everywhere inside it.\n * The PSFs are incoherent intensity PSFs.\n *\n * This is the CPU float64 reference for the optics-as-convolution step (#277). The `webgpu-f32`\n * path runs the same equation; see docs/WEBGPU.md.\n *\n * Equation (indexed by source site; an energy-conserving scatter written as a gather):\n *   out(x, y) = Σ_{v=-r..r} Σ_{u=-r..r} in(x − u, y − v) · K_{p(x − u, y − v)}(u, v)\n * where `p(s)` is the patch containing the **source** site and `K_p(u, v)` is the fraction of the source\n * site's quantity that lands at displacement (+u right, +v down), in grid sites.\n *\n * Incoherent imaging is a linear sum of intensity PSFs: the output intensity is the input intensity\n * convolved with the intensity PSF (MIT OpenCourseWare 2.71 Optics, Spring 2009, lecture 22, free).\n * A PSF that varies across the field is approximated as constant inside each patch, as in methods\n * that divide the image into isoplanatic subsections (van de Ketterij et al., arXiv:2011.01738,\n * sec. 1, free). Also see: Goodman, Introduction to Fourier Optics, 4th ed., 2017, §7.1. Camera\n * simulators such as ISETCam use the same piecewise shift-invariant method\n * (Farrell et al., Proc. SPIE 8299, 2012). The caller chooses the patch size and kernels, and so\n * owns the patch approximation.\n *\n * Edges: the grid holds every source. Sources outside the grid are not represented. Anything\n * scattered outside the grid is lost; compare `inputSum` with `outputSum` to see how much. There is no\n * padding, reflection or renormalization.\n */\nimport { InvalidConfigurationError } from \"../core/configuration-error.js\";\nimport { requireAllowlistedRecord } from \"../core/record-validation.js\";\nimport { FULL_RESOLUTION_CEILING_PIXELS, resolveRaisedLimit } from \"../core/pixel-limits.js\";\n\nexport const IRRADIANCE_CONVOLUTION_SCHEMA_VERSION = \"0.1.0\" as const;\n/** `maximumPixels` is the default grid limit; a request may raise it with its own `maximumPixels`, up to `ceilingPixels`. */\nexport const IRRADIANCE_CONVOLUTION_LIMITS = Object.freeze({\n  maximumPixels: 24_000_000,\n  ceilingPixels: FULL_RESOLUTION_CEILING_PIXELS,\n  maximumDimension: 16_384,\n  maximumKernelRadius: 32,\n  maximumPatchCount: 4_096,\n  maximumKernelWeights: 4_194_304,\n  /** A kernel may lose energy (vignetting, transmission) but may not create it, beyond rounding slack. */\n  maximumKernelSum: 1 + 2 ** -20,\n  /** Inputs larger than this are rejected so that float32 execution cannot overflow. */\n  maximumSampleMagnitude: 2 ** 100\n});\n\nexport interface IrradianceConvolutionInput {\n  schemaVersion: typeof IRRADIANCE_CONVOLUTION_SCHEMA_VERSION;\n  width: number;\n  height: number;\n  /** Row-major samples. Finite, non-negative, linear and additive (for example irradiance, or expected photons per site). */\n  samples: Float32Array | Float64Array;\n  /** Kernel half-width in grid sites. Each kernel is (2r+1)×(2r+1). */\n  kernelRadius: number;\n  /** Patch size in grid sites. The last row and column of patches may be partial. */\n  patchWidth: number;\n  patchHeight: number;\n  /**\n   * Finite non-negative weights, `[patchRow][patchColumn][v + r][u + r]`, one kernel per patch in\n   * row-major patch order. Each kernel sums to at most 1 (+2^-20 rounding slack).\n   */\n  kernels: Float32Array | Float64Array;\n  boundary: \"zero-outside-grid\";\n  /**\n   * Optional opt-in grid limit for large images: an integer from 1 to `IRRADIANCE_CONVOLUTION_LIMITS.ceilingPixels`.\n   * Default `IRRADIANCE_CONVOLUTION_LIMITS.maximumPixels` (24 000 000). It does not change the result.\n   */\n  maximumPixels?: number;\n}\n\n/** Validated, immutable geometry of a convolution request. Sample arrays are borrowed, not copied. */\nexport interface IrradianceConvolutionPlan {\n  width: number;\n  height: number;\n  kernelRadius: number;\n  kernelSize: number;\n  tapCount: number;\n  patchWidth: number;\n  patchHeight: number;\n  patchColumns: number;\n  patchRows: number;\n  /** Largest kernel sum and largest input magnitude. The float32 error bound uses them. */\n  maximumKernelSum: number;\n  maximumSampleMagnitude: number;\n}\n\nfunction positiveInteger(value: unknown, maximum: number): number {\n  if (typeof value !== \"number\" || !Number.isSafeInteger(value) || value <= 0 || value > maximum) {\n    throw new InvalidConfigurationError(\"Invalid irradiance convolution geometry.\");\n  }\n  return value;\n}\n\n/** Validates the request and returns its plan. Every sample and weight is checked once. */\nexport function calculateIrradianceConvolutionPlan(input: IrradianceConvolutionInput): IrradianceConvolutionPlan {\n  const r = requireAllowlistedRecord(input, [\"schemaVersion\", \"width\", \"height\", \"samples\", \"kernelRadius\", \"patchWidth\",\n    \"patchHeight\", \"kernels\", \"boundary\", \"maximumPixels\"], \"Invalid irradiance convolution input.\");\n  const L = IRRADIANCE_CONVOLUTION_LIMITS;\n  if (r.schemaVersion !== IRRADIANCE_CONVOLUTION_SCHEMA_VERSION || r.boundary !== \"zero-outside-grid\") {\n    throw new InvalidConfigurationError(\"Unsupported irradiance convolution schema or boundary.\");\n  }\n  const width = positiveInteger(r.width, L.maximumDimension), height = positiveInteger(r.height, L.maximumDimension);\n  const maximumPixels = resolveRaisedLimit(r.maximumPixels, L.maximumPixels, L.ceilingPixels,\n    \"Irradiance convolution maximumPixels must be an integer from 1 to \" + L.ceilingPixels + \".\");\n  if (width * height > maximumPixels) throw new InvalidConfigurationError(\"Irradiance convolution grid exceeds the pixel limit.\");\n  const radius = r.kernelRadius;\n  if (typeof radius !== \"number\" || !Number.isSafeInteger(radius) || radius < 0 || radius > L.maximumKernelRadius) {\n    throw new InvalidConfigurationError(\"Invalid irradiance convolution kernel radius.\");\n  }\n  const patchWidth = positiveInteger(r.patchWidth, width), patchHeight = positiveInteger(r.patchHeight, height);\n  const patchColumns = Math.ceil(width / patchWidth), patchRows = Math.ceil(height / patchHeight);\n  const kernelSize = 2 * radius + 1, tapCount = kernelSize * kernelSize, patchCount = patchColumns * patchRows;\n  if (patchCount > L.maximumPatchCount || patchCount * tapCount > L.maximumKernelWeights) {\n    throw new InvalidConfigurationError(\"Irradiance convolution kernels exceed the declared limits.\");\n  }\n  const samples = r.samples, kernels = r.kernels;\n  if (!(samples instanceof Float32Array || samples instanceof Float64Array) || samples.length !== width * height ||\n      !(kernels instanceof Float32Array || kernels instanceof Float64Array) || kernels.length !== patchCount * tapCount) {\n    throw new InvalidConfigurationError(\"Irradiance convolution arrays do not match the declared geometry.\");\n  }\n  let maximumSampleMagnitude = 0;\n  for (let i = 0; i < samples.length; i++) {\n    const v = samples[i]!;\n    // Negated comparison also rejects NaN.\n    if (!(v >= 0 && v <= L.maximumSampleMagnitude)) {\n      throw new InvalidConfigurationError(\"Irradiance samples must be finite, non-negative and within the declared magnitude.\");\n    }\n    if (v > maximumSampleMagnitude) maximumSampleMagnitude = v;\n  }\n  let maximumKernelSum = 0;\n  for (let p = 0; p < patchCount; p++) {\n    let sum = 0;\n    for (let t = p * tapCount, end = t + tapCount; t < end; t++) {\n      const w = kernels[t]!;\n      if (!(w >= 0 && w <= 1)) throw new InvalidConfigurationError(\"PSF kernel weights must be finite and within [0, 1].\");\n      sum += w;\n    }\n    if (sum > L.maximumKernelSum) throw new InvalidConfigurationError(\"A PSF kernel may not create energy (sum exceeds 1).\");\n    if (sum > maximumKernelSum) maximumKernelSum = sum;\n  }\n  return Object.freeze({ width, height, kernelRadius: radius, kernelSize, tapCount, patchWidth, patchHeight, patchColumns, patchRows,\n    maximumKernelSum, maximumSampleMagnitude });\n}\n\n/**\n * Reference gather for output rows `[y0, y1)`, written into `output` at `(y - y0) * width + x`.\n * The summation order (v outer, u inner, both ascending) is part of the contract shared with the\n * WGSL (WebGPU shader) kernel.\n * @internal\n */\nexport function convolveIrradianceRows(plan: IrradianceConvolutionPlan, samples: Float32Array | Float64Array,\n  kernels: Float32Array | Float64Array, y0: number, y1: number, output: Float64Array): void {\n  const { width, height, kernelRadius: r, kernelSize, tapCount, patchWidth, patchHeight, patchColumns, patchRows } = plan;\n  for (let y = y0; y < y1; y++) {\n    for (let x = 0; x < width; x++) {\n      let sum = 0;\n      for (let v = -r; v <= r; v++) {\n        const sy = y - v;\n        if (sy < 0 || sy >= height) continue;\n        const patchRowBase = Math.min(Math.floor(sy / patchHeight), patchRows - 1) * patchColumns;\n        const kernelRow = (v + r) * kernelSize + r, sampleRow = sy * width;\n        for (let u = -r; u <= r; u++) {\n          const sx = x - u;\n          if (sx < 0 || sx >= width) continue;\n          const patch = patchRowBase + Math.min(Math.floor(sx / patchWidth), patchColumns - 1);\n          sum += samples[sampleRow + sx]! * kernels[patch * tapCount + kernelRow + u]!;\n        }\n      }\n      output[(y - y0) * width + x] = sum;\n    }\n  }\n}\n\n/** Synchronous float64 reference for bounded grids. Use the image-kernel accelerator for large grids. */\nexport function calculateIrradianceConvolution(input: IrradianceConvolutionInput): { plan: IrradianceConvolutionPlan; values: Float64Array } {\n  const plan = calculateIrradianceConvolutionPlan(input);\n  const values = new Float64Array(plan.width * plan.height);\n  convolveIrradianceRows(plan, input.samples, input.kernels, 0, plan.height, values);\n  return { plan, values };\n}\n"]}
+{"version":3,"file":"irradiance-convolution.js","sourceRoot":"","sources":["../../src/optics/irradiance-convolution.ts"],"names":[],"mappings":"AAAA,sCAAsC;AAEtC;;;;;;;;;;;;;;;;;;;;;;;;;;GA0BG;AACH,OAAO,EAAE,yBAAyB,EAAE,MAAM,gCAAgC,CAAC;AAC3E,OAAO,EAAE,wBAAwB,EAAE,MAAM,8BAA8B,CAAC;AACxE,OAAO,EAAE,8BAA8B,EAAE,kBAAkB,EAAE,MAAM,yBAAyB,CAAC;AAE7F,MAAM,CAAC,MAAM,qCAAqC,GAAG,OAAgB,CAAC;AACtE,6HAA6H;AAC7H,MAAM,CAAC,MAAM,6BAA6B,GAAG,MAAM,CAAC,MAAM,CAAC;IACzD,aAAa,EAAE,UAAU;IACzB,aAAa,EAAE,8BAA8B;IAC7C,gBAAgB,EAAE,MAAM;IACxB,mBAAmB,EAAE,EAAE;IACvB,iBAAiB,EAAE,KAAK;IACxB,oBAAoB,EAAE,SAAS;IAC/B,wGAAwG;IACxG,gBAAgB,EAAE,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE;IAC9B,sFAAsF;IACtF,sBAAsB,EAAE,CAAC,IAAI,GAAG;CACjC,CAAC,CAAC;AA0CH,SAAS,eAAe,CAAC,KAAc,EAAE,OAAe;IACtD,IAAI,OAAO,KAAK,KAAK,QAAQ,IAAI,CAAC,MAAM,CAAC,aAAa,CAAC,KAAK,CAAC,IAAI,KAAK,IAAI,CAAC,IAAI,KAAK,GAAG,OAAO,EAAE,CAAC;QAC/F,MAAM,IAAI,yBAAyB,CAAC,0CAA0C,CAAC,CAAC;IAClF,CAAC;IACD,OAAO,KAAK,CAAC;AACf,CAAC;AAED,yGAAyG;AACzG,SAAS,6BAA6B,CAAC,OAAoC;IACzE,IAAI,sBAAsB,GAAG,CAAC,CAAC;IAC/B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACxC,MAAM,CAAC,GAAG,OAAO,CAAC,CAAC,CAAE,CAAC;QACtB,uCAAuC;QACvC,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,6BAA6B,CAAC,sBAAsB,CAAC,EAAE,CAAC;YAC3E,MAAM,IAAI,yBAAyB,CAAC,oFAAoF,CAAC,CAAC;QAC5H,CAAC;QACD,IAAI,CAAC,GAAG,sBAAsB;YAAE,sBAAsB,GAAG,CAAC,CAAC;IAC7D,CAAC;IACD,OAAO,sBAAsB,CAAC;AAChC,CAAC;AAED,wHAAwH;AACxH,SAAS,uBAAuB,CAAC,OAAoC,EAAE,UAAkB,EAAE,QAAgB;IACzG,IAAI,gBAAgB,GAAG,CAAC,CAAC;IACzB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,UAAU,EAAE,CAAC,EAAE,EAAE,CAAC;QACpC,IAAI,GAAG,GAAG,CAAC,CAAC;QACZ,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,QAAQ,EAAE,GAAG,GAAG,CAAC,GAAG,QAAQ,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;YAC5D,MAAM,CAAC,GAAG,OAAO,CAAC,CAAC,CAAE,CAAC;YACtB,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;gBAAE,MAAM,IAAI,yBAAyB,CAAC,sDAAsD,CAAC,CAAC;YACrH,GAAG,IAAI,CAAC,CAAC;QACX,CAAC;QACD,IAAI,GAAG,GAAG,6BAA6B,CAAC,gBAAgB;YAAE,MAAM,IAAI,yBAAyB,CAAC,qDAAqD,CAAC,CAAC;QACrJ,IAAI,GAAG,GAAG,gBAAgB;YAAE,gBAAgB,GAAG,GAAG,CAAC;IACrD,CAAC;IACD,OAAO,gBAAgB,CAAC;AAC1B,CAAC;AAED,2FAA2F;AAC3F,MAAM,UAAU,kCAAkC,CAAC,KAAiC;IAClF,MAAM,CAAC,GAAG,wBAAwB,CAAC,KAAK,EAAE,CAAC,eAAe,EAAE,OAAO,EAAE,QAAQ,EAAE,SAAS,EAAE,cAAc,EAAE,YAAY;QACpH,aAAa,EAAE,SAAS,EAAE,UAAU,EAAE,eAAe,CAAC,EAAE,uCAAuC,CAAC,CAAC;IACnG,MAAM,CAAC,GAAG,6BAA6B,CAAC;IACxC,IAAI,CAAC,CAAC,aAAa,KAAK,qCAAqC,IAAI,CAAC,CAAC,QAAQ,KAAK,mBAAmB,EAAE,CAAC;QACpG,MAAM,IAAI,yBAAyB,CAAC,wDAAwD,CAAC,CAAC;IAChG,CAAC;IACD,MAAM,KAAK,GAAG,eAAe,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,gBAAgB,CAAC,EAAE,MAAM,GAAG,eAAe,CAAC,CAAC,CAAC,MAAM,EAAE,CAAC,CAAC,gBAAgB,CAAC,CAAC;IACnH,MAAM,aAAa,GAAG,kBAAkB,CAAC,CAAC,CAAC,aAAa,EAAE,CAAC,CAAC,aAAa,EAAE,CAAC,CAAC,aAAa,EACxF,oEAAoE,GAAG,CAAC,CAAC,aAAa,GAAG,GAAG,CAAC,CAAC;IAChG,IAAI,KAAK,GAAG,MAAM,GAAG,aAAa;QAAE,MAAM,IAAI,yBAAyB,CAAC,sDAAsD,CAAC,CAAC;IAChI,MAAM,MAAM,GAAG,CAAC,CAAC,YAAY,CAAC;IAC9B,IAAI,OAAO,MAAM,KAAK,QAAQ,IAAI,CAAC,MAAM,CAAC,aAAa,CAAC,MAAM,CAAC,IAAI,MAAM,GAAG,CAAC,IAAI,MAAM,GAAG,CAAC,CAAC,mBAAmB,EAAE,CAAC;QAChH,MAAM,IAAI,yBAAyB,CAAC,+CAA+C,CAAC,CAAC;IACvF,CAAC;IACD,MAAM,UAAU,GAAG,eAAe,CAAC,CAAC,CAAC,UAAU,EAAE,KAAK,CAAC,EAAE,WAAW,GAAG,eAAe,CAAC,CAAC,CAAC,WAAW,EAAE,MAAM,CAAC,CAAC;IAC9G,MAAM,YAAY,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,GAAG,UAAU,CAAC,EAAE,SAAS,GAAG,IAAI,CAAC,IAAI,CAAC,MAAM,GAAG,WAAW,CAAC,CAAC;IAChG,MAAM,UAAU,GAAG,CAAC,GAAG,MAAM,GAAG,CAAC,EAAE,QAAQ,GAAG,UAAU,GAAG,UAAU,EAAE,UAAU,GAAG,YAAY,GAAG,SAAS,CAAC;IAC7G,IAAI,UAAU,GAAG,CAAC,CAAC,iBAAiB,IAAI,UAAU,GAAG,QAAQ,GAAG,CAAC,CAAC,oBAAoB,EAAE,CAAC;QACvF,MAAM,IAAI,yBAAyB,CAAC,4DAA4D,CAAC,CAAC;IACpG,CAAC;IACD,MAAM,OAAO,GAAG,CAAC,CAAC,OAAO,EAAE,OAAO,GAAG,CAAC,CAAC,OAAO,CAAC;IAC/C,IAAI,CAAC,CAAC,OAAO,YAAY,YAAY,IAAI,OAAO,YAAY,YAAY,CAAC,IAAI,OAAO,CAAC,MAAM,KAAK,KAAK,GAAG,MAAM;QAC1G,CAAC,CAAC,OAAO,YAAY,YAAY,IAAI,OAAO,YAAY,YAAY,CAAC,IAAI,OAAO,CAAC,MAAM,KAAK,UAAU,GAAG,QAAQ,EAAE,CAAC;QACtH,MAAM,IAAI,yBAAyB,CAAC,mEAAmE,CAAC,CAAC;IAC3G,CAAC;IACD,MAAM,sBAAsB,GAAG,6BAA6B,CAAC,OAAO,CAAC,CAAC;IACtE,MAAM,gBAAgB,GAAG,uBAAuB,CAAC,OAAO,EAAE,UAAU,EAAE,QAAQ,CAAC,CAAC;IAChF,OAAO,MAAM,CAAC,MAAM,CAAC,EAAE,KAAK,EAAE,MAAM,EAAE,YAAY,EAAE,MAAM,EAAE,UAAU,EAAE,QAAQ,EAAE,UAAU,EAAE,WAAW,EAAE,YAAY,EAAE,SAAS;QAChI,gBAAgB,EAAE,sBAAsB,EAAE,CAAC,CAAC;AAChD,CAAC;AAED;;;;;GAKG;AACH,MAAM,UAAU,sBAAsB,CAAC,IAA+B,EAAE,OAAoC,EAC1G,OAAoC,EAAE,EAAU,EAAE,EAAU,EAAE,MAAoB;IAClF,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,YAAY,EAAE,CAAC,EAAE,UAAU,EAAE,QAAQ,EAAE,UAAU,EAAE,WAAW,EAAE,YAAY,EAAE,SAAS,EAAE,GAAG,IAAI,CAAC;IACxH,KAAK,IAAI,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;QAC7B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;YAC/B,IAAI,GAAG,GAAG,CAAC,CAAC;YACZ,KAAK,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC7B,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;gBACjB,IAAI,EAAE,GAAG,CAAC,IAAI,EAAE,IAAI,MAAM;oBAAE,SAAS;gBACrC,MAAM,YAAY,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,GAAG,WAAW,CAAC,EAAE,SAAS,GAAG,CAAC,CAAC,GAAG,YAAY,CAAC;gBAC1F,MAAM,SAAS,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,UAAU,GAAG,CAAC,EAAE,SAAS,GAAG,EAAE,GAAG,KAAK,CAAC;gBACnE,KAAK,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC7B,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;oBACjB,IAAI,EAAE,GAAG,CAAC,IAAI,EAAE,IAAI,KAAK;wBAAE,SAAS;oBACpC,MAAM,KAAK,GAAG,YAAY,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,GAAG,UAAU,CAAC,EAAE,YAAY,GAAG,CAAC,CAAC,CAAC;oBACrF,GAAG,IAAI,OAAO,CAAC,SAAS,GAAG,EAAE,CAAE,GAAG,OAAO,CAAC,KAAK,GAAG,QAAQ,GAAG,SAAS,GAAG,CAAC,CAAE,CAAC;gBAC/E,CAAC;YACH,CAAC;YACD,MAAM,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC;QACrC,CAAC;IACH,CAAC;AACH,CAAC;AAED,yGAAyG;AACzG,MAAM,UAAU,8BAA8B,CAAC,KAAiC;IAC9E,MAAM,IAAI,GAAG,kCAAkC,CAAC,KAAK,CAAC,CAAC;IACvD,MAAM,MAAM,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC;IAC1D,sBAAsB,CAAC,IAAI,EAAE,KAAK,CAAC,OAAO,EAAE,KAAK,CAAC,OAAO,EAAE,CAAC,EAAE,IAAI,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IACnF,OAAO,EAAE,IAAI,EAAE,MAAM,EAAE,CAAC;AAC1B,CAAC;AAED;;;;;;;;;GASG;AACH,MAAM,UAAU,iCAAiC,CAAC,QAAiE;IAEjH,IAAI,IAA2C,CAAC;IAChD,OAAO,CAAC,OAAO,EAAE,OAAO,EAAE,EAAE;QAC1B,IAAI,IAAI,KAAK,SAAS,EAAE,CAAC;YACvB,MAAM,KAAK,GAAG,8BAA8B,CAAC,EAAE,GAAG,QAAQ,EAAE,OAAO,EAAE,OAAO,EAAE,CAAC,CAAC;YAChF,IAAI,GAAG,KAAK,CAAC,IAAI,CAAC;YAClB,OAAO,KAAK,CAAC,MAAM,CAAC;QACtB,CAAC;QACD,MAAM,UAAU,GAAG,IAAI,CAAC,YAAY,GAAG,IAAI,CAAC,SAAS,CAAC;QACtD,IAAI,CAAC,CAAC,OAAO,YAAY,YAAY,IAAI,OAAO,YAAY,YAAY,CAAC,IAAI,OAAO,CAAC,MAAM,KAAK,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,MAAM;YACpH,CAAC,CAAC,OAAO,YAAY,YAAY,IAAI,OAAO,YAAY,YAAY,CAAC,IAAI,OAAO,CAAC,MAAM,KAAK,UAAU,GAAG,IAAI,CAAC,QAAQ,EAAE,CAAC;YAC3H,MAAM,IAAI,yBAAyB,CAAC,mEAAmE,CAAC,CAAC;QAC3G,CAAC;QACD,6BAA6B,CAAC,OAAO,CAAC,CAAC;QACvC,uBAAuB,CAAC,OAAO,EAAE,UAAU,EAAE,IAAI,CAAC,QAAQ,CAAC,CAAC;QAC5D,MAAM,MAAM,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC;QAC1D,sBAAsB,CAAC,IAAI,EAAE,OAAO,EAAE,OAAO,EAAE,CAAC,EAAE,IAAI,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;QACvE,OAAO,MAAM,CAAC;IAChB,CAAC,CAAC;AACJ,CAAC","sourcesContent":["// SPDX-License-Identifier: Apache-2.0\n\n/**\n * Blurs a sampled image-plane quantity with caller-supplied PSF (point-spread function) kernels, one\n * kernel per patch of the image. The quantity must be linear and additive (for example irradiance).\n * Each patch is treated as shift-invariant (\"isoplanatic\"): the same PSF applies everywhere inside it.\n * The PSFs are incoherent intensity PSFs.\n *\n * This is the CPU float64 reference for the optics-as-convolution step (#277). The `webgpu-f32`\n * path runs the same equation; see docs/WEBGPU.md.\n *\n * Equation (indexed by source site; an energy-conserving scatter written as a gather):\n *   out(x, y) = Σ_{v=-r..r} Σ_{u=-r..r} in(x − u, y − v) · K_{p(x − u, y − v)}(u, v)\n * where `p(s)` is the patch containing the **source** site and `K_p(u, v)` is the fraction of the source\n * site's quantity that lands at displacement (+u right, +v down), in grid sites.\n *\n * Incoherent imaging is a linear sum of intensity PSFs: the output intensity is the input intensity\n * convolved with the intensity PSF (MIT OpenCourseWare 2.71 Optics, Spring 2009, lecture 22, free).\n * A PSF that varies across the field is approximated as constant inside each patch, as in methods\n * that divide the image into isoplanatic subsections (van de Ketterij et al., arXiv:2011.01738,\n * sec. 1, free). Also see: Goodman, Introduction to Fourier Optics, 4th ed., 2017, §7.1. Camera\n * simulators such as ISETCam use the same piecewise shift-invariant method\n * (Farrell et al., Proc. SPIE 8299, 2012). The caller chooses the patch size and kernels, and so\n * owns the patch approximation.\n *\n * Edges: the grid holds every source. Sources outside the grid are not represented. Anything\n * scattered outside the grid is lost; compare `inputSum` with `outputSum` to see how much. There is no\n * padding, reflection or renormalization.\n */\nimport { InvalidConfigurationError } from \"../core/configuration-error.js\";\nimport { requireAllowlistedRecord } from \"../core/record-validation.js\";\nimport { FULL_RESOLUTION_CEILING_PIXELS, resolveRaisedLimit } from \"../core/pixel-limits.js\";\n\nexport const IRRADIANCE_CONVOLUTION_SCHEMA_VERSION = \"0.1.0\" as const;\n/** `maximumPixels` is the default grid limit; a request may raise it with its own `maximumPixels`, up to `ceilingPixels`. */\nexport const IRRADIANCE_CONVOLUTION_LIMITS = Object.freeze({\n  maximumPixels: 24_000_000,\n  ceilingPixels: FULL_RESOLUTION_CEILING_PIXELS,\n  maximumDimension: 16_384,\n  maximumKernelRadius: 32,\n  maximumPatchCount: 4_096,\n  maximumKernelWeights: 4_194_304,\n  /** A kernel may lose energy (vignetting, transmission) but may not create it, beyond rounding slack. */\n  maximumKernelSum: 1 + 2 ** -20,\n  /** Inputs larger than this are rejected so that float32 execution cannot overflow. */\n  maximumSampleMagnitude: 2 ** 100\n});\n\nexport interface IrradianceConvolutionInput {\n  schemaVersion: typeof IRRADIANCE_CONVOLUTION_SCHEMA_VERSION;\n  width: number;\n  height: number;\n  /** Row-major samples. Finite, non-negative, linear and additive (for example irradiance, or expected photons per site). */\n  samples: Float32Array | Float64Array;\n  /** Kernel half-width in grid sites. Each kernel is (2r+1)×(2r+1). */\n  kernelRadius: number;\n  /** Patch size in grid sites. The last row and column of patches may be partial. */\n  patchWidth: number;\n  patchHeight: number;\n  /**\n   * Finite non-negative weights, `[patchRow][patchColumn][v + r][u + r]`, one kernel per patch in\n   * row-major patch order. Each kernel sums to at most 1 (+2^-20 rounding slack).\n   */\n  kernels: Float32Array | Float64Array;\n  boundary: \"zero-outside-grid\";\n  /**\n   * Optional opt-in grid limit for large images: an integer from 1 to `IRRADIANCE_CONVOLUTION_LIMITS.ceilingPixels`.\n   * Default `IRRADIANCE_CONVOLUTION_LIMITS.maximumPixels` (24 000 000). It does not change the result.\n   */\n  maximumPixels?: number;\n}\n\n/** Validated, immutable geometry of a convolution request. Sample arrays are borrowed, not copied. */\nexport interface IrradianceConvolutionPlan {\n  width: number;\n  height: number;\n  kernelRadius: number;\n  kernelSize: number;\n  tapCount: number;\n  patchWidth: number;\n  patchHeight: number;\n  patchColumns: number;\n  patchRows: number;\n  /** Largest kernel sum and largest input magnitude. The float32 error bound uses them. */\n  maximumKernelSum: number;\n  maximumSampleMagnitude: number;\n}\n\nfunction positiveInteger(value: unknown, maximum: number): number {\n  if (typeof value !== \"number\" || !Number.isSafeInteger(value) || value <= 0 || value > maximum) {\n    throw new InvalidConfigurationError(\"Invalid irradiance convolution geometry.\");\n  }\n  return value;\n}\n\n/** Checks every sample (finite, non-negative, within the declared magnitude) and returns the largest. */\nfunction checkedMaximumSampleMagnitude(samples: Float32Array | Float64Array): number {\n  let maximumSampleMagnitude = 0;\n  for (let i = 0; i < samples.length; i++) {\n    const v = samples[i]!;\n    // Negated comparison also rejects NaN.\n    if (!(v >= 0 && v <= IRRADIANCE_CONVOLUTION_LIMITS.maximumSampleMagnitude)) {\n      throw new InvalidConfigurationError(\"Irradiance samples must be finite, non-negative and within the declared magnitude.\");\n    }\n    if (v > maximumSampleMagnitude) maximumSampleMagnitude = v;\n  }\n  return maximumSampleMagnitude;\n}\n\n/** Checks every kernel weight (finite, in [0, 1]) and kernel sum (at most 1 plus slack) and returns the largest sum. */\nfunction checkedMaximumKernelSum(kernels: Float32Array | Float64Array, patchCount: number, tapCount: number): number {\n  let maximumKernelSum = 0;\n  for (let p = 0; p < patchCount; p++) {\n    let sum = 0;\n    for (let t = p * tapCount, end = t + tapCount; t < end; t++) {\n      const w = kernels[t]!;\n      if (!(w >= 0 && w <= 1)) throw new InvalidConfigurationError(\"PSF kernel weights must be finite and within [0, 1].\");\n      sum += w;\n    }\n    if (sum > IRRADIANCE_CONVOLUTION_LIMITS.maximumKernelSum) throw new InvalidConfigurationError(\"A PSF kernel may not create energy (sum exceeds 1).\");\n    if (sum > maximumKernelSum) maximumKernelSum = sum;\n  }\n  return maximumKernelSum;\n}\n\n/** Validates the request and returns its plan. Every sample and weight is checked once. */\nexport function calculateIrradianceConvolutionPlan(input: IrradianceConvolutionInput): IrradianceConvolutionPlan {\n  const r = requireAllowlistedRecord(input, [\"schemaVersion\", \"width\", \"height\", \"samples\", \"kernelRadius\", \"patchWidth\",\n    \"patchHeight\", \"kernels\", \"boundary\", \"maximumPixels\"], \"Invalid irradiance convolution input.\");\n  const L = IRRADIANCE_CONVOLUTION_LIMITS;\n  if (r.schemaVersion !== IRRADIANCE_CONVOLUTION_SCHEMA_VERSION || r.boundary !== \"zero-outside-grid\") {\n    throw new InvalidConfigurationError(\"Unsupported irradiance convolution schema or boundary.\");\n  }\n  const width = positiveInteger(r.width, L.maximumDimension), height = positiveInteger(r.height, L.maximumDimension);\n  const maximumPixels = resolveRaisedLimit(r.maximumPixels, L.maximumPixels, L.ceilingPixels,\n    \"Irradiance convolution maximumPixels must be an integer from 1 to \" + L.ceilingPixels + \".\");\n  if (width * height > maximumPixels) throw new InvalidConfigurationError(\"Irradiance convolution grid exceeds the pixel limit.\");\n  const radius = r.kernelRadius;\n  if (typeof radius !== \"number\" || !Number.isSafeInteger(radius) || radius < 0 || radius > L.maximumKernelRadius) {\n    throw new InvalidConfigurationError(\"Invalid irradiance convolution kernel radius.\");\n  }\n  const patchWidth = positiveInteger(r.patchWidth, width), patchHeight = positiveInteger(r.patchHeight, height);\n  const patchColumns = Math.ceil(width / patchWidth), patchRows = Math.ceil(height / patchHeight);\n  const kernelSize = 2 * radius + 1, tapCount = kernelSize * kernelSize, patchCount = patchColumns * patchRows;\n  if (patchCount > L.maximumPatchCount || patchCount * tapCount > L.maximumKernelWeights) {\n    throw new InvalidConfigurationError(\"Irradiance convolution kernels exceed the declared limits.\");\n  }\n  const samples = r.samples, kernels = r.kernels;\n  if (!(samples instanceof Float32Array || samples instanceof Float64Array) || samples.length !== width * height ||\n      !(kernels instanceof Float32Array || kernels instanceof Float64Array) || kernels.length !== patchCount * tapCount) {\n    throw new InvalidConfigurationError(\"Irradiance convolution arrays do not match the declared geometry.\");\n  }\n  const maximumSampleMagnitude = checkedMaximumSampleMagnitude(samples);\n  const maximumKernelSum = checkedMaximumKernelSum(kernels, patchCount, tapCount);\n  return Object.freeze({ width, height, kernelRadius: radius, kernelSize, tapCount, patchWidth, patchHeight, patchColumns, patchRows,\n    maximumKernelSum, maximumSampleMagnitude });\n}\n\n/**\n * Reference gather for output rows `[y0, y1)`, written into `output` at `(y - y0) * width + x`.\n * The summation order (v outer, u inner, both ascending) is part of the contract shared with the\n * WGSL (WebGPU shader) kernel.\n * @internal\n */\nexport function convolveIrradianceRows(plan: IrradianceConvolutionPlan, samples: Float32Array | Float64Array,\n  kernels: Float32Array | Float64Array, y0: number, y1: number, output: Float64Array): void {\n  const { width, height, kernelRadius: r, kernelSize, tapCount, patchWidth, patchHeight, patchColumns, patchRows } = plan;\n  for (let y = y0; y < y1; y++) {\n    for (let x = 0; x < width; x++) {\n      let sum = 0;\n      for (let v = -r; v <= r; v++) {\n        const sy = y - v;\n        if (sy < 0 || sy >= height) continue;\n        const patchRowBase = Math.min(Math.floor(sy / patchHeight), patchRows - 1) * patchColumns;\n        const kernelRow = (v + r) * kernelSize + r, sampleRow = sy * width;\n        for (let u = -r; u <= r; u++) {\n          const sx = x - u;\n          if (sx < 0 || sx >= width) continue;\n          const patch = patchRowBase + Math.min(Math.floor(sx / patchWidth), patchColumns - 1);\n          sum += samples[sampleRow + sx]! * kernels[patch * tapCount + kernelRow + u]!;\n        }\n      }\n      output[(y - y0) * width + x] = sum;\n    }\n  }\n}\n\n/** Synchronous float64 reference for bounded grids. Use the image-kernel accelerator for large grids. */\nexport function calculateIrradianceConvolution(input: IrradianceConvolutionInput): { plan: IrradianceConvolutionPlan; values: Float64Array } {\n  const plan = calculateIrradianceConvolutionPlan(input);\n  const values = new Float64Array(plan.width * plan.height);\n  convolveIrradianceRows(plan, input.samples, input.kernels, 0, plan.height, values);\n  return { plan, values };\n}\n\n/**\n * Internal (#453; not exported from the package): runs several convolutions that share one geometry (width, height,\n * kernel radius, patch size, boundary, pixel limit) with one validation of that geometry. The first call is\n * `calculateIrradianceConvolution` itself, with every check. Each later call checks its own arrays exactly as the plan\n * does (array types and lengths against the geometry, every sample, every kernel weight and kernel sum, with the same\n * messages) and then runs the same gather (`convolveIrradianceRows`) with the first call's plan, whose geometry fields are\n * the same numbers a full validation of this call would give. The values are therefore identical to one\n * `calculateIrradianceConvolution` per call; only the repeated record and geometry checks are skipped.\n * @internal\n */\nexport function createIrradianceConvolutionRunner(geometry: Omit<IrradianceConvolutionInput, \"samples\" | \"kernels\">):\n  (samples: Float32Array | Float64Array, kernels: Float32Array | Float64Array) => Float64Array {\n  let plan: IrradianceConvolutionPlan | undefined;\n  return (samples, kernels) => {\n    if (plan === undefined) {\n      const first = calculateIrradianceConvolution({ ...geometry, samples, kernels });\n      plan = first.plan;\n      return first.values;\n    }\n    const patchCount = plan.patchColumns * plan.patchRows;\n    if (!(samples instanceof Float32Array || samples instanceof Float64Array) || samples.length !== plan.width * plan.height ||\n        !(kernels instanceof Float32Array || kernels instanceof Float64Array) || kernels.length !== patchCount * plan.tapCount) {\n      throw new InvalidConfigurationError(\"Irradiance convolution arrays do not match the declared geometry.\");\n    }\n    checkedMaximumSampleMagnitude(samples);\n    checkedMaximumKernelSum(kernels, patchCount, plan.tapCount);\n    const values = new Float64Array(plan.width * plan.height);\n    convolveIrradianceRows(plan, samples, kernels, 0, plan.height, values);\n    return values;\n  };\n}\n"]}

package/dist/optics/scene-to-sensor-irradiance.d.ts

@@ -212,6 +212,19 @@ export declare function parseSceneToSensorIrradianceProfile(value: unknown): Sce
  * factor 1/(1 + m)² from it, so the two cannot drift apart.
  */
 export declare function thinLensWorkingFNumberFactor(focalLengthMm: number, objectDistanceM: number): number;
+/**
+ * Internal (#453; not exported from the package): the sensor-plane spectral irradiance that
+ * {@link calculateSceneRadianceToSensorIrradiance} returns, for repeated queries of one optical state. This call runs the
+ * request-independent checks once (the public function's own code, in its order); each call of the returned function then
+ * runs the public function's per-query steps on an engine-built request and a result already parsed by the caller: the identity
+ * match, the wavelength-basis check, the transmission or T-stop path, the front-of-lens filters, the field throughput and
+ * the finite, non-negative irradiance, with the same operations, order and errors. It skips only the re-validation of the
+ * unchanged optical inputs, the re-parse of the request and result, and the result envelope.
+ */
+export declare function prepareSceneRadianceToSensorIrradianceValue(optics: Pick<CalculateSceneRadianceToSensorIrradianceInput, "profile" | "focalLengthMm" | "nominalFNumber" | "focus" | "frontOfLensFilters">): (request: SceneRadianceEvaluationRequest, result: SceneRadianceEvaluationResult, imagePointMm: {
+    x: number;
+    y: number;
+}, fieldThroughput: OpticalBridgeFieldThroughput) => number;
 /**
  * Applies the declared pupil acceptance and exactly one transmission or field-throughput path to
  * spectral scene radiance. Returns the irradiance before the sensor stack, with conservative

package/dist/optics/scene-to-sensor-irradiance.js

@@ -484,30 +484,7 @@ function validateSceneIdentity(request, result) {
         throw new InvalidScientificInputError("Scene-radiance request/result identity must match exactly.");
     }
 }
-/**
- * Applies the declared pupil acceptance and exactly one transmission or field-throughput path to
- * spectral scene radiance. Returns the irradiance before the sensor stack, with conservative
- * evidence.
- *
- * The optical bridge turns outgoing scene radiance (W/m²/sr/nm) into irradiance (W/m²/nm) arriving
- * at the front of the sensor stack.
- * - Circular paraxial acceptance is pi/(4*Nworking²). At finite focus with unit pupil magnification,
- *   the working f-number is N*(1+m).
- * - Spectral transmission and an effective working T-stop are alternative paths; use one.
- * - Field throughput is applied exactly once.
- * - Sensor response and PSF (point-spread function) blur happen later in the pipeline.
- * @param input - CalculateSceneRadianceToSensorIrradianceInput. See that type for coordinates, units and profile binding.
- * @returns CalculationResult<
-  SceneToSensorIrradianceResult
->. 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 calculateSceneRadianceToSensorIrradiance(input) {
-    requireInputKeys(input, "calculateSceneRadianceToSensorIrradiance", ["sceneRadianceRequest", "sceneRadianceResult", "profile", "focalLengthMm", "nominalFNumber", "focus", "imagePointMm", "fieldThroughput", "frontOfLensFilters"]);
-    const request = parseSceneRadianceEvaluationRequest(input.sceneRadianceRequest);
-    const sceneResult = parseSceneRadianceEvaluationResult(input.sceneRadianceResult);
-    validateSceneIdentity(request, sceneResult);
+function prepareSceneRadianceToSensorIrradiance(input) {
     const profile = parseSceneToSensorIrradianceProfile(input.profile);
     const focalLengthMm = requirePositiveFiniteScientific(input.focalLengthMm, "focalLengthMm");
     const nominalFNumber = requirePositiveFiniteScientific(input.nominalFNumber, "nominalFNumber");
@@ -535,6 +512,15 @@ export function calculateSceneRadianceToSensorIrradiance(input) {
         throw new InvalidScientificInputError("Paraxial geometric acceptance must remain finite and greater than zero.");
     }
     const transmission = profile.transmission;
+    return { profile, focalLengthMm, nominalFNumber, working, workingFNumberValue, geometricAcceptance, transmission };
+}
+/**
+ * The per-query part of {@link calculateSceneRadianceToSensorIrradiance}, after the request and result are parsed and their
+ * identities matched: the wavelength-basis check, the transmission or T-stop path, the front-of-lens filters, the field
+ * throughput and the irradiance, with the same operations in the same order.
+ */
+function applySceneRadianceToSensorIrradiance(prepared, input, request, sceneResult) {
+    const { workingFNumberValue, geometricAcceptance, transmission } = prepared;
     if (request.wavelengthBasis !==
         transmission.wavelengthBasis) {
         throw new InvalidScientificInputError("Scene-radiance and optical-transmission wavelength bases must match; implicit air/vacuum conversion is not performed.");
@@ -616,6 +602,54 @@ export function calculateSceneRadianceToSensorIrradiance(input) {
         irradiance < 0) {
         throw new InvalidScientificInputError("Sensor-plane spectral irradiance must remain finite and non-negative.");
     }
+    return { spectralTransmissionFactor, effectiveWorkingTStop, tStopEquivalentTransmissionFactor, effectiveAcceptance, transmissionEvidence,
+        frontOfLensFilterCount, frontOfLensFilterPolarizationModeled, frontOfLensFilterTransmissionFactor, frontOfLensFilterAttenuationStops,
+        frontOfLensFilterEvidence, frontOfLensFilterLimitations, fieldThroughputFactor, sceneRadiance, irradiance };
+}
+/**
+ * Internal (#453; not exported from the package): the sensor-plane spectral irradiance that
+ * {@link calculateSceneRadianceToSensorIrradiance} returns, for repeated queries of one optical state. This call runs the
+ * request-independent checks once (the public function's own code, in its order); each call of the returned function then
+ * runs the public function's per-query steps on an engine-built request and a result already parsed by the caller: the identity
+ * match, the wavelength-basis check, the transmission or T-stop path, the front-of-lens filters, the field throughput and
+ * the finite, non-negative irradiance, with the same operations, order and errors. It skips only the re-validation of the
+ * unchanged optical inputs, the re-parse of the request and result, and the result envelope.
+ */
+export function prepareSceneRadianceToSensorIrradianceValue(optics) {
+    const prepared = prepareSceneRadianceToSensorIrradiance(optics);
+    const frontOfLensFilters = optics.frontOfLensFilters;
+    return (request, result, imagePointMm, fieldThroughput) => {
+        validateSceneIdentity(request, result);
+        return applySceneRadianceToSensorIrradiance(prepared, { imagePointMm, fieldThroughput, frontOfLensFilters }, request, result).irradiance;
+    };
+}
+/**
+ * Applies the declared pupil acceptance and exactly one transmission or field-throughput path to
+ * spectral scene radiance. Returns the irradiance before the sensor stack, with conservative
+ * evidence.
+ *
+ * The optical bridge turns outgoing scene radiance (W/m²/sr/nm) into irradiance (W/m²/nm) arriving
+ * at the front of the sensor stack.
+ * - Circular paraxial acceptance is pi/(4*Nworking²). At finite focus with unit pupil magnification,
+ *   the working f-number is N*(1+m).
+ * - Spectral transmission and an effective working T-stop are alternative paths; use one.
+ * - Field throughput is applied exactly once.
+ * - Sensor response and PSF (point-spread function) blur happen later in the pipeline.
+ * @param input - CalculateSceneRadianceToSensorIrradianceInput. See that type for coordinates, units and profile binding.
+ * @returns CalculationResult<
+  SceneToSensorIrradianceResult
+>. 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 calculateSceneRadianceToSensorIrradiance(input) {
+    requireInputKeys(input, "calculateSceneRadianceToSensorIrradiance", ["sceneRadianceRequest", "sceneRadianceResult", "profile", "focalLengthMm", "nominalFNumber", "focus", "imagePointMm", "fieldThroughput", "frontOfLensFilters"]);
+    const request = parseSceneRadianceEvaluationRequest(input.sceneRadianceRequest);
+    const sceneResult = parseSceneRadianceEvaluationResult(input.sceneRadianceResult);
+    validateSceneIdentity(request, sceneResult);
+    const prepared = prepareSceneRadianceToSensorIrradiance(input);
+    const { profile, focalLengthMm, nominalFNumber, working, workingFNumberValue, geometricAcceptance, transmission } = prepared;
+    const { spectralTransmissionFactor, effectiveWorkingTStop, tStopEquivalentTransmissionFactor, effectiveAcceptance, transmissionEvidence, frontOfLensFilterCount, frontOfLensFilterPolarizationModeled, frontOfLensFilterTransmissionFactor, frontOfLensFilterAttenuationStops, frontOfLensFilterEvidence, frontOfLensFilterLimitations, fieldThroughputFactor, sceneRadiance, irradiance } = applySceneRadianceToSensorIrradiance(prepared, input, request, sceneResult);
     const limitations = [
         ...profile.limitations,
         ...frontOfLensFilterLimitations,

package/dist/optics/scene-to-sensor-irradiance.js.map

@@ -1 +1 @@
-{"version":3,"file":"scene-to-sensor-irradiance.js","sourceRoot":"","sources":["../../src/optics/scene-to-sensor-irradiance.ts"],"names":[],"mappings":"AAAA,sCAAsC;AAEtC;;;;;;;;;;;;;;;;;;;;GAoBG;AAEH,OAAO,EACL,mBAAmB,EAEpB,MAAM,+BAA+B,CAAC;AACvC,OAAO,EAAE,yBAAyB,EAAE,MAAM,gCAAgC,CAAC;AAC3E,OAAO,EACL,iBAAiB,EAGlB,MAAM,gCAAgC,CAAC;AACxC,OAAO,EACL,4BAA4B,EAE7B,MAAM,qBAAqB,CAAC;AAC7B,OAAO,EAAE,2BAA2B,EAAE,MAAM,uBAAuB,CAAC;AACpE,OAAO,EACL,oCAAoC,EAErC,MAAM,2BAA2B,CAAC;AAEnC,OAAO,EAAE,8BAA8B,EAAE,MAAM,gBAAgB,CAAC;AAChE,OAAO,EACL,mCAAmC,EACnC,kCAAkC,EAGnC,MAAM,6BAA6B,CAAC;AACrC,OAAO,EAAE,gBAAgB,EAAE,mBAAmB,EAAE,0BAA0B,EAAE,0BAA0B,EAAE,MAAM,8BAA8B,CAAC;AAI7I,MAAM,CAAC,MAAM,iDAAiD,GAC5D,OAAgB,CAAC;AAoQnB,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,GAAG,8BAA8B,CACtC,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,eAAe,CACtB,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,iCAAiC,CACpC,CAAC;IACJ,CAAC;IACD,OAAO,MAAM,CAAC;AAChB,CAAC;AAED,SAAS,UAAU,CACjB,KAAc,EACd,IAAY;IAEZ,MAAM,MAAM,GACV,aAAa,CAAC,KAAK,EAAE,IAAI,EAAE,CAAC,SAAS,EAAE,SAAS,CAAC,CAAC,CAAC;IACrD,MAAM,OAAO,GACX,qBAAqB,CACnB,MAAM,CAAC,OAAO,EACd,IAAI,GAAG,UAAU,CAClB,CAAC;IACJ,MAAM,OAAO,GACX,qBAAqB,CACnB,MAAM,CAAC,OAAO,EACd,IAAI,GAAG,UAAU,CAClB,CAAC;IACJ,IAAI,OAAO,GAAG,OAAO,EAAE,CAAC;QACtB,MAAM,IAAI,yBAAyB,CACjC,IAAI;YACF,iDAAiD,CACpD,CAAC;IACJ,CAAC;IACD,OAAO;QACL,OAAO;QACP,OAAO;KACR,CAAC;AACJ,CAAC;AAED,SAAS,qBAAqB,CAC5B,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;IAE7E,IAAI,MAAM,CAAC,IAAI,KAAK,UAAU,EAAE,CAAC;QAC/B,OAAO;YACL,IAAI,EAAE,UAAU;YAChB,QAAQ,EACN,eAAe,CACb,MAAM,CAAC,QAAQ,EACf,IAAI,GAAG,WAAW,CACnB;YACH,KAAK,EACH,qBAAqB,CACnB,MAAM,CAAC,KAAK,EACZ,IAAI,GAAG,QAAQ,CAChB;SACJ,CAAC;IACJ,CAAC;IAED,IACE,MAAM,CAAC,IAAI;QACX,gBAAgB,EAChB,CAAC;QACD,OAAO;YACL,IAAI,EAAE,gBAAgB;YACtB,UAAU,EACR,qBAAqB,CACnB,MAAM,CAAC,UAAU,EACjB,IAAI,GAAG,aAAa,CACrB;SACJ,CAAC;IACJ,CAAC;IAED,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,GAAG,KAAK,CAAC,GAAG,CACtB,CAAC,KAAK,EAAE,KAAK,EAAE,EAAE,CACf,qBAAqB,CACnB,KAAK,EACL,IAAI;QACF,GAAG;QACH,KAAK;QACL,GAAG,CACN,CACJ,CAAC;IACF,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,4BAA4B,CACnC,KAAc,EACd,IAAY;IAEZ,MAAM,QAAQ,GACZ,iBAAiB,CAAC,KAAK,EAAE,IAAI,CAAC,CAAC;IAEjC,IACE,QAAQ,CAAC,IAAI,CACX,CAAC,KAAK,EAAE,EAAE,CACR,KAAK,CAAC,WAAW;QACjB,wBAAwB,CAC3B,EACD,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,IAAI;YACF,0EAA0E,CAC7E,CAAC;IACJ,CAAC;IACD,OAAO,QAAQ,CAAC;AAClB,CAAC;AAED,SAAS,oBAAoB,CAC3B,KAAc,EACd,IAAY;IAEZ,IACE,CAAC,KAAK,CAAC,OAAO,CAAC,KAAK,CAAC;QACrB,KAAK,CAAC,MAAM,GAAG,CAAC,EAChB,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,IAAI;YACF,qCAAqC,CACxC,CAAC;IACJ,CAAC;IAED,IAAI,QAAQ,GACV,MAAM,CAAC,iBAAiB,CAAC;IAC3B,OAAO,KAAK,CAAC,GAAG,CACd,CAAC,KAAK,EAAE,KAAK,EAAE,EAAE;QACf,MAAM,UAAU,GACd,IAAI;YACJ,GAAG;YACH,KAAK;YACL,GAAG,CAAC;QACN,MAAM,MAAM,GACV,aAAa,CACX,KAAK,EACL,UAAU,EAAE,CAAC,sBAAsB,EAAE,0BAA0B,CAAC,CACjE,CAAC;QACJ,MAAM,oBAAoB,GACxB,qBAAqB,CACnB,MAAM,CAAC,oBAAoB,EAC3B,UAAU;YACR,uBAAuB,CAC1B,CAAC;QACJ,IACE,oBAAoB;YACpB,QAAQ,EACR,CAAC;YACD,MAAM,IAAI,yBAAyB,CACjC,IAAI;gBACF,2CAA2C,CAC9C,CAAC;QACJ,CAAC;QACD,QAAQ;YACN,oBAAoB,CAAC;QACvB,OAAO;YACL,oBAAoB;YACpB,wBAAwB,EACtB,eAAe,CACb,MAAM,CAAC,wBAAwB,EAC/B,UAAU;gBACR,2BAA2B,CAC9B;SACJ,CAAC;IACJ,CAAC,CACF,CAAC;AACJ,CAAC;AAED,SAAS,iBAAiB,CACxB,KAAc,EACd,IAAY;IAEZ,MAAM,MAAM,GACV,aAAa,CAAC,KAAK,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;IACnC,0BAA0B,CAAC,MAAM,EAAE,MAAM,EAAE;QACzC,iIAAiI;QACjI,uBAAuB,EAAE,CAAC,MAAM,EAAE,kBAAkB,EAAE,iBAAiB,EAAE,SAAS,EAAE,aAAa,EAAE,cAAc,EAAE,2BAA2B,CAAC;QAC/I,wCAAwC,EAAE,CAAC,MAAM,EAAE,kBAAkB,EAAE,cAAc,EAAE,iBAAiB,EAAE,2BAA2B;YACnI,aAAa,EAAE,YAAY,EAAE,SAAS,CAAC;KAC1C,EAAE,IAAI,CAAC,CAAC;IAET,IACE,MAAM,CAAC,IAAI;QACX,uBAAuB,EACvB,CAAC;QACD,IACE,MAAM,CAAC,YAAY;YACnB,SAAS;YACT,MAAM,CAAC,yBAAyB;gBAChC,SAAS,EACT,CAAC;YACD,MAAM,IAAI,yBAAyB,CACjC,IAAI;gBACF,6DAA6D,CAChE,CAAC;QACJ,CAAC;QAED,MAAM,MAAM,GACV,qBAAqB,CACnB,MAAM,CAAC,gBAAgB,EACvB,IAAI;YACF,mBAAmB,CACtB,CAAC;QACJ,MAAM,WAAW,GACf,aAAa,CACX,MAAM,CAAC,OAAO,EACd,IAAI,GAAG,UAAU,EAAE,CAAC,OAAO,EAAE,UAAU,CAAC,CACzC,CAAC;QACJ,MAAM,WAAW,GACf,gBAAgB,CACd,MAAM,CAAC,WAAW,EAClB,IAAI,GAAG,cAAc,CACtB,CAAC;QAEJ,IACE,MAAM,KAAK,YAAY;YACvB,WAAW,CAAC,IAAI;gBACd,UAAU,EACZ,CAAC;YACD,MAAM,IAAI,yBAAyB,CACjC,IAAI;gBACF,6EAA6E,CAChF,CAAC;QACJ,CAAC;QAED,OAAO;YACL,IAAI,EACF,uBAAuB;YACzB,gBAAgB,EAAE,MAAM;YACxB,eAAe,EACb,kBAAkB,CAChB,MAAM,CAAC,eAAe,EACtB,IAAI;gBACF,kBAAkB,CACrB;YACH,OAAO,EAAE;gBACP,KAAK,EACH,oBAAoB,CAClB,WAAW,CAAC,KAAK,EACjB,IAAI;oBACF,gBAAgB,CACnB;gBACH,QAAQ,EACN,4BAA4B,CAC1B,WAAW,CAAC,QAAQ,EACpB,IAAI;oBACF,mBAAmB,CACtB;aACJ;YACD,WAAW;SACZ,CAAC;IACJ,CAAC;IAED,IACE,MAAM,CAAC,IAAI;QACX,wCAAwC,EACxC,CAAC;QACD,IACE,MAAM,CAAC,OAAO;YACd,SAAS,EACT,CAAC;YACD,MAAM,IAAI,yBAAyB,CACjC,IAAI;gBACF,uEAAuE,CAC1E,CAAC;QACJ,CAAC;QACD,IACE,MAAM,CAAC,gBAAgB;YACvB,eAAe,EACf,CAAC;YACD,MAAM,IAAI,yBAAyB,CACjC,IAAI;gBACF,mDAAmD,CACtD,CAAC;QACJ,CAAC;QACD,MAAM,gBAAgB,GACpB,aAAa,CACX,MAAM,CAAC,YAAY,EACnB,IAAI;YACF,eAAe,EAAE,CAAC,OAAO,EAAE,UAAU,CAAC,CACzC,CAAC;QACJ,OAAO;YACL,IAAI,EACF,wCAAwC;YAC1C,gBAAgB,EACd,eAAe;YACjB,YAAY,EAAE;gBACZ,KAAK,EACH,qBAAqB,CACnB,gBAAgB,CAAC,KAAK,EACtB,IAAI;oBACF,qBAAqB,CACxB;gBACH,QAAQ,EACN,iBAAiB,CACf,gBAAgB,CAAC,QAAQ,EACzB,IAAI;oBACF,wBAAwB,CAC3B;aACJ;YACD,eAAe,EACb,kBAAkB,CAChB,MAAM,CAAC,eAAe,EACtB,IAAI;gBACF,kBAAkB,CACrB;YACH,yBAAyB,EACvB,UAAU,CACR,MAAM,CAAC,yBAAyB,EAChC,IAAI;gBACF,4BAA4B,CAC/B;YACH,WAAW,EACT,gBAAgB,CACd,MAAM,CAAC,WAAW,EAClB,IAAI,GAAG,cAAc,CACtB;YACH,UAAU,EACR,qBAAqB,CACnB,MAAM,CAAC,UAAU,EACjB,IAAI;gBACF,aAAa,CAChB;SACJ,CAAC;IACJ,CAAC;IAED,MAAM,IAAI,yBAAyB,CACjC,IAAI,GAAG,mBAAmB,CAC3B,CAAC;AACJ,CAAC;AAED,SAAS,kBAAkB,CACzB,KAAc,EACd,IAAY;IAKZ,MAAM,MAAM,GACV,4BAA4B,CAC1B,KAAK,EACL,IAAI,CACL,CAAC;IACJ,IAAI,MAAM,KAAK,aAAa,EAAE,CAAC;QAC7B,MAAM,IAAI,yBAAyB,CACjC,IAAI;YACF,qCAAqC,CACxC,CAAC;IACJ,CAAC;IACD,OAAO,MAAM,CAAC;AAChB,CAAC;AAED,SAAS,uBAAuB,CAC9B,KAAc,EACd,IAAY;IAEZ,MAAM,MAAM,GACV,aAAa,CAAC,KAAK,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;IACnC,0BAA0B,CAAC,MAAM,EAAE,MAAM,EAAE,EAAE,GAAG,EAAE,CAAC,MAAM,CAAC,EAAE,uBAAuB,EAAE,CAAC,MAAM,EAAE,iBAAiB,CAAC,EAAE,EAAE,IAAI,CAAC,CAAC;IAE1H,IAAI,MAAM,CAAC,IAAI,KAAK,KAAK,EAAE,CAAC;QAC1B,OAAO;YACL,IAAI,EAAE,KAAK;SACZ,CAAC;IACJ,CAAC;IACD,IACE,MAAM,CAAC,IAAI;QACX,uBAAuB,EACvB,CAAC;QACD,OAAO;YACL,IAAI,EACF,uBAAuB;YACzB,eAAe,EACb,UAAU,CACR,MAAM,CAAC,eAAe,EACtB,IAAI;gBACF,kBAAkB,CACrB;SACJ,CAAC;IACJ,CAAC;IACD,MAAM,IAAI,yBAAyB,CACjC,IAAI,GAAG,mBAAmB,CAC3B,CAAC;AACJ,CAAC;AAED;;;;;;;;;;;;;;;GAeG;AACH,MAAM,UAAU,mCAAmC,CACjD,KAAc;IAEd,MAAM,MAAM,GACV,aAAa,CACX,KAAK,EACL,gCAAgC,EAAE,CAAC,eAAe,EAAE,eAAe,EAAE,kBAAkB,EAAE,cAAc,EAAE,gCAAgC,EAAE,4BAA4B,EAAE,4BAA4B,EAAE,oBAAoB,EAAE,qBAAqB,EAAE,mCAAmC,EAAE,6BAA6B,EAAE,WAAW,EAAE,gBAAgB,EAAE,eAAe,EAAE,UAAU,EAAE,aAAa,CAAC,CACnY,CAAC;IACJ,IACE,MAAM,CAAC,aAAa;QACpB,iDAAiD,EACjD,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,wDAAwD;YACtD,iDAAiD;YACjD,IAAI,CACP,CAAC;IACJ,CAAC;IAED,MAAM,aAAa,GACjB,aAAa,CACX,MAAM,CAAC,aAAa,EACpB,8CAA8C,EAAE,CAAC,eAAe,EAAE,gBAAgB,EAAE,OAAO,CAAC,CAC7F,CAAC;IACJ,MAAM,MAAM,GACV,qBAAqB,CACnB,MAAM,CAAC,gBAAgB,EACvB,iDAAiD,CAClD,CAAC;IACJ,MAAM,YAAY,GAChB,iBAAiB,CACf,MAAM,CAAC,YAAY,EACnB,6CAA6C,CAC9C,CAAC;IAEJ,IACE,MAAM,KAAK,YAAY;QACvB,YAAY,CAAC,gBAAgB;YAC3B,YAAY,EACd,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,gFAAgF,CACjF,CAAC;IACJ,CAAC;IAED,IACE,MAAM,CAAC,8BAA8B;QACnC,uBAAuB;QACzB,MAAM,CAAC,0BAA0B;YAC/B,6CAA6C;QAC/C,MAAM,CAAC,0BAA0B;YAC/B,KAAK;QACP,MAAM,CAAC,kBAAkB;YACvB,KAAK;QACP,MAAM,CAAC,mBAAmB;YACxB,KAAK;QACP,MAAM,CAAC,iCAAiC;YACtC,KAAK;QACP,MAAM,CAAC,2BAA2B;YAChC,KAAK,EACP,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,2FAA2F,CAC5F,CAAC;IACJ,CAAC;IAED,OAAO;QACL,aAAa,EACX,iDAAiD;QACnD,SAAS,EACP,qBAAqB,CACnB,MAAM,CAAC,SAAS,EAChB,0CAA0C,CAC3C;QACH,cAAc,EACZ,qBAAqB,CACnB,MAAM,CAAC,cAAc,EACrB,+CAA+C,CAChD;QACH,aAAa,EACX,qBAAqB,CACnB,MAAM,CAAC,aAAa,EACpB,8CAA8C,CAC/C;QACH,gBAAgB,EAAE,MAAM;QACxB,aAAa,EAAE;YACb,aAAa,EACX,UAAU,CACR,aAAa,CAAC,aAAa,EAC3B,4DAA4D,CAC7D;YACH,cAAc,EACZ,UAAU,CACR,aAAa,CAAC,cAAc,EAC5B,6DAA6D,CAC9D;YACH,KAAK,EACH,uBAAuB,CACrB,aAAa,CAAC,KAAK,EACnB,oDAAoD,CACrD;SACJ;QACD,YAAY;QACZ,8BAA8B,EAC5B,uBAAuB;QACzB,0BAA0B,EACxB,6CAA6C;QAC/C,0BAA0B,EACxB,KAAK;QACP,kBAAkB,EAAE,KAAK;QACzB,mBAAmB,EAAE,KAAK;QAC1B,iCAAiC,EAC/B,KAAK;QACP,2BAA2B,EACzB,KAAK;QACP,QAAQ,EACN,iBAAiB,CACf,MAAM,CAAC,QAAQ,EACf,yCAAyC,CAC1C;QACH,WAAW,EACT,gBAAgB,CACd,MAAM,CAAC,WAAW,EAClB,4CAA4C,CAC7C;KACJ,CAAC;AACJ,CAAC;AAED,SAAS,YAAY,CACnB,KAAa,EACb,KAAmB,EACnB,IAAY;IAEZ,IACE,KAAK,GAAG,KAAK,CAAC,OAAO;QACrB,KAAK,GAAG,KAAK,CAAC,OAAO,EACrB,CAAC;QACD,MAAM,IAAI,2BAA2B,CACnC,IAAI;YACF,wDAAwD,CAC3D,CAAC;IACJ,CAAC;AACH,CAAC;AAED;;;;;GAKG;AACH,MAAM,UAAU,4BAA4B,CAAC,aAAqB,EAAE,eAAuB;IACzF,OAAO,CAAC,GAAG,8BAA8B,CAAC,EAAE,aAAa,EAAE,eAAe,EAAE,CAAC,CAAC,KAAK,CAAC,aAAa,CAAC;AACpG,CAAC;AAED,SAAS,cAAc,CACrB,KAC2B,EAC3B,aAAqB,EACrB,cAAsB;IAMtB,IACE,KAAK,CAAC,IAAI,KAAK,gBAAgB,EAC/B,CAAC;QACD,OAAO;YACL,KAAK,EAAE,cAAc;YACrB,QAAQ,EAAE,EAAE;SACb,CAAC;IACJ,CAAC;IAED,IACE,KAAK,CAAC,IAAI;QACV,2BAA2B,EAC3B,CAAC;QACD,IACE,KAAK;aACF,4BAA4B;YAC/B,OAAO,EACP,CAAC;YACD,MAAM,IAAI,2BAA2B,CACnC,gFAAgF,CACjF,CAAC;QACJ,CAAC;QACD,OAAO;YACL,KAAK,EACH,cAAc;gBACd,4BAA4B,CAAC,aAAa,EAAE,KAAK,CAAC,eAAe,CAAC;YACpE,QAAQ,EAAE,EAAE;SACb,CAAC;IACJ,CAAC;IAED,IACE,KAAK,CAAC,IAAI;QACV,2BAA2B,EAC3B,CAAC;QACD,MAAM,KAAK,GACT,+BAA+B,CAC7B,KAAK,CAAC,cAAc,CAAC,KAAK,EAC1B,4BAA4B,CAC7B,CAAC;QACJ,OAAO;YACL,KAAK;YACL,QAAQ,EACN,iBAAiB,CACf,KAAK,CAAC,cAAc,CAAC,QAAQ,EAC7B,+BAA+B,CAChC;SACJ,CAAC;IACJ,CAAC;IAED,MAAM,IAAI,yBAAyB,CACjC,wBAAwB,CACzB,CAAC;AACJ,CAAC;AAED,SAAS,+BAA+B,CACtC,KAAc,EACd,IAAY;IAEZ,IACE,OAAO,KAAK,KAAK,QAAQ;QACzB,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC;QACvB,KAAK,IAAI,CAAC,EACV,CAAC;QACD,MAAM,IAAI,2BAA2B,CACnC,IAAI;YACF,wCAAwC,CAC3C,CAAC;IACJ,CAAC;IACD,OAAO,KAAK,CAAC;AACf,CAAC;AAED,SAAS,mBAAmB,CAC1B,KAC2B;IAE3B,IACE,KAAK,CAAC,IAAI;QACV,2BAA2B,EAC3B,CAAC;QACD,OAAO,+BAA+B,CACpC,KAAK,CAAC,eAAe,EACrB,uBAAuB,CACxB,CAAC;IACJ,CAAC;IACD,IACE,KAAK,CAAC,IAAI;QACR,2BAA2B;QAC7B,KAAK,CAAC,KAAK,CAAC,IAAI,KAAK,QAAQ,EAC7B,CAAC;QACD,OAAO,+BAA+B,CACpC,KAAK,CAAC,KAAK,CAAC,eAAe,EAC3B,6BAA6B,CAC9B,CAAC;IACJ,CAAC;IACD,OAAO,IAAI,CAAC;AACd,CAAC;AAED,SAAS,uBAAuB,CAC9B,KAIG,EACH,oBAA4B;IAE5B,MAAM,OAAO,GACX,KAAK,CAAC,OAAO,CAAC,KAAK,CAAC;IACtB,MAAM,KAAK,GACT,OAAO,CAAC,CAAC,CAAE,CAAC;IACd,MAAM,IAAI,GACR,OAAO,CACL,OAAO,CAAC,MAAM,GAAG,CAAC,CAClB,CAAC;IAEL,IACE,oBAAoB;QAClB,KAAK,CAAC,oBAAoB;QAC5B,oBAAoB;YAClB,IAAI,CAAC,oBAAoB,EAC3B,CAAC;QACD,MAAM,IAAI,2BAA2B,CACnC,uEAAuE,CACxE,CAAC;IACJ,CAAC;IAED,IACE,oBAAoB;QACpB,KAAK,CAAC,oBAAoB,EAC1B,CAAC;QACD,OAAO,KAAK;aACT,wBAAwB,CAAC;IAC9B,CAAC;IAED,KACE,IAAI,KAAK,GAAG,CAAC,EACb,KAAK,GAAG,OAAO,CAAC,MAAM,EACtB,KAAK,IAAI,CAAC,EACV,CAAC;QACD,MAAM,KAAK,GACT,OAAO,CAAC,KAAK,CAAE,CAAC;QAClB,IACE,oBAAoB;YACpB,KAAK,CAAC,oBAAoB,EAC1B,CAAC;YACD,MAAM,IAAI,GACR,OAAO,CAAC,KAAK,GAAG,CAAC,CAAE,CAAC;YACtB,IACE,oBAAoB;gBACpB,KAAK,CAAC,oBAAoB,EAC1B,CAAC;gBACD,OAAO,KAAK;qBACT,wBAAwB,CAAC;YAC9B,CAAC;YACD,MAAM,KAAK,GACT,CAAC,oBAAoB;gBACnB,IAAI,CAAC,oBAAoB,CAAC;gBAC5B,CAAC,KAAK,CAAC,oBAAoB;oBACzB,IAAI,CAAC,oBAAoB,CAAC,CAAC;YAC/B,OAAO,CACL,IAAI;iBACD,wBAAwB;gBAC3B,KAAK;oBACH,CAAC,KAAK;yBACD,wBAAwB;wBAC3B,IAAI;6BACD,wBAAwB,CAAC,CACjC,CAAC;QACJ,CAAC;IACH,CAAC;IAED,OAAO,IAAI;SACR,wBAAwB,CAAC;AAC9B,CAAC;AAED,SAAS,WAAW,CAClB,KAC8B,EAC9B,YAGC;IAED,IAAI,KAAK,CAAC,IAAI,KAAK,OAAO,EAAE,CAAC;QAC3B,OAAO,CAAC,CAAC;IACX,CAAC;IACD,IACE,KAAK,CAAC,IAAI;QACV,gCAAgC,EAChC,CAAC;QACD,MAAM,MAAM,GACV,KAAK,CAAC,MAAM,CAAC;QACf,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CACpB,CAAC,EACD,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC,CAAC,EACxB,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC,CAAC,EACxB,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,YAAY,CAAC,CAAC,CAAC,EAC/B,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,YAAY,CAAC,CAAC,CAAC,CAChC,CAAC;QACF,MAAM,SAAS,GACb,MAAM,CAAC,OAAO;YACd,EAAE;YACF,KAAK,CAAC;QACR,IACE,IAAI,CAAC,GAAG,CACN,MAAM,CAAC,YAAY,CAAC,CAAC;YACnB,YAAY,CAAC,CAAC,CACjB,GAAG,SAAS;YACb,IAAI,CAAC,GAAG,CACN,MAAM,CAAC,YAAY,CAAC,CAAC;gBACnB,YAAY,CAAC,CAAC,CACjB,GAAG,SAAS,EACb,CAAC;YACD,MAAM,IAAI,2BAA2B,CACnC,yFAAyF,CAC1F,CAAC;QACJ,CAAC;QACD,MAAM,MAAM,GACV,MAAM,CAAC,sBAAsB,CAAC;QAChC,IACE,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC;YACxB,MAAM,IAAI,CAAC;YACX,MAAM,GAAG,CAAC,EACV,CAAC;YACD,MAAM,IAAI,2BAA2B,CACnC,uDAAuD,CACxD,CAAC;QACJ,CAAC;QACD,OAAO,MAAM,CAAC;IAChB,CAAC;IACD,MAAM,IAAI,yBAAyB,CACjC,kCAAkC,CACnC,CAAC;AACJ,CAAC;AAED,SAAS,qBAAqB,CAC5B,OACgC,EAChC,MAC+B;IAE/B,IACE,OAAO,CAAC,QAAQ;QACd,MAAM,CAAC,QAAQ;QACjB,OAAO,CAAC,iBAAiB;YACvB,MAAM,CAAC,iBAAiB;QAC1B,OAAO,CAAC,OAAO;YACb,MAAM,CAAC,OAAO;QAChB,OAAO,CAAC,oBAAoB;YAC1B,MAAM,CAAC,oBAAoB;QAC7B,OAAO,CAAC,eAAe;YACrB,MAAM,CAAC,eAAe,EACxB,CAAC;QACD,MAAM,IAAI,2BAA2B,CACnC,4DAA4D,CAC7D,CAAC;IACJ,CAAC;AACH,CAAC;AAED;;;;;;;;;;;;;;;;;;GAkBG;AACH,MAAM,UAAU,wCAAwC,CACtD,KAC+C;IAI/C,gBAAgB,CAAC,KAAK,EAAE,0CAA0C,EAAE,CAAC,sBAAsB,EAAE,qBAAqB,EAAE,SAAS,EAAE,eAAe,EAAE,gBAAgB,EAAE,OAAO,EAAE,cAAc,EAAE,iBAAiB,EAAE,oBAAoB,CAAC,CAAC,CAAC;IACrO,MAAM,OAAO,GACX,mCAAmC,CACjC,KAAK,CAAC,oBAAoB,CAC3B,CAAC;IACJ,MAAM,WAAW,GACf,kCAAkC,CAChC,KAAK,CAAC,mBAAmB,CAC1B,CAAC;IACJ,qBAAqB,CACnB,OAAO,EACP,WAAW,CACZ,CAAC;IAEF,MAAM,OAAO,GACX,mCAAmC,CACjC,KAAK,CAAC,OAAO,CACd,CAAC;IACJ,MAAM,aAAa,GACjB,+BAA+B,CAC7B,KAAK,CAAC,aAAa,EACnB,eAAe,CAChB,CAAC;IACJ,MAAM,cAAc,GAClB,+BAA+B,CAC7B,KAAK,CAAC,cAAc,EACpB,gBAAgB,CACjB,CAAC;IAEJ,YAAY,CACV,aAAa,EACb,OAAO,CAAC,aAAa;SAClB,aAAa,EAChB,eAAe,CAChB,CAAC;IACF,YAAY,CACV,cAAc,EACd,OAAO,CAAC,aAAa;SAClB,cAAc,EACjB,gBAAgB,CACjB,CAAC;IAEF,MAAM,cAAc,GAClB,mBAAmB,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC;IACnC,IACE,OAAO,CAAC,aAAa,CAAC,KAAK,CAAC,IAAI;QAC9B,uBAAuB,EACzB,CAAC;QACD,IAAI,cAAc,KAAK,IAAI,EAAE,CAAC;YAC5B,MAAM,IAAI,2BAA2B,CACnC,mDAAmD,CACpD,CAAC;QACJ,CAAC;QACD,YAAY,CACV,cAAc,EACd,OAAO,CAAC,aAAa;aAClB,KAAK,CAAC,eAAe,EACxB,uBAAuB,CACxB,CAAC;IACJ,CAAC;IAED,MAAM,OAAO,GACX,cAAc,CACZ,KAAK,CAAC,KAAK,EACX,aAAa,EACb,cAAc,CACf,CAAC;IACJ,MAAM,mBAAmB,GACvB,OAAO,CAAC,KAAK,CAAC;IAChB,MAAM,mBAAmB,GACvB,IAAI,CAAC,EAAE;QACP,CAAC,CAAC;YACA,mBAAmB;YACnB,mBAAmB,CAAC,CAAC;IAEzB,IACE,CAAC,MAAM,CAAC,QAAQ,CACd,mBAAmB,CACpB;QACD,mBAAmB,IAAI,CAAC,EACxB,CAAC;QACD,MAAM,IAAI,2BAA2B,CACnC,yEAAyE,CAC1E,CAAC;IACJ,CAAC;IAED,MAAM,YAAY,GAChB,OAAO,CAAC,YAAY,CAAC;IAEvB,IACE,OAAO,CAAC,eAAe;QACvB,YAAY,CAAC,eAAe,EAC5B,CAAC;QACD,MAAM,IAAI,2BAA2B,CACnC,uHAAuH,CACxH,CAAC;IACJ,CAAC;IAED,IAAI,0BAA0B,GACZ,IAAI,CAAC;IACvB,IAAI,qBAAqB,GACP,IAAI,CAAC;IACvB,IAAI,iCAAiC,GACnB,IAAI,CAAC;IACvB,IAAI,mBAAmB,GACrB,mBAAmB,CAAC;IACtB,IAAI,oBAC2B,CAAC;IAEhC,IACE,YAAY,CAAC,IAAI;QACjB,uBAAuB,EACvB,CAAC;QACD,0BAA0B;YACxB,uBAAuB,CACrB,YAAY,EACZ,OAAO,CAAC,oBAAoB,CAC7B,CAAC;QACJ,mBAAmB;YACjB,0BAA0B,CAAC;QAC7B,oBAAoB;YAClB,YAAY,CAAC,OAAO,CAAC,QAAQ,CAAC;IAClC,CAAC;SAAM,CAAC;QACN,YAAY,CACV,OAAO,CAAC,oBAAoB,EAC5B,YAAY;aACT,yBAAyB,EAC5B,kBAAkB,CACnB,CAAC;QACF,qBAAqB;YACnB,YAAY,CAAC,YAAY,CAAC,KAAK,CAAC;QAClC,IACE,qBAAqB;YACrB,mBAAmB,EACnB,CAAC;YACD,MAAM,IAAI,2BAA2B,CACnC,uHAAuH,CACxH,CAAC;QACJ,CAAC;QACD,iCAAiC;YAC/B,CAAC,mBAAmB;gBAClB,qBAAqB,CAAC;gBACxB,CAAC,CAAC;QACJ,mBAAmB;YACjB,IAAI,CAAC,EAAE;gBACP,CAAC,CAAC;oBACA,qBAAqB;oBACrB,qBAAqB,CAAC,CAAC;QAC3B,oBAAoB;YAClB,YAAY;iBACT,YAAY,CAAC,QAAQ,CAAC;IAC7B,CAAC;IAED,IAAI,sBAAsB,GAAG,CAAC,CAAC;IAC/B,IAAI,oCAAoC,GAAG,KAAK,CAAC;IACjD,IAAI,mCAAmC,GAAG,CAAC,CAAC;IAC5C,IAAI,iCAAiC,GAAG,CAAC,CAAC;IAC1C,IAAI,yBAAyB,GACK,EAAE,CAAC;IACrC,IAAI,4BAA4B,GACV,EAAE,CAAC;IAEzB,IACE,KAAK,CAAC,kBAAkB;QACxB,SAAS,EACT,CAAC;QACD,MAAM,QAAQ,GACZ,oCAAoC,CAAC;YACnC,OAAO,EACL,KAAK,CAAC,kBAAkB;YAC1B,oBAAoB,EAClB,OAAO,CAAC,oBAAoB;YAC9B,eAAe,EACb,OAAO,CAAC,eAAe;SAC1B,CAAC,CAAC;QACL,sBAAsB;YACpB,QAAQ,CAAC,kBAAkB,CAAC;QAC9B,oCAAoC;YAClC,QAAQ,CAAC,mBAAmB,CAAC;QAC/B,mCAAmC;YACjC,QAAQ;iBACL,gCAAgC,CAAC;QACtC,iCAAiC;YAC/B,QAAQ;iBACL,wBAAwB,CAAC;QAC9B,yBAAyB;YACvB,QAAQ,CAAC,UAAU,CAAC,OAAO,CACzB,CAAC,SAAS,EAAE,EAAE,CACZ,SAAS,CAAC,QAAQ,CACrB,CAAC;QACJ,4BAA4B;YAC1B,QAAQ,CAAC,UAAU,CAAC,OAAO,CACzB,CAAC,SAAS,EAAE,EAAE,CACZ,SAAS,CAAC,WAAW,CACxB,CAAC;QACJ,mBAAmB;YACjB,mCAAmC,CAAC;IACxC,CAAC;IAED,MAAM,qBAAqB,GACzB,WAAW,CACT,KAAK,CAAC,eAAe,EACrB,KAAK,CAAC,YAAY,CACnB,CAAC;IACJ,mBAAmB;QACjB,qBAAqB,CAAC;IAExB,MAAM,aAAa,GACjB,WAAW;SACR,qDAAqD,CAAC;IAC3D,MAAM,UAAU,GACd,aAAa;QACb,mBAAmB,CAAC;IAEtB,IACE,CAAC,MAAM,CAAC,QAAQ,CAAC,UAAU,CAAC;QAC5B,UAAU,GAAG,CAAC,EACd,CAAC;QACD,MAAM,IAAI,2BAA2B,CACnC,uEAAuE,CACxE,CAAC;IACJ,CAAC;IAED,MAAM,WAAW,GAAG;QAClB,GAAG,OAAO,CAAC,WAAW;QACtB,GAAG,4BAA4B;QAC/B,4IAA4I;QAC5I,wGAAwG;QACxG,mEAAmE;QACnE,sFAAsF;QACtF,wIAAwI;QACxI,6FAA6F;QAC7F,oCAAoC;YAClC,CAAC,CAAC,0OAA0O;YAC5O,CAAC,CAAC,uHAAuH;KAC5H,CAAC;IAEF,OAAO,mBAAmB,CACxB;QACE,aAAa,EACX,iDAAiD;QACnD,SAAS,EAAE,OAAO,CAAC,SAAS;QAC5B,cAAc,EACZ,OAAO,CAAC,cAAc;QACxB,aAAa,EACX,OAAO,CAAC,aAAa;QACvB,QAAQ,EAAE,OAAO,CAAC,QAAQ;QAC1B,iBAAiB,EACf,OAAO,CAAC,iBAAiB;QAC3B,OAAO,EAAE,OAAO,CAAC,OAAO;QACxB,4BAA4B,EAC1B,OAAO;aACJ,4BAA4B;QACjC,oBAAoB,EAClB,OAAO,CAAC,oBAAoB;QAC9B,eAAe,EACb,OAAO,CAAC,eAAe;QACzB,aAAa,EACX,4BAA4B;QAC9B,SAAS,EACP,aAAa;QACf,cAAc,EACZ,kCAAkC;QACpC,UAAU,EAAE,UAAU;QACtB,0DAA0D,EACxD,aAAa;QACf,aAAa;QACb,cAAc;QACd,cAAc,EACZ,mBAAmB;QACrB,UAAU,EAAE,KAAK,CAAC,KAAK,CAAC,IAAI;QAC5B,+CAA+C,EAC7C,mBAAmB;QACrB,gBAAgB,EACd,YAAY,CAAC,IAAI;QACnB,0BAA0B;QAC1B,qBAAqB;QACrB,iCAAiC;QACjC,sBAAsB;QACtB,mCAAmC;QACnC,iCAAiC;QACjC,+CAA+C,EAC7C,IAAI;QACN,oCAAoC;QACpC,qBAAqB;QACrB,sDAAsD,EACpD,mBAAmB;QACrB,yDAAyD,EACvD,UAAU;QACZ,gBAAgB,EACd,eAAe;QACjB,6BAA6B,EAC3B,WAAW,CAAC,gBAAgB;QAC9B,8BAA8B,EAC5B,OAAO,CAAC,gBAAgB;QAC1B,wBAAwB,EACtB,WAAW,CAAC,WAAW;QACzB,kBAAkB,EAChB,YAAY,CAAC,WAAW;QAC1B,4BAA4B,EAC1B,KAAK;QACP,iCAAiC,EAC/B,IAAI;QACN,2BAA2B,EACzB,KAAK;QACP,+BAA+B,EAC7B,KAAK;QACP,wBAAwB,EACtB,KAAK;QACP,yBAAyB,EACvB,KAAK;QACP,UAAU,EAAE,KAAK;QACjB,wBAAwB,EAAE,KAAK;QAC/B,2BAA2B,EACzB,KAAK;QACP,iBAAiB,EAAE,KAAK;QACxB,mBAAmB,EAAE,KAAK;QAC1B,iCAAiC,EAC/B,KAAK;QACP,2BAA2B,EACzB,KAAK;QACP,8CAA8C,EAC5C,KAAK;QACP,iBAAiB,EAAE;YACjB,aAAa,EACX,WAAW,CAAC,QAAQ;YACtB,cAAc,EACZ,OAAO,CAAC,QAAQ;YAClB,YAAY,EACV,oBAAoB;YACtB,cAAc,EACZ,OAAO,CAAC,QAAQ;YAClB,kBAAkB,EAChB,yBAAyB;SAC5B;QACD,WAAW;KACZ,EACD,qCAAqC,EACrC,OAAO,EACP,WAAW,CACZ,CAAC;AACJ,CAAC","sourcesContent":["// SPDX-License-Identifier: Apache-2.0\n\n/**\n * Scene-to-sensor optical bridge: from spectral scene radiance to the irradiance that reaches the\n * sensor.\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 here.\n *\n * The calculation applies the declared pupil acceptance and exactly one transmission or\n * field-throughput path. It returns the irradiance before the sensor stack, with conservative\n * evidence.\n *\n * The optical bridge turns outgoing scene radiance (W/m²/sr/nm) into irradiance (W/m²/nm) arriving\n * at the front of the sensor stack.\n * - Circular paraxial acceptance is pi/(4*Nworking²). At finite focus with unit pupil magnification,\n *   the working f-number is N*(1+m).\n * - Spectral transmission and an effective working T-stop are alternative paths; use one.\n * - Field throughput is applied exactly once.\n * - Sensor response and PSF (point-spread function) blur happen later in the pipeline.\n * @see docs/SCENE_SENSOR_QUADRATURE.md for the equations, units, coordinate conventions, known\n * blockers and 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 EvidenceBackedFact,\n  type EvidenceProvenance\n} from \"../core/evidence-provenance.js\";\nimport {\n  parseSpectralWavelengthBasis,\n  type SpectralWavelengthBasis\n} from \"../core/spectral.js\";\nimport { InvalidScientificInputError } from \"../core/validation.js\";\nimport {\n  composeFrontOfLensFilterTransmission,\n  type FrontOfLensFilterProfile\n} from \"./front-of-lens-filter.js\";\nimport type { IlluminationVignetting } from \"./illumination-vignetting.js\";\nimport { calculateThinLensImageDistance } from \"./thin-lens.js\";\nimport {\n  parseSceneRadianceEvaluationRequest,\n  parseSceneRadianceEvaluationResult,\n  type SceneRadianceEvaluationRequest,\n  type SceneRadianceEvaluationResult\n} from \"../schema/scene-radiance.js\";\nimport { requireInputKeys, rejectUnknownFields, rejectUnknownVariantFields, UNCERTAINTY_VARIANT_FIELDS } from \"../core/record-validation.js\";\n\ntype UnknownRecord = Record<string, unknown>;\n\nexport const SCENE_TO_SENSOR_IRRADIANCE_PROFILE_SCHEMA_VERSION =\n  \"0.1.0\" as const;\n\nexport type OpticalBridgeScientificStatus =\n  | \"calibrated\"\n  | \"approximation\";\n\nexport type OpticalBridgeUncertainty =\n  | {\n      kind: \"relative\";\n      fraction: number;\n      basis: string;\n    }\n  | {\n      kind: \"not-quantified\";\n      limitation: string;\n    };\n\nexport interface NumericRange {\n  minimum: number;\n  maximum: number;\n}\n\nexport type OpticalBridgeFocusApplicability =\n  | {\n      kind: \"any\";\n    }\n  | {\n      kind: \"finite-distance-range\";\n      objectDistanceM: NumericRange;\n    };\n\nexport interface SpectralTransmissionSample {\n  wavelengthNanometers: number;\n  linearTransmissionFactor: number;\n}\n\nexport type OpticalTransmissionModel =\n  | {\n      kind: \"spectral-transmission\";\n      scientificStatus:\n        OpticalBridgeScientificStatus;\n      wavelengthBasis:\n        Exclude<\n          SpectralWavelengthBasis,\n          \"unspecified\"\n        >;\n      samples:\n        EvidenceBackedFact<\n          readonly SpectralTransmissionSample[]\n        >;\n      uncertainty:\n        OpticalBridgeUncertainty;\n    }\n  | {\n      kind:\n        \"effective-working-t-stop-approximation\";\n      scientificStatus: \"approximation\";\n      workingTStop:\n        EvidenceBackedFact<number>;\n      wavelengthBasis:\n        Exclude<\n          SpectralWavelengthBasis,\n          \"unspecified\"\n        >;\n      wavelengthRangeNanometers:\n        NumericRange;\n      uncertainty:\n        OpticalBridgeUncertainty;\n      limitation: string;\n    };\n\nexport interface SceneToSensorIrradianceProfile {\n  schemaVersion:\n    typeof SCENE_TO_SENSOR_IRRADIANCE_PROFILE_SCHEMA_VERSION;\n  profileId: string;\n  profileVersion: string;\n  lensProfileId: string;\n  scientificStatus:\n    OpticalBridgeScientificStatus;\n  applicability: {\n    focalLengthMm: NumericRange;\n    nominalFNumber: NumericRange;\n    focus:\n      OpticalBridgeFocusApplicability;\n  };\n  transmission:\n    OpticalTransmissionModel;\n  distortionAreaMappingOwnership:\n    \"not-applied-by-bridge\";\n  psfRedistributionOwnership:\n    \"downstream-normalized-energy-redistribution\";\n  sensorOpticalStackIncluded: false;\n  strayLightIncluded: false;\n  polarizationModeled: false;\n  wavelengthChangingBehaviorModeled: false;\n  volumetricScatteringModeled: false;\n  evidence: readonly EvidenceProvenance[];\n  limitations: readonly string[];\n}\n\nexport type OpticalBridgeFocusContext =\n  | {\n      kind: \"infinity-focus\";\n    }\n  | {\n      kind:\n        \"ideal-symmetric-thin-lens\";\n      objectDistanceM: number;\n      pupilMagnificationAssumption:\n        \"unity\";\n    }\n  | {\n      kind:\n        \"supplied-working-f-number\";\n      focus:\n        | {\n            kind: \"infinity\";\n          }\n        | {\n            kind: \"finite\";\n            objectDistanceM: number;\n          };\n      workingFNumber:\n        EvidenceBackedFact<number>;\n      basis: string;\n    };\n\nexport type OpticalBridgeFieldThroughput =\n  | {\n      kind: \"unity\";\n    }\n  | {\n      kind:\n        \"illumination-vignetting-result\";\n      result:\n        IlluminationVignetting;\n    };\n\nexport interface CalculateSceneRadianceToSensorIrradianceInput {\n  sceneRadianceRequest:\n    SceneRadianceEvaluationRequest;\n  sceneRadianceResult:\n    SceneRadianceEvaluationResult;\n  profile:\n    SceneToSensorIrradianceProfile;\n  focalLengthMm: number;\n  nominalFNumber: number;\n  focus:\n    OpticalBridgeFocusContext;\n  imagePointMm: {\n    x: number;\n    y: number;\n  };\n  fieldThroughput:\n    OpticalBridgeFieldThroughput;\n  frontOfLensFilters?:\n    readonly FrontOfLensFilterProfile[];\n}\n\nexport interface SceneToSensorIrradianceResult {\n  schemaVersion:\n    typeof SCENE_TO_SENSOR_IRRADIANCE_PROFILE_SCHEMA_VERSION;\n  profileId: string;\n  profileVersion: string;\n  lensProfileId: string;\n  sampleId: string;\n  providerProfileId: string;\n  sceneId: string;\n  timeSecondsFromExposureStart: number;\n  wavelengthNanometers: number;\n  wavelengthBasis:\n    Exclude<\n      SpectralWavelengthBasis,\n      \"unspecified\"\n    >;\n  inputQuantity:\n    \"outgoing-spectral-radiance\";\n  inputUnit:\n    \"W/m^2/sr/nm\";\n  outputQuantity:\n    \"sensor-plane-spectral-irradiance\";\n  outputUnit:\n    \"W/m^2/nm\";\n  sceneSpectralRadianceWattsPerSquareMeterSteradianNanometer:\n    number;\n  focalLengthMm: number;\n  nominalFNumber: number;\n  workingFNumber: number;\n  focusModel:\n    OpticalBridgeFocusContext[\"kind\"];\n  paraxialGeometricAcceptanceSolidAngleSteradians:\n    number;\n  transmissionPath:\n    OpticalTransmissionModel[\"kind\"];\n  spectralTransmissionFactor:\n    number | null;\n  effectiveWorkingTStop:\n    number | null;\n  tStopEquivalentTransmissionFactor:\n    number | null;\n  frontOfLensFilterCount: number;\n  frontOfLensFilterTransmissionFactor:\n    number;\n  frontOfLensFilterAttenuationStops:\n    number;\n  frontOfLensFilterTransmissionAppliedExactlyOnce:\n    true;\n  /**\n   * True when a front-of-lens linear polarizer was applied. Scene radiance is unpolarized in this\n   * bridge, so a polarizer transmits (k1 + k2) / 2.\n   */\n  frontOfLensFilterPolarizationModeled:\n    boolean;\n  fieldThroughputFactor: number;\n  effectiveAcceptanceAfterTransmissionAndFieldSteradians:\n    number;\n  sensorPlaneSpectralIrradianceWattsPerSquareMeterNanometer:\n    number;\n  scientificStatus: \"approximation\";\n  sceneRadianceScientificStatus:\n    SceneRadianceEvaluationResult[\"scientificStatus\"];\n  opticalProfileScientificStatus:\n    OpticalBridgeScientificStatus;\n  sceneRadianceUncertainty:\n    SceneRadianceEvaluationResult[\"uncertainty\"];\n  opticalUncertainty:\n    OpticalBridgeUncertainty;\n  radianceAmplificationApplied: false;\n  fieldThroughputAppliedExactlyOnce:\n    true;\n  universalCos4FalloffApplied: false;\n  distortionAreaCorrectionApplied: false;\n  psfRedistributionApplied: false;\n  sensorOpticalStackApplied: false;\n  cfaApplied: false;\n  quantumEfficiencyApplied: false;\n  spectralResponsivityApplied: false;\n  strayLightApplied: false;\n  polarizationModeled: false;\n  wavelengthChangingBehaviorModeled:\n    false;\n  volumetricScatteringModeled: false;\n  calibratedSensorPlaneIrradianceClaimAuthorized:\n    false;\n  componentEvidence: {\n    sceneRadiance:\n      readonly EvidenceProvenance[];\n    opticalProfile:\n      readonly EvidenceProvenance[];\n    transmission:\n      readonly EvidenceProvenance[];\n    workingFNumber:\n      readonly EvidenceProvenance[];\n    frontOfLensFilters:\n      readonly EvidenceProvenance[];\n  };\n  limitations:\n    readonly string[];\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 + \" 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 requireFraction(\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 be from zero through one.\"\n    );\n  }\n  return parsed;\n}\n\nfunction parseRange(\n  value: unknown,\n  path: string\n): NumericRange {\n  const record =\n    requireRecord(value, path, [\"minimum\", \"maximum\"]);\n  const minimum =\n    requirePositiveFinite(\n      record.minimum,\n      path + \".minimum\"\n    );\n  const maximum =\n    requirePositiveFinite(\n      record.maximum,\n      path + \".maximum\"\n    );\n  if (minimum > maximum) {\n    throw new InvalidConfigurationError(\n      path +\n        \".minimum must be less than or equal to maximum.\"\n    );\n  }\n  return {\n    minimum,\n    maximum\n  };\n}\n\nfunction parseScientificStatus(\n  value: unknown,\n  path: string\n): OpticalBridgeScientificStatus {\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): OpticalBridgeUncertainty {\n  const record =\n    requireRecord(value, path, null);\n  rejectUnknownVariantFields(record, \"kind\", UNCERTAINTY_VARIANT_FIELDS, path);\n\n  if (record.kind === \"relative\") {\n    return {\n      kind: \"relative\",\n      fraction:\n        requireFraction(\n          record.fraction,\n          path + \".fraction\"\n        ),\n      basis:\n        requireNonEmptyString(\n          record.basis,\n          path + \".basis\"\n        )\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 + \".limitation\"\n        )\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 = 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 parseReusableNumericEvidence(\n  value: unknown,\n  path: string\n): readonly EvidenceProvenance[] {\n  const evidence =\n    parseEvidenceList(value, path);\n\n  if (\n    evidence.some(\n      (entry) =>\n        entry.reuseStatus ===\n        \"factual-reference-only\"\n    )\n  ) {\n    throw new InvalidConfigurationError(\n      path +\n        \" for numeric optical curve data must be reusable-data or photivra-owned.\"\n    );\n  }\n  return evidence;\n}\n\nfunction parseSpectralSamples(\n  value: unknown,\n  path: string\n): readonly SpectralTransmissionSample[] {\n  if (\n    !Array.isArray(value) ||\n    value.length < 2\n  ) {\n    throw new InvalidConfigurationError(\n      path +\n        \" must contain at least two samples.\"\n    );\n  }\n\n  let previous =\n    Number.NEGATIVE_INFINITY;\n  return value.map(\n    (entry, index) => {\n      const samplePath =\n        path +\n        \"[\" +\n        index +\n        \"]\";\n      const record =\n        requireRecord(\n          entry,\n          samplePath, [\"wavelengthNanometers\", \"linearTransmissionFactor\"]\n        );\n      const wavelengthNanometers =\n        requirePositiveFinite(\n          record.wavelengthNanometers,\n          samplePath +\n            \".wavelengthNanometers\"\n        );\n      if (\n        wavelengthNanometers <=\n        previous\n      ) {\n        throw new InvalidConfigurationError(\n          path +\n            \" wavelengths must be strictly increasing.\"\n        );\n      }\n      previous =\n        wavelengthNanometers;\n      return {\n        wavelengthNanometers,\n        linearTransmissionFactor:\n          requireFraction(\n            record.linearTransmissionFactor,\n            samplePath +\n              \".linearTransmissionFactor\"\n          )\n      };\n    }\n  );\n}\n\nfunction parseTransmission(\n  value: unknown,\n  path: string\n): OpticalTransmissionModel {\n  const record =\n    requireRecord(value, path, null);\n  rejectUnknownVariantFields(record, \"kind\", {\n    // workingTStop and wavelengthRangeNanometers (spectral) and samples (working T-stop) are rejected below with their own messages.\n    \"spectral-transmission\": [\"kind\", \"scientificStatus\", \"wavelengthBasis\", \"samples\", \"uncertainty\", \"workingTStop\", \"wavelengthRangeNanometers\"],\n    \"effective-working-t-stop-approximation\": [\"kind\", \"scientificStatus\", \"workingTStop\", \"wavelengthBasis\", \"wavelengthRangeNanometers\",\n      \"uncertainty\", \"limitation\", \"samples\"]\n  }, path);\n\n  if (\n    record.kind ===\n    \"spectral-transmission\"\n  ) {\n    if (\n      record.workingTStop !==\n      undefined ||\n      record.wavelengthRangeNanometers !==\n      undefined\n    ) {\n      throw new InvalidConfigurationError(\n        path +\n          \" spectral-transmission must not include T-stop-only fields.\"\n      );\n    }\n\n    const status =\n      parseScientificStatus(\n        record.scientificStatus,\n        path +\n          \".scientificStatus\"\n      );\n    const samplesFact =\n      requireRecord(\n        record.samples,\n        path + \".samples\", [\"value\", \"evidence\"]\n      );\n    const uncertainty =\n      parseUncertainty(\n        record.uncertainty,\n        path + \".uncertainty\"\n      );\n\n    if (\n      status === \"calibrated\" &&\n      uncertainty.kind !==\n        \"relative\"\n    ) {\n      throw new InvalidConfigurationError(\n        path +\n          \" calibrated spectral transmission requires quantified relative uncertainty.\"\n      );\n    }\n\n    return {\n      kind:\n        \"spectral-transmission\",\n      scientificStatus: status,\n      wavelengthBasis:\n        parseResolvedBasis(\n          record.wavelengthBasis,\n          path +\n            \".wavelengthBasis\"\n        ),\n      samples: {\n        value:\n          parseSpectralSamples(\n            samplesFact.value,\n            path +\n              \".samples.value\"\n          ),\n        evidence:\n          parseReusableNumericEvidence(\n            samplesFact.evidence,\n            path +\n              \".samples.evidence\"\n          )\n      },\n      uncertainty\n    };\n  }\n\n  if (\n    record.kind ===\n    \"effective-working-t-stop-approximation\"\n  ) {\n    if (\n      record.samples !==\n      undefined\n    ) {\n      throw new InvalidConfigurationError(\n        path +\n          \" T-stop approximation must not include spectral transmission samples.\"\n      );\n    }\n    if (\n      record.scientificStatus !==\n      \"approximation\"\n    ) {\n      throw new InvalidConfigurationError(\n        path +\n          ' T-stop scientificStatus must be \"approximation\".'\n      );\n    }\n    const workingTStopFact =\n      requireRecord(\n        record.workingTStop,\n        path +\n          \".workingTStop\", [\"value\", \"evidence\"]\n      );\n    return {\n      kind:\n        \"effective-working-t-stop-approximation\",\n      scientificStatus:\n        \"approximation\",\n      workingTStop: {\n        value:\n          requirePositiveFinite(\n            workingTStopFact.value,\n            path +\n              \".workingTStop.value\"\n          ),\n        evidence:\n          parseEvidenceList(\n            workingTStopFact.evidence,\n            path +\n              \".workingTStop.evidence\"\n          )\n      },\n      wavelengthBasis:\n        parseResolvedBasis(\n          record.wavelengthBasis,\n          path +\n            \".wavelengthBasis\"\n        ),\n      wavelengthRangeNanometers:\n        parseRange(\n          record.wavelengthRangeNanometers,\n          path +\n            \".wavelengthRangeNanometers\"\n        ),\n      uncertainty:\n        parseUncertainty(\n          record.uncertainty,\n          path + \".uncertainty\"\n        ),\n      limitation:\n        requireNonEmptyString(\n          record.limitation,\n          path +\n            \".limitation\"\n        )\n    };\n  }\n\n  throw new InvalidConfigurationError(\n    path + \".kind is invalid.\"\n  );\n}\n\nfunction parseResolvedBasis(\n  value: unknown,\n  path: string\n): Exclude<\n  SpectralWavelengthBasis,\n  \"unspecified\"\n> {\n  const parsed =\n    parseSpectralWavelengthBasis(\n      value,\n      path\n    );\n  if (parsed === \"unspecified\") {\n    throw new InvalidConfigurationError(\n      path +\n        ' must resolve to \"air\" or \"vacuum\".'\n    );\n  }\n  return parsed;\n}\n\nfunction parseFocusApplicability(\n  value: unknown,\n  path: string\n): OpticalBridgeFocusApplicability {\n  const record =\n    requireRecord(value, path, null);\n  rejectUnknownVariantFields(record, \"kind\", { any: [\"kind\"], \"finite-distance-range\": [\"kind\", \"objectDistanceM\"] }, path);\n\n  if (record.kind === \"any\") {\n    return {\n      kind: \"any\"\n    };\n  }\n  if (\n    record.kind ===\n    \"finite-distance-range\"\n  ) {\n    return {\n      kind:\n        \"finite-distance-range\",\n      objectDistanceM:\n        parseRange(\n          record.objectDistanceM,\n          path +\n            \".objectDistanceM\"\n        )\n    };\n  }\n  throw new InvalidConfigurationError(\n    path + \".kind is invalid.\"\n  );\n}\n\n/**\n * Validates an untrusted scene-to-sensor profile and returns the normalized typed contract. Unknown\n * enum values, missing required fields and incompatible scientific data are rejected here.\n *\n * The optical bridge turns outgoing scene radiance (W/m²/sr/nm) into irradiance (W/m²/nm) arriving\n * at the front of the sensor stack.\n * - Circular paraxial acceptance is pi/(4*Nworking²). At finite focus with unit pupil magnification,\n *   the working f-number is N*(1+m).\n * - Spectral transmission and an effective working T-stop are alternative paths; use one.\n * - Field throughput is applied exactly once.\n * - Sensor response and PSF (point-spread function) blur happen later in the pipeline.\n * @param value - unknown. Untrusted input; TypeScript types alone are not validation.\n * @returns SceneToSensorIrradianceProfile. 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 parseSceneToSensorIrradianceProfile(\n  value: unknown\n): SceneToSensorIrradianceProfile {\n  const record =\n    requireRecord(\n      value,\n      \"sceneToSensorIrradianceProfile\", [\"schemaVersion\", \"applicability\", \"scientificStatus\", \"transmission\", \"distortionAreaMappingOwnership\", \"psfRedistributionOwnership\", \"sensorOpticalStackIncluded\", \"strayLightIncluded\", \"polarizationModeled\", \"wavelengthChangingBehaviorModeled\", \"volumetricScatteringModeled\", \"profileId\", \"profileVersion\", \"lensProfileId\", \"evidence\", \"limitations\"]\n    );\n  if (\n    record.schemaVersion !==\n    SCENE_TO_SENSOR_IRRADIANCE_PROFILE_SCHEMA_VERSION\n  ) {\n    throw new InvalidConfigurationError(\n      'sceneToSensorIrradianceProfile.schemaVersion must be \"' +\n        SCENE_TO_SENSOR_IRRADIANCE_PROFILE_SCHEMA_VERSION +\n        '\".'\n    );\n  }\n\n  const applicability =\n    requireRecord(\n      record.applicability,\n      \"sceneToSensorIrradianceProfile.applicability\", [\"focalLengthMm\", \"nominalFNumber\", \"focus\"]\n    );\n  const status =\n    parseScientificStatus(\n      record.scientificStatus,\n      \"sceneToSensorIrradianceProfile.scientificStatus\"\n    );\n  const transmission =\n    parseTransmission(\n      record.transmission,\n      \"sceneToSensorIrradianceProfile.transmission\"\n    );\n\n  if (\n    status === \"calibrated\" &&\n    transmission.scientificStatus !==\n      \"calibrated\"\n  ) {\n    throw new InvalidConfigurationError(\n      \"A calibrated optical bridge profile requires calibrated spectral transmission.\"\n    );\n  }\n\n  if (\n    record.distortionAreaMappingOwnership !==\n      \"not-applied-by-bridge\" ||\n    record.psfRedistributionOwnership !==\n      \"downstream-normalized-energy-redistribution\" ||\n    record.sensorOpticalStackIncluded !==\n      false ||\n    record.strayLightIncluded !==\n      false ||\n    record.polarizationModeled !==\n      false ||\n    record.wavelengthChangingBehaviorModeled !==\n      false ||\n    record.volumetricScatteringModeled !==\n      false\n  ) {\n    throw new InvalidConfigurationError(\n      \"sceneToSensorIrradianceProfile ownership/completeness flags are invalid for schema 0.1.0.\"\n    );\n  }\n\n  return {\n    schemaVersion:\n      SCENE_TO_SENSOR_IRRADIANCE_PROFILE_SCHEMA_VERSION,\n    profileId:\n      requireNonEmptyString(\n        record.profileId,\n        \"sceneToSensorIrradianceProfile.profileId\"\n      ),\n    profileVersion:\n      requireNonEmptyString(\n        record.profileVersion,\n        \"sceneToSensorIrradianceProfile.profileVersion\"\n      ),\n    lensProfileId:\n      requireNonEmptyString(\n        record.lensProfileId,\n        \"sceneToSensorIrradianceProfile.lensProfileId\"\n      ),\n    scientificStatus: status,\n    applicability: {\n      focalLengthMm:\n        parseRange(\n          applicability.focalLengthMm,\n          \"sceneToSensorIrradianceProfile.applicability.focalLengthMm\"\n        ),\n      nominalFNumber:\n        parseRange(\n          applicability.nominalFNumber,\n          \"sceneToSensorIrradianceProfile.applicability.nominalFNumber\"\n        ),\n      focus:\n        parseFocusApplicability(\n          applicability.focus,\n          \"sceneToSensorIrradianceProfile.applicability.focus\"\n        )\n    },\n    transmission,\n    distortionAreaMappingOwnership:\n      \"not-applied-by-bridge\",\n    psfRedistributionOwnership:\n      \"downstream-normalized-energy-redistribution\",\n    sensorOpticalStackIncluded:\n      false,\n    strayLightIncluded: false,\n    polarizationModeled: false,\n    wavelengthChangingBehaviorModeled:\n      false,\n    volumetricScatteringModeled:\n      false,\n    evidence:\n      parseEvidenceList(\n        record.evidence,\n        \"sceneToSensorIrradianceProfile.evidence\"\n      ),\n    limitations:\n      parseLimitations(\n        record.limitations,\n        \"sceneToSensorIrradianceProfile.limitations\"\n      )\n  };\n}\n\nfunction ensureWithin(\n  value: number,\n  range: NumericRange,\n  path: string\n): void {\n  if (\n    value < range.minimum ||\n    value > range.maximum\n  ) {\n    throw new InvalidScientificInputError(\n      path +\n        \" lies outside the optical profile applicability range.\"\n    );\n  }\n}\n\n/**\n * Internal (2.6; not exported from the package): N_w / N = 1 + m for an ideal symmetric thin lens with unity pupil\n * magnification focused at `objectDistanceM`, m = v/s from `calculateThinLensImageDistance`. The irradiance bridge's\n * `ideal-symmetric-thin-lens` focus uses exactly this factor, and the framing reports derive their close-up light\n * factor 1/(1 + m)² from it, so the two cannot drift apart.\n */\nexport function thinLensWorkingFNumberFactor(focalLengthMm: number, objectDistanceM: number): number {\n  return 1 + calculateThinLensImageDistance({ focalLengthMm, objectDistanceM }).value.magnification;\n}\n\nfunction workingFNumber(\n  focus:\n    OpticalBridgeFocusContext,\n  focalLengthMm: number,\n  nominalFNumber: number\n): {\n  value: number;\n  evidence:\n    readonly EvidenceProvenance[];\n} {\n  if (\n    focus.kind === \"infinity-focus\"\n  ) {\n    return {\n      value: nominalFNumber,\n      evidence: []\n    };\n  }\n\n  if (\n    focus.kind ===\n    \"ideal-symmetric-thin-lens\"\n  ) {\n    if (\n      focus\n        .pupilMagnificationAssumption !==\n      \"unity\"\n    ) {\n      throw new InvalidScientificInputError(\n        'Ideal symmetric thin-lens focus requires pupilMagnificationAssumption \"unity\".'\n      );\n    }\n    return {\n      value:\n        nominalFNumber *\n        thinLensWorkingFNumberFactor(focalLengthMm, focus.objectDistanceM),\n      evidence: []\n    };\n  }\n\n  if (\n    focus.kind ===\n    \"supplied-working-f-number\"\n  ) {\n    const value =\n      requirePositiveFiniteScientific(\n        focus.workingFNumber.value,\n        \"focus.workingFNumber.value\"\n      );\n    return {\n      value,\n      evidence:\n        parseEvidenceList(\n          focus.workingFNumber.evidence,\n          \"focus.workingFNumber.evidence\"\n        )\n    };\n  }\n\n  throw new InvalidConfigurationError(\n    \"focus.kind is invalid.\"\n  );\n}\n\nfunction requirePositiveFiniteScientific(\n  value: unknown,\n  path: string\n): number {\n  if (\n    typeof value !== \"number\" ||\n    !Number.isFinite(value) ||\n    value <= 0\n  ) {\n    throw new InvalidScientificInputError(\n      path +\n        \" must be finite and greater than zero.\"\n    );\n  }\n  return value;\n}\n\nfunction finiteFocusDistance(\n  focus:\n    OpticalBridgeFocusContext\n): number | null {\n  if (\n    focus.kind ===\n    \"ideal-symmetric-thin-lens\"\n  ) {\n    return requirePositiveFiniteScientific(\n      focus.objectDistanceM,\n      \"focus.objectDistanceM\"\n    );\n  }\n  if (\n    focus.kind ===\n      \"supplied-working-f-number\" &&\n    focus.focus.kind === \"finite\"\n  ) {\n    return requirePositiveFiniteScientific(\n      focus.focus.objectDistanceM,\n      \"focus.focus.objectDistanceM\"\n    );\n  }\n  return null;\n}\n\nfunction interpolateTransmission(\n  model:\n    Extract<\n      OpticalTransmissionModel,\n      { kind: \"spectral-transmission\" }\n    >,\n  wavelengthNanometers: number\n): number {\n  const samples =\n    model.samples.value;\n  const first =\n    samples[0]!;\n  const last =\n    samples[\n      samples.length - 1\n    ]!;\n\n  if (\n    wavelengthNanometers <\n      first.wavelengthNanometers ||\n    wavelengthNanometers >\n      last.wavelengthNanometers\n  ) {\n    throw new InvalidScientificInputError(\n      \"Scene wavelength lies outside the spectral transmission sample range.\"\n    );\n  }\n\n  if (\n    wavelengthNanometers ===\n    first.wavelengthNanometers\n  ) {\n    return first\n      .linearTransmissionFactor;\n  }\n\n  for (\n    let index = 1;\n    index < samples.length;\n    index += 1\n  ) {\n    const right =\n      samples[index]!;\n    if (\n      wavelengthNanometers <=\n      right.wavelengthNanometers\n    ) {\n      const left =\n        samples[index - 1]!;\n      if (\n        wavelengthNanometers ===\n        right.wavelengthNanometers\n      ) {\n        return right\n          .linearTransmissionFactor;\n      }\n      const phase =\n        (wavelengthNanometers -\n          left.wavelengthNanometers) /\n        (right.wavelengthNanometers -\n          left.wavelengthNanometers);\n      return (\n        left\n          .linearTransmissionFactor +\n        phase *\n          (right\n              .linearTransmissionFactor -\n            left\n              .linearTransmissionFactor)\n      );\n    }\n  }\n\n  return last\n    .linearTransmissionFactor;\n}\n\nfunction fieldFactor(\n  field:\n    OpticalBridgeFieldThroughput,\n  imagePointMm: {\n    x: number;\n    y: number;\n  }\n): number {\n  if (field.kind === \"unity\") {\n    return 1;\n  }\n  if (\n    field.kind ===\n    \"illumination-vignetting-result\"\n  ) {\n    const result =\n      field.result;\n    const scale = Math.max(\n      1,\n      Math.abs(imagePointMm.x),\n      Math.abs(imagePointMm.y),\n      Math.abs(result.imagePointMm.x),\n      Math.abs(result.imagePointMm.y)\n    );\n    const tolerance =\n      Number.EPSILON *\n      64 *\n      scale;\n    if (\n      Math.abs(\n        result.imagePointMm.x -\n          imagePointMm.x\n      ) > tolerance ||\n      Math.abs(\n        result.imagePointMm.y -\n          imagePointMm.y\n      ) > tolerance\n    ) {\n      throw new InvalidScientificInputError(\n        \"Illumination-vignetting result imagePointMm must match the optical bridge imagePointMm.\"\n      );\n    }\n    const factor =\n      result.linearThroughputFactor;\n    if (\n      !Number.isFinite(factor) ||\n      factor <= 0 ||\n      factor > 1\n    ) {\n      throw new InvalidScientificInputError(\n        \"Field throughput factor must be finite and in (0, 1].\"\n      );\n    }\n    return factor;\n  }\n  throw new InvalidConfigurationError(\n    \"fieldThroughput.kind is invalid.\"\n  );\n}\n\nfunction validateSceneIdentity(\n  request:\n    SceneRadianceEvaluationRequest,\n  result:\n    SceneRadianceEvaluationResult\n): void {\n  if (\n    request.sampleId !==\n      result.sampleId ||\n    request.providerProfileId !==\n      result.providerProfileId ||\n    request.sceneId !==\n      result.sceneId ||\n    request.wavelengthNanometers !==\n      result.wavelengthNanometers ||\n    request.wavelengthBasis !==\n      result.wavelengthBasis\n  ) {\n    throw new InvalidScientificInputError(\n      \"Scene-radiance request/result identity must match exactly.\"\n    );\n  }\n}\n\n/**\n * Applies the declared pupil acceptance and exactly one transmission or field-throughput path to\n * spectral scene radiance. Returns the irradiance before the sensor stack, with conservative\n * evidence.\n *\n * The optical bridge turns outgoing scene radiance (W/m²/sr/nm) into irradiance (W/m²/nm) arriving\n * at the front of the sensor stack.\n * - Circular paraxial acceptance is pi/(4*Nworking²). At finite focus with unit pupil magnification,\n *   the working f-number is N*(1+m).\n * - Spectral transmission and an effective working T-stop are alternative paths; use one.\n * - Field throughput is applied exactly once.\n * - Sensor response and PSF (point-spread function) blur happen later in the pipeline.\n * @param input - CalculateSceneRadianceToSensorIrradianceInput. See that type for coordinates, units and profile binding.\n * @returns CalculationResult<\n  SceneToSensorIrradianceResult\n>. 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 calculateSceneRadianceToSensorIrradiance(\n  input:\n    CalculateSceneRadianceToSensorIrradianceInput\n): CalculationResult<\n  SceneToSensorIrradianceResult\n> {\n  requireInputKeys(input, \"calculateSceneRadianceToSensorIrradiance\", [\"sceneRadianceRequest\", \"sceneRadianceResult\", \"profile\", \"focalLengthMm\", \"nominalFNumber\", \"focus\", \"imagePointMm\", \"fieldThroughput\", \"frontOfLensFilters\"]);\n  const request =\n    parseSceneRadianceEvaluationRequest(\n      input.sceneRadianceRequest\n    );\n  const sceneResult =\n    parseSceneRadianceEvaluationResult(\n      input.sceneRadianceResult\n    );\n  validateSceneIdentity(\n    request,\n    sceneResult\n  );\n\n  const profile =\n    parseSceneToSensorIrradianceProfile(\n      input.profile\n    );\n  const focalLengthMm =\n    requirePositiveFiniteScientific(\n      input.focalLengthMm,\n      \"focalLengthMm\"\n    );\n  const nominalFNumber =\n    requirePositiveFiniteScientific(\n      input.nominalFNumber,\n      \"nominalFNumber\"\n    );\n\n  ensureWithin(\n    focalLengthMm,\n    profile.applicability\n      .focalLengthMm,\n    \"focalLengthMm\"\n  );\n  ensureWithin(\n    nominalFNumber,\n    profile.applicability\n      .nominalFNumber,\n    \"nominalFNumber\"\n  );\n\n  const focusDistanceM =\n    finiteFocusDistance(input.focus);\n  if (\n    profile.applicability.focus.kind ===\n      \"finite-distance-range\"\n  ) {\n    if (focusDistanceM === null) {\n      throw new InvalidScientificInputError(\n        \"Optical profile requires a finite focus distance.\"\n      );\n    }\n    ensureWithin(\n      focusDistanceM,\n      profile.applicability\n        .focus.objectDistanceM,\n      \"focus object distance\"\n    );\n  }\n\n  const working =\n    workingFNumber(\n      input.focus,\n      focalLengthMm,\n      nominalFNumber\n    );\n  const workingFNumberValue =\n    working.value;\n  const geometricAcceptance =\n    Math.PI /\n    (4 *\n      workingFNumberValue *\n      workingFNumberValue);\n\n  if (\n    !Number.isFinite(\n      geometricAcceptance\n    ) ||\n    geometricAcceptance <= 0\n  ) {\n    throw new InvalidScientificInputError(\n      \"Paraxial geometric acceptance must remain finite and greater than zero.\"\n    );\n  }\n\n  const transmission =\n    profile.transmission;\n\n  if (\n    request.wavelengthBasis !==\n    transmission.wavelengthBasis\n  ) {\n    throw new InvalidScientificInputError(\n      \"Scene-radiance and optical-transmission wavelength bases must match; implicit air/vacuum conversion is not performed.\"\n    );\n  }\n\n  let spectralTransmissionFactor:\n    number | null = null;\n  let effectiveWorkingTStop:\n    number | null = null;\n  let tStopEquivalentTransmissionFactor:\n    number | null = null;\n  let effectiveAcceptance =\n    geometricAcceptance;\n  let transmissionEvidence:\n    readonly EvidenceProvenance[];\n\n  if (\n    transmission.kind ===\n    \"spectral-transmission\"\n  ) {\n    spectralTransmissionFactor =\n      interpolateTransmission(\n        transmission,\n        request.wavelengthNanometers\n      );\n    effectiveAcceptance *=\n      spectralTransmissionFactor;\n    transmissionEvidence =\n      transmission.samples.evidence;\n  } else {\n    ensureWithin(\n      request.wavelengthNanometers,\n      transmission\n        .wavelengthRangeNanometers,\n      \"scene wavelength\"\n    );\n    effectiveWorkingTStop =\n      transmission.workingTStop.value;\n    if (\n      effectiveWorkingTStop <\n      workingFNumberValue\n    ) {\n      throw new InvalidScientificInputError(\n        \"effective working T-stop must not be smaller than the geometric working f-number for a passive optical approximation.\"\n      );\n    }\n    tStopEquivalentTransmissionFactor =\n      (workingFNumberValue /\n        effectiveWorkingTStop) **\n      2;\n    effectiveAcceptance =\n      Math.PI /\n      (4 *\n        effectiveWorkingTStop *\n        effectiveWorkingTStop);\n    transmissionEvidence =\n      transmission\n        .workingTStop.evidence;\n  }\n\n  let frontOfLensFilterCount = 0;\n  let frontOfLensFilterPolarizationModeled = false;\n  let frontOfLensFilterTransmissionFactor = 1;\n  let frontOfLensFilterAttenuationStops = 0;\n  let frontOfLensFilterEvidence:\n    readonly EvidenceProvenance[] = [];\n  let frontOfLensFilterLimitations:\n    readonly string[] = [];\n\n  if (\n    input.frontOfLensFilters !==\n    undefined\n  ) {\n    const composed =\n      composeFrontOfLensFilterTransmission({\n        filters:\n          input.frontOfLensFilters,\n        wavelengthNanometers:\n          request.wavelengthNanometers,\n        wavelengthBasis:\n          request.wavelengthBasis\n      });\n    frontOfLensFilterCount =\n      composed.appliedFilterCount;\n    frontOfLensFilterPolarizationModeled =\n      composed.polarizationModeled;\n    frontOfLensFilterTransmissionFactor =\n      composed\n        .combinedLinearTransmissionFactor;\n    frontOfLensFilterAttenuationStops =\n      composed\n        .combinedAttenuationStops;\n    frontOfLensFilterEvidence =\n      composed.components.flatMap(\n        (component) =>\n          component.evidence\n      );\n    frontOfLensFilterLimitations =\n      composed.components.flatMap(\n        (component) =>\n          component.limitations\n      );\n    effectiveAcceptance *=\n      frontOfLensFilterTransmissionFactor;\n  }\n\n  const fieldThroughputFactor =\n    fieldFactor(\n      input.fieldThroughput,\n      input.imagePointMm\n    );\n  effectiveAcceptance *=\n    fieldThroughputFactor;\n\n  const sceneRadiance =\n    sceneResult\n      .spectralRadianceWattsPerSquareMeterSteradianNanometer;\n  const irradiance =\n    sceneRadiance *\n    effectiveAcceptance;\n\n  if (\n    !Number.isFinite(irradiance) ||\n    irradiance < 0\n  ) {\n    throw new InvalidScientificInputError(\n      \"Sensor-plane spectral irradiance must remain finite and non-negative.\"\n    );\n  }\n\n  const limitations = [\n    ...profile.limitations,\n    ...frontOfLensFilterLimitations,\n    \"Front-of-lens filter transmission is applied outside the lens throughput path exactly once and does not change geometric working f-number.\",\n    \"Paraxial circular-pupil acceptance uses pi/(4*N_working^2); no universal cos^4 field falloff is added.\",\n    \"Field throughput is supplied separately and applied exactly once.\",\n    \"Geometric distortion area-density/Jacobian correction is not applied by this bridge.\",\n    \"PSF/diffraction energy redistribution is downstream and must preserve normalized energy apart from separately owned throughput losses.\",\n    \"Sensor optical stack, microlens, CFA, QE/responsivity and charge conversion are downstream.\",\n    frontOfLensFilterPolarizationModeled\n      ? \"Flare, ghosting, veiling glare, scene-radiance polarization (scene light is treated as unpolarized; the front-of-lens linear polarizer transmits (k1 + k2) / 2), volumetric scattering and wavelength-changing behavior are not modeled.\"\n      : \"Flare, ghosting, veiling glare, polarization, volumetric scattering and wavelength-changing behavior are not modeled.\"\n  ];\n\n  return approximationResult(\n    {\n      schemaVersion:\n        SCENE_TO_SENSOR_IRRADIANCE_PROFILE_SCHEMA_VERSION,\n      profileId: profile.profileId,\n      profileVersion:\n        profile.profileVersion,\n      lensProfileId:\n        profile.lensProfileId,\n      sampleId: request.sampleId,\n      providerProfileId:\n        request.providerProfileId,\n      sceneId: request.sceneId,\n      timeSecondsFromExposureStart:\n        request\n          .timeSecondsFromExposureStart,\n      wavelengthNanometers:\n        request.wavelengthNanometers,\n      wavelengthBasis:\n        request.wavelengthBasis,\n      inputQuantity:\n        \"outgoing-spectral-radiance\",\n      inputUnit:\n        \"W/m^2/sr/nm\",\n      outputQuantity:\n        \"sensor-plane-spectral-irradiance\",\n      outputUnit: \"W/m^2/nm\",\n      sceneSpectralRadianceWattsPerSquareMeterSteradianNanometer:\n        sceneRadiance,\n      focalLengthMm,\n      nominalFNumber,\n      workingFNumber:\n        workingFNumberValue,\n      focusModel: input.focus.kind,\n      paraxialGeometricAcceptanceSolidAngleSteradians:\n        geometricAcceptance,\n      transmissionPath:\n        transmission.kind,\n      spectralTransmissionFactor,\n      effectiveWorkingTStop,\n      tStopEquivalentTransmissionFactor,\n      frontOfLensFilterCount,\n      frontOfLensFilterTransmissionFactor,\n      frontOfLensFilterAttenuationStops,\n      frontOfLensFilterTransmissionAppliedExactlyOnce:\n        true,\n      frontOfLensFilterPolarizationModeled,\n      fieldThroughputFactor,\n      effectiveAcceptanceAfterTransmissionAndFieldSteradians:\n        effectiveAcceptance,\n      sensorPlaneSpectralIrradianceWattsPerSquareMeterNanometer:\n        irradiance,\n      scientificStatus:\n        \"approximation\",\n      sceneRadianceScientificStatus:\n        sceneResult.scientificStatus,\n      opticalProfileScientificStatus:\n        profile.scientificStatus,\n      sceneRadianceUncertainty:\n        sceneResult.uncertainty,\n      opticalUncertainty:\n        transmission.uncertainty,\n      radianceAmplificationApplied:\n        false,\n      fieldThroughputAppliedExactlyOnce:\n        true,\n      universalCos4FalloffApplied:\n        false,\n      distortionAreaCorrectionApplied:\n        false,\n      psfRedistributionApplied:\n        false,\n      sensorOpticalStackApplied:\n        false,\n      cfaApplied: false,\n      quantumEfficiencyApplied: false,\n      spectralResponsivityApplied:\n        false,\n      strayLightApplied: false,\n      polarizationModeled: false,\n      wavelengthChangingBehaviorModeled:\n        false,\n      volumetricScatteringModeled:\n        false,\n      calibratedSensorPlaneIrradianceClaimAuthorized:\n        false,\n      componentEvidence: {\n        sceneRadiance:\n          sceneResult.evidence,\n        opticalProfile:\n          profile.evidence,\n        transmission:\n          transmissionEvidence,\n        workingFNumber:\n          working.evidence,\n        frontOfLensFilters:\n          frontOfLensFilterEvidence\n      },\n      limitations\n    },\n    \"scene-radiance-to-sensor-irradiance\",\n    \"1.1.0\",\n    limitations\n  );\n}\n"]}
+{"version":3,"file":"scene-to-sensor-irradiance.js","sourceRoot":"","sources":["../../src/optics/scene-to-sensor-irradiance.ts"],"names":[],"mappings":"AAAA,sCAAsC;AAEtC;;;;;;;;;;;;;;;;;;;;GAoBG;AAEH,OAAO,EACL,mBAAmB,EAEpB,MAAM,+BAA+B,CAAC;AACvC,OAAO,EAAE,yBAAyB,EAAE,MAAM,gCAAgC,CAAC;AAC3E,OAAO,EACL,iBAAiB,EAGlB,MAAM,gCAAgC,CAAC;AACxC,OAAO,EACL,4BAA4B,EAE7B,MAAM,qBAAqB,CAAC;AAC7B,OAAO,EAAE,2BAA2B,EAAE,MAAM,uBAAuB,CAAC;AACpE,OAAO,EACL,oCAAoC,EAErC,MAAM,2BAA2B,CAAC;AAEnC,OAAO,EAAE,8BAA8B,EAAE,MAAM,gBAAgB,CAAC;AAChE,OAAO,EACL,mCAAmC,EACnC,kCAAkC,EAGnC,MAAM,6BAA6B,CAAC;AACrC,OAAO,EAAE,gBAAgB,EAAE,mBAAmB,EAAE,0BAA0B,EAAE,0BAA0B,EAAE,MAAM,8BAA8B,CAAC;AAI7I,MAAM,CAAC,MAAM,iDAAiD,GAC5D,OAAgB,CAAC;AAoQnB,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,GAAG,8BAA8B,CACtC,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,eAAe,CACtB,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,iCAAiC,CACpC,CAAC;IACJ,CAAC;IACD,OAAO,MAAM,CAAC;AAChB,CAAC;AAED,SAAS,UAAU,CACjB,KAAc,EACd,IAAY;IAEZ,MAAM,MAAM,GACV,aAAa,CAAC,KAAK,EAAE,IAAI,EAAE,CAAC,SAAS,EAAE,SAAS,CAAC,CAAC,CAAC;IACrD,MAAM,OAAO,GACX,qBAAqB,CACnB,MAAM,CAAC,OAAO,EACd,IAAI,GAAG,UAAU,CAClB,CAAC;IACJ,MAAM,OAAO,GACX,qBAAqB,CACnB,MAAM,CAAC,OAAO,EACd,IAAI,GAAG,UAAU,CAClB,CAAC;IACJ,IAAI,OAAO,GAAG,OAAO,EAAE,CAAC;QACtB,MAAM,IAAI,yBAAyB,CACjC,IAAI;YACF,iDAAiD,CACpD,CAAC;IACJ,CAAC;IACD,OAAO;QACL,OAAO;QACP,OAAO;KACR,CAAC;AACJ,CAAC;AAED,SAAS,qBAAqB,CAC5B,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;IAE7E,IAAI,MAAM,CAAC,IAAI,KAAK,UAAU,EAAE,CAAC;QAC/B,OAAO;YACL,IAAI,EAAE,UAAU;YAChB,QAAQ,EACN,eAAe,CACb,MAAM,CAAC,QAAQ,EACf,IAAI,GAAG,WAAW,CACnB;YACH,KAAK,EACH,qBAAqB,CACnB,MAAM,CAAC,KAAK,EACZ,IAAI,GAAG,QAAQ,CAChB;SACJ,CAAC;IACJ,CAAC;IAED,IACE,MAAM,CAAC,IAAI;QACX,gBAAgB,EAChB,CAAC;QACD,OAAO;YACL,IAAI,EAAE,gBAAgB;YACtB,UAAU,EACR,qBAAqB,CACnB,MAAM,CAAC,UAAU,EACjB,IAAI,GAAG,aAAa,CACrB;SACJ,CAAC;IACJ,CAAC;IAED,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,GAAG,KAAK,CAAC,GAAG,CACtB,CAAC,KAAK,EAAE,KAAK,EAAE,EAAE,CACf,qBAAqB,CACnB,KAAK,EACL,IAAI;QACF,GAAG;QACH,KAAK;QACL,GAAG,CACN,CACJ,CAAC;IACF,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,4BAA4B,CACnC,KAAc,EACd,IAAY;IAEZ,MAAM,QAAQ,GACZ,iBAAiB,CAAC,KAAK,EAAE,IAAI,CAAC,CAAC;IAEjC,IACE,QAAQ,CAAC,IAAI,CACX,CAAC,KAAK,EAAE,EAAE,CACR,KAAK,CAAC,WAAW;QACjB,wBAAwB,CAC3B,EACD,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,IAAI;YACF,0EAA0E,CAC7E,CAAC;IACJ,CAAC;IACD,OAAO,QAAQ,CAAC;AAClB,CAAC;AAED,SAAS,oBAAoB,CAC3B,KAAc,EACd,IAAY;IAEZ,IACE,CAAC,KAAK,CAAC,OAAO,CAAC,KAAK,CAAC;QACrB,KAAK,CAAC,MAAM,GAAG,CAAC,EAChB,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,IAAI;YACF,qCAAqC,CACxC,CAAC;IACJ,CAAC;IAED,IAAI,QAAQ,GACV,MAAM,CAAC,iBAAiB,CAAC;IAC3B,OAAO,KAAK,CAAC,GAAG,CACd,CAAC,KAAK,EAAE,KAAK,EAAE,EAAE;QACf,MAAM,UAAU,GACd,IAAI;YACJ,GAAG;YACH,KAAK;YACL,GAAG,CAAC;QACN,MAAM,MAAM,GACV,aAAa,CACX,KAAK,EACL,UAAU,EAAE,CAAC,sBAAsB,EAAE,0BAA0B,CAAC,CACjE,CAAC;QACJ,MAAM,oBAAoB,GACxB,qBAAqB,CACnB,MAAM,CAAC,oBAAoB,EAC3B,UAAU;YACR,uBAAuB,CAC1B,CAAC;QACJ,IACE,oBAAoB;YACpB,QAAQ,EACR,CAAC;YACD,MAAM,IAAI,yBAAyB,CACjC,IAAI;gBACF,2CAA2C,CAC9C,CAAC;QACJ,CAAC;QACD,QAAQ;YACN,oBAAoB,CAAC;QACvB,OAAO;YACL,oBAAoB;YACpB,wBAAwB,EACtB,eAAe,CACb,MAAM,CAAC,wBAAwB,EAC/B,UAAU;gBACR,2BAA2B,CAC9B;SACJ,CAAC;IACJ,CAAC,CACF,CAAC;AACJ,CAAC;AAED,SAAS,iBAAiB,CACxB,KAAc,EACd,IAAY;IAEZ,MAAM,MAAM,GACV,aAAa,CAAC,KAAK,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;IACnC,0BAA0B,CAAC,MAAM,EAAE,MAAM,EAAE;QACzC,iIAAiI;QACjI,uBAAuB,EAAE,CAAC,MAAM,EAAE,kBAAkB,EAAE,iBAAiB,EAAE,SAAS,EAAE,aAAa,EAAE,cAAc,EAAE,2BAA2B,CAAC;QAC/I,wCAAwC,EAAE,CAAC,MAAM,EAAE,kBAAkB,EAAE,cAAc,EAAE,iBAAiB,EAAE,2BAA2B;YACnI,aAAa,EAAE,YAAY,EAAE,SAAS,CAAC;KAC1C,EAAE,IAAI,CAAC,CAAC;IAET,IACE,MAAM,CAAC,IAAI;QACX,uBAAuB,EACvB,CAAC;QACD,IACE,MAAM,CAAC,YAAY;YACnB,SAAS;YACT,MAAM,CAAC,yBAAyB;gBAChC,SAAS,EACT,CAAC;YACD,MAAM,IAAI,yBAAyB,CACjC,IAAI;gBACF,6DAA6D,CAChE,CAAC;QACJ,CAAC;QAED,MAAM,MAAM,GACV,qBAAqB,CACnB,MAAM,CAAC,gBAAgB,EACvB,IAAI;YACF,mBAAmB,CACtB,CAAC;QACJ,MAAM,WAAW,GACf,aAAa,CACX,MAAM,CAAC,OAAO,EACd,IAAI,GAAG,UAAU,EAAE,CAAC,OAAO,EAAE,UAAU,CAAC,CACzC,CAAC;QACJ,MAAM,WAAW,GACf,gBAAgB,CACd,MAAM,CAAC,WAAW,EAClB,IAAI,GAAG,cAAc,CACtB,CAAC;QAEJ,IACE,MAAM,KAAK,YAAY;YACvB,WAAW,CAAC,IAAI;gBACd,UAAU,EACZ,CAAC;YACD,MAAM,IAAI,yBAAyB,CACjC,IAAI;gBACF,6EAA6E,CAChF,CAAC;QACJ,CAAC;QAED,OAAO;YACL,IAAI,EACF,uBAAuB;YACzB,gBAAgB,EAAE,MAAM;YACxB,eAAe,EACb,kBAAkB,CAChB,MAAM,CAAC,eAAe,EACtB,IAAI;gBACF,kBAAkB,CACrB;YACH,OAAO,EAAE;gBACP,KAAK,EACH,oBAAoB,CAClB,WAAW,CAAC,KAAK,EACjB,IAAI;oBACF,gBAAgB,CACnB;gBACH,QAAQ,EACN,4BAA4B,CAC1B,WAAW,CAAC,QAAQ,EACpB,IAAI;oBACF,mBAAmB,CACtB;aACJ;YACD,WAAW;SACZ,CAAC;IACJ,CAAC;IAED,IACE,MAAM,CAAC,IAAI;QACX,wCAAwC,EACxC,CAAC;QACD,IACE,MAAM,CAAC,OAAO;YACd,SAAS,EACT,CAAC;YACD,MAAM,IAAI,yBAAyB,CACjC,IAAI;gBACF,uEAAuE,CAC1E,CAAC;QACJ,CAAC;QACD,IACE,MAAM,CAAC,gBAAgB;YACvB,eAAe,EACf,CAAC;YACD,MAAM,IAAI,yBAAyB,CACjC,IAAI;gBACF,mDAAmD,CACtD,CAAC;QACJ,CAAC;QACD,MAAM,gBAAgB,GACpB,aAAa,CACX,MAAM,CAAC,YAAY,EACnB,IAAI;YACF,eAAe,EAAE,CAAC,OAAO,EAAE,UAAU,CAAC,CACzC,CAAC;QACJ,OAAO;YACL,IAAI,EACF,wCAAwC;YAC1C,gBAAgB,EACd,eAAe;YACjB,YAAY,EAAE;gBACZ,KAAK,EACH,qBAAqB,CACnB,gBAAgB,CAAC,KAAK,EACtB,IAAI;oBACF,qBAAqB,CACxB;gBACH,QAAQ,EACN,iBAAiB,CACf,gBAAgB,CAAC,QAAQ,EACzB,IAAI;oBACF,wBAAwB,CAC3B;aACJ;YACD,eAAe,EACb,kBAAkB,CAChB,MAAM,CAAC,eAAe,EACtB,IAAI;gBACF,kBAAkB,CACrB;YACH,yBAAyB,EACvB,UAAU,CACR,MAAM,CAAC,yBAAyB,EAChC,IAAI;gBACF,4BAA4B,CAC/B;YACH,WAAW,EACT,gBAAgB,CACd,MAAM,CAAC,WAAW,EAClB,IAAI,GAAG,cAAc,CACtB;YACH,UAAU,EACR,qBAAqB,CACnB,MAAM,CAAC,UAAU,EACjB,IAAI;gBACF,aAAa,CAChB;SACJ,CAAC;IACJ,CAAC;IAED,MAAM,IAAI,yBAAyB,CACjC,IAAI,GAAG,mBAAmB,CAC3B,CAAC;AACJ,CAAC;AAED,SAAS,kBAAkB,CACzB,KAAc,EACd,IAAY;IAKZ,MAAM,MAAM,GACV,4BAA4B,CAC1B,KAAK,EACL,IAAI,CACL,CAAC;IACJ,IAAI,MAAM,KAAK,aAAa,EAAE,CAAC;QAC7B,MAAM,IAAI,yBAAyB,CACjC,IAAI;YACF,qCAAqC,CACxC,CAAC;IACJ,CAAC;IACD,OAAO,MAAM,CAAC;AAChB,CAAC;AAED,SAAS,uBAAuB,CAC9B,KAAc,EACd,IAAY;IAEZ,MAAM,MAAM,GACV,aAAa,CAAC,KAAK,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;IACnC,0BAA0B,CAAC,MAAM,EAAE,MAAM,EAAE,EAAE,GAAG,EAAE,CAAC,MAAM,CAAC,EAAE,uBAAuB,EAAE,CAAC,MAAM,EAAE,iBAAiB,CAAC,EAAE,EAAE,IAAI,CAAC,CAAC;IAE1H,IAAI,MAAM,CAAC,IAAI,KAAK,KAAK,EAAE,CAAC;QAC1B,OAAO;YACL,IAAI,EAAE,KAAK;SACZ,CAAC;IACJ,CAAC;IACD,IACE,MAAM,CAAC,IAAI;QACX,uBAAuB,EACvB,CAAC;QACD,OAAO;YACL,IAAI,EACF,uBAAuB;YACzB,eAAe,EACb,UAAU,CACR,MAAM,CAAC,eAAe,EACtB,IAAI;gBACF,kBAAkB,CACrB;SACJ,CAAC;IACJ,CAAC;IACD,MAAM,IAAI,yBAAyB,CACjC,IAAI,GAAG,mBAAmB,CAC3B,CAAC;AACJ,CAAC;AAED;;;;;;;;;;;;;;;GAeG;AACH,MAAM,UAAU,mCAAmC,CACjD,KAAc;IAEd,MAAM,MAAM,GACV,aAAa,CACX,KAAK,EACL,gCAAgC,EAAE,CAAC,eAAe,EAAE,eAAe,EAAE,kBAAkB,EAAE,cAAc,EAAE,gCAAgC,EAAE,4BAA4B,EAAE,4BAA4B,EAAE,oBAAoB,EAAE,qBAAqB,EAAE,mCAAmC,EAAE,6BAA6B,EAAE,WAAW,EAAE,gBAAgB,EAAE,eAAe,EAAE,UAAU,EAAE,aAAa,CAAC,CACnY,CAAC;IACJ,IACE,MAAM,CAAC,aAAa;QACpB,iDAAiD,EACjD,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,wDAAwD;YACtD,iDAAiD;YACjD,IAAI,CACP,CAAC;IACJ,CAAC;IAED,MAAM,aAAa,GACjB,aAAa,CACX,MAAM,CAAC,aAAa,EACpB,8CAA8C,EAAE,CAAC,eAAe,EAAE,gBAAgB,EAAE,OAAO,CAAC,CAC7F,CAAC;IACJ,MAAM,MAAM,GACV,qBAAqB,CACnB,MAAM,CAAC,gBAAgB,EACvB,iDAAiD,CAClD,CAAC;IACJ,MAAM,YAAY,GAChB,iBAAiB,CACf,MAAM,CAAC,YAAY,EACnB,6CAA6C,CAC9C,CAAC;IAEJ,IACE,MAAM,KAAK,YAAY;QACvB,YAAY,CAAC,gBAAgB;YAC3B,YAAY,EACd,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,gFAAgF,CACjF,CAAC;IACJ,CAAC;IAED,IACE,MAAM,CAAC,8BAA8B;QACnC,uBAAuB;QACzB,MAAM,CAAC,0BAA0B;YAC/B,6CAA6C;QAC/C,MAAM,CAAC,0BAA0B;YAC/B,KAAK;QACP,MAAM,CAAC,kBAAkB;YACvB,KAAK;QACP,MAAM,CAAC,mBAAmB;YACxB,KAAK;QACP,MAAM,CAAC,iCAAiC;YACtC,KAAK;QACP,MAAM,CAAC,2BAA2B;YAChC,KAAK,EACP,CAAC;QACD,MAAM,IAAI,yBAAyB,CACjC,2FAA2F,CAC5F,CAAC;IACJ,CAAC;IAED,OAAO;QACL,aAAa,EACX,iDAAiD;QACnD,SAAS,EACP,qBAAqB,CACnB,MAAM,CAAC,SAAS,EAChB,0CAA0C,CAC3C;QACH,cAAc,EACZ,qBAAqB,CACnB,MAAM,CAAC,cAAc,EACrB,+CAA+C,CAChD;QACH,aAAa,EACX,qBAAqB,CACnB,MAAM,CAAC,aAAa,EACpB,8CAA8C,CAC/C;QACH,gBAAgB,EAAE,MAAM;QACxB,aAAa,EAAE;YACb,aAAa,EACX,UAAU,CACR,aAAa,CAAC,aAAa,EAC3B,4DAA4D,CAC7D;YACH,cAAc,EACZ,UAAU,CACR,aAAa,CAAC,cAAc,EAC5B,6DAA6D,CAC9D;YACH,KAAK,EACH,uBAAuB,CACrB,aAAa,CAAC,KAAK,EACnB,oDAAoD,CACrD;SACJ;QACD,YAAY;QACZ,8BAA8B,EAC5B,uBAAuB;QACzB,0BAA0B,EACxB,6CAA6C;QAC/C,0BAA0B,EACxB,KAAK;QACP,kBAAkB,EAAE,KAAK;QACzB,mBAAmB,EAAE,KAAK;QAC1B,iCAAiC,EAC/B,KAAK;QACP,2BAA2B,EACzB,KAAK;QACP,QAAQ,EACN,iBAAiB,CACf,MAAM,CAAC,QAAQ,EACf,yCAAyC,CAC1C;QACH,WAAW,EACT,gBAAgB,CACd,MAAM,CAAC,WAAW,EAClB,4CAA4C,CAC7C;KACJ,CAAC;AACJ,CAAC;AAED,SAAS,YAAY,CACnB,KAAa,EACb,KAAmB,EACnB,IAAY;IAEZ,IACE,KAAK,GAAG,KAAK,CAAC,OAAO;QACrB,KAAK,GAAG,KAAK,CAAC,OAAO,EACrB,CAAC;QACD,MAAM,IAAI,2BAA2B,CACnC,IAAI;YACF,wDAAwD,CAC3D,CAAC;IACJ,CAAC;AACH,CAAC;AAED;;;;;GAKG;AACH,MAAM,UAAU,4BAA4B,CAAC,aAAqB,EAAE,eAAuB;IACzF,OAAO,CAAC,GAAG,8BAA8B,CAAC,EAAE,aAAa,EAAE,eAAe,EAAE,CAAC,CAAC,KAAK,CAAC,aAAa,CAAC;AACpG,CAAC;AAED,SAAS,cAAc,CACrB,KAC2B,EAC3B,aAAqB,EACrB,cAAsB;IAMtB,IACE,KAAK,CAAC,IAAI,KAAK,gBAAgB,EAC/B,CAAC;QACD,OAAO;YACL,KAAK,EAAE,cAAc;YACrB,QAAQ,EAAE,EAAE;SACb,CAAC;IACJ,CAAC;IAED,IACE,KAAK,CAAC,IAAI;QACV,2BAA2B,EAC3B,CAAC;QACD,IACE,KAAK;aACF,4BAA4B;YAC/B,OAAO,EACP,CAAC;YACD,MAAM,IAAI,2BAA2B,CACnC,gFAAgF,CACjF,CAAC;QACJ,CAAC;QACD,OAAO;YACL,KAAK,EACH,cAAc;gBACd,4BAA4B,CAAC,aAAa,EAAE,KAAK,CAAC,eAAe,CAAC;YACpE,QAAQ,EAAE,EAAE;SACb,CAAC;IACJ,CAAC;IAED,IACE,KAAK,CAAC,IAAI;QACV,2BAA2B,EAC3B,CAAC;QACD,MAAM,KAAK,GACT,+BAA+B,CAC7B,KAAK,CAAC,cAAc,CAAC,KAAK,EAC1B,4BAA4B,CAC7B,CAAC;QACJ,OAAO;YACL,KAAK;YACL,QAAQ,EACN,iBAAiB,CACf,KAAK,CAAC,cAAc,CAAC,QAAQ,EAC7B,+BAA+B,CAChC;SACJ,CAAC;IACJ,CAAC;IAED,MAAM,IAAI,yBAAyB,CACjC,wBAAwB,CACzB,CAAC;AACJ,CAAC;AAED,SAAS,+BAA+B,CACtC,KAAc,EACd,IAAY;IAEZ,IACE,OAAO,KAAK,KAAK,QAAQ;QACzB,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC;QACvB,KAAK,IAAI,CAAC,EACV,CAAC;QACD,MAAM,IAAI,2BAA2B,CACnC,IAAI;YACF,wCAAwC,CAC3C,CAAC;IACJ,CAAC;IACD,OAAO,KAAK,CAAC;AACf,CAAC;AAED,SAAS,mBAAmB,CAC1B,KAC2B;IAE3B,IACE,KAAK,CAAC,IAAI;QACV,2BAA2B,EAC3B,CAAC;QACD,OAAO,+BAA+B,CACpC,KAAK,CAAC,eAAe,EACrB,uBAAuB,CACxB,CAAC;IACJ,CAAC;IACD,IACE,KAAK,CAAC,IAAI;QACR,2BAA2B;QAC7B,KAAK,CAAC,KAAK,CAAC,IAAI,KAAK,QAAQ,EAC7B,CAAC;QACD,OAAO,+BAA+B,CACpC,KAAK,CAAC,KAAK,CAAC,eAAe,EAC3B,6BAA6B,CAC9B,CAAC;IACJ,CAAC;IACD,OAAO,IAAI,CAAC;AACd,CAAC;AAED,SAAS,uBAAuB,CAC9B,KAIG,EACH,oBAA4B;IAE5B,MAAM,OAAO,GACX,KAAK,CAAC,OAAO,CAAC,KAAK,CAAC;IACtB,MAAM,KAAK,GACT,OAAO,CAAC,CAAC,CAAE,CAAC;IACd,MAAM,IAAI,GACR,OAAO,CACL,OAAO,CAAC,MAAM,GAAG,CAAC,CAClB,CAAC;IAEL,IACE,oBAAoB;QAClB,KAAK,CAAC,oBAAoB;QAC5B,oBAAoB;YAClB,IAAI,CAAC,oBAAoB,EAC3B,CAAC;QACD,MAAM,IAAI,2BAA2B,CACnC,uEAAuE,CACxE,CAAC;IACJ,CAAC;IAED,IACE,oBAAoB;QACpB,KAAK,CAAC,oBAAoB,EAC1B,CAAC;QACD,OAAO,KAAK;aACT,wBAAwB,CAAC;IAC9B,CAAC;IAED,KACE,IAAI,KAAK,GAAG,CAAC,EACb,KAAK,GAAG,OAAO,CAAC,MAAM,EACtB,KAAK,IAAI,CAAC,EACV,CAAC;QACD,MAAM,KAAK,GACT,OAAO,CAAC,KAAK,CAAE,CAAC;QAClB,IACE,oBAAoB;YACpB,KAAK,CAAC,oBAAoB,EAC1B,CAAC;YACD,MAAM,IAAI,GACR,OAAO,CAAC,KAAK,GAAG,CAAC,CAAE,CAAC;YACtB,IACE,oBAAoB;gBACpB,KAAK,CAAC,oBAAoB,EAC1B,CAAC;gBACD,OAAO,KAAK;qBACT,wBAAwB,CAAC;YAC9B,CAAC;YACD,MAAM,KAAK,GACT,CAAC,oBAAoB;gBACnB,IAAI,CAAC,oBAAoB,CAAC;gBAC5B,CAAC,KAAK,CAAC,oBAAoB;oBACzB,IAAI,CAAC,oBAAoB,CAAC,CAAC;YAC/B,OAAO,CACL,IAAI;iBACD,wBAAwB;gBAC3B,KAAK;oBACH,CAAC,KAAK;yBACD,wBAAwB;wBAC3B,IAAI;6BACD,wBAAwB,CAAC,CACjC,CAAC;QACJ,CAAC;IACH,CAAC;IAED,OAAO,IAAI;SACR,wBAAwB,CAAC;AAC9B,CAAC;AAED,SAAS,WAAW,CAClB,KAC8B,EAC9B,YAGC;IAED,IAAI,KAAK,CAAC,IAAI,KAAK,OAAO,EAAE,CAAC;QAC3B,OAAO,CAAC,CAAC;IACX,CAAC;IACD,IACE,KAAK,CAAC,IAAI;QACV,gCAAgC,EAChC,CAAC;QACD,MAAM,MAAM,GACV,KAAK,CAAC,MAAM,CAAC;QACf,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CACpB,CAAC,EACD,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC,CAAC,EACxB,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC,CAAC,EACxB,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,YAAY,CAAC,CAAC,CAAC,EAC/B,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,YAAY,CAAC,CAAC,CAAC,CAChC,CAAC;QACF,MAAM,SAAS,GACb,MAAM,CAAC,OAAO;YACd,EAAE;YACF,KAAK,CAAC;QACR,IACE,IAAI,CAAC,GAAG,CACN,MAAM,CAAC,YAAY,CAAC,CAAC;YACnB,YAAY,CAAC,CAAC,CACjB,GAAG,SAAS;YACb,IAAI,CAAC,GAAG,CACN,MAAM,CAAC,YAAY,CAAC,CAAC;gBACnB,YAAY,CAAC,CAAC,CACjB,GAAG,SAAS,EACb,CAAC;YACD,MAAM,IAAI,2BAA2B,CACnC,yFAAyF,CAC1F,CAAC;QACJ,CAAC;QACD,MAAM,MAAM,GACV,MAAM,CAAC,sBAAsB,CAAC;QAChC,IACE,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC;YACxB,MAAM,IAAI,CAAC;YACX,MAAM,GAAG,CAAC,EACV,CAAC;YACD,MAAM,IAAI,2BAA2B,CACnC,uDAAuD,CACxD,CAAC;QACJ,CAAC;QACD,OAAO,MAAM,CAAC;IAChB,CAAC;IACD,MAAM,IAAI,yBAAyB,CACjC,kCAAkC,CACnC,CAAC;AACJ,CAAC;AAED,SAAS,qBAAqB,CAC5B,OACgC,EAChC,MAC+B;IAE/B,IACE,OAAO,CAAC,QAAQ;QACd,MAAM,CAAC,QAAQ;QACjB,OAAO,CAAC,iBAAiB;YACvB,MAAM,CAAC,iBAAiB;QAC1B,OAAO,CAAC,OAAO;YACb,MAAM,CAAC,OAAO;QAChB,OAAO,CAAC,oBAAoB;YAC1B,MAAM,CAAC,oBAAoB;QAC7B,OAAO,CAAC,eAAe;YACrB,MAAM,CAAC,eAAe,EACxB,CAAC;QACD,MAAM,IAAI,2BAA2B,CACnC,4DAA4D,CAC7D,CAAC;IACJ,CAAC;AACH,CAAC;AAgBD,SAAS,sCAAsC,CAC7C,KAAoH;IAEpH,MAAM,OAAO,GACX,mCAAmC,CACjC,KAAK,CAAC,OAAO,CACd,CAAC;IACJ,MAAM,aAAa,GACjB,+BAA+B,CAC7B,KAAK,CAAC,aAAa,EACnB,eAAe,CAChB,CAAC;IACJ,MAAM,cAAc,GAClB,+BAA+B,CAC7B,KAAK,CAAC,cAAc,EACpB,gBAAgB,CACjB,CAAC;IAEJ,YAAY,CACV,aAAa,EACb,OAAO,CAAC,aAAa;SAClB,aAAa,EAChB,eAAe,CAChB,CAAC;IACF,YAAY,CACV,cAAc,EACd,OAAO,CAAC,aAAa;SAClB,cAAc,EACjB,gBAAgB,CACjB,CAAC;IAEF,MAAM,cAAc,GAClB,mBAAmB,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC;IACnC,IACE,OAAO,CAAC,aAAa,CAAC,KAAK,CAAC,IAAI;QAC9B,uBAAuB,EACzB,CAAC;QACD,IAAI,cAAc,KAAK,IAAI,EAAE,CAAC;YAC5B,MAAM,IAAI,2BAA2B,CACnC,mDAAmD,CACpD,CAAC;QACJ,CAAC;QACD,YAAY,CACV,cAAc,EACd,OAAO,CAAC,aAAa;aAClB,KAAK,CAAC,eAAe,EACxB,uBAAuB,CACxB,CAAC;IACJ,CAAC;IAED,MAAM,OAAO,GACX,cAAc,CACZ,KAAK,CAAC,KAAK,EACX,aAAa,EACb,cAAc,CACf,CAAC;IACJ,MAAM,mBAAmB,GACvB,OAAO,CAAC,KAAK,CAAC;IAChB,MAAM,mBAAmB,GACvB,IAAI,CAAC,EAAE;QACP,CAAC,CAAC;YACA,mBAAmB;YACnB,mBAAmB,CAAC,CAAC;IAEzB,IACE,CAAC,MAAM,CAAC,QAAQ,CACd,mBAAmB,CACpB;QACD,mBAAmB,IAAI,CAAC,EACxB,CAAC;QACD,MAAM,IAAI,2BAA2B,CACnC,yEAAyE,CAC1E,CAAC;IACJ,CAAC;IAED,MAAM,YAAY,GAChB,OAAO,CAAC,YAAY,CAAC;IACvB,OAAO,EAAE,OAAO,EAAE,aAAa,EAAE,cAAc,EAAE,OAAO,EAAE,mBAAmB,EAAE,mBAAmB,EAAE,YAAY,EAAE,CAAC;AACrH,CAAC;AAmBD;;;;GAIG;AACH,SAAS,oCAAoC,CAC3C,QAAiD,EACjD,KAC0E,EAC1E,OAAuC,EACvC,WAA0C;IAE1C,MAAM,EAAE,mBAAmB,EAAE,mBAAmB,EAAE,YAAY,EAAE,GAAG,QAAQ,CAAC;IAC5E,IACE,OAAO,CAAC,eAAe;QACvB,YAAY,CAAC,eAAe,EAC5B,CAAC;QACD,MAAM,IAAI,2BAA2B,CACnC,uHAAuH,CACxH,CAAC;IACJ,CAAC;IAED,IAAI,0BAA0B,GACZ,IAAI,CAAC;IACvB,IAAI,qBAAqB,GACP,IAAI,CAAC;IACvB,IAAI,iCAAiC,GACnB,IAAI,CAAC;IACvB,IAAI,mBAAmB,GACrB,mBAAmB,CAAC;IACtB,IAAI,oBAC2B,CAAC;IAEhC,IACE,YAAY,CAAC,IAAI;QACjB,uBAAuB,EACvB,CAAC;QACD,0BAA0B;YACxB,uBAAuB,CACrB,YAAY,EACZ,OAAO,CAAC,oBAAoB,CAC7B,CAAC;QACJ,mBAAmB;YACjB,0BAA0B,CAAC;QAC7B,oBAAoB;YAClB,YAAY,CAAC,OAAO,CAAC,QAAQ,CAAC;IAClC,CAAC;SAAM,CAAC;QACN,YAAY,CACV,OAAO,CAAC,oBAAoB,EAC5B,YAAY;aACT,yBAAyB,EAC5B,kBAAkB,CACnB,CAAC;QACF,qBAAqB;YACnB,YAAY,CAAC,YAAY,CAAC,KAAK,CAAC;QAClC,IACE,qBAAqB;YACrB,mBAAmB,EACnB,CAAC;YACD,MAAM,IAAI,2BAA2B,CACnC,uHAAuH,CACxH,CAAC;QACJ,CAAC;QACD,iCAAiC;YAC/B,CAAC,mBAAmB;gBAClB,qBAAqB,CAAC;gBACxB,CAAC,CAAC;QACJ,mBAAmB;YACjB,IAAI,CAAC,EAAE;gBACP,CAAC,CAAC;oBACA,qBAAqB;oBACrB,qBAAqB,CAAC,CAAC;QAC3B,oBAAoB;YAClB,YAAY;iBACT,YAAY,CAAC,QAAQ,CAAC;IAC7B,CAAC;IAED,IAAI,sBAAsB,GAAG,CAAC,CAAC;IAC/B,IAAI,oCAAoC,GAAG,KAAK,CAAC;IACjD,IAAI,mCAAmC,GAAG,CAAC,CAAC;IAC5C,IAAI,iCAAiC,GAAG,CAAC,CAAC;IAC1C,IAAI,yBAAyB,GACK,EAAE,CAAC;IACrC,IAAI,4BAA4B,GACV,EAAE,CAAC;IAEzB,IACE,KAAK,CAAC,kBAAkB;QACxB,SAAS,EACT,CAAC;QACD,MAAM,QAAQ,GACZ,oCAAoC,CAAC;YACnC,OAAO,EACL,KAAK,CAAC,kBAAkB;YAC1B,oBAAoB,EAClB,OAAO,CAAC,oBAAoB;YAC9B,eAAe,EACb,OAAO,CAAC,eAAe;SAC1B,CAAC,CAAC;QACL,sBAAsB;YACpB,QAAQ,CAAC,kBAAkB,CAAC;QAC9B,oCAAoC;YAClC,QAAQ,CAAC,mBAAmB,CAAC;QAC/B,mCAAmC;YACjC,QAAQ;iBACL,gCAAgC,CAAC;QACtC,iCAAiC;YAC/B,QAAQ;iBACL,wBAAwB,CAAC;QAC9B,yBAAyB;YACvB,QAAQ,CAAC,UAAU,CAAC,OAAO,CACzB,CAAC,SAAS,EAAE,EAAE,CACZ,SAAS,CAAC,QAAQ,CACrB,CAAC;QACJ,4BAA4B;YAC1B,QAAQ,CAAC,UAAU,CAAC,OAAO,CACzB,CAAC,SAAS,EAAE,EAAE,CACZ,SAAS,CAAC,WAAW,CACxB,CAAC;QACJ,mBAAmB;YACjB,mCAAmC,CAAC;IACxC,CAAC;IAED,MAAM,qBAAqB,GACzB,WAAW,CACT,KAAK,CAAC,eAAe,EACrB,KAAK,CAAC,YAAY,CACnB,CAAC;IACJ,mBAAmB;QACjB,qBAAqB,CAAC;IAExB,MAAM,aAAa,GACjB,WAAW;SACR,qDAAqD,CAAC;IAC3D,MAAM,UAAU,GACd,aAAa;QACb,mBAAmB,CAAC;IAEtB,IACE,CAAC,MAAM,CAAC,QAAQ,CAAC,UAAU,CAAC;QAC5B,UAAU,GAAG,CAAC,EACd,CAAC;QACD,MAAM,IAAI,2BAA2B,CACnC,uEAAuE,CACxE,CAAC;IACJ,CAAC;IAED,OAAO,EAAE,0BAA0B,EAAE,qBAAqB,EAAE,iCAAiC,EAAE,mBAAmB,EAAE,oBAAoB;QACtI,sBAAsB,EAAE,oCAAoC,EAAE,mCAAmC,EAAE,iCAAiC;QACpI,yBAAyB,EAAE,4BAA4B,EAAE,qBAAqB,EAAE,aAAa,EAAE,UAAU,EAAE,CAAC;AAChH,CAAC;AAED;;;;;;;;GAQG;AACH,MAAM,UAAU,2CAA2C,CACzD,MAA4I;IAG5I,MAAM,QAAQ,GAAG,sCAAsC,CAAC,MAAM,CAAC,CAAC;IAChE,MAAM,kBAAkB,GAAG,MAAM,CAAC,kBAAkB,CAAC;IACrD,OAAO,CAAC,OAAO,EAAE,MAAM,EAAE,YAAY,EAAE,eAAe,EAAE,EAAE;QACxD,qBAAqB,CAAC,OAAO,EAAE,MAAM,CAAC,CAAC;QACvC,OAAO,oCAAoC,CAAC,QAAQ,EAAE,EAAE,YAAY,EAAE,eAAe,EAAE,kBAAkB,EAAE,EAAE,OAAO,EAAE,MAAM,CAAC,CAAC,UAAU,CAAC;IAC3I,CAAC,CAAC;AACJ,CAAC;AAED;;;;;;;;;;;;;;;;;;GAkBG;AACH,MAAM,UAAU,wCAAwC,CACtD,KAC+C;IAI/C,gBAAgB,CAAC,KAAK,EAAE,0CAA0C,EAAE,CAAC,sBAAsB,EAAE,qBAAqB,EAAE,SAAS,EAAE,eAAe,EAAE,gBAAgB,EAAE,OAAO,EAAE,cAAc,EAAE,iBAAiB,EAAE,oBAAoB,CAAC,CAAC,CAAC;IACrO,MAAM,OAAO,GACX,mCAAmC,CACjC,KAAK,CAAC,oBAAoB,CAC3B,CAAC;IACJ,MAAM,WAAW,GACf,kCAAkC,CAChC,KAAK,CAAC,mBAAmB,CAC1B,CAAC;IACJ,qBAAqB,CACnB,OAAO,EACP,WAAW,CACZ,CAAC;IAEF,MAAM,QAAQ,GAAG,sCAAsC,CAAC,KAAK,CAAC,CAAC;IAC/D,MAAM,EAAE,OAAO,EAAE,aAAa,EAAE,cAAc,EAAE,OAAO,EAAE,mBAAmB,EAAE,mBAAmB,EAAE,YAAY,EAAE,GAAG,QAAQ,CAAC;IAC7H,MAAM,EAAE,0BAA0B,EAAE,qBAAqB,EAAE,iCAAiC,EAAE,mBAAmB,EAAE,oBAAoB,EACrI,sBAAsB,EAAE,oCAAoC,EAAE,mCAAmC,EAAE,iCAAiC,EACpI,yBAAyB,EAAE,4BAA4B,EAAE,qBAAqB,EAAE,aAAa,EAAE,UAAU,EAAE,GAC3G,oCAAoC,CAAC,QAAQ,EAAE,KAAK,EAAE,OAAO,EAAE,WAAW,CAAC,CAAC;IAE9E,MAAM,WAAW,GAAG;QAClB,GAAG,OAAO,CAAC,WAAW;QACtB,GAAG,4BAA4B;QAC/B,4IAA4I;QAC5I,wGAAwG;QACxG,mEAAmE;QACnE,sFAAsF;QACtF,wIAAwI;QACxI,6FAA6F;QAC7F,oCAAoC;YAClC,CAAC,CAAC,0OAA0O;YAC5O,CAAC,CAAC,uHAAuH;KAC5H,CAAC;IAEF,OAAO,mBAAmB,CACxB;QACE,aAAa,EACX,iDAAiD;QACnD,SAAS,EAAE,OAAO,CAAC,SAAS;QAC5B,cAAc,EACZ,OAAO,CAAC,cAAc;QACxB,aAAa,EACX,OAAO,CAAC,aAAa;QACvB,QAAQ,EAAE,OAAO,CAAC,QAAQ;QAC1B,iBAAiB,EACf,OAAO,CAAC,iBAAiB;QAC3B,OAAO,EAAE,OAAO,CAAC,OAAO;QACxB,4BAA4B,EAC1B,OAAO;aACJ,4BAA4B;QACjC,oBAAoB,EAClB,OAAO,CAAC,oBAAoB;QAC9B,eAAe,EACb,OAAO,CAAC,eAAe;QACzB,aAAa,EACX,4BAA4B;QAC9B,SAAS,EACP,aAAa;QACf,cAAc,EACZ,kCAAkC;QACpC,UAAU,EAAE,UAAU;QACtB,0DAA0D,EACxD,aAAa;QACf,aAAa;QACb,cAAc;QACd,cAAc,EACZ,mBAAmB;QACrB,UAAU,EAAE,KAAK,CAAC,KAAK,CAAC,IAAI;QAC5B,+CAA+C,EAC7C,mBAAmB;QACrB,gBAAgB,EACd,YAAY,CAAC,IAAI;QACnB,0BAA0B;QAC1B,qBAAqB;QACrB,iCAAiC;QACjC,sBAAsB;QACtB,mCAAmC;QACnC,iCAAiC;QACjC,+CAA+C,EAC7C,IAAI;QACN,oCAAoC;QACpC,qBAAqB;QACrB,sDAAsD,EACpD,mBAAmB;QACrB,yDAAyD,EACvD,UAAU;QACZ,gBAAgB,EACd,eAAe;QACjB,6BAA6B,EAC3B,WAAW,CAAC,gBAAgB;QAC9B,8BAA8B,EAC5B,OAAO,CAAC,gBAAgB;QAC1B,wBAAwB,EACtB,WAAW,CAAC,WAAW;QACzB,kBAAkB,EAChB,YAAY,CAAC,WAAW;QAC1B,4BAA4B,EAC1B,KAAK;QACP,iCAAiC,EAC/B,IAAI;QACN,2BAA2B,EACzB,KAAK;QACP,+BAA+B,EAC7B,KAAK;QACP,wBAAwB,EACtB,KAAK;QACP,yBAAyB,EACvB,KAAK;QACP,UAAU,EAAE,KAAK;QACjB,wBAAwB,EAAE,KAAK;QAC/B,2BAA2B,EACzB,KAAK;QACP,iBAAiB,EAAE,KAAK;QACxB,mBAAmB,EAAE,KAAK;QAC1B,iCAAiC,EAC/B,KAAK;QACP,2BAA2B,EACzB,KAAK;QACP,8CAA8C,EAC5C,KAAK;QACP,iBAAiB,EAAE;YACjB,aAAa,EACX,WAAW,CAAC,QAAQ;YACtB,cAAc,EACZ,OAAO,CAAC,QAAQ;YAClB,YAAY,EACV,oBAAoB;YACtB,cAAc,EACZ,OAAO,CAAC,QAAQ;YAClB,kBAAkB,EAChB,yBAAyB;SAC5B;QACD,WAAW;KACZ,EACD,qCAAqC,EACrC,OAAO,EACP,WAAW,CACZ,CAAC;AACJ,CAAC","sourcesContent":["// SPDX-License-Identifier: Apache-2.0\n\n/**\n * Scene-to-sensor optical bridge: from spectral scene radiance to the irradiance that reaches the\n * sensor.\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 here.\n *\n * The calculation applies the declared pupil acceptance and exactly one transmission or\n * field-throughput path. It returns the irradiance before the sensor stack, with conservative\n * evidence.\n *\n * The optical bridge turns outgoing scene radiance (W/m²/sr/nm) into irradiance (W/m²/nm) arriving\n * at the front of the sensor stack.\n * - Circular paraxial acceptance is pi/(4*Nworking²). At finite focus with unit pupil magnification,\n *   the working f-number is N*(1+m).\n * - Spectral transmission and an effective working T-stop are alternative paths; use one.\n * - Field throughput is applied exactly once.\n * - Sensor response and PSF (point-spread function) blur happen later in the pipeline.\n * @see docs/SCENE_SENSOR_QUADRATURE.md for the equations, units, coordinate conventions, known\n * blockers and 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 EvidenceBackedFact,\n  type EvidenceProvenance\n} from \"../core/evidence-provenance.js\";\nimport {\n  parseSpectralWavelengthBasis,\n  type SpectralWavelengthBasis\n} from \"../core/spectral.js\";\nimport { InvalidScientificInputError } from \"../core/validation.js\";\nimport {\n  composeFrontOfLensFilterTransmission,\n  type FrontOfLensFilterProfile\n} from \"./front-of-lens-filter.js\";\nimport type { IlluminationVignetting } from \"./illumination-vignetting.js\";\nimport { calculateThinLensImageDistance } from \"./thin-lens.js\";\nimport {\n  parseSceneRadianceEvaluationRequest,\n  parseSceneRadianceEvaluationResult,\n  type SceneRadianceEvaluationRequest,\n  type SceneRadianceEvaluationResult\n} from \"../schema/scene-radiance.js\";\nimport { requireInputKeys, rejectUnknownFields, rejectUnknownVariantFields, UNCERTAINTY_VARIANT_FIELDS } from \"../core/record-validation.js\";\n\ntype UnknownRecord = Record<string, unknown>;\n\nexport const SCENE_TO_SENSOR_IRRADIANCE_PROFILE_SCHEMA_VERSION =\n  \"0.1.0\" as const;\n\nexport type OpticalBridgeScientificStatus =\n  | \"calibrated\"\n  | \"approximation\";\n\nexport type OpticalBridgeUncertainty =\n  | {\n      kind: \"relative\";\n      fraction: number;\n      basis: string;\n    }\n  | {\n      kind: \"not-quantified\";\n      limitation: string;\n    };\n\nexport interface NumericRange {\n  minimum: number;\n  maximum: number;\n}\n\nexport type OpticalBridgeFocusApplicability =\n  | {\n      kind: \"any\";\n    }\n  | {\n      kind: \"finite-distance-range\";\n      objectDistanceM: NumericRange;\n    };\n\nexport interface SpectralTransmissionSample {\n  wavelengthNanometers: number;\n  linearTransmissionFactor: number;\n}\n\nexport type OpticalTransmissionModel =\n  | {\n      kind: \"spectral-transmission\";\n      scientificStatus:\n        OpticalBridgeScientificStatus;\n      wavelengthBasis:\n        Exclude<\n          SpectralWavelengthBasis,\n          \"unspecified\"\n        >;\n      samples:\n        EvidenceBackedFact<\n          readonly SpectralTransmissionSample[]\n        >;\n      uncertainty:\n        OpticalBridgeUncertainty;\n    }\n  | {\n      kind:\n        \"effective-working-t-stop-approximation\";\n      scientificStatus: \"approximation\";\n      workingTStop:\n        EvidenceBackedFact<number>;\n      wavelengthBasis:\n        Exclude<\n          SpectralWavelengthBasis,\n          \"unspecified\"\n        >;\n      wavelengthRangeNanometers:\n        NumericRange;\n      uncertainty:\n        OpticalBridgeUncertainty;\n      limitation: string;\n    };\n\nexport interface SceneToSensorIrradianceProfile {\n  schemaVersion:\n    typeof SCENE_TO_SENSOR_IRRADIANCE_PROFILE_SCHEMA_VERSION;\n  profileId: string;\n  profileVersion: string;\n  lensProfileId: string;\n  scientificStatus:\n    OpticalBridgeScientificStatus;\n  applicability: {\n    focalLengthMm: NumericRange;\n    nominalFNumber: NumericRange;\n    focus:\n      OpticalBridgeFocusApplicability;\n  };\n  transmission:\n    OpticalTransmissionModel;\n  distortionAreaMappingOwnership:\n    \"not-applied-by-bridge\";\n  psfRedistributionOwnership:\n    \"downstream-normalized-energy-redistribution\";\n  sensorOpticalStackIncluded: false;\n  strayLightIncluded: false;\n  polarizationModeled: false;\n  wavelengthChangingBehaviorModeled: false;\n  volumetricScatteringModeled: false;\n  evidence: readonly EvidenceProvenance[];\n  limitations: readonly string[];\n}\n\nexport type OpticalBridgeFocusContext =\n  | {\n      kind: \"infinity-focus\";\n    }\n  | {\n      kind:\n        \"ideal-symmetric-thin-lens\";\n      objectDistanceM: number;\n      pupilMagnificationAssumption:\n        \"unity\";\n    }\n  | {\n      kind:\n        \"supplied-working-f-number\";\n      focus:\n        | {\n            kind: \"infinity\";\n          }\n        | {\n            kind: \"finite\";\n            objectDistanceM: number;\n          };\n      workingFNumber:\n        EvidenceBackedFact<number>;\n      basis: string;\n    };\n\nexport type OpticalBridgeFieldThroughput =\n  | {\n      kind: \"unity\";\n    }\n  | {\n      kind:\n        \"illumination-vignetting-result\";\n      result:\n        IlluminationVignetting;\n    };\n\nexport interface CalculateSceneRadianceToSensorIrradianceInput {\n  sceneRadianceRequest:\n    SceneRadianceEvaluationRequest;\n  sceneRadianceResult:\n    SceneRadianceEvaluationResult;\n  profile:\n    SceneToSensorIrradianceProfile;\n  focalLengthMm: number;\n  nominalFNumber: number;\n  focus:\n    OpticalBridgeFocusContext;\n  imagePointMm: {\n    x: number;\n    y: number;\n  };\n  fieldThroughput:\n    OpticalBridgeFieldThroughput;\n  frontOfLensFilters?:\n    readonly FrontOfLensFilterProfile[];\n}\n\nexport interface SceneToSensorIrradianceResult {\n  schemaVersion:\n    typeof SCENE_TO_SENSOR_IRRADIANCE_PROFILE_SCHEMA_VERSION;\n  profileId: string;\n  profileVersion: string;\n  lensProfileId: string;\n  sampleId: string;\n  providerProfileId: string;\n  sceneId: string;\n  timeSecondsFromExposureStart: number;\n  wavelengthNanometers: number;\n  wavelengthBasis:\n    Exclude<\n      SpectralWavelengthBasis,\n      \"unspecified\"\n    >;\n  inputQuantity:\n    \"outgoing-spectral-radiance\";\n  inputUnit:\n    \"W/m^2/sr/nm\";\n  outputQuantity:\n    \"sensor-plane-spectral-irradiance\";\n  outputUnit:\n    \"W/m^2/nm\";\n  sceneSpectralRadianceWattsPerSquareMeterSteradianNanometer:\n    number;\n  focalLengthMm: number;\n  nominalFNumber: number;\n  workingFNumber: number;\n  focusModel:\n    OpticalBridgeFocusContext[\"kind\"];\n  paraxialGeometricAcceptanceSolidAngleSteradians:\n    number;\n  transmissionPath:\n    OpticalTransmissionModel[\"kind\"];\n  spectralTransmissionFactor:\n    number | null;\n  effectiveWorkingTStop:\n    number | null;\n  tStopEquivalentTransmissionFactor:\n    number | null;\n  frontOfLensFilterCount: number;\n  frontOfLensFilterTransmissionFactor:\n    number;\n  frontOfLensFilterAttenuationStops:\n    number;\n  frontOfLensFilterTransmissionAppliedExactlyOnce:\n    true;\n  /**\n   * True when a front-of-lens linear polarizer was applied. Scene radiance is unpolarized in this\n   * bridge, so a polarizer transmits (k1 + k2) / 2.\n   */\n  frontOfLensFilterPolarizationModeled:\n    boolean;\n  fieldThroughputFactor: number;\n  effectiveAcceptanceAfterTransmissionAndFieldSteradians:\n    number;\n  sensorPlaneSpectralIrradianceWattsPerSquareMeterNanometer:\n    number;\n  scientificStatus: \"approximation\";\n  sceneRadianceScientificStatus:\n    SceneRadianceEvaluationResult[\"scientificStatus\"];\n  opticalProfileScientificStatus:\n    OpticalBridgeScientificStatus;\n  sceneRadianceUncertainty:\n    SceneRadianceEvaluationResult[\"uncertainty\"];\n  opticalUncertainty:\n    OpticalBridgeUncertainty;\n  radianceAmplificationApplied: false;\n  fieldThroughputAppliedExactlyOnce:\n    true;\n  universalCos4FalloffApplied: false;\n  distortionAreaCorrectionApplied: false;\n  psfRedistributionApplied: false;\n  sensorOpticalStackApplied: false;\n  cfaApplied: false;\n  quantumEfficiencyApplied: false;\n  spectralResponsivityApplied: false;\n  strayLightApplied: false;\n  polarizationModeled: false;\n  wavelengthChangingBehaviorModeled:\n    false;\n  volumetricScatteringModeled: false;\n  calibratedSensorPlaneIrradianceClaimAuthorized:\n    false;\n  componentEvidence: {\n    sceneRadiance:\n      readonly EvidenceProvenance[];\n    opticalProfile:\n      readonly EvidenceProvenance[];\n    transmission:\n      readonly EvidenceProvenance[];\n    workingFNumber:\n      readonly EvidenceProvenance[];\n    frontOfLensFilters:\n      readonly EvidenceProvenance[];\n  };\n  limitations:\n    readonly string[];\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 + \" 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 requireFraction(\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 be from zero through one.\"\n    );\n  }\n  return parsed;\n}\n\nfunction parseRange(\n  value: unknown,\n  path: string\n): NumericRange {\n  const record =\n    requireRecord(value, path, [\"minimum\", \"maximum\"]);\n  const minimum =\n    requirePositiveFinite(\n      record.minimum,\n      path + \".minimum\"\n    );\n  const maximum =\n    requirePositiveFinite(\n      record.maximum,\n      path + \".maximum\"\n    );\n  if (minimum > maximum) {\n    throw new InvalidConfigurationError(\n      path +\n        \".minimum must be less than or equal to maximum.\"\n    );\n  }\n  return {\n    minimum,\n    maximum\n  };\n}\n\nfunction parseScientificStatus(\n  value: unknown,\n  path: string\n): OpticalBridgeScientificStatus {\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): OpticalBridgeUncertainty {\n  const record =\n    requireRecord(value, path, null);\n  rejectUnknownVariantFields(record, \"kind\", UNCERTAINTY_VARIANT_FIELDS, path);\n\n  if (record.kind === \"relative\") {\n    return {\n      kind: \"relative\",\n      fraction:\n        requireFraction(\n          record.fraction,\n          path + \".fraction\"\n        ),\n      basis:\n        requireNonEmptyString(\n          record.basis,\n          path + \".basis\"\n        )\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 + \".limitation\"\n        )\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 = 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 parseReusableNumericEvidence(\n  value: unknown,\n  path: string\n): readonly EvidenceProvenance[] {\n  const evidence =\n    parseEvidenceList(value, path);\n\n  if (\n    evidence.some(\n      (entry) =>\n        entry.reuseStatus ===\n        \"factual-reference-only\"\n    )\n  ) {\n    throw new InvalidConfigurationError(\n      path +\n        \" for numeric optical curve data must be reusable-data or photivra-owned.\"\n    );\n  }\n  return evidence;\n}\n\nfunction parseSpectralSamples(\n  value: unknown,\n  path: string\n): readonly SpectralTransmissionSample[] {\n  if (\n    !Array.isArray(value) ||\n    value.length < 2\n  ) {\n    throw new InvalidConfigurationError(\n      path +\n        \" must contain at least two samples.\"\n    );\n  }\n\n  let previous =\n    Number.NEGATIVE_INFINITY;\n  return value.map(\n    (entry, index) => {\n      const samplePath =\n        path +\n        \"[\" +\n        index +\n        \"]\";\n      const record =\n        requireRecord(\n          entry,\n          samplePath, [\"wavelengthNanometers\", \"linearTransmissionFactor\"]\n        );\n      const wavelengthNanometers =\n        requirePositiveFinite(\n          record.wavelengthNanometers,\n          samplePath +\n            \".wavelengthNanometers\"\n        );\n      if (\n        wavelengthNanometers <=\n        previous\n      ) {\n        throw new InvalidConfigurationError(\n          path +\n            \" wavelengths must be strictly increasing.\"\n        );\n      }\n      previous =\n        wavelengthNanometers;\n      return {\n        wavelengthNanometers,\n        linearTransmissionFactor:\n          requireFraction(\n            record.linearTransmissionFactor,\n            samplePath +\n              \".linearTransmissionFactor\"\n          )\n      };\n    }\n  );\n}\n\nfunction parseTransmission(\n  value: unknown,\n  path: string\n): OpticalTransmissionModel {\n  const record =\n    requireRecord(value, path, null);\n  rejectUnknownVariantFields(record, \"kind\", {\n    // workingTStop and wavelengthRangeNanometers (spectral) and samples (working T-stop) are rejected below with their own messages.\n    \"spectral-transmission\": [\"kind\", \"scientificStatus\", \"wavelengthBasis\", \"samples\", \"uncertainty\", \"workingTStop\", \"wavelengthRangeNanometers\"],\n    \"effective-working-t-stop-approximation\": [\"kind\", \"scientificStatus\", \"workingTStop\", \"wavelengthBasis\", \"wavelengthRangeNanometers\",\n      \"uncertainty\", \"limitation\", \"samples\"]\n  }, path);\n\n  if (\n    record.kind ===\n    \"spectral-transmission\"\n  ) {\n    if (\n      record.workingTStop !==\n      undefined ||\n      record.wavelengthRangeNanometers !==\n      undefined\n    ) {\n      throw new InvalidConfigurationError(\n        path +\n          \" spectral-transmission must not include T-stop-only fields.\"\n      );\n    }\n\n    const status =\n      parseScientificStatus(\n        record.scientificStatus,\n        path +\n          \".scientificStatus\"\n      );\n    const samplesFact =\n      requireRecord(\n        record.samples,\n        path + \".samples\", [\"value\", \"evidence\"]\n      );\n    const uncertainty =\n      parseUncertainty(\n        record.uncertainty,\n        path + \".uncertainty\"\n      );\n\n    if (\n      status === \"calibrated\" &&\n      uncertainty.kind !==\n        \"relative\"\n    ) {\n      throw new InvalidConfigurationError(\n        path +\n          \" calibrated spectral transmission requires quantified relative uncertainty.\"\n      );\n    }\n\n    return {\n      kind:\n        \"spectral-transmission\",\n      scientificStatus: status,\n      wavelengthBasis:\n        parseResolvedBasis(\n          record.wavelengthBasis,\n          path +\n            \".wavelengthBasis\"\n        ),\n      samples: {\n        value:\n          parseSpectralSamples(\n            samplesFact.value,\n            path +\n              \".samples.value\"\n          ),\n        evidence:\n          parseReusableNumericEvidence(\n            samplesFact.evidence,\n            path +\n              \".samples.evidence\"\n          )\n      },\n      uncertainty\n    };\n  }\n\n  if (\n    record.kind ===\n    \"effective-working-t-stop-approximation\"\n  ) {\n    if (\n      record.samples !==\n      undefined\n    ) {\n      throw new InvalidConfigurationError(\n        path +\n          \" T-stop approximation must not include spectral transmission samples.\"\n      );\n    }\n    if (\n      record.scientificStatus !==\n      \"approximation\"\n    ) {\n      throw new InvalidConfigurationError(\n        path +\n          ' T-stop scientificStatus must be \"approximation\".'\n      );\n    }\n    const workingTStopFact =\n      requireRecord(\n        record.workingTStop,\n        path +\n          \".workingTStop\", [\"value\", \"evidence\"]\n      );\n    return {\n      kind:\n        \"effective-working-t-stop-approximation\",\n      scientificStatus:\n        \"approximation\",\n      workingTStop: {\n        value:\n          requirePositiveFinite(\n            workingTStopFact.value,\n            path +\n              \".workingTStop.value\"\n          ),\n        evidence:\n          parseEvidenceList(\n            workingTStopFact.evidence,\n            path +\n              \".workingTStop.evidence\"\n          )\n      },\n      wavelengthBasis:\n        parseResolvedBasis(\n          record.wavelengthBasis,\n          path +\n            \".wavelengthBasis\"\n        ),\n      wavelengthRangeNanometers:\n        parseRange(\n          record.wavelengthRangeNanometers,\n          path +\n            \".wavelengthRangeNanometers\"\n        ),\n      uncertainty:\n        parseUncertainty(\n          record.uncertainty,\n          path + \".uncertainty\"\n        ),\n      limitation:\n        requireNonEmptyString(\n          record.limitation,\n          path +\n            \".limitation\"\n        )\n    };\n  }\n\n  throw new InvalidConfigurationError(\n    path + \".kind is invalid.\"\n  );\n}\n\nfunction parseResolvedBasis(\n  value: unknown,\n  path: string\n): Exclude<\n  SpectralWavelengthBasis,\n  \"unspecified\"\n> {\n  const parsed =\n    parseSpectralWavelengthBasis(\n      value,\n      path\n    );\n  if (parsed === \"unspecified\") {\n    throw new InvalidConfigurationError(\n      path +\n        ' must resolve to \"air\" or \"vacuum\".'\n    );\n  }\n  return parsed;\n}\n\nfunction parseFocusApplicability(\n  value: unknown,\n  path: string\n): OpticalBridgeFocusApplicability {\n  const record =\n    requireRecord(value, path, null);\n  rejectUnknownVariantFields(record, \"kind\", { any: [\"kind\"], \"finite-distance-range\": [\"kind\", \"objectDistanceM\"] }, path);\n\n  if (record.kind === \"any\") {\n    return {\n      kind: \"any\"\n    };\n  }\n  if (\n    record.kind ===\n    \"finite-distance-range\"\n  ) {\n    return {\n      kind:\n        \"finite-distance-range\",\n      objectDistanceM:\n        parseRange(\n          record.objectDistanceM,\n          path +\n            \".objectDistanceM\"\n        )\n    };\n  }\n  throw new InvalidConfigurationError(\n    path + \".kind is invalid.\"\n  );\n}\n\n/**\n * Validates an untrusted scene-to-sensor profile and returns the normalized typed contract. Unknown\n * enum values, missing required fields and incompatible scientific data are rejected here.\n *\n * The optical bridge turns outgoing scene radiance (W/m²/sr/nm) into irradiance (W/m²/nm) arriving\n * at the front of the sensor stack.\n * - Circular paraxial acceptance is pi/(4*Nworking²). At finite focus with unit pupil magnification,\n *   the working f-number is N*(1+m).\n * - Spectral transmission and an effective working T-stop are alternative paths; use one.\n * - Field throughput is applied exactly once.\n * - Sensor response and PSF (point-spread function) blur happen later in the pipeline.\n * @param value - unknown. Untrusted input; TypeScript types alone are not validation.\n * @returns SceneToSensorIrradianceProfile. 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 parseSceneToSensorIrradianceProfile(\n  value: unknown\n): SceneToSensorIrradianceProfile {\n  const record =\n    requireRecord(\n      value,\n      \"sceneToSensorIrradianceProfile\", [\"schemaVersion\", \"applicability\", \"scientificStatus\", \"transmission\", \"distortionAreaMappingOwnership\", \"psfRedistributionOwnership\", \"sensorOpticalStackIncluded\", \"strayLightIncluded\", \"polarizationModeled\", \"wavelengthChangingBehaviorModeled\", \"volumetricScatteringModeled\", \"profileId\", \"profileVersion\", \"lensProfileId\", \"evidence\", \"limitations\"]\n    );\n  if (\n    record.schemaVersion !==\n    SCENE_TO_SENSOR_IRRADIANCE_PROFILE_SCHEMA_VERSION\n  ) {\n    throw new InvalidConfigurationError(\n      'sceneToSensorIrradianceProfile.schemaVersion must be \"' +\n        SCENE_TO_SENSOR_IRRADIANCE_PROFILE_SCHEMA_VERSION +\n        '\".'\n    );\n  }\n\n  const applicability =\n    requireRecord(\n      record.applicability,\n      \"sceneToSensorIrradianceProfile.applicability\", [\"focalLengthMm\", \"nominalFNumber\", \"focus\"]\n    );\n  const status =\n    parseScientificStatus(\n      record.scientificStatus,\n      \"sceneToSensorIrradianceProfile.scientificStatus\"\n    );\n  const transmission =\n    parseTransmission(\n      record.transmission,\n      \"sceneToSensorIrradianceProfile.transmission\"\n    );\n\n  if (\n    status === \"calibrated\" &&\n    transmission.scientificStatus !==\n      \"calibrated\"\n  ) {\n    throw new InvalidConfigurationError(\n      \"A calibrated optical bridge profile requires calibrated spectral transmission.\"\n    );\n  }\n\n  if (\n    record.distortionAreaMappingOwnership !==\n      \"not-applied-by-bridge\" ||\n    record.psfRedistributionOwnership !==\n      \"downstream-normalized-energy-redistribution\" ||\n    record.sensorOpticalStackIncluded !==\n      false ||\n    record.strayLightIncluded !==\n      false ||\n    record.polarizationModeled !==\n      false ||\n    record.wavelengthChangingBehaviorModeled !==\n      false ||\n    record.volumetricScatteringModeled !==\n      false\n  ) {\n    throw new InvalidConfigurationError(\n      \"sceneToSensorIrradianceProfile ownership/completeness flags are invalid for schema 0.1.0.\"\n    );\n  }\n\n  return {\n    schemaVersion:\n      SCENE_TO_SENSOR_IRRADIANCE_PROFILE_SCHEMA_VERSION,\n    profileId:\n      requireNonEmptyString(\n        record.profileId,\n        \"sceneToSensorIrradianceProfile.profileId\"\n      ),\n    profileVersion:\n      requireNonEmptyString(\n        record.profileVersion,\n        \"sceneToSensorIrradianceProfile.profileVersion\"\n      ),\n    lensProfileId:\n      requireNonEmptyString(\n        record.lensProfileId,\n        \"sceneToSensorIrradianceProfile.lensProfileId\"\n      ),\n    scientificStatus: status,\n    applicability: {\n      focalLengthMm:\n        parseRange(\n          applicability.focalLengthMm,\n          \"sceneToSensorIrradianceProfile.applicability.focalLengthMm\"\n        ),\n      nominalFNumber:\n        parseRange(\n          applicability.nominalFNumber,\n          \"sceneToSensorIrradianceProfile.applicability.nominalFNumber\"\n        ),\n      focus:\n        parseFocusApplicability(\n          applicability.focus,\n          \"sceneToSensorIrradianceProfile.applicability.focus\"\n        )\n    },\n    transmission,\n    distortionAreaMappingOwnership:\n      \"not-applied-by-bridge\",\n    psfRedistributionOwnership:\n      \"downstream-normalized-energy-redistribution\",\n    sensorOpticalStackIncluded:\n      false,\n    strayLightIncluded: false,\n    polarizationModeled: false,\n    wavelengthChangingBehaviorModeled:\n      false,\n    volumetricScatteringModeled:\n      false,\n    evidence:\n      parseEvidenceList(\n        record.evidence,\n        \"sceneToSensorIrradianceProfile.evidence\"\n      ),\n    limitations:\n      parseLimitations(\n        record.limitations,\n        \"sceneToSensorIrradianceProfile.limitations\"\n      )\n  };\n}\n\nfunction ensureWithin(\n  value: number,\n  range: NumericRange,\n  path: string\n): void {\n  if (\n    value < range.minimum ||\n    value > range.maximum\n  ) {\n    throw new InvalidScientificInputError(\n      path +\n        \" lies outside the optical profile applicability range.\"\n    );\n  }\n}\n\n/**\n * Internal (2.6; not exported from the package): N_w / N = 1 + m for an ideal symmetric thin lens with unity pupil\n * magnification focused at `objectDistanceM`, m = v/s from `calculateThinLensImageDistance`. The irradiance bridge's\n * `ideal-symmetric-thin-lens` focus uses exactly this factor, and the framing reports derive their close-up light\n * factor 1/(1 + m)² from it, so the two cannot drift apart.\n */\nexport function thinLensWorkingFNumberFactor(focalLengthMm: number, objectDistanceM: number): number {\n  return 1 + calculateThinLensImageDistance({ focalLengthMm, objectDistanceM }).value.magnification;\n}\n\nfunction workingFNumber(\n  focus:\n    OpticalBridgeFocusContext,\n  focalLengthMm: number,\n  nominalFNumber: number\n): {\n  value: number;\n  evidence:\n    readonly EvidenceProvenance[];\n} {\n  if (\n    focus.kind === \"infinity-focus\"\n  ) {\n    return {\n      value: nominalFNumber,\n      evidence: []\n    };\n  }\n\n  if (\n    focus.kind ===\n    \"ideal-symmetric-thin-lens\"\n  ) {\n    if (\n      focus\n        .pupilMagnificationAssumption !==\n      \"unity\"\n    ) {\n      throw new InvalidScientificInputError(\n        'Ideal symmetric thin-lens focus requires pupilMagnificationAssumption \"unity\".'\n      );\n    }\n    return {\n      value:\n        nominalFNumber *\n        thinLensWorkingFNumberFactor(focalLengthMm, focus.objectDistanceM),\n      evidence: []\n    };\n  }\n\n  if (\n    focus.kind ===\n    \"supplied-working-f-number\"\n  ) {\n    const value =\n      requirePositiveFiniteScientific(\n        focus.workingFNumber.value,\n        \"focus.workingFNumber.value\"\n      );\n    return {\n      value,\n      evidence:\n        parseEvidenceList(\n          focus.workingFNumber.evidence,\n          \"focus.workingFNumber.evidence\"\n        )\n    };\n  }\n\n  throw new InvalidConfigurationError(\n    \"focus.kind is invalid.\"\n  );\n}\n\nfunction requirePositiveFiniteScientific(\n  value: unknown,\n  path: string\n): number {\n  if (\n    typeof value !== \"number\" ||\n    !Number.isFinite(value) ||\n    value <= 0\n  ) {\n    throw new InvalidScientificInputError(\n      path +\n        \" must be finite and greater than zero.\"\n    );\n  }\n  return value;\n}\n\nfunction finiteFocusDistance(\n  focus:\n    OpticalBridgeFocusContext\n): number | null {\n  if (\n    focus.kind ===\n    \"ideal-symmetric-thin-lens\"\n  ) {\n    return requirePositiveFiniteScientific(\n      focus.objectDistanceM,\n      \"focus.objectDistanceM\"\n    );\n  }\n  if (\n    focus.kind ===\n      \"supplied-working-f-number\" &&\n    focus.focus.kind === \"finite\"\n  ) {\n    return requirePositiveFiniteScientific(\n      focus.focus.objectDistanceM,\n      \"focus.focus.objectDistanceM\"\n    );\n  }\n  return null;\n}\n\nfunction interpolateTransmission(\n  model:\n    Extract<\n      OpticalTransmissionModel,\n      { kind: \"spectral-transmission\" }\n    >,\n  wavelengthNanometers: number\n): number {\n  const samples =\n    model.samples.value;\n  const first =\n    samples[0]!;\n  const last =\n    samples[\n      samples.length - 1\n    ]!;\n\n  if (\n    wavelengthNanometers <\n      first.wavelengthNanometers ||\n    wavelengthNanometers >\n      last.wavelengthNanometers\n  ) {\n    throw new InvalidScientificInputError(\n      \"Scene wavelength lies outside the spectral transmission sample range.\"\n    );\n  }\n\n  if (\n    wavelengthNanometers ===\n    first.wavelengthNanometers\n  ) {\n    return first\n      .linearTransmissionFactor;\n  }\n\n  for (\n    let index = 1;\n    index < samples.length;\n    index += 1\n  ) {\n    const right =\n      samples[index]!;\n    if (\n      wavelengthNanometers <=\n      right.wavelengthNanometers\n    ) {\n      const left =\n        samples[index - 1]!;\n      if (\n        wavelengthNanometers ===\n        right.wavelengthNanometers\n      ) {\n        return right\n          .linearTransmissionFactor;\n      }\n      const phase =\n        (wavelengthNanometers -\n          left.wavelengthNanometers) /\n        (right.wavelengthNanometers -\n          left.wavelengthNanometers);\n      return (\n        left\n          .linearTransmissionFactor +\n        phase *\n          (right\n              .linearTransmissionFactor -\n            left\n              .linearTransmissionFactor)\n      );\n    }\n  }\n\n  return last\n    .linearTransmissionFactor;\n}\n\nfunction fieldFactor(\n  field:\n    OpticalBridgeFieldThroughput,\n  imagePointMm: {\n    x: number;\n    y: number;\n  }\n): number {\n  if (field.kind === \"unity\") {\n    return 1;\n  }\n  if (\n    field.kind ===\n    \"illumination-vignetting-result\"\n  ) {\n    const result =\n      field.result;\n    const scale = Math.max(\n      1,\n      Math.abs(imagePointMm.x),\n      Math.abs(imagePointMm.y),\n      Math.abs(result.imagePointMm.x),\n      Math.abs(result.imagePointMm.y)\n    );\n    const tolerance =\n      Number.EPSILON *\n      64 *\n      scale;\n    if (\n      Math.abs(\n        result.imagePointMm.x -\n          imagePointMm.x\n      ) > tolerance ||\n      Math.abs(\n        result.imagePointMm.y -\n          imagePointMm.y\n      ) > tolerance\n    ) {\n      throw new InvalidScientificInputError(\n        \"Illumination-vignetting result imagePointMm must match the optical bridge imagePointMm.\"\n      );\n    }\n    const factor =\n      result.linearThroughputFactor;\n    if (\n      !Number.isFinite(factor) ||\n      factor <= 0 ||\n      factor > 1\n    ) {\n      throw new InvalidScientificInputError(\n        \"Field throughput factor must be finite and in (0, 1].\"\n      );\n    }\n    return factor;\n  }\n  throw new InvalidConfigurationError(\n    \"fieldThroughput.kind is invalid.\"\n  );\n}\n\nfunction validateSceneIdentity(\n  request:\n    SceneRadianceEvaluationRequest,\n  result:\n    SceneRadianceEvaluationResult\n): void {\n  if (\n    request.sampleId !==\n      result.sampleId ||\n    request.providerProfileId !==\n      result.providerProfileId ||\n    request.sceneId !==\n      result.sceneId ||\n    request.wavelengthNanometers !==\n      result.wavelengthNanometers ||\n    request.wavelengthBasis !==\n      result.wavelengthBasis\n  ) {\n    throw new InvalidScientificInputError(\n      \"Scene-radiance request/result identity must match exactly.\"\n    );\n  }\n}\n\n/**\n * The request-independent part of {@link calculateSceneRadianceToSensorIrradiance}: the profile, focal length, f-number,\n * applicability, focus and working f-number checks and the paraxial acceptance, in the public function's order.\n */\ninterface PreparedSceneRadianceToSensorIrradiance {\n  profile: ReturnType<typeof parseSceneToSensorIrradianceProfile>;\n  focalLengthMm: number;\n  nominalFNumber: number;\n  working: ReturnType<typeof workingFNumber>;\n  workingFNumberValue: number;\n  geometricAcceptance: number;\n  transmission: ReturnType<typeof parseSceneToSensorIrradianceProfile>[\"transmission\"];\n}\n\nfunction prepareSceneRadianceToSensorIrradiance(\n  input: Pick<CalculateSceneRadianceToSensorIrradianceInput, \"profile\" | \"focalLengthMm\" | \"nominalFNumber\" | \"focus\">\n): PreparedSceneRadianceToSensorIrradiance {\n  const profile =\n    parseSceneToSensorIrradianceProfile(\n      input.profile\n    );\n  const focalLengthMm =\n    requirePositiveFiniteScientific(\n      input.focalLengthMm,\n      \"focalLengthMm\"\n    );\n  const nominalFNumber =\n    requirePositiveFiniteScientific(\n      input.nominalFNumber,\n      \"nominalFNumber\"\n    );\n\n  ensureWithin(\n    focalLengthMm,\n    profile.applicability\n      .focalLengthMm,\n    \"focalLengthMm\"\n  );\n  ensureWithin(\n    nominalFNumber,\n    profile.applicability\n      .nominalFNumber,\n    \"nominalFNumber\"\n  );\n\n  const focusDistanceM =\n    finiteFocusDistance(input.focus);\n  if (\n    profile.applicability.focus.kind ===\n      \"finite-distance-range\"\n  ) {\n    if (focusDistanceM === null) {\n      throw new InvalidScientificInputError(\n        \"Optical profile requires a finite focus distance.\"\n      );\n    }\n    ensureWithin(\n      focusDistanceM,\n      profile.applicability\n        .focus.objectDistanceM,\n      \"focus object distance\"\n    );\n  }\n\n  const working =\n    workingFNumber(\n      input.focus,\n      focalLengthMm,\n      nominalFNumber\n    );\n  const workingFNumberValue =\n    working.value;\n  const geometricAcceptance =\n    Math.PI /\n    (4 *\n      workingFNumberValue *\n      workingFNumberValue);\n\n  if (\n    !Number.isFinite(\n      geometricAcceptance\n    ) ||\n    geometricAcceptance <= 0\n  ) {\n    throw new InvalidScientificInputError(\n      \"Paraxial geometric acceptance must remain finite and greater than zero.\"\n    );\n  }\n\n  const transmission =\n    profile.transmission;\n  return { profile, focalLengthMm, nominalFNumber, working, workingFNumberValue, geometricAcceptance, transmission };\n}\n\ninterface AppliedSceneRadianceToSensorIrradiance {\n  spectralTransmissionFactor: number | null;\n  effectiveWorkingTStop: number | null;\n  tStopEquivalentTransmissionFactor: number | null;\n  effectiveAcceptance: number;\n  transmissionEvidence: readonly EvidenceProvenance[];\n  frontOfLensFilterCount: number;\n  frontOfLensFilterPolarizationModeled: boolean;\n  frontOfLensFilterTransmissionFactor: number;\n  frontOfLensFilterAttenuationStops: number;\n  frontOfLensFilterEvidence: readonly EvidenceProvenance[];\n  frontOfLensFilterLimitations: readonly string[];\n  fieldThroughputFactor: number;\n  sceneRadiance: number;\n  irradiance: number;\n}\n\n/**\n * The per-query part of {@link calculateSceneRadianceToSensorIrradiance}, after the request and result are parsed and their\n * identities matched: the wavelength-basis check, the transmission or T-stop path, the front-of-lens filters, the field\n * throughput and the irradiance, with the same operations in the same order.\n */\nfunction applySceneRadianceToSensorIrradiance(\n  prepared: PreparedSceneRadianceToSensorIrradiance,\n  input: Pick<CalculateSceneRadianceToSensorIrradianceInput, \"imagePointMm\" | \"fieldThroughput\"> &\n    { frontOfLensFilters?: readonly FrontOfLensFilterProfile[] | undefined },\n  request: SceneRadianceEvaluationRequest,\n  sceneResult: SceneRadianceEvaluationResult\n): AppliedSceneRadianceToSensorIrradiance {\n  const { workingFNumberValue, geometricAcceptance, transmission } = prepared;\n  if (\n    request.wavelengthBasis !==\n    transmission.wavelengthBasis\n  ) {\n    throw new InvalidScientificInputError(\n      \"Scene-radiance and optical-transmission wavelength bases must match; implicit air/vacuum conversion is not performed.\"\n    );\n  }\n\n  let spectralTransmissionFactor:\n    number | null = null;\n  let effectiveWorkingTStop:\n    number | null = null;\n  let tStopEquivalentTransmissionFactor:\n    number | null = null;\n  let effectiveAcceptance =\n    geometricAcceptance;\n  let transmissionEvidence:\n    readonly EvidenceProvenance[];\n\n  if (\n    transmission.kind ===\n    \"spectral-transmission\"\n  ) {\n    spectralTransmissionFactor =\n      interpolateTransmission(\n        transmission,\n        request.wavelengthNanometers\n      );\n    effectiveAcceptance *=\n      spectralTransmissionFactor;\n    transmissionEvidence =\n      transmission.samples.evidence;\n  } else {\n    ensureWithin(\n      request.wavelengthNanometers,\n      transmission\n        .wavelengthRangeNanometers,\n      \"scene wavelength\"\n    );\n    effectiveWorkingTStop =\n      transmission.workingTStop.value;\n    if (\n      effectiveWorkingTStop <\n      workingFNumberValue\n    ) {\n      throw new InvalidScientificInputError(\n        \"effective working T-stop must not be smaller than the geometric working f-number for a passive optical approximation.\"\n      );\n    }\n    tStopEquivalentTransmissionFactor =\n      (workingFNumberValue /\n        effectiveWorkingTStop) **\n      2;\n    effectiveAcceptance =\n      Math.PI /\n      (4 *\n        effectiveWorkingTStop *\n        effectiveWorkingTStop);\n    transmissionEvidence =\n      transmission\n        .workingTStop.evidence;\n  }\n\n  let frontOfLensFilterCount = 0;\n  let frontOfLensFilterPolarizationModeled = false;\n  let frontOfLensFilterTransmissionFactor = 1;\n  let frontOfLensFilterAttenuationStops = 0;\n  let frontOfLensFilterEvidence:\n    readonly EvidenceProvenance[] = [];\n  let frontOfLensFilterLimitations:\n    readonly string[] = [];\n\n  if (\n    input.frontOfLensFilters !==\n    undefined\n  ) {\n    const composed =\n      composeFrontOfLensFilterTransmission({\n        filters:\n          input.frontOfLensFilters,\n        wavelengthNanometers:\n          request.wavelengthNanometers,\n        wavelengthBasis:\n          request.wavelengthBasis\n      });\n    frontOfLensFilterCount =\n      composed.appliedFilterCount;\n    frontOfLensFilterPolarizationModeled =\n      composed.polarizationModeled;\n    frontOfLensFilterTransmissionFactor =\n      composed\n        .combinedLinearTransmissionFactor;\n    frontOfLensFilterAttenuationStops =\n      composed\n        .combinedAttenuationStops;\n    frontOfLensFilterEvidence =\n      composed.components.flatMap(\n        (component) =>\n          component.evidence\n      );\n    frontOfLensFilterLimitations =\n      composed.components.flatMap(\n        (component) =>\n          component.limitations\n      );\n    effectiveAcceptance *=\n      frontOfLensFilterTransmissionFactor;\n  }\n\n  const fieldThroughputFactor =\n    fieldFactor(\n      input.fieldThroughput,\n      input.imagePointMm\n    );\n  effectiveAcceptance *=\n    fieldThroughputFactor;\n\n  const sceneRadiance =\n    sceneResult\n      .spectralRadianceWattsPerSquareMeterSteradianNanometer;\n  const irradiance =\n    sceneRadiance *\n    effectiveAcceptance;\n\n  if (\n    !Number.isFinite(irradiance) ||\n    irradiance < 0\n  ) {\n    throw new InvalidScientificInputError(\n      \"Sensor-plane spectral irradiance must remain finite and non-negative.\"\n    );\n  }\n\n  return { spectralTransmissionFactor, effectiveWorkingTStop, tStopEquivalentTransmissionFactor, effectiveAcceptance, transmissionEvidence,\n    frontOfLensFilterCount, frontOfLensFilterPolarizationModeled, frontOfLensFilterTransmissionFactor, frontOfLensFilterAttenuationStops,\n    frontOfLensFilterEvidence, frontOfLensFilterLimitations, fieldThroughputFactor, sceneRadiance, irradiance };\n}\n\n/**\n * Internal (#453; not exported from the package): the sensor-plane spectral irradiance that\n * {@link calculateSceneRadianceToSensorIrradiance} returns, for repeated queries of one optical state. This call runs the\n * request-independent checks once (the public function's own code, in its order); each call of the returned function then\n * runs the public function's per-query steps on an engine-built request and a result already parsed by the caller: the identity\n * match, the wavelength-basis check, the transmission or T-stop path, the front-of-lens filters, the field throughput and\n * the finite, non-negative irradiance, with the same operations, order and errors. It skips only the re-validation of the\n * unchanged optical inputs, the re-parse of the request and result, and the result envelope.\n */\nexport function prepareSceneRadianceToSensorIrradianceValue(\n  optics: Pick<CalculateSceneRadianceToSensorIrradianceInput, \"profile\" | \"focalLengthMm\" | \"nominalFNumber\" | \"focus\" | \"frontOfLensFilters\">\n): (request: SceneRadianceEvaluationRequest, result: SceneRadianceEvaluationResult, imagePointMm: { x: number; y: number },\n    fieldThroughput: OpticalBridgeFieldThroughput) => number {\n  const prepared = prepareSceneRadianceToSensorIrradiance(optics);\n  const frontOfLensFilters = optics.frontOfLensFilters;\n  return (request, result, imagePointMm, fieldThroughput) => {\n    validateSceneIdentity(request, result);\n    return applySceneRadianceToSensorIrradiance(prepared, { imagePointMm, fieldThroughput, frontOfLensFilters }, request, result).irradiance;\n  };\n}\n\n/**\n * Applies the declared pupil acceptance and exactly one transmission or field-throughput path to\n * spectral scene radiance. Returns the irradiance before the sensor stack, with conservative\n * evidence.\n *\n * The optical bridge turns outgoing scene radiance (W/m²/sr/nm) into irradiance (W/m²/nm) arriving\n * at the front of the sensor stack.\n * - Circular paraxial acceptance is pi/(4*Nworking²). At finite focus with unit pupil magnification,\n *   the working f-number is N*(1+m).\n * - Spectral transmission and an effective working T-stop are alternative paths; use one.\n * - Field throughput is applied exactly once.\n * - Sensor response and PSF (point-spread function) blur happen later in the pipeline.\n * @param input - CalculateSceneRadianceToSensorIrradianceInput. See that type for coordinates, units and profile binding.\n * @returns CalculationResult<\n  SceneToSensorIrradianceResult\n>. 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 calculateSceneRadianceToSensorIrradiance(\n  input:\n    CalculateSceneRadianceToSensorIrradianceInput\n): CalculationResult<\n  SceneToSensorIrradianceResult\n> {\n  requireInputKeys(input, \"calculateSceneRadianceToSensorIrradiance\", [\"sceneRadianceRequest\", \"sceneRadianceResult\", \"profile\", \"focalLengthMm\", \"nominalFNumber\", \"focus\", \"imagePointMm\", \"fieldThroughput\", \"frontOfLensFilters\"]);\n  const request =\n    parseSceneRadianceEvaluationRequest(\n      input.sceneRadianceRequest\n    );\n  const sceneResult =\n    parseSceneRadianceEvaluationResult(\n      input.sceneRadianceResult\n    );\n  validateSceneIdentity(\n    request,\n    sceneResult\n  );\n\n  const prepared = prepareSceneRadianceToSensorIrradiance(input);\n  const { profile, focalLengthMm, nominalFNumber, working, workingFNumberValue, geometricAcceptance, transmission } = prepared;\n  const { spectralTransmissionFactor, effectiveWorkingTStop, tStopEquivalentTransmissionFactor, effectiveAcceptance, transmissionEvidence,\n    frontOfLensFilterCount, frontOfLensFilterPolarizationModeled, frontOfLensFilterTransmissionFactor, frontOfLensFilterAttenuationStops,\n    frontOfLensFilterEvidence, frontOfLensFilterLimitations, fieldThroughputFactor, sceneRadiance, irradiance } =\n    applySceneRadianceToSensorIrradiance(prepared, input, request, sceneResult);\n\n  const limitations = [\n    ...profile.limitations,\n    ...frontOfLensFilterLimitations,\n    \"Front-of-lens filter transmission is applied outside the lens throughput path exactly once and does not change geometric working f-number.\",\n    \"Paraxial circular-pupil acceptance uses pi/(4*N_working^2); no universal cos^4 field falloff is added.\",\n    \"Field throughput is supplied separately and applied exactly once.\",\n    \"Geometric distortion area-density/Jacobian correction is not applied by this bridge.\",\n    \"PSF/diffraction energy redistribution is downstream and must preserve normalized energy apart from separately owned throughput losses.\",\n    \"Sensor optical stack, microlens, CFA, QE/responsivity and charge conversion are downstream.\",\n    frontOfLensFilterPolarizationModeled\n      ? \"Flare, ghosting, veiling glare, scene-radiance polarization (scene light is treated as unpolarized; the front-of-lens linear polarizer transmits (k1 + k2) / 2), volumetric scattering and wavelength-changing behavior are not modeled.\"\n      : \"Flare, ghosting, veiling glare, polarization, volumetric scattering and wavelength-changing behavior are not modeled.\"\n  ];\n\n  return approximationResult(\n    {\n      schemaVersion:\n        SCENE_TO_SENSOR_IRRADIANCE_PROFILE_SCHEMA_VERSION,\n      profileId: profile.profileId,\n      profileVersion:\n        profile.profileVersion,\n      lensProfileId:\n        profile.lensProfileId,\n      sampleId: request.sampleId,\n      providerProfileId:\n        request.providerProfileId,\n      sceneId: request.sceneId,\n      timeSecondsFromExposureStart:\n        request\n          .timeSecondsFromExposureStart,\n      wavelengthNanometers:\n        request.wavelengthNanometers,\n      wavelengthBasis:\n        request.wavelengthBasis,\n      inputQuantity:\n        \"outgoing-spectral-radiance\",\n      inputUnit:\n        \"W/m^2/sr/nm\",\n      outputQuantity:\n        \"sensor-plane-spectral-irradiance\",\n      outputUnit: \"W/m^2/nm\",\n      sceneSpectralRadianceWattsPerSquareMeterSteradianNanometer:\n        sceneRadiance,\n      focalLengthMm,\n      nominalFNumber,\n      workingFNumber:\n        workingFNumberValue,\n      focusModel: input.focus.kind,\n      paraxialGeometricAcceptanceSolidAngleSteradians:\n        geometricAcceptance,\n      transmissionPath:\n        transmission.kind,\n      spectralTransmissionFactor,\n      effectiveWorkingTStop,\n      tStopEquivalentTransmissionFactor,\n      frontOfLensFilterCount,\n      frontOfLensFilterTransmissionFactor,\n      frontOfLensFilterAttenuationStops,\n      frontOfLensFilterTransmissionAppliedExactlyOnce:\n        true,\n      frontOfLensFilterPolarizationModeled,\n      fieldThroughputFactor,\n      effectiveAcceptanceAfterTransmissionAndFieldSteradians:\n        effectiveAcceptance,\n      sensorPlaneSpectralIrradianceWattsPerSquareMeterNanometer:\n        irradiance,\n      scientificStatus:\n        \"approximation\",\n      sceneRadianceScientificStatus:\n        sceneResult.scientificStatus,\n      opticalProfileScientificStatus:\n        profile.scientificStatus,\n      sceneRadianceUncertainty:\n        sceneResult.uncertainty,\n      opticalUncertainty:\n        transmission.uncertainty,\n      radianceAmplificationApplied:\n        false,\n      fieldThroughputAppliedExactlyOnce:\n        true,\n      universalCos4FalloffApplied:\n        false,\n      distortionAreaCorrectionApplied:\n        false,\n      psfRedistributionApplied:\n        false,\n      sensorOpticalStackApplied:\n        false,\n      cfaApplied: false,\n      quantumEfficiencyApplied: false,\n      spectralResponsivityApplied:\n        false,\n      strayLightApplied: false,\n      polarizationModeled: false,\n      wavelengthChangingBehaviorModeled:\n        false,\n      volumetricScatteringModeled:\n        false,\n      calibratedSensorPlaneIrradianceClaimAuthorized:\n        false,\n      componentEvidence: {\n        sceneRadiance:\n          sceneResult.evidence,\n        opticalProfile:\n          profile.evidence,\n        transmission:\n          transmissionEvidence,\n        workingFNumber:\n          working.evidence,\n        frontOfLensFilters:\n          frontOfLensFilterEvidence\n      },\n      limitations\n    },\n    \"scene-radiance-to-sensor-irradiance\",\n    \"1.1.0\",\n    limitations\n  );\n}\n"]}

package/dist/optics/sensor-aperture-rays.js

@@ -9,6 +9,7 @@ import { InvalidScientificInputError, requirePositiveFinite } from "../core/vali
 import { calculateFocusPlaneImageDistance, parseFocusPlane } from "./focus-state.js";
 import { applyPhysicalCameraRotation, calculateInverseCameraRotationImageMapping } from "../motion/camera-rotation.js";
 import { requireInputKeys } from "../core/record-validation.js";
+import { hypot3 } from "../core/portable-math.js";
 /**
  * Polar midpoint rule with equal-area cells. Weights are normalized as an average, so they never
  * apply the f-number a second time.
@@ -59,7 +60,7 @@ export function sensorApertureRayValues(owned, pupil, focus) {
         const x = owned.sourcePointNativeSensorMm.x / v, y = -owned.sourcePointNativeSensorMm.y / v;
         const direction = focus.kind === "finite" ? { x: x * focus.distanceM - origin.x, y: y * focus.distanceM - origin.y, z: focus.distanceM }
             : { x, y, z: 1 };
-        const length = Math.hypot(direction.x, direction.y, direction.z);
+        const length = hypot3(direction.x, direction.y, direction.z);
         const unit = { x: direction.x / length, y: direction.y / length, z: direction.z / length };
         // Same physical right-hand camera rotation as the inverse mapping above.
         const rotatedOrigin = applyPhysicalCameraRotation(origin, displacement, "camera-to-reference"), rotatedDirection = applyPhysicalCameraRotation(unit, displacement, "camera-to-reference");