sg het (.py)

Computation

The computation scripts

Each fermion-sector paper carries a computation script that performs the scans behind its result: it builds the prism and the metric edge operator from scratch, forms the pieces of a chain as vectors on the edges, and searches the assignments the paper's rules allow. The scripts are the source of the numbers; the verification scripts then check the reported values against them.

Running them

Each script is one file, with no sibling imports and no external data, and runs in a single cell. Requirements are Python 3.8+ with NumPy, and Matplotlib only for the figure options. Run with no argument for the solution of the paper; the scan options reproduce the search tables and take longer.

papercomputation scriptno argumentoptions
sg hetcompute_het_sg.pythe prism quantities, the summary of results, the three neutrino Yukawa couplings and the selection over cell factors—

Table 1. The computation scripts. Each is one file with no sibling imports and no external data, and runs in a single cell.


What they share

It builds these objects before anything else, so each stands on its own:

stepwhat is built
the prismthe 16 vertices v = (x, y, z, τ) and the 32 edges, numbered by the lower vertex read as a binary number and, within a vertex, in the order x, y, z, τ
the metric operatorthe weights ge, the phase S, the gradient A = S|g|1∕2d0, the Laplacian Lω = A†A, the three edge adjacencies Iz, Dxy, Iτ, the bond operator ME and its inverse GE
the pieceseach edge, xy loop, τ face and staircase as a vector on the 32 edges, signed by traversal
the linksℓ1 = ⟨b1|GEA|vo⟩, ℓk+1 = ⟨bk+1|GE|bk⟩, ℓn+1 = ⟨vs|PHA†GE|bn⟩, and the amplitude as their product
the referencesthe neutrino masses from the NuFIT splittings with the DESI sum, and the MZ reading surface from the self-consistent root of MZ = (v∕2)√(gL2 + gY2)

Table 2. The construction every script builds before its own search.