Characteristic polynomial: (4*lambda**3 - 23*lambda**2 + 38*lambda - 17)/4
H =
 [[ 1.   -0.5   0.  ]
 [-0.5   1.75 -0.5 ]
 [ 0.   -0.5   3.  ]]
spec(H) = [0.726377103195317 1.834876148944763 3.18874674785992 ]
vacuum vector = [0.872386059638559 0.477409601540651 0.104988738944263]
P_vac =
 [[0.761057437051691 0.416485481121663 0.091590712274007]
 [0.416485481121663 0.227919927643203 0.050122632025636]
 [0.091590712274007 0.050122632025636 0.011022635305107]]
restricted electric eigenvalues = [0.576371368409391 1.980643415247995]
s_IR = [-0.482104205973867  0.875798927039865  0.023489869691604]
u_UV = [ 0.080734735588944  0.071107747486603 -0.994195951871041]

Verified constants
E0                   = 0.726377103195317
E1                   = 1.83487614894476
E2                   = 3.18874674785992
mu_IR                = 0.576371368409391
mu_UV                = 1.980643415248
sigma_vac            = 0.878680870599103
tau_capture          = 0.0711077474866026
J                    = 0.875798927039865
G                    = 0.767023760604179
M_phys               = 1.40165394364111
E_phys               = 0.576371368409391
L_tilde              = 2.60748115341905
b_max                = 0.53755093945864
b                    = 0.5
D_mag                = 0.0979133669315835
R_cal                = -0.052092367854342
r_sec                = 0.052092367854342
c_IR                 = 0.947907632145658
c_UV                 = 1.25426631205268
epsilon_mix          = 0.416309692765243
det_margin           = 1.01561484964763
kappa                = 0.657490679870655
actual_gap           = 1.10849904574945
certified_fraction   = 0.593135990862429
