# Alignment test of the closure constant n* = 1/2 + 2*rho_mass/log(3)
# Exact-arithmetic residuals Delta = n*_fit - (1/2 + 2*rho_mass/log 3)
# Executed with mpmath at 50 digits on 2026-10-31.
from mpmath import mp, mpf, log
mp.dps = 50
ln3 = log(3)
nstar_fit = [mpf("5.636"), mpf("5.680"), mpf("5.683"), mpf("5.686")]
rho_mass  = [mpf("2.830"), mpf("2.841"), mpf("2.845"), mpf("2.847")]
rho_11    = mpf("2.8476")
for a, b in zip(nstar_fit, rho_mass):
    print("Delta = %+.6f" % (a - (mpf('0.5') + 2*b/ln3)))
print("last point with rho_mass(11) = 2.8476: Delta = %+.6f"
      % (nstar_fit[3] - (mpf('0.5') + 2*rho_11/ln3)))