==================== rs03_push.out ====================

W1 sharp step, 400 random cases (89 beyond the turning point): max |r - Gamma| = 1.7e-11  PASS
W2 tanh: th10: 7.5e-04 (order 1.01)  th25: 1.1e-03 (order 1.01)  th29: 1.4e-03 (order 1.00)  th31: 1.8e-03 (order 0.98)  th60: 9.9e-04 (order 0.82)
W2 erf : th10: 5.9e-04 (order 1.01)  th25: 8.3e-04 (order 1.00)  th29: 1.1e-03 (order 1.00)  th31: 1.4e-03 (order 0.98)  th60: 8.0e-04 (order 0.87)
W2 ramp: th10: 4.7e-04 (order 1.01)  th25: 6.7e-04 (order 1.00)  th29: 9.0e-04 (order 1.00)  th31: 1.2e-03 (order 0.98)  th60: 6.5e-04 (order 0.90)
W2 Richardson to w = 0, all profiles and angles: max |r0 - Gamma| = 7.4e-04  KILL
W3 reciprocity: |r21 + r12| step 8.5e-15; ||r21| - |r12|| smooth 1.3e-13  PASS
W4 flux |r|^2 + (Y2/Y1)|t|^2 = 1: max error 6.4e-14  PASS
W5 opposite: thick layer (omega w = 3), b2 x4 changes |r| by up to 0.160  PASS (g_perp matters off the step)
   omega w = 10: half-reflection at 29.9960 deg
   omega w = 30: half-reflection at 29.9995 deg
   omega w = 100: half-reflection at 29.9999 deg
W6 |r| = 1 beyond theta_c = 30.0000: max dev 2.2e-16; half-reflection at omega w = 10, 30, 100: 29.996, 30.000, 30.000  PASS
W7 gap family Gamma = tanh(ln(D2/D1)/2), 200 pairs: max error 2.9e-13  PASS

==================== rs03b.out ====================
W2b tanh: th10: 7.5e-04 (order 1.05)  th25: 1.1e-03 (order 1.02)  th29: 1.4e-03 (order 1.01)  th31: 1.8e-03 (order 0.95)  th60: 9.9e-04 (order 0.57)
W2b erf : th10: 5.9e-04 (order 1.03)  th25: 8.3e-04 (order 1.01)  th29: 1.1e-03 (order 1.01)  th31: 1.4e-03 (order 0.96)  th60: 8.0e-04 (order 0.70)
W2b ramp: th10: 4.7e-04 (order 1.02)  th25: 6.7e-04 (order 1.01)  th29: 9.0e-04 (order 1.00)  th31: 1.2e-03 (order 0.97)  th60: 6.5e-04 (order 0.77)
W2b Richardson to w = 0: max |r0 - Gamma| = 7.4e-04  KILL

==================== rs03c.out ====================
theta 10: |r - Gamma|/(omega w) = 0.0747, 0.0746, 0.0746, 0.0746  at omega w = 0.01, 0.005, 0.0025, 0.00125
theta 25: |r - Gamma|/(omega w) = 0.1054, 0.1053, 0.1053, 0.1053  at omega w = 0.01, 0.005, 0.0025, 0.00125
theta 31: |r - Gamma|/(omega w) = 0.1807, 0.1818, 0.1823, 0.1826  at omega w = 0.01, 0.005, 0.0025, 0.00125
theta 60: |r - Gamma|/(omega w) = 0.0992, 0.1029, 0.1048, 0.1057  at omega w = 0.01, 0.005, 0.0025, 0.00125
RS03C: PASS (first-order convergence to Gamma)