PyStratify against treams: oblique incidence

← Verification · Cylinders · Plane wave

Reference: treams, oblique incidence
Kind: code: another group's code
Calculated: cross widths, coefficients
Agreement: 1e-12

A cylinder of radius 100 nm, n = 1.5 + 0.01i, at 60° to its axis: cross widths against treams Oblique incidence couples the two polarizations at the surface; both codes average the two incident polarizations. treams' cylindrical T-matrix is per axial wavenumber, which a plane wave fixes.

Q_ext, PyStratify Q_sca, PyStratify Q_ext, treams Q_sca, treams

Hover to read values.

largest relative difference: 1.25×10⁻¹⁵
PyStratify the lines, the reference the circles

Run it yourself: unified_checks.py · test_unified.py · figures.py · UNIFIED_CHECKS.md on GitHub, PyStratify 97e52a9