[reduction in_plane_lc_pillar normal] gauge=-1j (shipped)    dR=5.024e-15 dT=2.098e-14 dJ=7.961e-15
[reduction in_plane_lc_pillar normal] gauge=+1j              dR=1.064e-02 dT=7.839e-02 dJ=5.036e-02
C:\tmp\lum_aniso\validation\probe_verify_staggered_oop\v2c_hgauge.py:76: UserWarning: PMM2DStackPure.solve: lossless energy closure violated for incident E_x (sum R+T = 6.484, off by +5.48; every permittivity in this stack is real or Hermitian, so R+T = 1 is exact).  The staggered basis is under-resolved or the solve sits inside a Rayleigh cutoff: raise n_modes (M), or detune the wavelength / use PMM2DStackHybrid near a cutoff.
  res = st.solve(jones=True, retain_internal=retain)
C:\tmp\lum_aniso\validation\probe_verify_staggered_oop\v2c_hgauge.py:76: UserWarning: PMM2DStackPure.solve: lossless energy closure violated for incident E_y (sum R+T = 3.367, off by +2.37; every permittivity in this stack is real or Hermitian, so R+T = 1 is exact).  The staggered basis is under-resolved or the solve sits inside a Rayleigh cutoff: raise n_modes (M), or detune the wavelength / use PMM2DStackHybrid near a cutoff.
  res = st.solve(jones=True, retain_internal=retain)
[reduction in_plane_lc_pillar normal] gauge=+1               dR=3.180e+00 dT=2.218e+00 dJ=1.421e+00
C:\tmp\lum_aniso\validation\probe_verify_staggered_oop\v2c_hgauge.py:76: UserWarning: PMM2DStackPure.solve: lossless energy closure violated for incident E_x (sum R+T = 21.28, off by +20.3; every permittivity in this stack is real or Hermitian, so R+T = 1 is exact).  The staggered basis is under-resolved or the solve sits inside a Rayleigh cutoff: raise n_modes (M), or detune the wavelength / use PMM2DStackHybrid near a cutoff.
  res = st.solve(jones=True, retain_internal=retain)
C:\tmp\lum_aniso\validation\probe_verify_staggered_oop\v2c_hgauge.py:76: UserWarning: PMM2DStackPure.solve: lossless energy closure violated for incident E_y (sum R+T = 8.375, off by +7.37; every permittivity in this stack is real or Hermitian, so R+T = 1 is exact).  The staggered basis is under-resolved or the solve sits inside a Rayleigh cutoff: raise n_modes (M), or detune the wavelength / use PMM2DStackHybrid near a cutoff.
  res = st.solve(jones=True, retain_internal=retain)
[reduction in_plane_lc_pillar normal] gauge=-1               dR=7.864e+00 dT=6.355e+00 dJ=1.998e+00
[reduction in_plane_lc_pillar normal] gauge=-2j (scale x2)   dR=3.681e-01 dT=3.747e-01 dJ=3.678e-01
[reduction in_plane_lc_pillar conical20_35] gauge=-1j (shipped)    dR=4.600e-15 dT=1.188e-14 dJ=1.200e-14
[reduction in_plane_lc_pillar conical20_35] gauge=+1j              dR=2.241e-03 dT=2.354e-02 dJ=3.429e-02
C:\tmp\lum_aniso\validation\probe_verify_staggered_oop\v2c_hgauge.py:76: UserWarning: PMM2DStackPure.solve: lossless energy closure violated for incident E_x (sum R+T = 4.604, off by +3.6; every permittivity in this stack is real or Hermitian, so R+T = 1 is exact).  The staggered basis is under-resolved or the solve sits inside a Rayleigh cutoff: raise n_modes (M), or detune the wavelength / use PMM2DStackHybrid near a cutoff.
  res = st.solve(jones=True, retain_internal=retain)
