photoelastic op• Data kinds: image2d × image2d × image2d → table
• Call: import fullseye as fs; fs.ledger.cubic_photoelastic_birefringence(sxx, syy, txy, pi11_minus_pi12, pi44, refractive_index, crystal_angle_deg=0.0) (implementation को सीधे call करना हो तो import photoelastic; photoelastic.cubic_photoelastic_birefringence(sxx, syy, txy, pi11_minus_pi12, pi44, refractive_index, crystal_angle_deg=0.0), registry से लेना हो तो opsoptics.get("cubic_photoelastic_birefringence"))
• GPU: नहीं (केवल CPU reference implementation)
{100} तल पर घनीय क्रिस्टल (जैसे Si वेफ़र) का प्रतिबल द्विअपवर्तन Δn और मंद अक्ष → dict।
> नीचे का विस्तृत विवरण मूल text में है — सारांश और headings का अनुवाद हो चुका है।
Light travels along the [100] normal of a {100} plate. crystal_angle_deg φ is the angle between the [010] direction and the x axis (Stoehr et al. 2020, Eq. 19). The 4th-order piezo-optic tensor of the cubic crystal is rotated, π'_ijkl = R_im R_jn R_ko R_lp π_mnop, the impermeability change is ΔB_ij = π'_ijkl σ_kl and the in-plane birefringence is Δn = (n³/2)(ΔB_II − ΔB_I) — which reproduces Eq. (26) of that paper::
Δn = (n³/4) √( 2((π11−π12)² + π44²)(4τ² + (σx−σy)²)
+ 2((π11−π12)² − π44²)[((σx−σy)² − 4τ²) cos 4φ
+ 4τ(σx−σy) sin 4φ] )
With π44 = π11 − π12 the material is photoelastically isotropic and Δn = (n³/2)|π11 − π12|(σ1 − σ2). In the sign of photoelastic_retardation (n_σ1 − n_σ2 = C(σ1 − σ2); a change of impermeability ΔB changes the index by −(n³/2)ΔB) that is C = −n³(π11−π12)/2, returned as isotropic_coefficient: π11 − π12 > 0 makes σ1 the fast axis, and the value can be passed to photoelastic_retardation as it is.
Coefficients are in Pa⁻¹ (silicon_piezo_optic_coefficients). Returns delta_n (≥ 0), slow_axis_deg (direction of the larger refractive index, (−90, 90]; 0 where Δn = 0 to 1e-12 of that pixel's own stress scale — the fast axis that optics.mueller_element and photoelastic_polarization_images take is slow_axis_deg + 90), effective_coefficient (|Δn|/(σ1−σ2) ≥ 0 [Pa⁻¹], 0 where σ1 = σ2), isotropic_coefficient (signed, as above) and isotropic (σ1 = σ2).
Raises ValueError: non-finite / empty / non-broadcasting stresses, |π| outside [1e-15, 1e-9] Pa⁻¹ (a table value in ×10⁻⁷ MPa⁻¹ or in MPa⁻¹ passed as Pa⁻¹), or a refractive index outside [1, 10].
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 · photoelastic_retardation · retardation_convert · 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.