research/lab/py/gasket-witness-weights
0 directories and 2 files in research/lab/py/gasket-witness-weights.
gasket-witness-weights
- Reparametrises the gasket residual
R(level)by the witnesszinstead of by the multiplier pair(s,t): every off-diagonal collinear pair ofG_levelis(sz, tz)withgcd(s,t) = 1andza positive integer vector, uniquely. - Checks the free-digit automaton
B(s,t)is a constrained tensor square of a single-coordinate carry automatonC(s,t):T = S (x) S - U (x) U - V (x) V + W (x) W, soBneeds no four-tuple state graph. - Checks the box construction
P_level(s,t) = #{z >= (1,1) : z_1 + z_2 <= (3^level-1)/(2 max(s,t)), sz, tz in G_level}againstB(s,t)and against the tensor iteration at multipliers where a forward build ofBblows up. - Checks the witness weight
w = z_1 + z_2is never below 4, that the weight layers scale exactly asR_{3w}(level) = R_w(level-1), and that the weight-four layer is counted by the no-adjacent-ones setF_levelwith|F_level| = Fib(level+1) - 1. - Checks every multiplier pair above
(3^level-1)/10contributes exactly 4 ordered collinear pairs. - Checks the golden ceiling
M_level(z) <= Fib(level+1) - 1on every coprime direction in a box, at every level to 40, withz = (1,3)the only attainer, and on six families chosen to favour a breach at every level to 45. - Carries the second moment
E(level)and the residualR(level)to level 17 by grouping all3^levelpoints by primitive direction, and checks bothR/3^levelandR/phi^(2 level)fall at every level through 17. - Rebuilds the ray automaton in the direction coordinate: a multiplier word is a word over the increments
{0, z_2, -z_1}summing to zero, soM_level(z) + 1counts the closed paths of a carry automaton whose states live in[-z_1/2, z_2/2]. - Checks the branch structure that turns the golden ceiling into a theorem: out-degree at most two, branch states in one residue class mod 3, and the two successors of a branch state differing by
q/3whereqis the unique member of{z_1, z_2, w}divisible by 3. - Checks the proved cases over the box and over the de-duplicated union of the six adversarial families, names the directions left over and runs them to level 60, and verifies nine explicit Fibonacci certificates in exact integer arithmetic.
- Checks the renewal criterion
Sum_{j>=2} f_j Fib(level+1-j) <= Fib(level-1)on the first-return counts of every occupied direction of the box. - Solves the golden potential in exact
Q(sqrt5)arithmetic:u(0) = 1andphi u(c) = sum of u over the successors of cfor every livec != 0, read off asU(z) = sum of u over the successors of the start state other than itself. - Checks the criterion
U(z) <= phi^-2, which impliesM_level(z) <= Fib(level+1) - 1at every level, on the box, the six adversarial families and a stress list of highv_3directions, and confirms every solution against both inequalities of the criterion rather than only against the linear system. - Checks the occupancy residue law: writing
q = 3^k q1for the coordinate divisible by 3 andpfor the other, a direction carrying mass at any level hasq1 = p mod 3, so the congruence alone empties a large share of the directions with3 | z1 z2. - Checks the short first returns: no first return has length between 2 and
v_3(q);f_2is nonzero only at{1,3}andf_3only at{1,9},{1,12},{3,10},{4,9}, each equal to 1. - Checks the degree potential:
pi(c) = 1where a live state has out-degree two,phi^-1where it has out-degree one,pi(0) = 1, is a super-solution of the golden criterion, and sweeping it under the same operator gives a decreasing chain of exact upper bounds onU; the sweep runs over consecutive depths, so the depth reported is the least one that works. - Checks the burst identity that blocks every valuation-graded potential:
Uisphi^-(k-1)times the sum ofuover the2^(k-1)burst-floor states, all of valuation 0, whileu(p) = phi^-1at the valuation-0 statep. - Checks the golden partition bound at
v_3(q) = 1:U <= phi^-1 (1 - phi^-max(t,2))witht = v_3(q1 - p), henceU <= phi^-2on the arithmetic classt <= 2. - Checks the two automaton-free restatements of the criterion:
Sum_level (M_level + 1) phi^-level <= phi^4, andSum_m phi^-l(m) <= phiover the multipliersmof the direction, wherel(m)is the number of base-3 digits of(z1 + z2) m. - Checks the route that dies: the state maximum
G(level) = max_c N(c,level)does not obeyG(level) <= G(level-1) + G(level-2).
RUN
uv run --with numpy python research/lab/py/gasket-witness-weights/witness.py- About 165 seconds on one core, peak under 3GB at level 17; prints one line per law and exits nonzero on the first failure.
- numpy is used only for the level 15 to 17 array sort; every other law is stdlib.
WITNESSES
- 473 coprime pairs below 40: the tensor square reproduces the
B(s,t)return counts at every level to 9, zero mismatches. - Carry states against reachable
Bstates:(365,1094)729 against 26931,(41,122)81 against 835,(25,52)38 against 393,(31,40)35 against 354. - 812 coprime pairs: the box construction agrees with
B(s,t)at level 9, zero mismatches; it agrees with the tensor iteration at(365,1094),(41,122),(122,123),(1,2460)and(2431,2458)at level 9 and level 12. R(level) = 20, 88, 432, 1624, 5512, 15896, 46064, 124928, 335704, 863848for level 4..13.R_{3w}(level) = R_w(level-1)on all 1869 weight layers divisible by 3 across level 5..13, zero failures.- No witness of weight below 4 at any level to 13.
R_4(level) = 12, 36, 108, 336, 988, 2596, 6672, 17480, 45720for level 4..12, each equal to twice the number of ordered coprime non-3-power-ratio pairs drawn fromF_level, and|F_level| = Fib(level+1) - 1at each of those levels.- Pairs above
(3^level-1)/10number 18, 57, 163, 402, 1019, 2702, 7060, 18607 at level 6..13; every one contributes exactly 4. M_level(1,3) = Fib(level+1) - 1forlevel <= 40; over 13158 coprime directions withz_1 <= 120andz_1 <= z_2 <= 240there is nolevel <= 40and nozwithM_level(z) > Fib(level+1) - 1, and(1,3)is the sole attainer at level 40.- Six adversarial families at every
level <= 45, coprime counts 16940, 253, 2998, 1499, 1199, 1199, total 24088, zero breaches: binary base-3 pairs below3^8, no-adjacent-ones pairs below3^7,(1,t)witht < 3000, consecutive below 1500,(s,3s-1)and(s,3s+1)withs < 1200. The shelf scriptgasket-ray-machine/scripts/verify.pyships the same six with the binary family at3^7, 11369 directions, to stay inside its time budget. E(level) = 4003372, 11679626, 34050692, 99800950, 292848756andR(level) = 863848, 2211960, 5549452, 14100688, 35354824for level 13..17.R(level)/3^levelpeaks at0.8401158at level 8 and falls at every level to0.2737709at level 17;R(level)/phi^(2 level)peaks at3.2378233at level 12 and falls at every level to2.7724831; the level ratioR(level+1)/R(level)reads2.5073119at level 17, belowphi^2 = 2.6180339.- The direction automaton reproduces the gasket-digit automaton on all 947 coprime pairs below 40 at every level to 20, zero mismatches.
- Of the 13158 coprime
zwithz_1 <= 120andz_1 <= z_2 <= 240, exactly 218 have a live automaton beyond the start state; every one has out-degree at most two, one branch class, and carries inside[-z_1/2, z_2/2], the largest live set 37 states. The branch argument settles 206 of them, 107 of which havev_3(q) = 1. The twelve left are(1,9),(1,27),(1,81),(1,90),(4,117),(9,73),(9,82),(9,235),(10,81),(13,108),(27,217),(27,226). - The first three of those twelve are shift rays, closed by
Fib(p+2) Fib(q+2) = Fib(p+q+3) - Fib(p+1) Fib(q+1); the other nine carry Fibonacci certificates of denominators 18, 40, 381, 18, 2013, 18, 40, 2013, 34. - The six adversarial families overlap: their 24088 coprime members are 23435 distinct directions, of which 717 lie in the box and 22718 are new. Of the 22718: 20945 have zero mass, 1693 fall to the branch argument, 3 are shift rays, and 77 are left to the enumeration. Those 77 hold at every level to 60, zero breaches, worst ratio to the ceiling below
0.1516. - The renewal criterion holds on all 218 occupied directions of the box at every level to 46;
f_1 = 1everywhere andf_2 = 1only at(1,3). M_level(z) <= D_level(w)on all 829 coprime directions withz_1 <= 30,z_1 <= z_2 <= 60at level 12; the bound is far weaker than the ceiling,D_24(w) = 4196351, 1683971, 613817, 228519atw = 4, 10, 28, 82againstFib(25) - 1 = 75024.- The state maximum at
(1,9)runs1, 1, 1, 2, 4, 6, 9, soG(4) = 4 > G(3) + G(2) = 3; the recursion fails for 8 directions of the box. - The golden potential is nonnegative and satisfies both criterion inequalities on all 218 occupied directions of the box; 214 have
U <= phi^-2and the four failures are the shift rays(1,3),(1,9),(1,27),(1,81), each withU = phi - 1exactly. Utakes 57 distinct values over the box, in orderphi - 1on the four shift rays,2 - phiat(1,12),(3,10),(4,9), then(9 phi - 14)/2at(1,90),(9,82),(10,81).- Over 36037 directions - the box, the six families and a stress list of high
v_3(q)directions - 1995 are occupied and 1987 satisfyU <= phi^-2; the 8 failures are exactly the shift rays(1,3^j)forj = 1..8, no direction hasUin the open interval(2 - phi, phi - 1), the maximum among the passing directions is2 - phiattained only at(1,12),(3,10),(4,9), and the largest live set is 256 states. - Of the 13158 coprime box directions, 6566 have
3dividing neither coordinate and carry no mass at any level; 6374 more have3dividing a coordinate and still carry none; 218 are occupied. - Adversarial cross-check on all 218:
f_1 = 1, the exact partial sumSum_{j=2}^{46} f_j phi^-jfrom the independent first-return enumeration is dominated byphi^-1 Ufrom the linear solve, andU <= phi^-1at every direction, with equality only on the shift rays. - Of the 36037 directions, 20193 have
3 | z1 z2and 15556 of those matchq1 = p mod 3; all 1995 occupied directions match, zero violations, so the congruence empties 4637 directions with no automaton built. - No occupied direction has a first return of length between 2 and
v_3(q), zero failures over 1995. - The burst identity
U = phi^-(k-1) Sum_m u(q1 m)over the2^(k-1)burst-floor states, andu(p) = phi^-1whenever2p <= q, hold exactly on all 1995. - The degree potential is a super-solution on 1968 of the 1995 occupied directions, 1902 of which have no double branching and 66 of which do, so the tool reaches strictly past the branch case.
- Sweeping it settles
U <= phi^-2on 1966 directions - least depth 1 on 1804, 3 on 101, 4 on 44, 5 on 11, 6 on 6 - and leaves 29: the 8 shift rays(1,3^j)and(1,756),(1,2196),(1,2214),(1,2268),(1,2430),(9,2188),(10,2187),(13,1080),(13,3267),(27,730),(27,2188),(28,729),(28,2187),(40,1053),(81,2188),(82,2187),(91,2214),(121,3159),(243,2188),(244,2187),(819,2539). - 757 occupied directions have
v_3(q) = 1; the boundU <= phi^-1 (1 - phi^-max(t,2))holds at every one, is attained exactly at(1,12)and(3,10), and 512 of them fall in the classt <= 2where it givesU <= phi^-2with zero failures. - The two automaton-free restatements agree with the linear solve on 111 directions with
z1 < 40,z2 < 120, checked in exactQ(sqrt5)arithmetic.
- README.md11.4 kB
- witness.py40.2 kB