C:\tmp\lum_aniso\validation\probe_verify_staggered_oop\v2c_hgauge.py:76: UserWarning: PMM2DStackPure.solve: lossless energy closure violated for incident E_y (sum R+T = 2.593, off by +1.59; every permittivity in this stack is real or Hermitian, so R+T = 1 is exact).  The staggered basis is under-resolved or the solve sits inside a Rayleigh cutoff: raise n_modes (M), or detune the wavelength / use PMM2DStackHybrid near a cutoff.
  res = st.solve(jones=True, retain_internal=retain)
[reduction in_plane_lc_pillar conical20_35] gauge=+1               dR=2.168e+00 dT=1.359e+00 dJ=1.123e+00
C:\tmp\lum_aniso\validation\probe_verify_staggered_oop\v2c_hgauge.py:76: UserWarning: PMM2DStackPure.solve: lossless energy closure violated for incident E_x (sum R+T = 60.46, off by +59.5; every permittivity in this stack is real or Hermitian, so R+T = 1 is exact).  The staggered basis is under-resolved or the solve sits inside a Rayleigh cutoff: raise n_modes (M), or detune the wavelength / use PMM2DStackHybrid near a cutoff.
  res = st.solve(jones=True, retain_internal=retain)
C:\tmp\lum_aniso\validation\probe_verify_staggered_oop\v2c_hgauge.py:76: UserWarning: PMM2DStackPure.solve: lossless energy closure violated for incident E_y (sum R+T = 22.94, off by +21.9; every permittivity in this stack is real or Hermitian, so R+T = 1 is exact).  The staggered basis is under-resolved or the solve sits inside a Rayleigh cutoff: raise n_modes (M), or detune the wavelength / use PMM2DStackHybrid near a cutoff.
  res = st.solve(jones=True, retain_internal=retain)
[reduction in_plane_lc_pillar conical20_35] gauge=-1               dR=1.025e+01 dT=1.247e+01 dJ=1.842e+00
[reduction in_plane_lc_pillar conical20_35] gauge=-2j (scale x2)   dR=3.641e-01 dT=3.669e-01 dJ=3.619e-01
[reduction in_plane_lc_uniform normal] gauge=-1j (shipped)    dR=4.857e-15 dT=1.654e-14 dJ=1.191e-14
[reduction in_plane_lc_uniform normal] gauge=+1j              dR=8.299e-15 dT=1.998e-14 dJ=5.508e-02
C:\tmp\lum_aniso\validation\probe_verify_staggered_oop\v2c_hgauge.py:76: UserWarning: PMM2DStackPure.solve: lossless energy closure violated for incident E_x (sum R+T = 10.27, off by +9.27; every permittivity in this stack is real or Hermitian, so R+T = 1 is exact).  The staggered basis is under-resolved or the solve sits inside a Rayleigh cutoff: raise n_modes (M), or detune the wavelength / use PMM2DStackHybrid near a cutoff.
  res = st.solve(jones=True, retain_internal=retain)
C:\tmp\lum_aniso\validation\probe_verify_staggered_oop\v2c_hgauge.py:76: UserWarning: PMM2DStackPure.solve: lossless energy closure violated for incident E_y (sum R+T = 4.541, off by +3.54; every permittivity in this stack is real or Hermitian, so R+T = 1 is exact).  The staggered basis is under-resolved or the solve sits inside a Rayleigh cutoff: raise n_modes (M), or detune the wavelength / use PMM2DStackHybrid near a cutoff.
  res = st.solve(jones=True, retain_internal=retain)
[reduction in_plane_lc_uniform normal] gauge=+1               dR=5.312e+00 dT=3.956e+00 dJ=1.839e+00
C:\tmp\lum_aniso\validation\probe_verify_staggered_oop\v2c_hgauge.py:76: UserWarning: PMM2DStackPure.solve: lossless energy closure violated for incident E_x (sum R+T = 11.13, off by +10.1; every permittivity in this stack is real or Hermitian, so R+T = 1 is exact).  The staggered basis is under-resolved or the solve sits inside a Rayleigh cutoff: raise n_modes (M), or detune the wavelength / use PMM2DStackHybrid near a cutoff.
  res = st.solve(jones=True, retain_internal=retain)
