Two-point Green's dyadic against extended precision

← Verification · Spheres · Point dipole
Reference
two-point Green's dyadic against extended precision
Kind
precision: the same mathematics in mpmath at 60–80 digits, written separately in src/pystratify/references.py and benchmarks/mp_reference.py: proves the numerics, not the formulation
Calculated
G(r, r′)
Agreement
1e-13
Points 56 and 60 nm from the centre of a gold nanoshell (50 nm silica core, 5 nm gold) in water, 614 nm: the two-point Green's dyadic per order against 80 digits The radial part of the Green's dyadic between two points, per order and polarisation, which every dipole–dipole coupling near the sphere is built from: PyStratify's normalized form against transfer matrices in 80-digit arithmetic.
2040608010012010⁻¹⁴10⁻¹²10⁻¹⁰10⁻⁸10⁻⁶10⁻⁴10⁻²1order ℓ|g_ℓ|
TM, PyStratify TE, PyStratify TM, 80 digits TE, 80 digits
Hover to read values.
  • Σ(2ℓ+1)|g − g₈₀| / Σ(2ℓ+1)|g₈₀|: 2.96×10⁻¹⁵
  • PyStratify the lines, the reference the circles

Run it yourself: sphere_emitter_checks.py · test_pairs.py · figures.py · SPHERE_EMITTER_CHECKS.md on GitHub, PyStratify 97e52a9