r2_onedim.py (verification run 2, AMR-022-4023)  seed=2030  numpy 2.5.2 sympy 1.14.0
(A) EXACT, n = 1, 2
  n = 1: right end 1 = 2*2^(-1): True
  n = 2: q(sqrt(1+d^2/4)) = 1: True; |q| at the critical point = d^2/4: True; b increasing in d, b(2) = sqrt 2: True
  (A) ALL PASSED
(B) EXHAUSTIVE GRIDS
  n = 3: grid 2563201 points (gaps in [0, 2.4], step 0.0015), feasible 1123169; best grid b = 1.5868394600; refined b = 1.587401051968; 2*2^(-1/n) = 1.587401051968; refined - Chebyshev = -1.80e-14; zeros [-1.374730e+00 -1.376677e-14  1.374730e+00]
  n = 4: grid 4173281 points (gaps in [0, 2.4], step 0.0150), feasible 700578; best grid b = 1.6771878017; refined b = 1.681646653120; 2*2^(-1/n) = 1.681792830507; refined - Chebyshev = -1.46e-04; zeros [-1.553748 -0.643555  0.643726  1.553576]
  (B) no value above 2*2^(-1/n)  (17.2 s)
(C) DIFFERENTIAL EVOLUTION, n = 3..10, two seeds
  n =  3  2*2^(-1/n) = 1.5874010520 | seed 2030: DE 1.5874010520, refined 1.587401051968 (crit 1.000e+00) | seed 2031: DE 1.5874010520, refined 1.587401051968 (crit 1.000e+00) | best - Chebyshev -2.86e-14
  n =  4  2*2^(-1/n) = 1.6817928305 | seed 2030: DE 1.6817928305, refined 1.681792830507 (crit 1.000e+00) | seed 2031: DE 1.6817928305, refined 1.681792830507 (crit 1.000e+00) | best - Chebyshev -2.69e-14
  n =  5  2*2^(-1/n) = 1.7411011266 | seed 2030: DE 1.7411011266, refined 1.741101126592 (crit 1.000e+00) | seed 2031: DE 1.7411011266, refined 1.741101126592 (crit 1.000e+00) | best - Chebyshev -1.32e-13
  n =  6  2*2^(-1/n) = 1.7817974363 | seed 2030: DE 1.7817974363, refined 1.781797436280 (crit 1.000e+00) | seed 2031: DE 1.7817974363, refined 1.781797436281 (crit 1.000e+00) | best - Chebyshev -8.90e-14
  n =  7  2*2^(-1/n) = 1.8114473285 | seed 2030: DE 1.8114473285, refined 1.811447328527 (crit 1.000e+00) | seed 2031: DE 1.8114473285, refined 1.811447328527 (crit 1.000e+00) | best - Chebyshev -7.21e-13
  n =  8  2*2^(-1/n) = 1.8340080864 | seed 2030: DE 1.8340080864, refined 1.834008086409 (crit 1.000e+00) | seed 2031: DE 1.8340080864, refined 1.834008086409 (crit 1.000e+00) | best - Chebyshev -3.04e-13
  n =  9  2*2^(-1/n) = 1.8517494246 | seed 2030: DE 1.8517494246, refined 1.851749424573 (crit 1.000e+00) | seed 2031: DE 1.8517494246, refined 1.851749424574 (crit 1.000e+00) | best - Chebyshev -4.80e-13
  n = 10  2*2^(-1/n) = 1.8660659831 | seed 2030: DE 1.8660659831, refined 1.866065983073 (crit 1.000e+00) | seed 2031: DE 1.8660659831, refined 1.866065983073 (crit 1.000e+00) | best - Chebyshev -8.56e-13
  (C) no value above 2*2^(-1/n)  (641.9 s)
(D) EMBEDDING into R^3 (Corollary 1.3, lower bound): ray casting from c over 3000 directions
  n = 3: 1-D value 1.5874010520, R^3 value 1.5874010520, difference +2.7e-15
  n = 4: 1-D value 1.6817928305, R^3 value 1.6817928305, difference -3.6e-15
  n = 5: 1-D value 1.7411011266, R^3 value 1.7411011266, difference -1.0e-14
  n = 6: 1-D value 1.7817974363, R^3 value 1.7817974363, difference +5.8e-15
  n = 7: 1-D value 1.8114473285, R^3 value 1.8114473285, difference +3.9e-13
  n = 8: 1-D value 1.8340080864, R^3 value 1.8340080864, difference -2.6e-12
  (D) ALL PASSED  (642.8 s)
TOTAL: ALL CHECKS PASSED  (642.8 s)