C:\tmp\lum_aniso\validation\probe_verify_staggered_oop\v2c_hgauge.py:76: UserWarning: PMM2DStackPure.solve: lossless energy closure violated for incident E_y (sum R+T = 4.8, off by +3.8; every permittivity in this stack is real or Hermitian, so R+T = 1 is exact).  The staggered basis is under-resolved or the solve sits inside a Rayleigh cutoff: raise n_modes (M), or detune the wavelength / use PMM2DStackHybrid near a cutoff.
  res = st.solve(jones=True, retain_internal=retain)
[reduction in_plane_lc_uniform normal] gauge=-1               dR=5.670e+00 dT=4.461e+00 dJ=1.664e+00
[reduction in_plane_lc_uniform normal] gauge=-2j (scale x2)   dR=4.027e-01 dT=4.027e-01 dJ=3.655e-01
[reduction in_plane_lc_uniform conical20_35] gauge=-1j (shipped)    dR=8.868e-15 dT=2.864e-14 dJ=1.787e-14
[reduction in_plane_lc_uniform conical20_35] gauge=+1j              dR=7.896e-11 dT=1.502e-11 dJ=4.827e-02
C:\tmp\lum_aniso\validation\probe_verify_staggered_oop\v2c_hgauge.py:76: UserWarning: PMM2DStackPure.solve: lossless energy closure violated for incident E_x (sum R+T = 6.346, off by +5.35; every permittivity in this stack is real or Hermitian, so R+T = 1 is exact).  The staggered basis is under-resolved or the solve sits inside a Rayleigh cutoff: raise n_modes (M), or detune the wavelength / use PMM2DStackHybrid near a cutoff.
  res = st.solve(jones=True, retain_internal=retain)
C:\tmp\lum_aniso\validation\probe_verify_staggered_oop\v2c_hgauge.py:76: UserWarning: PMM2DStackPure.solve: lossless energy closure violated for incident E_y (sum R+T = 3.158, off by +2.16; every permittivity in this stack is real or Hermitian, so R+T = 1 is exact).  The staggered basis is under-resolved or the solve sits inside a Rayleigh cutoff: raise n_modes (M), or detune the wavelength / use PMM2DStackHybrid near a cutoff.
  res = st.solve(jones=True, retain_internal=retain)
[reduction in_plane_lc_uniform conical20_35] gauge=+1               dR=3.250e+00 dT=2.096e+00 dJ=1.368e+00
C:\tmp\lum_aniso\validation\probe_verify_staggered_oop\v2c_hgauge.py:76: UserWarning: PMM2DStackPure.solve: lossless energy closure violated for incident E_x (sum R+T = 31.48, off by +30.5; every permittivity in this stack is real or Hermitian, so R+T = 1 is exact).  The staggered basis is under-resolved or the solve sits inside a Rayleigh cutoff: raise n_modes (M), or detune the wavelength / use PMM2DStackHybrid near a cutoff.
  res = st.solve(jones=True, retain_internal=retain)
C:\tmp\lum_aniso\validation\probe_verify_staggered_oop\v2c_hgauge.py:76: UserWarning: PMM2DStackPure.solve: lossless energy closure violated for incident E_y (sum R+T = 12.2, off by +11.2; every permittivity in this stack is real or Hermitian, so R+T = 1 is exact).  The staggered basis is under-resolved or the solve sits inside a Rayleigh cutoff: raise n_modes (M), or detune the wavelength / use PMM2DStackHybrid near a cutoff.
  res = st.solve(jones=True, retain_internal=retain)
