research/lab/rs/spin-census
1 directory and 2 files in research/lab/rs/spin-census.
Spin Census
- Runs the four unrun questions of spin plus the Gaussian Farey, and prints only.
- The spin mass
M(r)is read twice: exactly, as the integer shell histogram of filled cells about a fixed point, and throughmrlynum::spin::profilewithmass_within, the two agreeing to1.7e-6on the total and to0.5%on partial radii. The hole about the raster centre is measured on cell rectangles in exact integer arithmetic, never on cell centres, which would returnhole + 1/2whatever the hole. - Every fixed point is the attractor of one filled digit
d, the pointd/2of the unit square, about which the design is exactly self-similar; the corner digit gives the widest window,r <= side, the centre digit onlyr <= side/2. - The ripple is the residual
ln M(r) - (log(fill)/log 3) ln rtaken exactly, never fitted, folded into 24 bins oflog_3 rover a whole number of periods so every bin carries the same sample count. Its drift bar is the same fold on the first half of the window against the second, computed from the code itself. - Level 6 and level 7 for the census over all 256 codes with a filled corner digit, de-duplicated to transpose classes because the corner ripple is a class function; the transpose control asks the same code and its mirror for the same ripple and gets it to the last bit. The window's whole periods are counted by repeated multiplication, never by a logarithm, which floors to the wrong period at level 8.
- The powder is the arithmetic ring average of
|F(k)|^2over 240 logarithmic bins, band3 pad/sidetopad/8in frequency index, at level 7 with pad 4096 and again with pad 8192 on three codes; the instrument spread is a three-period fit window slid a quarter period at a time across the band, not a split of the band in two, which understates it by an order. - The spin spectrum is
mrlynum::spin::harmonicsat levels 1 and 2 over all 511 nonempty codes, ordersm = 0..12, compared against the 101 orbits of the square group. - The shape reading reduces the spin spectrum to exact integers. At
leveleveryP_mis a quadratic form in the indicators of the render's cells whose Gram matrix commutes with the raster's symmetry group, soP_mis a linear functional of the pair censusPhi_level: the number of filled cell pairs in each orbit of that group on pairs. Equal pair censuses force equalP_mat every order, every ring count and every truncation, in exact arithmetic and with no transform run. - The pair census is the whole instrument at level 1: base 3 plane has 11 classes, base 5 plane 55, base 3
dim 324 under the order-48 cube group; at level 2 the render is the Kronecker square and the counts are 461 and 24805. - The control is the full 511. Codes de-duplicate to orbits by canonical form and the census reports the number of distinct
Phi_1andPhi_2values against the Burnside orbit count, which the same pass computes from the cycle index and never from the sweep. - The reduction is checked against the instrument rather than assumed: the level-1 coefficient matrix is solved from independent censuses by elimination, then predicts
mrlynum::spin::harmonicsat 1024 rings andm = 0..12on all 511 codes, printing the worst residual and the rank of the 13 orders against the 11 classes. It holds for a constant-valued0/1render on a raster of sidebase^level; the spectrum count it is compared against is a greedy first-match bucketing at tolerance1e-9. - Base 5 and
dim 3runPhi_1in full over all2^25and2^27codes, storing canonical representatives only, and runPhi_2over the colliding groups inside a budget-capped weight window the generator prints; the cap did not bind and the window covered every group. - The sponge shadow counts lattice lines in direction
(a,b,c)meeting the sponge atlevel, as classes of filled cells underx -> x cross v, with the solid cube in the same direction as the exact ceiling. - The Gaussian Farey counts radii new at scale
nthree ways: the direct union over reduced squared radii, the square-free rule, and the Mobius identity over the radical.
RUN
CARGO_BUILD_JOBS=4 cargo run --release -p spin-census- About thirty seconds; prints only, writes nothing. Peak memory is the pad-8192 transform, about 1.1 GB.
CARGO_BUILD_JOBS=4 cargo run --release -p spin-census -- shaperuns the shape reading, the one verb group the default pass leaves out.- About two minutes and 0.4 GB, the peak being the base-5 sweep holding 4.2 million packed censuses; the
dim 3sweep visits all2^27codes and keeps 2.85 million.
WITNESSES
- spin.md the spin dimension rows: slopes
1.465054, 1.649432, 1.783588, 1.761814, 1.897854, 1.879522, 2.000100about the corner at level 6 against the exactlog(fill)/log 3. - spin.md the centre hole: codes 239 and 495 have first occupied radius
41.000againstside/6 = 40.5at level 5, mass zero inside. - spin.md the acid test:
127against239ripple gap0.11984on drift bar0.04126;255against495gap0.12042on bar0.01461; the solid 511 ripple flat at0.00277under its own bar0.00578. - spin.md the ripple census: 256 codes read at levels 6 and 7, transpose control
0.00e0, 13 equal-fill transpose class pairs inside their own bar at level 6 and 6 at level 7 with ratios0.71to0.95, closest7and273at0.01371with swings0.01114and0.01217. - spin.md the powder: pad 4096 slopes
-1.37986, -1.51762, -1.73012, -1.83071, -1.97886, -2.00433, slide spreads0.1652to0.4405, pad 8192 moving127, 255, 495to-1.81607, -2.03225, -2.12289; the fractal slide floor-2.27308against the solid's-2.75781. - spin.md the centre hole: four times the squared distance to the nearest filled cell
6561 = (side/3)^2for239and495at level 5,12802for79. - spin.md the spectrum: 0 isospectral pairs, 101 distinct spectra against 101 nonempty square classes.
- spin.md the shadow:
(0,0,1)sponge8, 64, 512, 4096, 32768against cube9, 81, 729, 6561;(1,1,1)sponge and cube both19, 217, 2107, 19441, 176419. - spin.md the Gaussian Farey:
2, 3, 9, 11, 22, 18, 40, 38, 55, 52, 91, 64, 123, 97, 128, 126, 199, 136, 243, 180, the Mobius identity ton = 64, the radical-six ratios0.56250to0.63801. - spin.md the shape reading at base 3: 11 pair classes at level 1 and 461 at level 2, 97 level-1 censuses against 101 orbits, the homometric pairs
45-105,61-121,78-102,94-118at level-1 spectrum gaps1.30e-16,1.03e-17,6.51e-17,1.64e-16and level-2 gaps0.151,0.0689,0.253,0.105;P_mread at level 1 alone, bucketed greedily at1e-9, gives 97 spectra, the count the census predicts. - spin.md the reduction: the level-1 coefficients solved from 11 independent censuses reproduce
mrlynum::spin::harmonicson all 511 codes at worst relative residual1.14e-14, coefficient rank 9 of 11, the six odd orders at the cap 3 that half-turning one member of a pair forces and the seven even orders at 6 of a possible 8. - spin.md the completeness away from base 3:
3993511level-1 censuses on4211743base-5 orbits with204856ties over423088orbits, and1461693on2852287cube orbits with757066ties over2147660; every tie broken at level 2, the window covering every group out to weights 21 and 24, with no pair surviving, and the canonical counts matching the Burnside averages4211744and2852288. mrlynum::spin::mass_withinand the crate testthe_spin_mass_scales_by_the_fill_about_a_filled_corner; the window regression testthe_window_counts_whole_periods_at_every_level.
- src/8 items
- Cargo.toml174 B
- README.md7.5 kB