research/lab/py/spectra-from-cuts
0 directories and 2 files in research/lab/py/spectra-from-cuts.
spectra-from-cuts
- Composes the two-tile claim of
notes/spectra.mdwith the carry automaton ofnotes/cuts.mdrun at an odd basebinstead of base 3. - The objects:
P_bis the digit polynomial of the solid with at most one odd coordinate,P_b(t) = E^2 (E + 3 O)withE = 1 + t^2 + .. + t^(b-1)andO = t + t^3 + .. + t^(b-2);g = 3(b-1)/2is the target digit; the automaton isM[c, c'] = P_b[c + g - b c']onabs(c) <= 1. blockrecomputes the even block[[M[0,0], M[0,1]], [M[1,0] + M[-1,0], M[1,1] + M[-1,1]]]at every oddb = 3..25, checks that no carry leavesabs(c) <= 1, that the block equals the closed form spectra prints for the class ofb mod 4, that(hexagons, triangles) = ((M^n)[0,0], (M^n)[1,0] + (M^n)[-1,0])obeys the block at levels0..8, and that a direct digit-sum census at heightsh_n - 1, h_n, h_n + 1,h_n = 3(b^n - 1)/2, returns the same pairs at low levels.formsderives the five coefficientsP_b[g], P_b[g+1], P_b[g+b], P_b[g+b-1], P_b[g+b+1]as polynomials inbwithin each class ofb mod 4, by inclusion-exclusion on((1 - x^n)/(1 - x))^3and(1 - x^(n-1))(1 - x^n)^2/(1 - x)^3withn = (b+1)/2, prints the base from which each extraction's active terms are stable, and checks every form against the digit polynomial at oddb = 3..101, which covers every base below the thresholds.splittakes the trace and determinant of the block in each class, and provesrho_b > fill/bforb = 3 mod 4andrho_b < fill/bforb = 1 mod 4as polynomial inequalities inn: real roots isolated exactly, all below the firstnof the class, and the sign read at thatn.classeschecks that the residue-class sums ofP_bmodulobare the row sums of the automaton, prints them as polynomials inb, and checksP_b(w) (1 + w)^3 + 2is divisible by1 + w + .. + w^(b-1), soP_b(w) = -2/(1 + w)^3at everyb-th root of unityw != 1.- Exact arithmetic throughout, integers and sympy; the whole run takes under one second.
RUN
uv run python research/lab/py/spectra-from-cuts/compose.py
uv run python research/lab/py/spectra-from-cuts/compose.py block
uv run python research/lab/py/spectra-from-cuts/compose.py forms
uv run python research/lab/py/spectra-from-cuts/compose.py split
uv run python research/lab/py/spectra-from-cuts/compose.py classes
WITNESSES
spectra.mdTHE CLAIM, the matrices and the closed form: verbblock, twelve odd bases, no mismatch.spectra.mdTHE COMPOSITION, the five coefficients and the block in each class: verbforms, run insidesplit.spectra.mdTHE COMPOSITION, the sign ofrho_b - fill/bin each class: verbsplit.spectra.mdTHE COMPOSITION, the residue-class sums and the root-of-unity value: verbclasses.cuts.mdThe ladder in every dimension, the base axis atdim = 3: verbblock.
- compose.py10.8 kB
- README.md2.8 kB