design op• Data kinds: table → table
• Call: import fullseye as fs; fs.ledger.seidel_coefficients(system, field=None) (implementation को सीधे call करना हो तो import raytrace; raytrace.seidel_coefficients(system, field=None), registry से लेना हो तो opsoptics.get("seidel_coefficients"))
हर surface पर और कुल third-order (Seidel) aberration sums (table)।
> नीचे का विस्तृत विवरण मूल text में है — सारांश और headings का अनुवाद हो चुका है।
Uses the paraxial marginal ray (through the stop edge) and chief ray (through the stop centre) at the system's default field (or *field*), with Welford's refraction-invariant form: A = n(yc + u), Δ(u/n), Lagrange invariant H = n(ūy − uȳ)::
S_I = −Σ A² y Δ(u/n) spherical
S_II = −Σ A Ā y Δ(u/n) coma
S_III = −Σ Ā² y Δ(u/n) astigmatism
S_IV = −Σ H² c Δ(1/n) Petzval (field curvature)
S_V = −Σ (Ā/A)[Ā² y Δ(u/n) + H² c Δ(1/n)] distortion
C_L = −Σ A y Δ(δn/n), C_T = −Σ Ā y Δ(δn/n) axial / lateral colour
plus the aspheric deformation contribution 8 G (n'−n) y⁴ with G = k c³/8 + A4 (conic and fourth-order coefficient) to S_I (and its ȳ/y powers to S_II, S_III, S_V); higher aspheric orders are fifth order and above and do not enter the Seidel sums. Sums are in millimetres of wavefront times 8:
the third-order wavefront at the pupil edge is W040 = S_I/8, W131 = S_II/2, W222 = S_III/2, W220 = (S_III + S_IV)/4, W311 = S_V/2 (Welford's normalisation). waves gives the same numbers divided by the wavelength. Checked against the exact ray-traced OPD at small aperture (tests/test_raytrace.py).
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।
• 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 में दिए गए हैं।
• lens_design_demo — py -3.11 examples/lens_design_demo.py
table को input के रूप में लेते हैं)abcd_matrix · wavefront_stats · sfr_from_edge · paraxial_trace · spot_stats · tolerance_analysis · wavefront_from_opd · spot_diagram
design)lens_system · thick_lens · glass · example_system · glass_catalog · sellmeier · paraxial_trace · spot_stats
*Provenance: raytrace.py — OPTICS operator registry. यह per-op note tools/opdocs.py md से अपने-आप generate होता है (हाथ से edit न करें)।*
© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.