[reduction in_plane_lc_uniform conical20_35] gauge=-1               dR=1.594e+01 dT=1.454e+01 dJ=3.045e+00
[reduction in_plane_lc_uniform conical20_35] gauge=-2j (scale x2)   dR=3.999e-01 dT=3.999e-01 dJ=3.568e-01
[reduction gyrotropic_pillar normal] gauge=-1j (shipped)    dR=2.498e-15 dT=7.438e-15 dJ=1.025e-14
[reduction gyrotropic_pillar normal] gauge=+1j              dR=5.240e-03 dT=2.207e-02 dJ=1.922e-02
C:\tmp\lum_aniso\validation\probe_verify_staggered_oop\v2c_hgauge.py:76: UserWarning: PMM2DStackPure.solve: lossless energy closure violated for incident E_x (sum R+T = 3.677, off by +2.68; every permittivity in this stack is real or Hermitian, so R+T = 1 is exact).  The staggered basis is under-resolved or the solve sits inside a Rayleigh cutoff: raise n_modes (M), or detune the wavelength / use PMM2DStackHybrid near a cutoff.
  res = st.solve(jones=True, retain_internal=retain)
C:\tmp\lum_aniso\validation\probe_verify_staggered_oop\v2c_hgauge.py:76: UserWarning: PMM2DStackPure.solve: lossless energy closure violated for incident E_y (sum R+T = 3.677, off by +2.68; every permittivity in this stack is real or Hermitian, so R+T = 1 is exact).  The staggered basis is under-resolved or the solve sits inside a Rayleigh cutoff: raise n_modes (M), or detune the wavelength / use PMM2DStackHybrid near a cutoff.
  res = st.solve(jones=True, retain_internal=retain)
[reduction gyrotropic_pillar normal] gauge=+1               dR=1.742e+00 dT=9.235e-01 dJ=1.069e+00
C:\tmp\lum_aniso\validation\probe_verify_staggered_oop\v2c_hgauge.py:76: UserWarning: PMM2DStackPure.solve: lossless energy closure violated for incident E_x (sum R+T = 538.6, off by +538; every permittivity in this stack is real or Hermitian, so R+T = 1 is exact).  The staggered basis is under-resolved or the solve sits inside a Rayleigh cutoff: raise n_modes (M), or detune the wavelength / use PMM2DStackHybrid near a cutoff.
  res = st.solve(jones=True, retain_internal=retain)
C:\tmp\lum_aniso\validation\probe_verify_staggered_oop\v2c_hgauge.py:76: UserWarning: PMM2DStackPure.solve: lossless energy closure violated for incident E_y (sum R+T = 538.6, off by +538; every permittivity in this stack is real or Hermitian, so R+T = 1 is exact).  The staggered basis is under-resolved or the solve sits inside a Rayleigh cutoff: raise n_modes (M), or detune the wavelength / use PMM2DStackHybrid near a cutoff.
  res = st.solve(jones=True, retain_internal=retain)
[reduction gyrotropic_pillar normal] gauge=-1               dR=2.274e+02 dT=2.715e+02 dJ=1.148e+01
[reduction gyrotropic_pillar normal] gauge=-2j (scale x2)   dR=3.849e-01 dT=3.877e-01 dJ=4.242e-01
[reduction gyrotropic_pillar conical20_35] gauge=-1j (shipped)    dR=1.209e-14 dT=2.331e-14 dJ=2.658e-14
[reduction gyrotropic_pillar conical20_35] gauge=+1j              dR=2.417e-03 dT=2.298e-02 dJ=1.804e-02
C:\tmp\lum_aniso\validation\probe_verify_staggered_oop\v2c_hgauge.py:76: UserWarning: PMM2DStackPure.solve: lossless energy closure violated for incident E_x (sum R+T = 2.956, off by +1.96; every permittivity in this stack is real or Hermitian, so R+T = 1 is exact).  The staggered basis is under-resolved or the solve sits inside a Rayleigh cutoff: raise n_modes (M), or detune the wavelength / use PMM2DStackHybrid near a cutoff.
  res = st.solve(jones=True, retain_internal=retain)
