photoelastic op• Data kinds: image2d → table
• Call: import fullseye as fs; fs.ledger.photoelastic_retardation(stress_difference_pa, thickness_m, stress_optic_coefficient, wavelength_m) (implementation को सीधे call करना हो तो import photoelastic; photoelastic.photoelastic_retardation(stress_difference_pa, thickness_m, stress_optic_coefficient, wavelength_m), registry से लेना हो तो opsoptics.get("photoelastic_retardation"))
• GPU: नहीं (केवल CPU reference implementation)
समदैशिक पदार्थ का प्रतिबल-प्रकाशिक नियम δ = 2π t C (σ1−σ2)/λ → dict (rad, nm, फ्रिंज क्रम और पदार्थ फ्रिंज मान)।
> नीचे का विस्तृत विवरण मूल text में है — सारांश और headings का अनुवाद हो चुका है।
stress_optic_coefficient C is in Pa⁻¹ (1 brewster = 1e-12 Pa⁻¹); the birefringence is Δn = C (σ1 − σ2) and the material fringe value f_σ = λ/C [N/m per fringe], so N = t (σ1 − σ2) / f_σ.
Sign: Δn = n_σ1 − n_σ2, so for C > 0 the σ1 direction (principal_stress angle_deg) is the slow axis. The retarder angle of optics.mueller_element and photoelastic_polarization_images is the fast axis (exp(−iωt), the fast axis leads), so pass axis_deg = angle_deg + 90 together with this delta_rad — for either sign of C (equivalently −delta_rad with angle_deg).
Returns delta_rad, retardation_nm, fringe_order, birefringence and fringe_value_n_per_m.
Raises ValueError: a non-finite/empty stress, |C| outside [1e-15, 1e-8] Pa⁻¹ (a value in brewster or in MPa⁻¹ passed as Pa⁻¹), and the thickness / wavelength checks of retardation_from_birefringence.
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 में दिए गए हैं।
• photoelastic_roundtrip — py -3.11 examples/photoelastic_roundtrip.py
table को input के रूप में लेते हैं)abcd_matrix · wavefront_stats · sfr_from_edge · paraxial_trace · seidel_coefficients · spot_stats · tolerance_analysis · wavefront_from_opd
photoelastic)brazilian_disc_stress · kirsch_plate_stress · principal_stress · retardation_from_birefringence · retardation_convert · cubic_photoelastic_birefringence · silicon_piezo_optic_coefficients · polariscope_intensity
*Provenance: photoelastic.py — OPTICS operator registry. यह per-op note tools/opdocs.py md से अपने-आप generate होता है (हाथ से edit न करें)।*
© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.