README.md
4.7 kB · markdown
interval-digits
- The numbers of the half interval on mobius: the digit set
F = {0..fill-1}at the odd basebase = 2 fill - 1, its transfer certificate, the walls it sets at the bars1/5and1/4, a falsification of every inequality behind it, readings of the true rate, and a sanity meter. wall: the constantc_infof the tail boundlambda/fill <= (2/pi) log base + c_inf + 3.4/base, proved from base101and checked directly at every odd base101..3001and at94939,200001,10^6 + 1,10^8 + 1, the least base where the tail bound clears each bar, the closed formlambda/fillcertified at 120 bits at every odd base from the wall to that point and failing at the base below,bar - alpha_1and the gap at the wall, the bars read as the exact intervals1/5and1/4, the lower boundmin G/fill >= (2/pi) log base - 1.31at every odd base9..9999with its tail from101, the density floor1 - alpha_base < 1/5,alpha_1at three larger bases, and the base where the constant of the Section 9 sketch of Maynard 2022 clears1/5on the same set, calibrated against its one-missing-digit crossing1520573.check: the three sums of the certificate's proof, thebP,bQand signedaPsums, against their bounds at every odd base3..401and at1001,10001,100001on201shifts; then at every odd base3..401on801shifts and at five larger bases on401shifts, the closed form ofS(t)against the direct sum, the bound onG(t), the inequalityT phi <= lambda phiand the lower bound onG(t), each asserted; the readingsG(0)/fill - (2/pi) log baseagainst1 + gamma'andmax G/fill - (2 sqrt2/pi) log base; the level sums at 40 digits at nine small bases, four shifts and every level up tobase^level <= 9261against(3/2) lambda^level; and the transfer identity at level 2.rate: power iteration of the transfer operator on a grid of cells with linear interpolation, readings of the true rate near the wall, and the one-step constantG(1/2)/fillagainst both bars.meter: the Mertens meter and the prime count on the half interval at the prime bases101,1009and10007below10^7, a sanity print far below every wall and never evidence for the theorem.
THE ROUNDING
wallevaluates every closed form inmpmathinterval arithmetic at 120 bits; lower bounds print their interval's lower end truncated down, upper bounds the upper end rounded up.checkprints every upper bound rounded up and every lower bound truncated down at its last printed digit; it runs in float64 and asserts strict inequalities whose margins sit far above float error except att = 0, where the bound onGis tight by design to1.5 * 10^-4of its size; the level sums run at 40 digits.rateprints readings: the interpolated power iteration and the one-step scan bound nothing.
RUN
uv run python research/lab/py/interval-digits/interval.py wallin under a second.uv run python research/lab/py/interval-digits/interval.py checkin 13 seconds.uv run python research/lab/py/interval-digits/interval.py ratein about 20 seconds.uv run python research/lab/py/interval-digits/interval.py meterin 2 seconds.- Prints only, writes nothing; every check raises if it fails.
WITNESSES
- The tail constant
c_inf <= 2.60043004, the direct check's smallest gap>= 1.3332 * 10^-7; the wall at the bar1/5is94939, the closed form certified at every odd base94939..94947and the tail bound from94946, failing at94937with gap<= -2.6733 * 10^-5,1/5 - alpha_1 >= 1.3678 * 10^-8at the wall. - The lower bound holds at every odd base from
9, its tail margin>= 0.016518from101, and the factor is negative at7; the density floor1 - alpha_base < 1/5holds at every odd base from27and fails at25. - The wall at the bar
1/4is3789, certified at3789..3793, failing at3787,1/4 - alpha_1 >= 7.9625 * 10^-6at the wall. alpha_1 <= 0.1993872at100003,0.1761232at1000003,0.1331636at10^9 + 7.- The Section 9 constant of Maynard 2022 at the consecutive half interval clears
1/5from7777884825. - The check: the
bPsum under its bound by at least3.99 * 10^-5, largestG/bound0.999857, largestT phi/(lambda phi)0.999866, smallestG/lower1.534638, largest level sum over(3/2) lambda^level0.557408, transfer identity to10^-15. - Readings:
G(0)/fill - (2/pi) log baseis1.962430at100001against1 + gamma' = 1.962523; the true rate reads(2/pi) log base + 2.26near7 * 10^4and crossesbase^(1/5)between70001and80001; the one-step constant reads(2 sqrt2/pi) log base + 1.19and clears1/4from7075and1/5from317063.