C:\tmp\lum_aniso\validation\probe_verify_staggered_oop\v2c_hgauge.py:76: UserWarning: PMM2DStackPure.solve: lossless energy closure violated for incident E_y (sum R+T = 2.918, off by +1.92; every permittivity in this stack is real or Hermitian, so R+T = 1 is exact).  The staggered basis is under-resolved or the solve sits inside a Rayleigh cutoff: raise n_modes (M), or detune the wavelength / use PMM2DStackHybrid near a cutoff.
  res = st.solve(jones=True, retain_internal=retain)
[reduction gyrotropic_pillar conical20_35] gauge=+1               dR=1.322e+00 dT=6.232e-01 dJ=9.124e-01
C:\tmp\lum_aniso\validation\probe_verify_staggered_oop\v2c_hgauge.py:76: UserWarning: PMM2DStackPure.solve: lossless energy closure violated for incident E_x (sum R+T = 49.31, off by +48.3; every permittivity in this stack is real or Hermitian, so R+T = 1 is exact).  The staggered basis is under-resolved or the solve sits inside a Rayleigh cutoff: raise n_modes (M), or detune the wavelength / use PMM2DStackHybrid near a cutoff.
  res = st.solve(jones=True, retain_internal=retain)
C:\tmp\lum_aniso\validation\probe_verify_staggered_oop\v2c_hgauge.py:76: UserWarning: PMM2DStackPure.solve: lossless energy closure violated for incident E_y (sum R+T = 51.4, off by +50.4; every permittivity in this stack is real or Hermitian, so R+T = 1 is exact).  The staggered basis is under-resolved or the solve sits inside a Rayleigh cutoff: raise n_modes (M), or detune the wavelength / use PMM2DStackHybrid near a cutoff.
  res = st.solve(jones=True, retain_internal=retain)
[reduction gyrotropic_pillar conical20_35] gauge=-1               dR=1.908e+01 dT=2.610e+01 dJ=3.776e+00
[reduction gyrotropic_pillar conical20_35] gauge=-2j (scale x2)   dR=3.597e-01 dT=3.614e-01 dJ=4.085e-01
[reduction gyrotropic_uniform normal] gauge=-1j (shipped)    dR=6.384e-15 dT=2.764e-14 dJ=1.432e-14
[reduction gyrotropic_uniform normal] gauge=+1j              dR=4.816e-15 dT=7.883e-15 dJ=2.667e-02
C:\tmp\lum_aniso\validation\probe_verify_staggered_oop\v2c_hgauge.py:76: UserWarning: PMM2DStackPure.solve: lossless energy closure violated for incident E_x (sum R+T = 4.423, off by +3.42; every permittivity in this stack is real or Hermitian, so R+T = 1 is exact).  The staggered basis is under-resolved or the solve sits inside a Rayleigh cutoff: raise n_modes (M), or detune the wavelength / use PMM2DStackHybrid near a cutoff.
  res = st.solve(jones=True, retain_internal=retain)
C:\tmp\lum_aniso\validation\probe_verify_staggered_oop\v2c_hgauge.py:76: UserWarning: PMM2DStackPure.solve: lossless energy closure violated for incident E_y (sum R+T = 4.423, off by +3.42; every permittivity in this stack is real or Hermitian, so R+T = 1 is exact).  The staggered basis is under-resolved or the solve sits inside a Rayleigh cutoff: raise n_modes (M), or detune the wavelength / use PMM2DStackHybrid near a cutoff.
  res = st.solve(jones=True, retain_internal=retain)
[reduction gyrotropic_uniform normal] gauge=+1               dR=2.198e+00 dT=1.225e+00 dJ=1.197e+00
C:\tmp\lum_aniso\validation\probe_verify_staggered_oop\v2c_hgauge.py:76: UserWarning: PMM2DStackPure.solve: lossless energy closure violated for incident E_x (sum R+T = 259.4, off by +258; every permittivity in this stack is real or Hermitian, so R+T = 1 is exact).  The staggered basis is under-resolved or the solve sits inside a Rayleigh cutoff: raise n_modes (M), or detune the wavelength / use PMM2DStackHybrid near a cutoff.
  res = st.solve(jones=True, retain_internal=retain)
