merit_evaluate — OPTICS lensmerit op

• Data kinds: table → table

• Call: import fullseye as fs; fs.ledger.merit_evaluate(system=None, operands=None, fields=None, wavelengths=None, field_weights=None, wavelength_weights=None) -> 'dict' (implementation को सीधे call करना हो तो import lensmerit; lensmerit.merit_evaluate(system=None, operands=None, fields=None, wavelengths=None, field_weights=None, wavelength_weights=None) -> 'dict', registry से लेना हो तो opsoptics.get("merit_evaluate"))

• GPU: नहीं (केवल CPU reference implementation)

उपयोग

operand table को prescription पर लागू करके हर पंक्ति का मान, लक्ष्य, भार, योगदान और उसका हिस्सा लौटाता है (table)।

> नीचे का विस्तृत विवरण मूल text में है — सारांश और headings का अनुवाद हो चुका है।

*operands* is a list of rows {"op": "EFFL", "target": 50.0, "weight": 1.0, ...arguments} (see merit_operands_list), or the dict merit_default returns. Each row contributes (weight (value - target))^2; a one-sided row (ETGT ...) contributes 0 while the bound holds. A group row (TRAN, OPDW) contributes the sum over its rays. *fields* (default: on axis, or 0 / 0.7 / 1.0 of the system field) and *wavelengths* (default the system wavelength; the primary is the one equal to it) are what the rows' field / wavelength indices point into; a row with "all" is expanded to one row per field and wavelength, its weight times sqrt(field_weight) and sqrt(wavelength_weight) (default 1).

Returns {"table", "merit", "n_operands", "n_residuals", "dominant", "fields", "wavelengths"}. Each table row has index op label value target weight sense margin contribution percent satisfied n_residuals source args (margin is the slack of a one-sided row, negative when violated). Fail-closed: an unknown operand, a missing or unknown argument, an index out of range, a lost ray, a value that is not finite and a table of zero weights are all ValueError.

Example::

t = merit_evaluate(lens, [{"op": "EFFL", "target": 50.0}, {"op": "ETGT", "surf": 0, "target": 1.0}])

t["merit"], t["dominant"]

पूरी family का input contract (fail-closed)

optics के सभी ops गणना से पहले input की जाँच करते हैं (कुछ भी चुपचाप पास नहीं होता):

• Unit argument के नाम में ही लिखी होती है — _mm / _um / _deg / _mrad। mm और µm की गड़बड़ी crash नहीं करती, बल्कि "सही दिखने वाला पर गलत उत्तर" देती है, इसलिए नाम से इसे रोका जाता है। मान के आकार से unit का अनुमान बिल्कुल नहीं लगाया जाता।

• String से ValueError — float('50') सफल हो जाता है, इसलिए बिना parse किया config value length के रूप में निकल जाता (मापा गया: thin_lens('50', '200') सही दिखने वाला 66.667 mm लौटा रहा था)। bool को भी True == 1 वाले implicit promotion के रूप में अस्वीकार किया जाता है।

• complex / masked array से ValueError (केवल real-valued slots; imaginary part को चुपचाप छोड़ना या mask हटाना अस्वीकार)। हर input पर NaN/Inf से ValueError।

• शून्य से भाग और उससे जुड़े मामलों को नाम लेकर अस्वीकार किया जाता है: focal length 0, radius of curvature 0, refractive index <= 0, पूरी तरह opaque aperture (सब 0, इसलिए normalisation 0/0), ऐसा PSF जिसका योग <= 0 हो, S0 = 0 वाला Stokes vector, और front focal point पर रखा object (image अनंत पर)।

• केवल 2 ops non-finite value लौटाते हैं, और वह भी contract के रूप में लिखा हुआ: depth_of_field hyperfocal distance से आगे far_mm = inf (hyperfocal distance की परिभाषा यही है) और gaussian_beam waist पर wavefront_radius_mm = inf (plane wavefront की radius of curvature)। दोनों साथ में एक finite साथी (far_is_infinite / curvature_per_mm) भी लौटाते हैं। इसके अलावा कोई भी चुपचाप आया NaN/Inf अंदर ही पकड़ा जाता है और ValueError देता है — "float64 overflow हुआ" और "उत्तर अनंत है" अलग-अलग दावे हैं, इसलिए पहले को दूसरे का चेहरा पहनाकर नहीं लौटाया जाता।

• Size की सीमा: generate होने वाले grid के लिए optics.MAX_GRID (4096), दिए गए field/PSF/aperture के लिए optics.MAX_FIELD_ELEMENTS (2^24), ABCD element chain के लिए optics.MAX_SYSTEM_ELEMENTS (1024), Zernike के लिए MAX_ZERNIKE_TERMS (512) / MAX_ZERNIKE_ORDER (40) / MAX_ZERNIKE_BASIS (2^25)। छोटे argument से बहुत बड़ा internal allocation होने वाले रास्तों (मापा गया: n_max=40 × 4096² पर 108 GB) को fail-closed ढंग से बंद किया गया है।

• भौतिक रूप से असंभव states भी अस्वीकार: degree of polarisation > 1 वाला Stokes vector, ऋणात्मक transmittance, ऋणात्मक intensity, और n-|m| odd जैसे अमान्य Zernike indices।

उपयोग की विस्तृत guide

• optics_imaging family guide

संदर्भ (sample data, साहित्य)

• Sample data catalog (DL URL / license) — 2-D के लिए skimage.data (BSD/public) + synthetic, 3-D के लिए real data sources (Stanford/PDS आदि) के DL URL।

• Operators का स्रोत और संदर्भ — इस op family के मूल research/methods के स्रोत।

• Algorithm का मानक स्रोत (लेखक, वर्ष) और उपयोग ऊपर की family usage guide में दिए गए हैं।

चलाए जा सकने वाले उदाहरण (verified samples जो सच में यह op call करते हैं)

• poc_lens_merit_operands — py -3.11 examples/poc_lens_merit_operands.py

वे अगले ops जिनसे type जुड़ता है (table को input के रूप में लेते हैं)

abcd_matrix · wavefront_stats · sfr_from_edge · paraxial_trace · seidel_coefficients · spot_stats · tolerance_analysis · wavefront_from_opd

इसी category के (lensmerit)

merit_operands_list · merit_default · merit_boundary_operands · optimize_lens_operands · merit_contributions


*Provenance: lensmerit.py — OPTICS operator registry. यह per-op note tools/opdocs.py md से अपने-आप generate होता है (हाथ से edit न करें)।*

© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.