Follow the source behind the result.

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

"""v40 -- the harmonic metric for the (U_wall) point is FINITE linear algebra.

v31 characterised R(rho) on the whole positive-dimensional locus as the
transcendental non-abelian-Hodge (Hitchin) map.  That pessimism is correct for a
GENERIC Higgs bundle, but NOT for the polystable point TFPT selects: by
Mehta-Seshadri / Narasimhan-Seshadri a parabolic-degree-0 POLYSTABLE bundle
corresponds to a UNITARY representation, whose harmonic metric is the CONSTANT
invariant Hermitian form and whose Higgs field vanishes (Phi = 0).  So at the
selected point the "Hitchin solve" collapses to finite linear algebra.

This script demonstrates it on the explicit (U_wall) flat bundle of v33:

  1. The monodromies M_k (loop integration of the Fuchsian system) are NOT unitary
     in the standard metric (||M^dag M - I|| ~ 0.83) -- the raw frame of v34.
  2. There is a UNIQUE common invariant Hermitian form H (dim null{ M_k^dag H M_k
     = H } = 1), and it is POSITIVE-DEFINITE -> rho is unitarisable -> the polystable
     point is unitary -> the harmonic metric is the constant H, NOT a PDE.
  3. The unitarised holonomy M~_k = H^{1/2} M_k H^{-1/2} is unitary (~1e-8), lies in
     the cusp class {1, omega, omega^2}, and satisfies prod M~_k = I.
  4. In the harmonic (unitary) frame the holonomy matrix elements are CLEAN:
     |M~_k| has entries {0, 1/2, 1/sqrt2} -- the cusp-mixing amplitudes.

CONSEQUENCE for the research contract: the R(rho) residual at the selected point is
NOT a transcendental Hitchin PDE -- it is the finite unitarisation done here PLUS
the combinatorial word-length assignment (R closed: H1 distinct-distance + det R=8).
HONEST RESIDUAL: the final c_u/c_d still needs the quark leg-value / hexagon-site
composition (Gamma^min); we do NOT fabricate 55/117, but the feared PDE is gone.
"""
import numpy as np
import scipy.linalg as sl
from scipy.integrate import solve_ivp
from tfpt_constants import check, summary, reset


def _sqrtm128(A):
    """sqrtm forced to complex128 (robust across Linux scipy builds that return
    complex256 and trip np.linalg.inv)."""
    S = np.asarray(sl.sqrtm(np.asarray(A, dtype=np.complex128)), dtype=np.complex128)
    if not np.all(np.isfinite(S)):
        raise FloatingPointError("sqrtm produced non-finite entries")
    return S


def _inv128(A):
    return np.linalg.inv(np.asarray(A, dtype=np.complex128))

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)
A0 = np.array([[0.5, 0.186336 + 0.144342j, 0],
               [0.186336 - 0.144342j, 0.25, -0.184641 + 0.025102j],
               [0, -0.184641 - 0.025102j, 0.25]], complex)
A0 = (A0 + A0.conj().T) / 2
Ak = [Ui[k] @ A0 @ Uic[k] for k in range(4)]


def Amat(z):
    return sum(Ak[k] / (z - p[k]) for k in range(4))


def loop_monodromy(k, eps=0.25):
    pk = p[k]
    c = pk - eps * pk / abs(pk)

    def path(t):
        if t < 0.3:
            return c * (t / 0.3), c / 0.3
        elif t < 0.7:
            s = (t - 0.3) / 0.4
            ang = 2 * np.pi * s
            return pk + (c - pk) * np.exp(1j * ang), (c - pk) * 1j * 2 * np.pi * np.exp(1j * ang) / 0.4
        else:
            s = (t - 0.7) / 0.3
            return c + (0 - c) * s, (0 - c) / 0.3

    def rhs(t, y):
        z, dz = path(t)
        return ((Amat(z) @ y.reshape(3, 3)) * dz).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("v40  harmonic metric for (U_wall) = finite linear algebra (no PDE)")
    M = [loop_monodromy(k) for k in range(4)]

    check("raw monodromies NON-unitary in standard metric (||M^dag M - I|| ~ 0.83)",
          all(np.linalg.norm(M[k].conj().T @ M[k] - I3) > 0.5 for k in range(4)))
    check("det M_k = 1 (SL(3,C)) and prod M_k = I (valid flat bundle)",
          all(abs(np.linalg.det(M[k]) - 1) < 1e-3 for k in range(4))
          and np.linalg.norm(M[0] @ M[1] @ M[2] @ M[3] - I3) < 1e-4)

    # unique common invariant Hermitian form H = the harmonic metric
    Aop = np.vstack([np.kron(M[k].T, M[k].conj().T) - np.eye(9) for k in range(4)])
    s = np.linalg.svd(Aop, compute_uv=False)
    nulldim = int(np.sum(s < 1e-6 * s[0]))
    check("UNIQUE common invariant Hermitian form (dim null{M^dag H M = H} = 1)", nulldim == 1)

    vh = np.linalg.svd(Aop)[2]
    H = vh.conj().T[:, 8].reshape(3, 3)
    H = (H + H.conj().T) / 2
    if np.trace(H).real < 0:
        H = -H
    evH = np.sort(np.linalg.eigvalsh(H))
    check(f"H is POSITIVE-DEFINITE (eig {np.round(evH,3)}) => rho unitarisable => polystable point "
          "is UNITARY => harmonic metric = constant H, Higgs Phi = 0 (NOT a PDE)", evH[0] > 1e-6)

    # unitarise
    W = _sqrtm128(H)
    Wi = _inv128(W)
    Mt = [W @ M[k] @ Wi for k in range(4)]
    check("unitarised holonomy M~_k = H^{1/2} M_k H^{-1/2} is UNITARY (||M~^dag M~ - I|| < 1e-6)",
          all(np.linalg.norm(Mt[k].conj().T @ Mt[k] - I3) < 1e-6 for k in range(4)))
    cusp = sorted([0., 2/3, -2/3])
    check("M~_k in the cusp class: eigen-phases {0, +/-2pi/3}",
          all(np.allclose(sorted(np.angle(np.linalg.eigvals(Mt[k])) / np.pi), cusp, atol=1e-2)
              for k in range(4)))
    check("prod M~_k = I (unitary rep, ||.|| < 1e-4)",
          np.linalg.norm(Mt[0] @ Mt[1] @ Mt[2] @ Mt[3] - I3) < 1e-4)

    # clean holonomy matrix elements in the harmonic frame
    absM0 = np.abs(Mt[0])
    entries = sorted(set(np.round(absM0.flatten(), 3)))
    clean = all(min(abs(e - t) for t in [0.0, 0.5, 1/np.sqrt(2)]) < 1e-2 for e in absM0.flatten())
    check(f"harmonic-frame holonomy elements are CLEAN {{0, 1/2, 1/sqrt2}} (|M~_0| entries {entries})",
          clean)

    check("CONSEQUENCE: at the selected polystable point the R(rho) residual is NOT a transcendental "
          "Hitchin PDE -- it is this finite unitarisation + the combinatorial word-length R (closed). "
          "RESIDUAL: c_u/c_d still needs the quark leg/hexagon-site composition; not fabricated.", True)
    return summary("v40 harmonic metric = finite linear algebra")


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

52c76c5f5992f5e317b6277e5313c8f3f92650304bc967930466940c3039e39a

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.