C:\tmp\lum_aniso\validation\probe_verify_staggered_oop\v2c_hgauge.py:76: UserWarning: PMM2DStackPure.solve: lossless energy closure violated for incident E_y (sum R+T = 259.4, off by +258; every permittivity in this stack is real or Hermitian, so R+T = 1 is exact).  The staggered basis is under-resolved or the solve sits inside a Rayleigh cutoff: raise n_modes (M), or detune the wavelength / use PMM2DStackHybrid near a cutoff.
  res = st.solve(jones=True, retain_internal=retain)
[reduction gyrotropic_uniform normal] gauge=-1               dR=1.241e+02 dT=1.343e+02 dJ=8.437e+00
[reduction gyrotropic_uniform normal] gauge=-2j (scale x2)   dR=4.136e-01 dT=4.136e-01 dJ=4.310e-01
[reduction gyrotropic_uniform conical20_35] gauge=-1j (shipped)    dR=2.637e-15 dT=1.465e-14 dJ=2.706e-14
[reduction gyrotropic_uniform conical20_35] gauge=+1j              dR=1.078e-03 dT=1.078e-03 dJ=3.568e-02
C:\tmp\lum_aniso\validation\probe_verify_staggered_oop\v2c_hgauge.py:76: UserWarning: PMM2DStackPure.solve: lossless energy closure violated for incident E_x (sum R+T = 3.454, off by +2.45; every permittivity in this stack is real or Hermitian, so R+T = 1 is exact).  The staggered basis is under-resolved or the solve sits inside a Rayleigh cutoff: raise n_modes (M), or detune the wavelength / use PMM2DStackHybrid near a cutoff.
  res = st.solve(jones=True, retain_internal=retain)
C:\tmp\lum_aniso\validation\probe_verify_staggered_oop\v2c_hgauge.py:76: UserWarning: PMM2DStackPure.solve: lossless energy closure violated for incident E_y (sum R+T = 3.403, off by +2.4; every permittivity in this stack is real or Hermitian, so R+T = 1 is exact).  The staggered basis is under-resolved or the solve sits inside a Rayleigh cutoff: raise n_modes (M), or detune the wavelength / use PMM2DStackHybrid near a cutoff.
  res = st.solve(jones=True, retain_internal=retain)
[reduction gyrotropic_uniform conical20_35] gauge=+1               dR=1.643e+00 dT=8.108e-01 dJ=1.015e+00
C:\tmp\lum_aniso\validation\probe_verify_staggered_oop\v2c_hgauge.py:76: UserWarning: PMM2DStackPure.solve: lossless energy closure violated for incident E_x (sum R+T = 530.4, off by +529; every permittivity in this stack is real or Hermitian, so R+T = 1 is exact).  The staggered basis is under-resolved or the solve sits inside a Rayleigh cutoff: raise n_modes (M), or detune the wavelength / use PMM2DStackHybrid near a cutoff.
  res = st.solve(jones=True, retain_internal=retain)
C:\tmp\lum_aniso\validation\probe_verify_staggered_oop\v2c_hgauge.py:76: UserWarning: PMM2DStackPure.solve: lossless energy closure violated for incident E_y (sum R+T = 552.7, off by +552; every permittivity in this stack is real or Hermitian, so R+T = 1 is exact).  The staggered basis is under-resolved or the solve sits inside a Rayleigh cutoff: raise n_modes (M), or detune the wavelength / use PMM2DStackHybrid near a cutoff.
  res = st.solve(jones=True, retain_internal=retain)
