# orbit_drag sweep  N=48  basis/constraints identical to weil_form.py
# lam_full: all n=p^k<e^2L | lam_primes: k=1 only | lam_kle2: k<=2 | drag = lam_primes - lam_full
#     L  #pp  #pr    lam_arch   lam_full  lam_primes   lam_kle2       drag  k>=2 orbits active
   0.60    2    2   -0.000385  +0.000078   +0.000078  +0.000078  +0.000000  []
   0.65    2    2   -0.081888  +0.000003   +0.000003  +0.000003  +0.000000  []
   0.70    3    2   -0.157622  +0.000000   +0.000000  +0.000000  +0.000000  [4]
   0.75    3    2   -0.228417  +0.000000   -0.055546  +0.000000  -0.055546  [4]
   0.80    3    2   -0.294947  +0.000000   -0.095220  +0.000000  -0.095220  [4]
   0.85    4    3   -0.357767  +0.000000   -0.143677  +0.000000  -0.143677  [4]
   0.90    4    3   -0.417339  +0.000000   -0.261406  +0.000000  -0.261406  [4]
   0.95    4    3   -0.474052  +0.000000   -0.315976  +0.000000  -0.315976  [4]
   1.00    5    4   -0.528239  -0.000000   -0.329369  -0.000000  -0.329369  [4]
   1.05    6    4   -0.580186  -0.000000   -0.337204  -0.000010  -0.337204  [4, 8]
   1.10    7    4   -0.630142  -0.000000   -0.347832  -0.101605  -0.347832  [4, 8, 9]
   1.15    7    4   -0.678326  -0.000000   -0.433717  -0.142180  -0.433717  [4, 8, 9]
   1.20    8    5   -0.724929  -0.000000   -0.476465  -0.174086  -0.476465  [4, 8, 9]
   1.25    8    5   -0.770121  -0.000000   -0.561638  -0.176786  -0.561638  [4, 8, 9]
   1.30    9    6   -0.814054  -0.000000   -0.635663  -0.188280  -0.635663  [4, 8, 9]
   1.35    9    6   -0.856865  -0.000000   -0.674543  -0.225266  -0.674543  [4, 8, 9]
   1.40   10    6   -0.898675  -0.000000   -0.700660  -0.242254  -0.700660  [4, 8, 9, 16]
   1.45   11    7   -0.939596  -0.000000   -0.711506  -0.287302  -0.711506  [4, 8, 9, 16]
   1.50   12    8   -0.979728  -0.000000   -0.714161  -0.296509  -0.714161  [4, 8, 9, 16]
   1.55   12    8   -1.019161  -0.000000   -0.731221  -0.325349  -0.731221  [4, 8, 9, 16]
   1.60   13    9   -1.057978  -0.000000   -0.747665  -0.335619  -0.747665  [4, 8, 9, 16]
   1.65   15    9   -1.096255  -0.000000   -0.752133  -0.339878  -0.752133  [4, 8, 9, 16, 25, 27]
   1.70   16   10   -1.134059  -0.000000   -0.756376  -0.349083  -0.756376  [4, 8, 9, 16, 25, 27]
   1.75   18   11   -1.171455  -0.000000   -0.764622  -0.370773  -0.764622  [4, 8, 9, 16, 25, 27, 32]
   1.80   18   11   -1.208497  -0.000000   -0.808425  -0.378733  -0.808425  [4, 8, 9, 16, 25, 27, 32]
   1.85   19   12   -1.245240  -0.000000   -0.906734  -0.386818  -0.906734  [4, 8, 9, 16, 25, 27, 32]
   1.90   21   14   -1.281728  -0.000000   -0.970274  -0.423458  -0.970274  [4, 8, 9, 16, 25, 27, 32]
   1.95   23   15   -1.318006  -0.000000   -1.005008  -0.445954  -1.005008  [4, 8, 9, 16, 25, 27, 32, 49]
   2.00   24   16   -1.354111  -0.000000   -1.011218  -0.465515  -1.011218  [4, 8, 9, 16, 25, 27, 32, 49]
# grid done in 49s

# lam_primes zero crossings (bisection to ~0.004 in L)
#   bisect L=0.7250  lam_primes=-0.015676
#   bisect L=0.7125  lam_primes=-0.000022
#   bisect L=0.7062  lam_primes=-0.000001
#   bisect L=0.7031  lam_primes=+0.000000
  crossing: lam_primes changes sign in L=[0.7031, 0.7062] (+0.000000 -> -0.000001)

# single-orbit ablations: lambda_min of the FULL sum minus one orbit
#     L        full   drop n=4   drop n=8   drop n=9  k<=2 only     primes
   0.80   +0.000000  -0.095220         --         --  +0.000000  -0.095220
   1.20   -0.000000  -0.345944  -0.174086  -0.185453  -0.174086  -0.476465
   1.60   -0.000000  -0.468974  -0.245011  -0.365249  -0.335619  -0.747665
   2.00   -0.000000  -0.490124  -0.245065  -0.366204  -0.465515  -1.011218

# N=64 convergence check (same quantities, larger basis)
#     L   lam_full48  lam_full64  lam_prim48  lam_prim64
   0.80    +0.000000   +0.000000   -0.095220   -0.095704
   1.40    -0.000000   -0.000003   -0.700660   -0.701190
   2.00    -0.000000   -0.000000   -1.011218   -1.011384
