# Black Hole Lifecycle: Singularity Elimination Results
# Generated by blackhole_lifecycle.py (stdout redirect)
#
# DESI w = -0.85
# K/V = 0.08108
# V0 = 5.973180e-27 kg/m^3 = 0.7 rho_crit
# rho_crit = 8.533115e-27 kg/m^3
# c_g is NOT bounded by MICROSCOPE. MICROSCOPE is a differential test
# between two test masses in the same field, so a coupling universal
# across materials cancels in the ratio. Only the composition-dependent
# part is bounded, and that part is zero by construction.
# c_g is therefore unbounded from above; values below survey the family.
# illustrative c_g for the integration = 1.0e+00
#
# TIME-GRADIENT FIELD  G = 1/(2x sqrt(1-1/x)),  x = r/r_s
# r/r_s	G
10000.000000	5.000250e-05
1000.000000	5.002502e-04
100.000000	5.025189e-03
10.000000	5.270463e-02
5.000000	1.118034e-01
2.000000	3.535534e-01
1.500000	5.773503e-01
1.100000	1.507557e+00
1.010000	4.975186e+00
1.001000	1.580349e+01
1.000100	4.999750e+01
#
# FREEZE-OUT  x = (1 + sqrt(1 + c_g^2))/2
# c_g	r_freeze/r_s	outside_by
1.000e-04	1.00000000	2.500000e-09
1.000e-03	1.00000025	2.499999e-07
1.000e-02	1.00002500	2.499938e-05
1.000e-01	1.00249378	2.493781e-03
5.000e-01	1.05901699	5.901699e-02
1.000e+00	1.20710678	2.071068e-01
2.000e+00	1.61803399	6.180340e-01
5.000e+00	3.04950976	2.049510e+00
#
# illustrative c_g = 1.0e+00 (family survey, not a measurement)
# r_freeze = 1.20710678 r_s  (2.071e-01 r_s outside horizon)
# G at freeze-out = 1.0000e+00
# x_freeze > 1 for ALL c_g > 0. No standoff distance is predicted.
# numerical free-fall stalled at r = 1.20710678 r_s
#
# BOUNCE  rho_G + 3P_G = -2 V(phi) < 0 for all V0 > 0
#
# REMNANT
# rho_rem = 2 V0 / c^2 = 1.329145e-43 kg/m^3 = 1.5576e-17 rho_crit
# m_eff = 2.577930e-24 kg = 1.184711e-16 M_Pl
# lambda_C = 1.364499e-19 m
# remnant mass: not computed. The remnant is the time-information
# content of what collapsed, preserved in full; a mass is a classical
# description of that content, not a quantity the remnant carries.
#
# CHECKS: 12/12 PASS
# [PASS] Background field from DESI energy split :: K/V = 0.0811, V0 = 0.7 rho_crit
# [PASS] Time-gradient field diverges at r_s :: G(1.0001 r_s) = 5.00e+01, unbounded at 1 r_s
# [PASS] Freeze-out solved in closed form :: x = (1 + sqrt(1 + c_g^2))/2
# [PASS] Freeze-out always outside horizon :: min margin 2.50e-09 r_s over c_g in [1e-4, 5]
# [PASS] Horizon never crossed in integration :: stalled at r = 1.20710678 r_s
# [PASS] Bounce condition satisfied for all V0 > 0 :: rho_G + 3P_G = -2V(phi) < 0
# [PASS] MICROSCOPE null predicted by construction :: composition-dependent coupling = 0 (buoyancy_origin.py); a universal c_g*G is common mode and cancels
# [PASS] Standoff distance not predicted :: c_g unbounded from above; standoff is an open input, not a result
# [PASS] Field built from clock observables :: dtau/dt from the local field matches the correction law at every x
# [PASS] Far field recovered, not assumed :: Gcal -> 1/(2x) -> 0 as x -> infinity
# [PASS] Bounce is a clock reading :: dtau/dt falls into the stall then returns to unity below it
# [PASS] Remnant density from the balance law :: rho = 2V0/c^2 = 1.329e-43 kg/m^3 = 1.56e-17 rho_crit; 1/c^2 restored, agrees with the paper's derivation
#
# HORIZON NEVER CROSSED FOR ANY c_g > 0. NO SINGULARITY.