[reduction gyrotropic_uniform conical20_35] gauge=-1               dR=2.500e+02 dT=3.018e+02 dJ=1.204e+01
[reduction gyrotropic_uniform conical20_35] gauge=-2j (scale x2)   dR=3.903e-01 dT=3.903e-01 dJ=4.174e-01
[reduction lossy_diagonal_pillar normal] gauge=-1j (shipped)    dR=9.673e-15 dT=3.042e-14 dJ=1.553e-14
[reduction lossy_diagonal_pillar normal] gauge=+1j              dR=3.040e-02 dT=2.854e-01 dJ=2.379e-01
[reduction lossy_diagonal_pillar normal] gauge=+1               dR=5.524e+01 dT=7.282e+01 dJ=7.597e+00
[reduction lossy_diagonal_pillar normal] gauge=-1               dR=3.379e+00 dT=2.185e+00 dJ=1.614e+00
[reduction lossy_diagonal_pillar normal] gauge=-2j (scale x2)   dR=3.890e-01 dT=2.972e-01 dJ=3.769e-01
[reduction lossy_diagonal_pillar conical20_35] gauge=-1j (shipped)    dR=1.638e-15 dT=2.698e-14 dJ=1.865e-14
[reduction lossy_diagonal_pillar conical20_35] gauge=+1j              dR=2.630e-02 dT=2.745e-01 dJ=2.458e-01
[reduction lossy_diagonal_pillar conical20_35] gauge=+1               dR=8.882e+01 dT=9.835e+01 dJ=9.346e+00
[reduction lossy_diagonal_pillar conical20_35] gauge=-1               dR=1.368e+01 dT=1.420e+01 dJ=1.296e+00
[reduction lossy_diagonal_pillar conical20_35] gauge=-2j (scale x2)   dR=3.876e-01 dT=3.045e-01 dJ=3.757e-01
[reduction lossy_diagonal_uniform normal] gauge=-1j (shipped)    dR=3.636e-15 dT=1.998e-14 dJ=5.372e-15
[reduction lossy_diagonal_uniform normal] gauge=+1j              dR=4.165e-02 dT=4.497e-01 dJ=2.539e-01
[reduction lossy_diagonal_uniform normal] gauge=+1               dR=1.302e+01 dT=1.776e+01 dJ=3.541e+00
[reduction lossy_diagonal_uniform normal] gauge=-1               dR=3.064e+00 dT=1.893e+00 dJ=1.447e+00
[reduction lossy_diagonal_uniform normal] gauge=-2j (scale x2)   dR=4.323e-01 dT=3.547e-01 dJ=3.830e-01
[reduction lossy_diagonal_uniform conical20_35] gauge=-1j (shipped)    dR=1.169e-14 dT=1.699e-14 dJ=2.131e-14
[reduction lossy_diagonal_uniform conical20_35] gauge=+1j              dR=4.025e-02 dT=4.504e-01 dJ=2.403e-01
[reduction lossy_diagonal_uniform conical20_35] gauge=+1               dR=3.592e+01 dT=5.224e+01 dJ=6.261e+00
[reduction lossy_diagonal_uniform conical20_35] gauge=-1               dR=4.172e+00 dT=2.844e+00 dJ=1.736e+00
[reduction lossy_diagonal_uniform conical20_35] gauge=-2j (scale x2)   dR=4.344e-01 dT=3.667e-01 dJ=3.862e-01
[absorption] gauge=-1j (shipped)    A=[0.23998689, -0.0, -0.0] budget_dev=3.132e-11
[absorption] gauge=+1j              A=[-0.30212362, -0.0, -0.0] budget_dev=2.843e-10
[absorption] gauge=+1               A=[-0.0256308, 0.0, 0.0] budget_dev=3.017e-10
[absorption] gauge=-1               A=[0.01861994, 0.0, -0.0] budget_dev=2.138e-10
[absorption] gauge=-2j (scale x2)   A=[0.25264586, 0.0, -0.0] budget_dev=2.825e-11
DONE
