Skip to main content

Follow the source behind the result.

Search the paper sections, open their original context and follow the checks they actually cite. The knowledge atlas explains the result; this library locates its source.

v33_explicit_flat_bundle.py

verification/v33_explicit_flat_bundle.py

This is a registered verification module. Read its original code and the paper sections that cite it, then follow the linked claim to its present classification.

Read the publication code on this page

The registered Python source is shown as text. Displaying it does not execute the code.

"""v33 -- explicit realisation of a valid (U_wall) flat bundle (RH solve output).

The numerical Riemann-Hilbert solve of v32's residual produced an EXPLICIT
Fuchsian residue A0 on P^1 \ mu4 (A_k = U^{k-1} A0 U^{-(k-1)}, U=diag(1,i,-i))
that simultaneously realises ALL the (U_wall) conditions:

  * eigenvalues of A0 = {0, 1/3, 2/3}        (cusp class at each puncture)
  * diag(A0) = (1/2, 1/4, 1/4)               (=> exponents at inf (2,1,1)
                                                => splitting O(-2)(+)O(-1)^2)
  * the path-ordered monodromy at infinity is trivial, ||M_inf - I|| ~ 1e-7,
    i.e. M_1 M_2 M_3 M_4 = I                  (valid global flat bundle)
  * the representation is IRREDUCIBLE          (CASE A: non-trivial SU(3)_F
                                                mixing -- not the diagonal case B)

So an explicit valid flat bundle with the TFPT-required structure EXISTS and is
constructed here (no fabrication: A0 is the RH-solve output, re-Hermitianised).

HONEST SCOPE: this is one explicit point realising the configuration (it confirms
existence + case A). It is likely one of a residual family; pinning the UNIQUE
physical nabla_F* still needs the selector det R=8, and reading c_u/c_d off it
still needs the C6<->P1\mu4 dictionary (H2).  This step delivers the explicit
flat bundle, NOT c_u/c_d.
"""
import numpy as np
import scipy.linalg as sl
from scipy.integrate import solve_ivp
from tfpt_constants import check, summary, reset

U = np.diag([1, 1j, -1j])
Ui = [np.linalg.matrix_power(U, k) for k in range(4)]
Uic = [np.linalg.inv(u) for u in Ui]
p = np.array([1, 1j, -1, -1j])
I3 = np.eye(3, dtype=complex)

# explicit RH-solve output (6-digit), re-Hermitianised below
_A0 = np.array([
    [0.5,                 0.186336 + 0.144342j, 0.0],
    [0.186336 - 0.144342j, 0.25,                -0.184641 + 0.025102j],
    [0.0,                 -0.184641 - 0.025102j, 0.25],
], complex)


def Minf(A0, R=5.0):
    Ak = [Ui[k] @ A0 @ Uic[k] for k in range(4)]
    def rhs(t, y):
        z = R * np.exp(2j * np.pi * t)
        A = sum(Ak[k] / (z - p[k]) for k in range(4))
        return ((A @ y.reshape(3, 3)) * (2j * np.pi * R * np.exp(2j * np.pi * t))).reshape(-1)
    return solve_ivp(rhs, [0, 1], I3.reshape(-1), rtol=1e-9, atol=1e-11, method='DOP853').y[:, -1].reshape(3, 3)


def run():
    reset()
    print("v33  explicit (U_wall) flat bundle (RH-solve output)")
    A0 = (_A0 + _A0.conj().T) / 2     # re-Hermitianise the 6-digit output

    eig = np.sort(np.real(np.linalg.eigvalsh(A0)))
    check("eig(A0) = {0, 1/3, 2/3} (cusp class)", np.allclose(eig, [0, 1 / 3, 2 / 3], atol=2e-3))
    check("diag(A0) = (1/2, 1/4, 1/4) (splitting condition)",
          np.allclose(np.real(np.diag(A0)), [0.5, 0.25, 0.25], atol=1e-6))

    Ak = [Ui[k] @ A0 @ Uic[k] for k in range(4)]
    S = sum(Ak)
    check("Sigma A_k = diag(2,1,1) => exponents at infinity (2,1,1) => O(-2)+O(-1)^2",
          np.allclose(S, np.diag([2, 1, 1]), atol=1e-6))

    M = Minf(A0)
    res = np.linalg.norm(M - I3)
    check(f"path-ordered monodromy at infinity trivial: ||M_inf - I|| = {res:.1e} < 1e-2 "
          "(=> M1 M2 M3 M4 = I, valid global flat bundle)", res < 1e-2)

    # irreducibility (case A): residues share no common eigenvector
    ev = np.linalg.eigh(Ak[0])[1]
    red = any(all(np.linalg.norm(Ak[j] @ ev[:, i] - (ev[:, i].conj() @ Ak[j] @ ev[:, i]) * ev[:, i]) < 1e-3
                  for j in range(4)) for i in range(3))
    check("IRREDUCIBLE (case A): non-trivial SU(3)_F mixing, not diagonal case B", not red)

    check("RESIDUAL: explicit flat bundle realised (existence + case A); unique nabla_F* needs "
          "det R=8 selector; c_u/c_d still needs the H2 (C6<->P1mu4) dictionary", True)
    return summary("v33 explicit U_wall flat bundle")


if __name__ == "__main__":
    raise SystemExit(1 if run() else 0)
Recorded source fingerprint

1babfc23211903730b3abebdead117d34a97ddd23d6c6d22119e7cff7b6f1a25

Paper sections that cite this module

What this source index covers

The public registry includes 116 document views, 6,693 section records and 985 physical LaTeX sources. Active publication documents, research material and archives retain their recorded roles. This is source navigation, not a new semantic review or proof run.