nDot.io physics
07-ilang-weight / 03-weight-of-mark
03weight of markverified

Determines the cost of a held body with its own mark and its weight in the field of a distant source, for two couplings and three field graphs, and whether the weight follows the cost for every level.

Code of version 1

check_weight_of_mark.py16 KB · raw
"""Check for 07-ilang-weight/03-weight-of-mark, derivation v1.

One field object (N_f = 1, n = 1/|X|) on a ring, a square torus and a cubic torus.
Held body: point term g_b P_x, g_b = sigma*gamma*M_k; held source: point term g_s P_xs.
Exact diagonalization is compared with the closed forms of the derivation:
cost (3.2), (3.4), (3.5); weight (3.6)-(3.8); ring form (7.6); weight against cost
(3.12)-(3.17); mirror (3.18). Outputs: output/check_log.txt, output/check_results.json.
"""
import json
from pathlib import Path

import numpy as np
import scipy.sparse as sp
from scipy.optimize import brentq
from scipy.sparse.linalg import eigsh

OUT = Path(__file__).parent / "output"
OUT.mkdir(parents=True, exist_ok=True)
np.seterr(divide="ignore", invalid="ignore")   # F = 0 only at the reference place x', which is masked
T_F = 1.0
M_LEV = (1.0, 1.5)                      # half-gaps of two clock levels
GHAT = (0.005, 0.01, 0.02, 0.04)        # targets for max |g_b| * ell_0 / t_f
SIZES = {1: (32, 64), 2: (12, 20), 3: (8, 12)}
LOG, RES = [], {}


def log(s=""):
    print(s)
    LOG.append(s)


def graph(D, L):
    n = L**D
    co = np.array(np.unravel_index(np.arange(n), (L,) * D)).T
    rows, cols = [], []
    for mu in range(D):
        sh = co.copy()
        sh[:, mu] = (sh[:, mu] + 1) % L
        nb = np.ravel_multi_index(tuple(sh.T), (L,) * D)
        rows += [np.arange(n), nb]
        cols += [nb, np.arange(n)]
    A = sp.csr_matrix((np.ones(2 * D * n), (np.concatenate(rows), np.concatenate(cols))), shape=(n, n))
    k = 2 * np.pi * np.arange(L) / L
    grids = np.meshgrid(*([k] * D), indexing="ij")
    om = sum(4 * np.sin(g / 2) ** 2 for g in grids)
    return co, (-T_F * A).tocsr(), om


def ellpow(om, p=1, shift=0.0):
    """(1/|X|) sum_{k!=0} e^{iky}/(omega-shift)^p as an array over displacements y."""
    inv = np.zeros_like(om)
    m = om > 1e-12
    inv[m] = 1.0 / (om[m] - shift) ** p
    return np.fft.ifftn(inv).real.ravel()


def lowest(H, vec=True):
    if H.shape[0] <= 2000:
        if vec:
            w, v = np.linalg.eigh(H.toarray())
            return w[0], v[:, 0]
        return np.linalg.eigvalsh(H.toarray())[0], None
    w, v = eigsh(H, k=1, which="SA", tol=0)
    return w[0], v[:, 0]


def point(n, i, g):
    return sp.csr_matrix(([g], ([i], [i])), shape=(n, n))


def exact_rho(om, g, X):
    """Lowest eigenpair of C_1 + g P_0 from (3.2)-(3.3): returns delta, eta, rho(y)."""
    wmin = om[om > 1e-12].min()
    f = lambda d: 1.0 / (X * d) - ellpow(om, 1, d / T_F)[0] / T_F - 1.0 / g
    d = brentq(f, 1e-14, T_F * wmin * (1 - 1e-12)) if g > 0 else brentq(f, -1e3, -1e-14)
    eta = X * d / T_F
    ld = ellpow(om, 1, d / T_F)
    l2d0 = np.sum(ld**2)
    return d, eta, (1 - eta * ld) ** 2 / (X + eta**2 * l2d0)


def slope(x, y):
    return float(np.polyfit(np.log(np.abs(x)), np.log(np.abs(y)), 1)[0])


def run_graph(D, L):
    tag = f"D{D}_L{L}"
    co, C1, om = graph(D, L)
    X = L**D
    n = 1.0 / X
    l1, l2, l3 = ellpow(om, 1), ellpow(om, 2), ellpow(om, 3)
    l0, l20 = l1[0], l2[0]
    ap = int(np.ravel_multi_index((L // 2,) * D, (L,) * D))   # antipode: r' of the pairs
    F = -2 * (l1 - l1[ap]) / T_F                              # F / g_s, (3.13)
    F2 = (l1**2 - l1[ap] ** 2 - 2 * (l2 - l2[ap]) / X) / (n * T_F**2)
    F3 = 2 * (l1 * l2 - l1[ap] * l2[ap] + l20 * (l1 - l1[ap]) - (l3 - l3[ap]) / X) / (n**2 * T_F**3 * X)
    dN = n * F                                                # [dN(x)-dN(x')]/g_s
    dTh = n**2 * F2 - n * l0 * F / T_F                        # [Theta_D(x)-Theta_D(x')]/g_s, (3.8)
    mask = np.arange(X) != ap
    eps00 = -2 * D * T_F
    r = {"L": L, "D": D, "X": X, "ell_0": l0, "ell2_0": l20}
    log(f"\n=== {tag}: |X|={X}, ell_0={l0:.6f}, ell2_0/|X|={l20 / X:.6e}")

    # (A) weight from eigenvalues with both point terms vs Hellmann-Feynman density
    gb = 0.02 / l0
    _, v = lowest(C1 + point(X, 0, gb))
    rho_hf = v**2
    errs = {}
    rng = np.random.default_rng(D * 100 + L)
    pl = [(int(x), int(xs)) for x, xs in rng.integers(0, X, size=(4, 2))]
    for gs in (1e-3, 5e-4):
        e = 0.0
        for x, xs in pl:
            ep = lowest(C1 + point(X, x, gb) + point(X, xs, gs), vec=False)[0]
            em = lowest(C1 + point(X, x, gb) + point(X, xs, -gs), vec=False)[0]
            ydisp = np.ravel_multi_index(tuple((co[xs] - co[x]) % L), (L,) * D)
            e = max(e, abs((ep - em) / (2 * gs) - rho_hf[ydisp]) / rho_hf.max())
        errs[gs] = e
    r["A_eig_vs_HF_relerr"] = errs
    log(f"(A) eigenvalue weight vs HF density, max rel err: g_s=1e-3: {errs[1e-3]:.2e}, g_s=5e-4: {errs[5e-4]:.2e}")

    # mirror (3.18): highest eigenvalue = - lowest of mirrored signs (L even, bipartite)
    gs = 1e-3
    Hp = (C1 + point(X, 0, gb) + point(X, 3 % X, gs))
    Hm = (C1 + point(X, 0, -gb) + point(X, 3 % X, -gs))
    top = -lowest(-Hp, vec=False)[0]
    r["mirror_err"] = abs(top + lowest(Hm, vec=False)[0])
    log(f"(M) highest state = -(lowest, mirrored signs): |diff| = {r['mirror_err']:.2e}")

    # scans over levels, branches and gamma
    gams = np.array(GHAT) * T_F / (l0 * max(M_LEV))
    rows = []
    for gam in gams:
        for M in M_LEV:
            for s in (1, -1):
                g = s * gam * M
                e, v = lowest(C1 + point(X, 0, g))
                d, eta, rho_ex = exact_rho(om, g, X)
                Cm = e - eps00                                   # mark cost, (3.4)
                w = v**2 - v[ap] ** 2                            # [W(x)-W(x')]/g_s, (3.6)
                rows.append(dict(gam=gam, M=M, s=s, g=g, Cm=Cm, w=w,
                                 cost_ex_err=abs(Cm - d), rho_ex_err=np.abs(v**2 - rho_ex).max() * X))
    r["B_rho_exact_maxrelerr"] = max(q["rho_ex_err"] for q in rows)
    r["C_cost_exact_maxerr"] = max(q["cost_ex_err"] for q in rows)
    log(f"(B) density vs exact form (3.3): max |X|*|diff| = {r['B_rho_exact_maxrelerr']:.2e}")
    log(f"(C) cost vs exact root (3.2): max |diff| = {r['C_cost_exact_maxerr']:.2e}")

    # (D),(E),(H) deviations from the expansions, scaling in gamma (level M=1.5, sigma=+1)
    sel = [q for q in rows if q["M"] == M_LEV[1] and q["s"] == 1]
    gg = np.array([q["g"] for q in sel])
    devC = [abs(q["Cm"] - (n * q["g"] - n * q["g"] ** 2 * l0 / T_F)) for q in sel]
    devW = [np.abs(q["w"] - (q["g"] * dN + q["g"] ** 2 * dTh))[mask].max() for q in sel]
    devU = [np.abs(q["w"] - (F * q["Cm"] + F2 * q["Cm"] ** 2 + F3 * q["Cm"] ** 3))[mask].max() for q in sel]
    r["D_cost_2nd_slope"] = slope(gg, devC)
    r["E_weight_2nd_slope"] = slope(gg, devW)
    r["H_W_of_C_3rd_slope"] = slope(gg, devU)
    log(f"(D) cost minus 2nd-order (3.5): slope in gamma = {r['D_cost_2nd_slope']:.3f}; devs {np.array(devC)}")
    log(f"(E) weight minus 2nd-order (3.8): slope = {r['E_weight_2nd_slope']:.3f}; devs {np.array(devW)}")
    log(f"(H) weight minus F C + F2 C^2 + F3 C^3 (3.12)-(3.13): slope = {r['H_W_of_C_3rd_slope']:.3f}")

    # (F) ring: Theta_D of (3.8) against Theta of (7.6)
    if D == 1:
        rr = np.arange(L)
        u = rr * (L - rr)
        th76 = n / (3 * L**2 * T_F**2) * (u**2 - (L**2 - 2) / 2 * u)
        diff = dTh - (th76 - th76[ap])
        r["F_ring_Theta_vs_7.6_maxdiff"] = float(np.abs(diff).max())
        log(f"(F) ring: Theta_D - Theta(7.6): max |diff| = {r['F_ring_Theta_vs_7.6_maxdiff']:.2e}")

    # (I) single-level break: B = dW/(Cost*F_c) - 1, both couplings
    ibr = []
    for q in rows:
        g, gam, M, s = q["g"], q["gam"], q["M"], q["s"]
        Bfg = q["w"][mask] / (q["Cm"] * F[mask]) - 1
        pfg = (q["Cm"] * F2 / F)[mask]
        cost_a = s * M + q["Cm"]
        Fa = gam * n * F / (1 + gam * n)
        Ba = q["w"][mask] / (cost_a * Fa[mask]) - 1
        pa = (g * dTh / dN)[mask]
        ibr.append((gam, np.abs(Bfg - pfg).max(), np.abs(Ba - pa).max(), np.abs(Bfg).max(), np.abs(Ba).max()))
    ib = np.array([t for t, q in zip(ibr, rows) if q["M"] == M_LEV[1] and q["s"] == 1])
    r["I_fullgrip_break_resid_slope"] = slope(ib[:, 0], ib[:, 1])
    r["I_a_break_resid_slope"] = slope(ib[:, 0], ib[:, 2])
    r["I_fullgrip_break_size_slope"] = slope(ib[:, 0], ib[:, 3])
    r["I_a_break_size_slope"] = slope(ib[:, 0], ib[:, 4])
    log(f"(I) single level: |B| slope fg {r['I_fullgrip_break_size_slope']:.3f}, (a) {r['I_a_break_size_slope']:.3f}; "
        f"|B - B_1st| slope fg {r['I_fullgrip_break_resid_slope']:.3f}, (a) {r['I_a_break_resid_slope']:.3f}")

    # (J) level and branch differences (3.17)
    jl = []
    for gam in gams:
        Q = {(q["M"], q["s"]): q for q in rows if q["gam"] == gam}
        Fa = gam * n * F / (1 + gam * n)
        for M in M_LEV:          # branch pair (M,+),(M,-)
            p, m_ = Q[(M, 1)], Q[(M, -1)]
            g = gam * M
            Rfg = (p["w"] - m_["w"])[mask] / ((p["Cm"] - m_["Cm"]) * F[mask]) - 1
            pfg = (g**2 * (-2 * n * l0 / T_F * F2 + n**2 * F3) / F)[mask]
            Ra = (p["w"] - m_["w"])[mask] / ((2 * M + p["Cm"] - m_["Cm"]) * Fa[mask]) - 1
            pa = pfg + g**2 * (l0**2 - l20 / X) / T_F**2
            jl.append(("branch", gam, M, np.abs(Rfg - pfg).max() / np.abs(pfg).max(),
                       np.abs(Ra - pa).max() / np.abs(pa).max()))
        p, m_ = Q[(M_LEV[1], 1)], Q[(M_LEV[0], 1)]  # level pair, sigma=+1
        Rfg = (p["w"] - m_["w"])[mask] / ((p["Cm"] - m_["Cm"]) * F[mask]) - 1
        pfg = ((p["Cm"] + m_["Cm"]) * F2 / F)[mask]
        Ra = (p["w"] - m_["w"])[mask] / ((M_LEV[1] - M_LEV[0] + p["Cm"] - m_["Cm"]) * Fa[mask]) - 1
        pa = (gam * (M_LEV[0] + M_LEV[1]) * dTh / dN)[mask]
        jl.append(("level", gam, None, np.abs(Rfg - pfg).max() / np.abs(pfg).max(),
                   np.abs(Ra - pa).max() / np.abs(pa).max()))
    r["J_rel_err_of_predicted_break_smallest_gamma"] = {
        f"{k}_{M}": [fg, a] for k, gm, M, fg, a in jl if gm == gams[0]}
    log("(J) relative error of the predicted difference breaks at the smallest gamma "
        "[full grip, (a)]: " + ", ".join(f"{k}{'' if M is None else M}: {fg:.1e}, {a:.1e}"
                                        for k, gm, M, fg, a in jl if gm == gams[0]))
    log("    ... and at the largest gamma: " + ", ".join(f"{k}{'' if M is None else M}: {fg:.1e}, {a:.1e}"
                                                         for k, gm, M, fg, a in jl if gm == gams[-1]))
    return r, (co, C1, om, l1, l2, ap)


def full_operator_check(D, L, gam=0.05, m=1.0, Nclock=2):
    """(R) reduction (3.1): spectrum of the full contract (body labels sigma, a; field) at tau=0
    equals the union over (k, sigma) of theta*sigma*M_k + spec(C_1 + sigma*gam*M_k P_x + g_s P_xs)."""
    co, C1, om = graph(D, L)
    X = L**D
    lam0 = 2 * np.pi                                       # gives M_k = m + k/2
    fk = np.exp(2j * np.pi * np.outer(np.arange(Nclock), np.arange(Nclock)) / Nclock) / np.sqrt(Nclock)
    Dc = (2 * np.pi / (Nclock * lam0)) * fk @ np.diag(np.arange(Nclock)) @ fk.conj().T
    body = np.kron(np.diag([1.0, -1.0]), m * np.eye(Nclock) + Dc)   # sigma_3 (x) (m 1 + D)
    x, xs, gs = 0, X // 3, 1e-3
    Px = np.zeros((X, X)); Px[x, x] = 1
    Ps = np.zeros((X, X)); Ps[xs, xs] = 1
    out = {}
    Mk = m + np.arange(Nclock) / 2
    for theta, name in ((1, "a"), (0, "full_grip")):
        H = theta * np.kron(body, np.eye(X)) + gam * np.kron(body, Px) \
            + np.kron(np.eye(2 * Nclock), C1.toarray() + gs * Ps)
        full = np.sort(np.linalg.eigvalsh(H))
        red = np.sort(np.concatenate([theta * s * M + np.linalg.eigvalsh(C1.toarray() + s * gam * M * Px + gs * Ps)
                                      for s in (1, -1) for M in Mk]))
        out[name] = float(np.abs(full - red).max())
    log(f"(R) {D}D L={L}: full contract vs reduction (3.1), max |diff|: (a) {out['a']:.2e}, full grip {out['full_grip']:.2e}")
    return out


def asymptotics():
    """(S) large-distance forms of ell = Lambda^+ used in (3.9)-(3.11) and (3.16)."""
    res = {}
    for L in (64, 128, 256, 512):
        k =2 * np.pi * np.arange(L) / L
        g1, g2 = np.meshgrid(k, k, indexing="ij")
        l1 = ellpow(4 * np.sin(g1 / 2) ** 2 + 4 * np.sin(g2 / 2) ** 2).reshape(L, L)
        res[f"2D_L{L}_ell0_minus_lnL_over_2pi"] = l1[0, 0] - np.log(L) / (2 * np.pi)
        if L == 512:   # ell(j,0) - [ell_0 - (ln j + gamma_E + 1.5 ln 2)/(2 pi) + j^2/(4 L^2)]
            res["2D_L512_ell(j)_minus_(3.10)_form"] = {
                j: l1[j, 0] - (l1[0, 0] - (np.log(j) + np.euler_gamma + 1.5 * np.log(2)) / (2 * np.pi)
                               + j**2 / (4 * L**2)) for j in (2, 4, 8, 16, 32, 64)}
    xi3 = 2.837297 / (4 * np.pi)
    for L in (32, 64, 96):
        k = 2 * np.pi * np.arange(L) / L
        gr = np.meshgrid(k, k, k, indexing="ij")
        l1 = ellpow(sum(4 * np.sin(g / 2) ** 2 for g in gr)).reshape(L, L, L)
        res[f"3D_L{L}_ell0"] = l1[0, 0, 0]
        res[f"3D_L{L}_ell0_plus_xi3_over_L"] = l1[0, 0, 0] + xi3 / L
        if L == 96:    # ell(j,0,0) - [1/(4 pi j) - xi3/L + j^2/(6 L^3)], relative to 1/(4 pi j)
            res["3D_L96_rel_dev_of_ell(j)_from_(3.11)_form"] = {
                j: (l1[j, 0, 0] - (1 / (4 * np.pi * j) - xi3 / L + j**2 / (6 * L**3))) * 4 * np.pi * j
                for j in (2, 4, 8, 16, 24)}
    log("\n(S) asymptotics of ell: " + json.dumps({k: (v if not isinstance(v, dict) else
                                                     {kk: round(float(vv), 5) for kk, vv in v.items()})
                                                 for k, v in res.items()}, default=float))
    return res


def distance_scan(D, L, data):
    """(K) break against distance along a lattice axis, r' = antipode; g_b*ell_0/t_f = +-0.05."""
    co, C1, om, l1, l2, ap = data
    X = L**D
    n = 1.0 / X
    l0 = l1[0]
    out = []
    log(f"\n(K) distance scan {D}D, L={L}: j, B_fg(num), B_fg(1st order), B_a(num), B_a(1st order), "
        f"-(eta/2)(ell(r)+ell(r')) [fg], -(g_b/t_f)ell_0 [(a) far limit]")
    for s in (1, -1):
        g = s * 0.05 / l0
        gam, M = abs(g), 1.0
        e, v = lowest(C1 + point(X, 0, g))
        Cm = e + 2 * D * T_F
        eta = Cm / (n * T_F)
        for j in range(1, L // 2 + 1):
            y = np.ravel_multi_index((j,) + (0,) * (D - 1), (L,) * D)
            if y == ap:
                continue
            Fy = -2 * (l1[y] - l1[ap]) / T_F
            F2y = (l1[y] ** 2 - l1[ap] ** 2 - 2 * (l2[y] - l2[ap]) / X) / (n * T_F**2)
            w = v[y] ** 2 - v[ap] ** 2
            Bfg = w / (Cm * Fy) - 1
            Fa = gam * n * Fy / (1 + gam * n)
            Ba = w / ((s * M + Cm) * Fa) - 1
            dThy = n**2 * F2y - n * l0 * Fy / T_F
            row = [s, j, Bfg, Cm * F2y / Fy, Ba, g * dThy / (n * Fy), -eta / 2 * (l1[y] + l1[ap]), -g * l0 / T_F]
            out.append(row)
            log("  s=%+d j=%2d  % .4e  % .4e  % .4e  % .4e  % .4e  % .4e" % tuple(row))
    return out


if __name__ == "__main__":
    data = {}
    for D, Ls in SIZES.items():
        for L in Ls:
            RES[f"D{D}_L{L}"], data[(D, L)] = run_graph(D, L)
    scans = {}
    for D, Ls in SIZES.items():
        scans[f"D{D}_L{Ls[-1]}"] = distance_scan(D, Ls[-1], data[(D, Ls[-1])])
    # larger cubic torus for the distance dependence (sparse Lanczos)
    co, C1, om = graph(3, 20)
    l1, l2 = ellpow(om, 1), ellpow(om, 2)
    ap = int(np.ravel_multi_index((10, 10, 10), (20,) * 3))
    RES["D3_L20"] = {"ell_0": l1[0]}
    log(f"\n=== D3_L20: ell_0={l1[0]:.6f}")
    scans["D3_L20"] = distance_scan(3, 20, (co, C1, om, l1, l2, ap))
    RES["distance_scans_columns"] = ["sigma", "j", "B_fg_num", "B_fg_1st", "B_a_num", "B_a_1st",
                                     "-(eta/2)(ell(r)+ell(r'))", "-(g_b/t_f)ell_0"]
    RES["distance_scans"] = scans
    log("")
    RES["R_reduction_full_operator"] = {f"D{D}_L{L}": full_operator_check(D, L)
                                        for D, L in ((1, 32), (2, 8), (3, 6))}
    RES["S_asymptotics"] = asymptotics()

    def clean(o):
        if isinstance(o, dict):
            return {str(k): clean(v) for k, v in o.items()}
        if isinstance(o, (list, tuple)):
            return [clean(v) for v in o]
        if isinstance(o, (np.floating, np.integer)):
            return o.item()
        return o
    (OUT / "check_results.json").write_text(json.dumps(clean(RES), indent=1), encoding="utf-8")
    (OUT / "check_log.txt").write_text("\n".join(LOG) + "\n", encoding="utf-8")
output/check_log.txt18 KB · raw

=== D1_L32: |X|=32, ell_0=2.664062, ell2_0/|X|=1.436096e+00
(A) eigenvalue weight vs HF density, max rel err: g_s=1e-3: 5.59e-06, g_s=5e-4: 1.40e-06
(M) highest state = -(lowest, mirrored signs): |diff| = 0.00e+00
(B) density vs exact form (3.3): max |X|*|diff| = 5.08e-11
(C) cost vs exact root (3.2): max |diff| = 2.95e-13
(D) cost minus 2nd-order (3.5): slope in gamma = 2.991; devs [1.16627836e-09 9.30385539e-09 7.40101657e-08 5.85392552e-07]
(E) weight minus 2nd-order (3.8): slope = 3.001; devs [6.91700019e-09 5.53543136e-08 4.43066650e-07 3.54637848e-06]
(H) weight minus F C + F2 C^2 + F3 C^3 (3.12)-(3.13): slope = 3.936
(F) ring: Theta_D - Theta(7.6): max |diff| = 2.11e-15
(I) single level: |B| slope fg 0.990, (a) 0.990; |B - B_1st| slope fg 1.966, (a) 1.981
(J) relative error of the predicted difference breaks at the smallest gamma [full grip, (a)]: branch1.0: 7.9e-06, 4.8e-05, branch1.5: 1.5e-05, 6.1e-05, level: 2.5e-03, 5.0e-03
    ... and at the largest gamma: branch1.0: 4.1e-04, 9.5e-04, branch1.5: 9.3e-04, 1.7e-03, level: 1.9e-02, 3.8e-02

=== D1_L64: |X|=64, ell_0=5.332031, ell2_0/|X|=5.702774e+00
(A) eigenvalue weight vs HF density, max rel err: g_s=1e-3: 2.25e-05, g_s=5e-4: 5.62e-06
(M) highest state = -(lowest, mirrored signs): |diff| = 0.00e+00
(B) density vs exact form (3.3): max |X|*|diff| = 3.42e-10
(C) cost vs exact root (3.2): max |diff| = 4.99e-13
(D) cost minus 2nd-order (3.5): slope in gamma = 2.991; devs [2.91991896e-10 2.32928109e-09 1.85281026e-08 1.46537716e-07]
(E) weight minus 2nd-order (3.8): slope = 3.000; devs [3.50192267e-09 2.80212992e-08 2.24235220e-07 1.79401083e-06]
(H) weight minus F C + F2 C^2 + F3 C^3 (3.12)-(3.13): slope = 3.937
(F) ring: Theta_D - Theta(7.6): max |diff| = 5.00e-16
(I) single level: |B| slope fg 0.990, (a) 0.990; |B - B_1st| slope fg 1.966, (a) 1.981
(J) relative error of the predicted difference breaks at the smallest gamma [full grip, (a)]: branch1.0: 6.4e-06, 2.0e-05, branch1.5: 1.5e-05, 3.2e-05, level: 2.5e-03, 5.0e-03
    ... and at the largest gamma: branch1.0: 4.1e-04, 7.1e-04, branch1.5: 9.2e-04, 1.5e-03, level: 1.9e-02, 3.8e-02

=== D2_L12: |X|=144, ell_0=0.444054, ell2_0/|X|=4.236694e-03
(A) eigenvalue weight vs HF density, max rel err: g_s=1e-3: 1.93e-07, g_s=5e-4: 4.85e-08
(M) highest state = -(lowest, mirrored signs): |diff| = 0.00e+00
(B) density vs exact form (3.3): max |X|*|diff| = 6.37e-11
(C) cost vs exact root (3.2): max |diff| = 4.97e-13
(D) cost minus 2nd-order (3.5): slope in gamma = 2.984; devs [1.90373019e-09 1.51576337e-08 1.20121630e-07 9.43235634e-07]
(E) weight minus 2nd-order (3.8): slope = 2.981; devs [3.41824170e-09 2.71900738e-08 2.15066569e-07 1.68249522e-06]
(H) weight minus F C + F2 C^2 + F3 C^3 (3.12)-(3.13): slope = 3.934
(I) single level: |B| slope fg 0.985, (a) 0.980; |B - B_1st| slope fg 1.965, (a) 1.981
(J) relative error of the predicted difference breaks at the smallest gamma [full grip, (a)]: branch1.0: 1.4e-05, 3.2e-05, branch1.5: 1.8e-05, 2.9e-05, level: 5.9e-04, 7.9e-03
    ... and at the largest gamma: branch1.0: 5.1e-04, 6.5e-04, branch1.5: 1.1e-03, 1.8e-03, level: 4.5e-03, 6.0e-02

=== D2_L20: |X|=400, ell_0=0.525479, ell2_0/|X|=4.025398e-03
(A) eigenvalue weight vs HF density, max rel err: g_s=1e-3: 2.73e-07, g_s=5e-4: 7.21e-08
(M) highest state = -(lowest, mirrored signs): |diff| = 0.00e+00
(B) density vs exact form (3.3): max |X|*|diff| = 1.29e-10
(C) cost vs exact root (3.2): max |diff| = 3.06e-13
(D) cost minus 2nd-order (3.5): slope in gamma = 2.984; devs [5.83192816e-10 4.64312159e-09 3.67910188e-08 2.88821749e-07]
(E) weight minus 2nd-order (3.8): slope = 2.981; devs [1.24641547e-09 9.91305574e-09 7.83874699e-08 6.12901570e-07]
(H) weight minus F C + F2 C^2 + F3 C^3 (3.12)-(3.13): slope = 3.891
(I) single level: |B| slope fg 0.985, (a) 0.979; |B - B_1st| slope fg 1.965, (a) 1.981
(J) relative error of the predicted difference breaks at the smallest gamma [full grip, (a)]: branch1.0: 1.5e-05, 8.5e-06, branch1.5: 3.3e-05, 2.4e-05, level: 3.9e-04, 8.0e-03
    ... and at the largest gamma: branch1.0: 5.7e-04, 8.7e-04, branch1.5: 1.3e-03, 2.1e-03, level: 3.0e-03, 6.1e-02

=== D3_L8: |X|=512, ell_0=0.224606, ell2_0/|X|=1.858469e-04
(A) eigenvalue weight vs HF density, max rel err: g_s=1e-3: 5.37e-08, g_s=5e-4: 1.68e-08
(M) highest state = -(lowest, mirrored signs): |diff| = 0.00e+00
(B) density vs exact form (3.3): max |X|*|diff| = 7.05e-11
(C) cost vs exact root (3.2): max |diff| = 3.09e-13
(D) cost minus 2nd-order (3.5): slope in gamma = 2.984; devs [1.07761861e-09 8.57860554e-09 6.79607966e-08 5.33302810e-07]
(E) weight minus 2nd-order (3.8): slope = 2.980; devs [9.78138329e-10 7.77746691e-09 6.14706875e-08 4.80186621e-07]
(H) weight minus F C + F2 C^2 + F3 C^3 (3.12)-(3.13): slope = 3.805
(I) single level: |B| slope fg 0.984, (a) 0.979; |B - B_1st| slope fg 1.964, (a) 1.980
(J) relative error of the predicted difference breaks at the smallest gamma [full grip, (a)]: branch1.0: 2.1e-05, 2.3e-05, branch1.5: 2.5e-05, 2.1e-05, level: 9.5e-05, 8.3e-03
    ... and at the largest gamma: branch1.0: 6.8e-04, 9.1e-04, branch1.5: 1.5e-03, 2.2e-03, level: 7.3e-04, 6.3e-02

=== D3_L12: |X|=1728, ell_0=0.233942, ell2_0/|X|=7.996514e-05
(A) eigenvalue weight vs HF density, max rel err: g_s=1e-3: 6.40e-08, g_s=5e-4: 6.45e-08
(M) highest state = -(lowest, mirrored signs): |diff| = 0.00e+00
(B) density vs exact form (3.3): max |X|*|diff| = 1.89e-10
(C) cost vs exact root (3.2): max |diff| = 2.42e-13
(D) cost minus 2nd-order (3.5): slope in gamma = 2.984; devs [3.07231407e-10 2.44570346e-09 1.93743548e-08 1.52022297e-07]
(E) weight minus 2nd-order (3.8): slope = 2.980; devs [2.89676337e-10 2.30318923e-09 1.82019437e-08 1.42160716e-07]
(H) weight minus F C + F2 C^2 + F3 C^3 (3.12)-(3.13): slope = 4.048
(I) single level: |B| slope fg 0.984, (a) 0.979; |B - B_1st| slope fg 1.963, (a) 1.980
(J) relative error of the predicted difference breaks at the smallest gamma [full grip, (a)]: branch1.0: 2.7e-05, 1.2e-05, branch1.5: 2.9e-05, 3.3e-05, level: 3.7e-05, 8.3e-03
    ... and at the largest gamma: branch1.0: 7.0e-04, 1.0e-03, branch1.5: 1.6e-03, 2.3e-03, level: 2.9e-04, 6.4e-02

(K) distance scan 1D, L=64: j, B_fg(num), B_fg(1st order), B_a(num), B_a(1st order), -(eta/2)(ell(r)+ell(r')) [fg], -(g_b/t_f)ell_0 [(a) far limit]
  s=+1 j= 1   2.5276e-03   2.9639e-03  -4.5645e-02  -4.6886e-02  -9.6938e-03  -5.0000e-02
  s=+1 j= 2   5.3999e-03   5.8000e-03  -4.2911e-02  -4.3907e-02  -7.5667e-03  -5.0000e-02
  s=+1 j= 3   8.1842e-03   8.5431e-03  -4.0260e-02  -4.1026e-02  -5.5094e-03  -5.0000e-02
  s=+1 j= 4   1.0880e-02   1.1193e-02  -3.7694e-02  -3.8242e-02  -3.5218e-03  -5.0000e-02
  s=+1 j= 5   1.3487e-02   1.3750e-02  -3.5212e-02  -3.5556e-02  -1.6040e-03  -5.0000e-02
  s=+1 j= 6   1.6004e-02   1.6214e-02  -3.2816e-02  -3.2967e-02   2.4409e-04  -5.0000e-02
  s=+1 j= 7   1.8430e-02   1.8586e-02  -3.0506e-02  -3.0476e-02   2.0224e-03  -5.0000e-02
  s=+1 j= 8   2.0766e-02   2.0864e-02  -2.8283e-02  -2.8083e-02   3.7311e-03  -5.0000e-02
  s=+1 j= 9   2.3011e-02   2.3049e-02  -2.6146e-02  -2.5788e-02   5.3699e-03  -5.0000e-02
  s=+1 j=10   2.5164e-02   2.5141e-02  -2.4097e-02  -2.3590e-02   6.9391e-03  -5.0000e-02
  s=+1 j=11   2.7224e-02   2.7140e-02  -2.2135e-02  -2.1490e-02   8.4385e-03  -5.0000e-02
  s=+1 j=12   2.9192e-02   2.9046e-02  -2.0262e-02  -1.9487e-02   9.8681e-03  -5.0000e-02
  s=+1 j=13   3.1066e-02   3.0860e-02  -1.8478e-02  -1.7582e-02   1.1228e-02  -5.0000e-02
  s=+1 j=14   3.2847e-02   3.2580e-02  -1.6782e-02  -1.5775e-02   1.2518e-02  -5.0000e-02
  s=+1 j=15   3.4534e-02   3.4207e-02  -1.5176e-02  -1.4066e-02   1.3739e-02  -5.0000e-02
  s=+1 j=16   3.6127e-02   3.5741e-02  -1.3660e-02  -1.2454e-02   1.4889e-02  -5.0000e-02
  s=+1 j=17   3.7624e-02   3.7183e-02  -1.2235e-02  -1.0940e-02   1.5970e-02  -5.0000e-02
  s=+1 j=18   3.9027e-02   3.8531e-02  -1.0899e-02  -9.5238e-03   1.6982e-02  -5.0000e-02
  s=+1 j=19   4.0335e-02   3.9786e-02  -9.6546e-03  -8.2051e-03   1.7923e-02  -5.0000e-02
  s=+1 j=20   4.1546e-02   4.0949e-02  -8.5012e-03  -6.9841e-03   1.8795e-02  -5.0000e-02
  s=+1 j=21   4.2662e-02   4.2018e-02  -7.4390e-03  -5.8608e-03   1.9597e-02  -5.0000e-02
  s=+1 j=22   4.3682e-02   4.2994e-02  -6.4685e-03  -4.8352e-03   2.0329e-02  -5.0000e-02
  s=+1 j=23   4.4605e-02   4.3878e-02  -5.5897e-03  -3.9072e-03   2.0992e-02  -5.0000e-02
  s=+1 j=24   4.5431e-02   4.4668e-02  -4.8029e-03  -3.0769e-03   2.1584e-02  -5.0000e-02
  s=+1 j=25   4.6161e-02   4.5365e-02  -4.1082e-03  -2.3443e-03   2.2107e-02  -5.0000e-02
  s=+1 j=26   4.6794e-02   4.5970e-02  -3.5059e-03  -1.7094e-03   2.2561e-02  -5.0000e-02
  s=+1 j=27   4.7329e-02   4.6481e-02  -2.9960e-03  -1.1722e-03   2.2944e-02  -5.0000e-02
  s=+1 j=28   4.7768e-02   4.6900e-02  -2.5786e-03  -7.3260e-04   2.3258e-02  -5.0000e-02
  s=+1 j=29   4.8109e-02   4.7225e-02  -2.2539e-03  -3.9072e-04   2.3502e-02  -5.0000e-02
  s=+1 j=30   4.8353e-02   4.7458e-02  -2.0219e-03  -1.4652e-04   2.3677e-02  -5.0000e-02
  s=+1 j=31   4.8499e-02   4.7597e-02  -1.8827e-03   1.2410e-15   2.3781e-02  -5.0000e-02
  s=-1 j= 1  -3.7745e-03  -3.2756e-03   4.8077e-02   4.6886e-02   1.0713e-02   5.0000e-02
  s=-1 j= 2  -6.8627e-03  -6.4100e-03   4.4828e-02   4.3907e-02   8.3625e-03   5.0000e-02
  s=-1 j= 3  -9.8424e-03  -9.4415e-03   4.1693e-02   4.1026e-02   6.0888e-03   5.0000e-02
  s=-1 j= 4  -1.2714e-02  -1.2370e-02   3.8672e-02   3.8242e-02   3.8922e-03   5.0000e-02
  s=-1 j= 5  -1.5479e-02  -1.5196e-02   3.5763e-02   3.5556e-02   1.7727e-03   5.0000e-02
  s=-1 j= 6  -1.8137e-02  -1.7920e-02   3.2967e-02   3.2967e-02  -2.6976e-04   5.0000e-02
  s=-1 j= 7  -2.0690e-02  -2.0540e-02   3.0281e-02   3.0476e-02  -2.2351e-03   5.0000e-02
  s=-1 j= 8  -2.3137e-02  -2.3058e-02   2.7706e-02   2.8083e-02  -4.1234e-03   5.0000e-02
  s=-1 j= 9  -2.5480e-02  -2.5473e-02   2.5241e-02   2.5788e-02  -5.9347e-03   5.0000e-02
  s=-1 j=10  -2.7719e-02  -2.7785e-02   2.2886e-02   2.3590e-02  -7.6688e-03   5.0000e-02
  s=-1 j=11  -2.9855e-02  -2.9995e-02   2.0639e-02   2.1490e-02  -9.3259e-03   5.0000e-02
  s=-1 j=12  -3.1888e-02  -3.2101e-02   1.8500e-02   1.9487e-02  -1.0906e-02   5.0000e-02
  s=-1 j=13  -3.3818e-02  -3.4105e-02   1.6469e-02   1.7582e-02  -1.2409e-02   5.0000e-02
  s=-1 j=14  -3.5647e-02  -3.6006e-02   1.4546e-02   1.5775e-02  -1.3835e-02   5.0000e-02
  s=-1 j=15  -3.7374e-02  -3.7805e-02   1.2728e-02   1.4066e-02  -1.5184e-02   5.0000e-02
  s=-1 j=16  -3.9001e-02  -3.9500e-02   1.1017e-02   1.2454e-02  -1.6455e-02   5.0000e-02
  s=-1 j=17  -4.0526e-02  -4.1093e-02   9.4123e-03   1.0940e-02  -1.7650e-02   5.0000e-02
  s=-1 j=18  -4.1952e-02  -4.2583e-02   7.9126e-03   9.5238e-03  -1.8767e-02   5.0000e-02
  s=-1 j=19  -4.3278e-02  -4.3971e-02   6.5178e-03   8.2051e-03  -1.9808e-02   5.0000e-02
  s=-1 j=20  -4.4504e-02  -4.5255e-02   5.2278e-03   6.9841e-03  -2.0771e-02   5.0000e-02
  s=-1 j=21  -4.5631e-02  -4.6437e-02   4.0422e-03   5.8608e-03  -2.1658e-02   5.0000e-02
  s=-1 j=22  -4.6659e-02  -4.7516e-02   2.9606e-03   4.8352e-03  -2.2467e-02   5.0000e-02
  s=-1 j=23  -4.7588e-02  -4.8492e-02   1.9829e-03   3.9072e-03  -2.3199e-02   5.0000e-02
  s=-1 j=24  -4.8419e-02  -4.9366e-02   1.1087e-03   3.0769e-03  -2.3854e-02   5.0000e-02
  s=-1 j=25  -4.9152e-02  -5.0136e-02   3.3783e-04   2.3443e-03  -2.4432e-02   5.0000e-02
  s=-1 j=26  -4.9787e-02  -5.0804e-02  -3.2984e-04   1.7094e-03  -2.4933e-02   5.0000e-02
  s=-1 j=27  -5.0323e-02  -5.1370e-02  -8.9452e-04   1.1722e-03  -2.5357e-02   5.0000e-02
  s=-1 j=28  -5.0762e-02  -5.1832e-02  -1.3563e-03   7.3260e-04  -2.5704e-02   5.0000e-02
  s=-1 j=29  -5.1104e-02  -5.2192e-02  -1.7154e-03   3.9072e-04  -2.5974e-02   5.0000e-02
  s=-1 j=30  -5.1347e-02  -5.2449e-02  -1.9718e-03   1.4652e-04  -2.6167e-02   5.0000e-02
  s=-1 j=31  -5.1493e-02  -5.2603e-02  -2.1257e-03  -1.2410e-15  -2.6282e-02   5.0000e-02

(K) distance scan 2D, L=20: j, B_fg(num), B_fg(1st order), B_a(num), B_a(1st order), -(eta/2)(ell(r)+ell(r')) [fg], -(g_b/t_f)ell_0 [(a) far limit]
  s=+1 j= 1  -8.4978e-03  -8.4456e-03  -5.5733e-02  -5.8868e-02  -9.9924e-03  -5.0000e-02
  s=+1 j= 2  -2.8846e-03  -2.8411e-03  -5.0387e-02  -5.2983e-02  -4.9406e-03  -5.0000e-02
  s=+1 j= 3   4.8955e-04   5.2406e-04  -4.7174e-02  -4.9450e-02  -2.0532e-03  -5.0000e-02
  s=+1 j= 4   2.8317e-03   2.8573e-03  -4.4943e-02  -4.7000e-02  -1.4103e-04  -5.0000e-02
  s=+1 j= 5   4.5826e-03   4.5997e-03  -4.3275e-02  -4.5170e-02   1.2180e-03  -5.0000e-02
  s=+1 j= 6   5.9393e-03   5.9485e-03  -4.1983e-02  -4.3754e-02   2.2061e-03  -5.0000e-02
  s=+1 j= 7   7.0020e-03   7.0042e-03  -4.0971e-02  -4.2645e-02   2.9146e-03  -5.0000e-02
  s=+1 j= 8   7.8122e-03   7.8084e-03  -4.0200e-02  -4.1801e-02   3.3937e-03  -5.0000e-02
  s=+1 j= 9   8.3500e-03   8.3420e-03  -3.9688e-02  -4.1241e-02   3.6715e-03  -5.0000e-02
  s=+1 j=10   8.5443e-03   8.5346e-03  -3.9503e-02  -4.1038e-02   3.7626e-03  -5.0000e-02
  s=-1 j= 1   9.2705e-03   9.3346e-03   6.2334e-02   5.8868e-02   1.1044e-02   5.0000e-02
  s=-1 j= 2   3.0870e-03   3.1401e-03   5.5825e-02   5.2983e-02   5.4607e-03   5.0000e-02
  s=-1 j= 3  -6.2103e-04  -5.7922e-04   5.1922e-02   4.9450e-02   2.2693e-03   5.0000e-02
  s=-1 j= 4  -3.1888e-03  -3.1580e-03   4.9220e-02   4.7000e-02   1.5587e-04   5.0000e-02
  s=-1 j= 5  -5.1041e-03  -5.0838e-03   4.7204e-02   4.5170e-02  -1.3462e-03   5.0000e-02
  s=-1 j= 6  -6.5852e-03  -6.5746e-03   4.5645e-02   4.3754e-02  -2.4384e-03   5.0000e-02
  s=-1 j= 7  -7.7434e-03  -7.7414e-03   4.4425e-02   4.2645e-02  -3.2214e-03   5.0000e-02
  s=-1 j= 8  -8.6250e-03  -8.6303e-03   4.3497e-02   4.1801e-02  -3.7509e-03   5.0000e-02
  s=-1 j= 9  -9.2096e-03  -9.2200e-03   4.2882e-02   4.1241e-02  -4.0580e-03   5.0000e-02
  s=-1 j=10  -9.4206e-03  -9.4330e-03   4.2660e-02   4.1038e-02  -4.1587e-03   5.0000e-02

(K) distance scan 3D, L=12: j, B_fg(num), B_fg(1st order), B_a(num), B_a(1st order), -(eta/2)(ell(r)+ell(r')) [fg], -(g_b/t_f)ell_0 [(a) far limit]
  s=+1 j= 1  -6.0983e-03  -6.0929e-03  -5.3424e-02  -5.6398e-02  -6.3037e-03  -5.0000e-02
  s=+1 j= 2  -1.5676e-03  -1.5633e-03  -4.9110e-02  -5.1641e-02  -1.9403e-03  -5.0000e-02
  s=+1 j= 3   1.0608e-04   1.0947e-04  -4.7516e-02  -4.9885e-02  -4.3436e-04  -5.0000e-02
  s=+1 j= 4   9.1017e-04   9.1284e-04  -4.6750e-02  -4.9042e-02   2.1651e-04  -5.0000e-02
  s=+1 j= 5   1.3332e-03   1.3354e-03  -4.6347e-02  -4.8598e-02   5.1213e-04  -5.0000e-02
  s=+1 j= 6   1.4763e-03   1.4782e-03  -4.6211e-02  -4.8448e-02   5.9926e-04  -5.0000e-02
  s=-1 j= 1   6.7276e-03   6.7342e-03   5.9702e-02   5.6398e-02   6.9672e-03   5.0000e-02
  s=-1 j= 2   1.7226e-03   1.7279e-03   5.4434e-02   5.1641e-02   2.1446e-03   5.0000e-02
  s=-1 j= 3  -1.2512e-04  -1.2099e-04   5.2489e-02   4.9885e-02   4.8009e-04   5.0000e-02
  s=-1 j= 4  -1.0122e-03  -1.0089e-03   5.1555e-02   4.9042e-02  -2.3930e-04   5.0000e-02
  s=-1 j= 5  -1.4786e-03  -1.4760e-03   5.1064e-02   4.8598e-02  -5.6604e-04   5.0000e-02
  s=-1 j= 6  -1.6362e-03  -1.6338e-03   5.0898e-02   4.8448e-02  -6.6234e-04   5.0000e-02

=== D3_L20: ell_0=0.241447

(K) distance scan 3D, L=20: j, B_fg(num), B_fg(1st order), B_a(num), B_a(1st order), -(eta/2)(ell(r)+ell(r')) [fg], -(g_b/t_f)ell_0 [(a) far limit]
  s=+1 j= 1  -6.9872e-03  -6.9853e-03  -5.4273e-02  -5.7335e-02  -7.0592e-03  -5.0000e-02
  s=+1 j= 2  -2.6715e-03  -2.6699e-03  -5.0163e-02  -5.2803e-02  -2.8079e-03  -5.0000e-02
  s=+1 j= 3  -1.1026e-03  -1.1012e-03  -4.8669e-02  -5.1156e-02  -1.3055e-03  -5.0000e-02
  s=+1 j= 4  -3.4119e-04  -3.3992e-04  -4.7944e-02  -5.0357e-02  -6.0565e-04  -5.0000e-02
  s=+1 j= 5   1.0698e-04   1.0808e-04  -4.7517e-02  -4.9887e-02  -2.1404e-04  -5.0000e-02
  s=+1 j= 6   4.0166e-04   4.0261e-04  -4.7236e-02  -4.9577e-02   2.7876e-05  -5.0000e-02
  s=+1 j= 7   6.0608e-04   6.0691e-04  -4.7042e-02  -4.9363e-02   1.8294e-04  -5.0000e-02
  s=+1 j= 8   7.4717e-04   7.4789e-04  -4.6907e-02  -4.9215e-02   2.8014e-04  -5.0000e-02
  s=+1 j= 9   8.3330e-04   8.3396e-04  -4.6825e-02  -4.9124e-02   3.3392e-04  -5.0000e-02
  s=+1 j=10   8.6280e-04   8.6343e-04  -4.6797e-02  -4.9093e-02   3.5117e-04  -5.0000e-02
  s=-1 j= 1   7.7183e-03   7.7206e-03   6.0753e-02   5.7335e-02   7.8022e-03   5.0000e-02
  s=-1 j= 2   2.9489e-03   2.9509e-03   5.5733e-02   5.2803e-02   3.1034e-03   5.0000e-02
  s=-1 j= 3   1.2153e-03   1.2171e-03   5.3908e-02   5.1156e-02   1.4429e-03   5.0000e-02
  s=-1 j= 4   3.7416e-04   3.7571e-04   5.3023e-02   5.0357e-02   6.6940e-04   5.0000e-02
  s=-1 j= 5  -1.2080e-04  -1.1946e-04   5.2502e-02   4.9887e-02   2.3657e-04   5.0000e-02
  s=-1 j= 6  -4.4615e-04  -4.4499e-04   5.2159e-02   4.9577e-02  -3.0810e-05   5.0000e-02
  s=-1 j= 7  -6.7180e-04  -6.7079e-04   5.1922e-02   4.9363e-02  -2.0219e-04   5.0000e-02
  s=-1 j= 8  -8.2750e-04  -8.2662e-04   5.1758e-02   4.9215e-02  -3.0963e-04   5.0000e-02
  s=-1 j= 9  -9.2255e-04  -9.2174e-04   5.1658e-02   4.9124e-02  -3.6907e-04   5.0000e-02
  s=-1 j=10  -9.5509e-04  -9.5432e-04   5.1623e-02   4.9093e-02  -3.8813e-04   5.0000e-02

(R) 1D L=32: full contract vs reduction (3.1), max |diff|: (a) 8.88e-15, full grip 3.55e-15
(R) 2D L=8: full contract vs reduction (3.1), max |diff|: (a) 1.24e-14, full grip 1.29e-14
(R) 3D L=6: full contract vs reduction (3.1), max |diff|: (a) 3.55e-14, full grip 1.64e-14

(S) asymptotics of ell: {"2D_L64_ell0_minus_lnL_over_2pi": 0.04875850142855509, "2D_L128_ell0_minus_lnL_over_2pi": 0.048763848880206684, "2D_L256_ell0_minus_lnL_over_2pi": 0.04876518701309496, "2D_L512_ell0_minus_lnL_over_2pi": 0.048765521625585895, "2D_L512_ell(j)_minus_(3.10)_form": {"2": 0.00428, "4": 0.00099, "8": 0.00021, "16": 5e-05, "32": 1e-05, "64": 3e-05}, "3D_L32_ell0": 0.2456765511865084, "3D_L32_ell0_plus_xi3_over_L": 0.2527323299766616, "3D_L64_ell0": 0.24920328389380492, "3D_L64_ell0_plus_xi3_over_L": 0.2527311732888815, "3D_L96_ell0": 0.25037913173616777, "3D_L96_ell0_plus_xi3_over_L": 0.2527310579995522, "3D_L96_rel_dev_of_ell(j)_from_(3.11)_form": {"2": 0.07793, "4": 0.01983, "8": 0.00413, "16": 0.0014, "24": 0.00353}}
output/check_results.json27 KB · raw
{
 "D1_L32": {
  "L": 32,
  "D": 1,
  "X": 32,
  "ell_0": 2.6640624999999956,
  "ell2_0": 45.9550781249998,
  "A_eig_vs_HF_relerr": {
   "0.001": 5.593997352368179e-06,
   "0.0005": 1.398521412502525e-06
  },
  "mirror_err": 0.0,
  "B_rho_exact_maxrelerr": 5.0804693785266863e-11,
  "C_cost_exact_maxerr": 2.951468821947578e-13,
  "D_cost_2nd_slope": 2.9905868110918195,
  "E_weight_2nd_slope": 3.000670899665059,
  "H_W_of_C_3rd_slope": 3.935874776880874,
  "F_ring_Theta_vs_7.6_maxdiff": 2.1094237467877974e-15,
  "I_fullgrip_break_resid_slope": 1.965790789376314,
  "I_a_break_resid_slope": 1.9808046223523343,
  "I_fullgrip_break_size_slope": 0.989885196782842,
  "I_a_break_size_slope": 0.9899200073495279,
  "J_rel_err_of_predicted_break_smallest_gamma": {
   "branch_1.0": [
    7.92286583059946e-06,
    4.8271872889614226e-05
   ],
   "branch_1.5": [
    1.4851634722096521e-05,
    6.125707876856205e-05
   ],
   "level_None": [
    0.002513606245519754,
    0.005031743612268052
   ]
  }
 },
 "D1_L64": {
  "L": 64,
  "D": 1,
  "X": 64,
  "ell_0": 5.33203125,
  "ell2_0": 364.97753906250057,
  "A_eig_vs_HF_relerr": {
   "0.001": 2.2498772426377858e-05,
   "0.0005": 5.624735272831155e-06
  },
  "mirror_err": 0.0,
  "B_rho_exact_maxrelerr": 3.415210336754626e-10,
  "C_cost_exact_maxerr": 4.990802591347838e-13,
  "D_cost_2nd_slope": 2.990514245403253,
  "E_weight_2nd_slope": 3.0002894123515746,
  "H_W_of_C_3rd_slope": 3.936973999291328,
  "F_ring_Theta_vs_7.6_maxdiff": 4.996003610813204e-16,
  "I_fullgrip_break_resid_slope": 1.9658742118016928,
  "I_a_break_resid_slope": 1.980826179823064,
  "I_fullgrip_break_size_slope": 0.9899014122970299,
  "I_a_break_size_slope": 0.9899101154249147,
  "J_rel_err_of_predicted_break_smallest_gamma": {
   "branch_1.0": [
    6.388140071694525e-06,
    1.9817788013936252e-05
   ],
   "branch_1.5": [
    1.4770131603217338e-05,
    3.1682750953346943e-05
   ],
   "level_None": [
    0.002519274992540249,
    0.005033163674340294
   ]
  }
 },
 "D2_L12": {
  "L": 12,
  "D": 2,
  "X": 144,
  "ell_0": 0.4440538049913048,
  "ell2_0": 0.6100839189036009,
  "A_eig_vs_HF_relerr": {
   "0.001": 1.9270573478764696e-07,
   "0.0005": 4.854699098408378e-08
  },
  "mirror_err": 0.0,
  "B_rho_exact_maxrelerr": 6.372177785429756e-11,
  "C_cost_exact_maxerr": 4.965628060647764e-13,
  "D_cost_2nd_slope": 2.984431567524043,
  "E_weight_2nd_slope": 2.9813028136666433,
  "H_W_of_C_3rd_slope": 3.9340053708821494,
  "I_fullgrip_break_resid_slope": 1.9649230428561733,
  "I_a_break_resid_slope": 1.9811867754348187,
  "I_fullgrip_break_size_slope": 0.9852116695179468,
  "I_a_break_size_slope": 0.9796468751185036,
  "J_rel_err_of_predicted_break_smallest_gamma": {
   "branch_1.0": [
    1.397384855593063e-05,
    3.2142580562300794e-05
   ],
   "branch_1.5": [
    1.797735611495196e-05,
    2.8618716241269706e-05
   ],
   "level_None": [
    0.0005853485210257128,
    0.007887520225396079
   ]
  }
 },
 "D2_L20": {
  "L": 20,
  "D": 2,
  "X": 400,
  "ell_0": 0.5254788958215179,
  "ell2_0": 1.6101590287737275,
  "A_eig_vs_HF_relerr": {
   "0.001": 2.733845906268316e-07,
   "0.0005": 7.209020735170592e-08
  },
  "mirror_err": 0.0,
  "B_rho_exact_maxrelerr": 1.2892638345807228e-10,
  "C_cost_exact_maxerr": 3.062704813882027e-13,
  "D_cost_2nd_slope": 2.9842158813322874,
  "E_weight_2nd_slope": 2.9808400960032637,
  "H_W_of_C_3rd_slope": 3.891436600811146,
  "I_fullgrip_break_resid_slope": 1.9646828367608753,
  "I_a_break_resid_slope": 1.9808060527474542,
  "I_fullgrip_break_size_slope": 0.9847190513743058,
  "I_a_break_size_slope": 0.9793698644972091,
  "J_rel_err_of_predicted_break_smallest_gamma": {
   "branch_1.0": [
    1.5310738485614688e-05,
    8.462004184107557e-06
   ],
   "branch_1.5": [
    3.302149205394912e-05,
    2.4309753573379333e-05
   ],
   "level_None": [
    0.0003864741837919538,
    0.008005991600455096
   ]
  }
 },
 "D3_L8": {
  "L": 8,
  "D": 3,
  "X": 512,
  "ell_0": 0.22460562539027706,
  "ell2_0": 0.09515363538880361,
  "A_eig_vs_HF_relerr": {
   "0.001": 5.37170971326829e-08,
   "0.0005": 1.6808438213811314e-08
  },
  "mirror_err": 0.0,
  "B_rho_exact_maxrelerr": 7.046230265927989e-11,
  "C_cost_exact_maxerr": 3.093339657773281e-13,
  "D_cost_2nd_slope": 2.9838781054221313,
  "E_weight_2nd_slope": 2.9800549323665346,
  "H_W_of_C_3rd_slope": 3.8048839104634387,
  "I_fullgrip_break_resid_slope": 1.9639749982714991,
  "I_a_break_resid_slope": 1.9799963309718023,
  "I_fullgrip_break_size_slope": 0.983995966467738,
  "I_a_break_size_slope": 0.978643186673062,
  "J_rel_err_of_predicted_break_smallest_gamma": {
   "branch_1.0": [
    2.1301191634222452e-05,
    2.275932529584844e-05
   ],
   "branch_1.5": [
    2.536110563448039e-05,
    2.1460439884637648e-05
   ],
   "level_None": [
    9.45152616173216e-05,
    0.00829329281608392
   ]
  }
 },
 "D3_L12": {
  "L": 12,
  "D": 3,
  "X": 1728,
  "ell_0": 0.23394212427744138,
  "ell2_0": 0.13817975408779593,
  "A_eig_vs_HF_relerr": {
   "0.001": 6.396102311823995e-08,
   "0.0005": 6.448775636097462e-08
  },
  "mirror_err": 0.0,
  "B_rho_exact_maxrelerr": 1.889639694718781e-10,
  "C_cost_exact_maxerr": 2.420098355656458e-13,
  "D_cost_2nd_slope": 2.9838051431478396,
  "E_weight_2nd_slope": 2.9798983553968057,
  "H_W_of_C_3rd_slope": 4.04826638642241,
  "I_fullgrip_break_resid_slope": 1.9627871359583995,
  "I_a_break_resid_slope": 1.9798818980933928,
  "I_fullgrip_break_size_slope": 0.9838540126643667,
  "I_a_break_size_slope": 0.9785498923284232,
  "J_rel_err_of_predicted_break_smallest_gamma": {
   "branch_1.0": [
    2.7096138434447814e-05,
    1.2334259632546264e-05
   ],
   "branch_1.5": [
    2.8760747194350076e-05,
    3.277840671627963e-05
   ],
   "level_None": [
    3.721847763348036e-05,
    0.008334470483390791
   ]
  }
 },
 "D3_L20": {
  "ell_0": 0.24144725262740885
 },
 "distance_scans_columns": [
  "sigma",
  "j",
  "B_fg_num",
  "B_fg_1st",
  "B_a_num",
  "B_a_1st",
  "-(eta/2)(ell(r)+ell(r'))",
  "-(g_b/t_f)ell_0"
 ],
 "distance_scans": {
  "D1_L64": [
   [
    1,
    1,
    0.002527589986057066,
    0.002963922581530643,
    -0.04564496756462,
    -0.04688644688644686,
    -0.009693770325476563,
    -0.05
   ],
   [
    1,
    2,
    0.005399874693970963,
    0.005799989679152086,
    -0.042910699308099653,
    -0.04390720390720387,
    -0.007566720002260478,
    -0.05
   ],
   [
    1,
    3,
    0.008184234217196718,
    0.008543070970294142,
    -0.04026013133406814,
    -0.04102564102564099,
    -0.005509409033903938,
    -0.05
   ],
   [
    1,
    4,
    0.010880054092164304,
    0.011193166454956804,
    -0.037693848580443,
    -0.03824175824175821,
    -0.003521837420406943,
    -0.05
   ],
   [
    1,
    5,
    0.013486738944518128,
    0.013750276133140067,
    -0.035212417813184516,
    -0.035555555555555535,
    -0.0016040051617694962,
    -0.05
   ],
   [
    1,
    6,
    0.016003712666155412,
    0.01621440000484395,
    -0.03281638745776494,
    -0.03296703296703294,
    0.000244087742008415,
    -0.05
   ],
   [
    1,
    7,
    0.018430418586177,
    0.018585538070068432,
    -0.030506287436432467,
    -0.030476190476190455,
    0.0020224412909267786,
    -0.05
   ],
   [
    1,
    8,
    0.02076631963629061,
    0.02086369032881352,
    -0.028282629010754734,
    -0.02808302808302807,
    0.003731055484985598,
    -0.05
   ],
   [
    1,
    9,
    0.02301089851055993,
    0.02304885678107923,
    -0.026145904629546468,
    -0.025787545787545767,
    0.005369930324184877,
    -0.05
   ],
   [
    1,
    10,
    0.025163657819314178,
    0.025141037426865537,
    -0.024096587782355905,
    -0.02358974358974357,
    0.006939065808524612,
    -0.05
   ],
   [
    1,
    11,
    0.027224120236902127,
    0.027140232266172456,
    -0.022135132858809636,
    -0.021489621489621472,
    0.008438461938004804,
    -0.05
   ],
   [
    1,
    12,
    0.02919182864401071,
    0.02904644129899998,
    -0.020261975013133315,
    -0.019487179487179474,
    0.00986811871262545,
    -0.05
   ],
   [
    1,
    13,
    0.03106634626377547,
    0.03085966452534812,
    -0.018477530034580414,
    -0.017582417582417572,
    0.011228036132386553,
    -0.05
   ],
   [
    1,
    14,
    0.03284725679161005,
    0.032579901945216865,
    -0.016782194223842528,
    -0.015775335775335766,
    0.012518214197288114,
    -0.05
   ],
   [
    1,
    15,
    0.0345341645195687,
    0.034207153558606214,
    -0.015176344274661746,
    -0.014065934065934054,
    0.013738652907330132,
    -0.05
   ],
   [
    1,
    16,
    0.03612669445416827,
    0.03574141936551618,
    -0.013660337161670144,
    -0.012454212454212441,
    0.014889352262512605,
    -0.05
   ],
   [
    1,
    17,
    0.03762449242793675,
    0.037182699365946745,
    -0.012234510034202062,
    -0.010940170940170931,
    0.015970312262835535,
    -0.05
   ],
   [
    1,
    18,
    0.039027225205491956,
    0.03853099355989792,
    -0.010899180115311324,
    -0.009523809523809521,
    0.01698153290829892,
    -0.05
   ],
   [
    1,
    19,
    0.04033458058252348,
    0.03978630194736971,
    -0.009654644607546614,
    -0.008205128205128205,
    0.017923014198902763,
    -0.05
   ],
   [
    1,
    20,
    0.04154626747914891,
    0.0409486245283621,
    -0.008501180604080338,
    -0.006984126984126983,
    0.018794756134647064,
    -0.05
   ],
   [
    1,
    21,
    0.04266201602701192,
    0.04201796130287511,
    -0.0074390450057968405,
    -0.005860805860805859,
    0.019596758715531817,
    -0.05
   ],
   [
    1,
    22,
    0.043681577649099745,
    0.04299431227090871,
    -0.006468474445308958,
    -0.004835164835164843,
    0.020329021941557026,
    -0.05
   ],
   [
    1,
    23,
    0.04460472513477609,
    0.04387767743246294,
    -0.005589685215530604,
    -0.003907203907203908,
    0.0209915458127227,
    -0.05
   ],
   [
    1,
    24,
    0.045431252707392344,
    0.04466805678753777,
    -0.004802873205315694,
    -0.003076923076923077,
    0.02158433032902882,
    -0.05
   ],
   [
    1,
    25,
    0.04616097608558878,
    0.045365450336133206,
    -0.004108213841100827,
    -0.0023443223443223443,
    0.022107375490475402,
    -0.05
   ],
   [
    1,
    26,
    0.04679373253758201,
    0.045969858078249264,
    -0.003505862035227847,
    -0.0017094017094017014,
    0.02256068129706244,
    -0.05
   ],
   [
    1,
    27,
    0.047329380932528586,
    0.04648128001388591,
    -0.0029959521370468423,
    -0.0011721611721611561,
    0.02294424774878993,
    -0.05
   ],
   [
    1,
    28,
    0.047767801779639685,
    0.04689971614304326,
    -0.0025785978956829325,
    -0.0007326007326006384,
    0.023258074845657877,
    -0.05
   ],
   [
    1,
    29,
    0.04810889726466061,
    0.04722516646572104,
    -0.002253892425307713,
    -0.0003907203907203907,
    0.023502162587666287,
    -0.05
   ],
   [
    1,
    30,
    0.048352591283290725,
    0.04745763098191948,
    -0.002021908173325926,
    -0.00014652014652018607,
    0.02367651097481514,
    -0.05
   ],
   [
    1,
    31,
    0.048498829448693215,
    0.04759710969163979,
    -0.0018826969132269555,
    1.2409772029686266e-15,
    0.023781120007104466,
    -0.05
   ],
   [
    -1,
    1,
    -0.003774491752387199,
    -0.0032756321772032066,
    0.04807701221935323,
    0.04688644688644686,
    0.010713244061911523,
    0.05
   ],
   [
    -1,
    2,
    -0.006862699830418761,
    -0.006409962574213296,
    0.0448280692152454,
    0.04390720390720387,
    0.008362496264153954,
    0.05
   ],
   [
    -1,
    3,
    -0.009842350224917307,
    -0.00944152804017388,
    0.04169333410048437,
    0.04102564102564099,
    0.006088822164683518,
    0.05
   ],
   [
    -1,
    4,
    -0.012714186649599224,
    -0.012370328575084945,
    0.03867202445257356,
    0.03824175824175821,
    0.003892221763500219,
    0.05
   ],
   [
    -1,
    5,
    -0.015478925348487937,
    -0.015196364178946495,
    0.03576338674845014,
    0.035555555555555535,
    0.001772695060604057,
    0.05
   ],
   [
    -1,
    6,
    -0.018137255334111813,
    -0.017919634851758546,
    0.03296669611389014,
    0.03296703296703294,
    -0.00026975794400498,
    0.05
   ],
   [
    -1,
    7,
    -0.02068983861656104,
    -0.020540140593521078,
    0.03028125608252985,
    0.030476190476190455,
    -0.0022351372503268793,
    0.05
   ],
   [
    -1,
    8,
    -0.023137310423588908,
    -0.023057881404234093,
    0.02770639836430777,
    0.02808302808302807,
    -0.004123442858361645,
    0.05
   ],
   [
    -1,
    9,
    -0.025480279411467466,
    -0.02547285728389761,
    0.025241482623636013,
    0.025787545787545767,
    -0.005934674768109281,
    0.05
   ],
   [
    -1,
    10,
    -0.02771932786737763,
    -0.02778506823251161,
    0.022885896266474415,
    0.02358974358974357,
    -0.007668832979569784,
    0.05
   ],
   [
    -1,
    11,
    -0.029855011902488515,
    -0.029994514250076095,
    0.02063905423720236,
    0.021489621489621472,
    -0.009325917492743151,
    0.05
   ],
   [
    -1,
    12,
    -0.03188786163649404,
    -0.032101195336591073,
    0.018500398824478292,
    0.019487179487179474,
    -0.010905928307629386,
    0.05
   ],
   [
    -1,
    13,
    -0.03381838137357995,
    -0.03410511149205654,
    0.016469399476113367,
    0.017582417582417572,
    -0.012408865424228487,
    0.05
   ],
   [
    -1,
    14,
    -0.03564704976933386,
    -0.0360062627164725,
    0.014545552623473457,
    0.015775335775335766,
    -0.013834728842540459,
    0.05
   ],
   [
    -1,
    15,
    -0.03737431998973684,
    -0.037804649009838945,
    0.012728381514213627,
    0.014065934065934054,
    -0.015183518562565295,
    0.05
   ],
   [
    -1,
    16,
    -0.03900061986085901,
    -0.03950027037215588,
    0.011017436054789886,
    0.012454212454212441,
    -0.016455234584302997,
    0.05
   ],
   [
    -1,
    17,
    -0.04052635201118182,
    -0.0410931268034233,
    0.00941229266073007,
    0.010940170940170931,
    -0.017649876907753565,
    0.05
   ],
   [
    -1,
    18,
    -0.04195189400444377,
    -0.0425832183036412,
    0.007912554116874304,
    0.009523809523809521,
    -0.018767445532917003,
    0.05
   ],
   [
    -1,
    19,
    -0.04327759846462009,
    -0.04397054487280961,
    0.006517849445889068,
    0.008205128205128205,
    -0.019807940459793304,
    0.05
   ],
   [
    -1,
    20,
    -0.044503793192851915,
    -0.04525510651092849,
    0.005227833785253377,
    0.006984126984126983,
    -0.020771361688382478,
    0.05
   ],
   [
    -1,
    21,
    -0.04563078127598974,
    -0.04643690321799788,
    0.004042188273064351,
    0.005860805860805859,
    -0.021657709218684512,
    0.05
   ],
   [
    -1,
    22,
    -0.04665884118610386,
    -0.04751593499401774,
    0.0029606199433491742,
    0.004835164835164843,
    -0.022466983050699413,
    0.05
   ],
   [
    -1,
    23,
    -0.04758822687358766,
    -0.04849220183898809,
    0.001982861628115229,
    0.003907203907203908,
    -0.023199183184427187,
    0.05
   ],
   [
    -1,
    24,
    -0.0484191678509045,
    -0.04936570375290894,
    0.0011086718692434605,
    0.003076923076923077,
    -0.02385430961986782,
    0.05
   ],
   [
    -1,
    25,
    -0.049151869267496484,
    -0.050136440735780274,
    0.00033783483968252703,
    0.0023443223443223443,
    -0.024432362357021325,
    0.05
   ],
   [
    -1,
    26,
    -0.04978651197872919,
    -0.0508044127876021,
    -0.0003298397290865118,
    0.0017094017094017014,
    -0.0249333413958877,
    0.05
   ],
   [
    -1,
    27,
    -0.050323252605525615,
    -0.05136961990837442,
    -0.0008945166097847901,
    0.0011721611721611561,
    -0.025357246736466933,
    0.05
   ],
   [
    -1,
    28,
    -0.05076222358595539,
    -0.0518320620980973,
    -0.001356335133613773,
    0.0007326007326006384,
    -0.025704078378759033,
    0.05
   ],
   [
    -1,
    29,
    -0.05110353321861438,
    -0.052191739356770496,
    -0.0017154092358904194,
    0.0003907203907203907,
    -0.025973836322764003,
    0.05
   ],
   [
    -1,
    30,
    -0.051347265697300526,
    -0.052448651684394235,
    -0.0019718274925272228,
    0.00014652014652018607,
    -0.026166520568481833,
    0.05
   ],
   [
    -1,
    31,
    -0.05149348113594443,
    -0.052602799080969855,
    -0.002125653146262896,
    -1.2409772029686266e-15,
    -0.02628213111591254,
    0.05
   ]
  ],
  "D2_L20": [
   [
    1,
    1,
    -0.008497837658679175,
    -0.008445643078092693,
    -0.0557327934305073,
    -0.05886821888200905,
    -0.00999235394824904,
    -0.05
   ],
   [
    1,
    2,
    -0.00288461865396461,
    -0.00284105126543199,
    -0.05038698700591204,
    -0.05298320261037477,
    -0.004940637314145103,
    -0.05
   ],
   [
    1,
    3,
    0.0004895511014693099,
    0.0005240642284487167,
    -0.047173561992363666,
    -0.04944971433872577,
    -0.002053196360715093,
    -0.05
   ],
   [
    1,
    4,
    0.0028317254530250136,
    0.002857303752008343,
    -0.04494296833735889,
    -0.04699973171363074,
    -0.00014102885832305707,
    -0.05
   ],
   [
    1,
    5,
    0.004582612270540354,
    0.004599686014085723,
    -0.04327549340186254,
    -0.04517016975684937,
    0.0012180387111356602,
    -0.05
   ],
   [
    1,
    6,
    0.005939268656292462,
    0.005948477159443807,
    -0.041983467842761524,
    -0.04375389215755347,
    0.0022061426221026905,
    -0.05
   ],
   [
    1,
    7,
    0.007001995826571683,
    0.007004153162835841,
    -0.040971368772744765,
    -0.04264539564876194,
    0.0029146454799306165,
    -0.05
   ],
   [
    1,
    8,
    0.007812211983358797,
    0.0078084130777393325,
    -0.040199751144317175,
    -0.041800894774471985,
    0.003393713178430251,
    -0.05
   ],
   [
    1,
    9,
    0.00835003838000059,
    0.00834199523029321,
    -0.039687546684795794,
    -0.04124061496195475,
    0.0036715333678174283,
    -0.05
   ],
   [
    1,
    10,
    0.008544281042796253,
    0.008534648389576236,
    -0.03950255770185562,
    -0.04103832244627123,
    0.0037626398233278112,
    -0.05
   ],
   [
    -1,
    1,
    0.009270521855733849,
    0.009334591036147241,
    0.06233392123407078,
    0.05886821888200905,
    0.011044101287832363,
    0.05
   ],
   [
    -1,
    2,
    0.0030869901048313775,
    0.003140086720492259,
    0.05582528416426569,
    0.05298320261037477,
    0.0054606651452157534,
    0.05
   ],
   [
    -1,
    3,
    -0.0006210330431644095,
    -0.0005792247202503816,
    0.05192230801908493,
    0.04944971433872577,
    0.002269305980251004,
    0.05
   ],
   [
    -1,
    4,
    -0.0031888146042367005,
    -0.0031580498659991127,
    0.049219522793939596,
    0.04699973171363074,
    0.00015587288079403248,
    0.05
   ],
   [
    -1,
    5,
    -0.0051041305351925415,
    -0.005083826943569209,
    0.04720350722698008,
    0.04517016975684937,
    -0.0013462436346783166,
    0.05
   ],
   [
    -1,
    6,
    -0.006585244699461201,
    -0.006574585387737049,
    0.04564452201553282,
    0.04375389215755347,
    -0.002438350632903235,
    0.05
   ],
   [
    -1,
    7,
    -0.007743409544690438,
    -0.007741376793343602,
    0.04442546550409765,
    0.04264539564876194,
    -0.003221427109686911,
    0.05
   ],
   [
    -1,
    8,
    -0.008625047066790681,
    -0.00863028925660698,
    0.043497475014255826,
    0.041800894774471985,
    -0.0037509191806604195,
    0.05
   ],
   [
    -1,
    9,
    -0.009209613413230233,
    -0.009220033712088089,
    0.0428821745119754,
    0.04124061496195475,
    -0.004057981393156846,
    0.05
   ],
   [
    -1,
    10,
    -0.009420606173775492,
    -0.009432964620616954,
    0.04266008859763626,
    0.04103832244627123,
    -0.004158677278014725,
    0.05
   ]
  ],
  "D3_L12": [
   [
    1,
    1,
    -0.006098251644688801,
    -0.006092883946521033,
    -0.05342446836999526,
    -0.05639754934460445,
    -0.006303688951855443,
    -0.05
   ],
   [
    1,
    2,
    -0.0015676078933928306,
    -0.0015633183383474887,
    -0.04910955844594522,
    -0.05164148969497667,
    -0.001940314726022822,
    -0.05
   ],
   [
    1,
    3,
    0.0001060791331743971,
    0.0001094667512354463,
    -0.0475155666961804,
    -0.049885059530303026,
    -0.00043436383887427235,
    -0.05
   ],
   [
    1,
    4,
    0.0009101667353155207,
    0.0009128416561389211,
    -0.046749767007495446,
    -0.04904151308473655,
    0.00021650696483286805,
    -0.05
   ],
   [
    1,
    5,
    0.0013332140852768681,
    0.0013353859756821606,
    -0.04634686373173713,
    -0.048597840078933814,
    0.0005121337730161345,
    -0.05
   ],
   [
    1,
    6,
    0.0014762516991735364,
    0.0014782276268491285,
    -0.04621063708192841,
    -0.04844785584817766,
    0.0005992607708214372,
    -0.05
   ],
   [
    -1,
    1,
    0.0067276317515891915,
    0.006734235218695972,
    0.05970211160661254,
    0.05639754934460445,
    0.006967230053926826,
    0.05
   ],
   [
    -1,
    2,
    0.0017226283670985776,
    0.001727876897793885,
    0.05443374259809031,
    0.05164148969497667,
    0.0021445568105393107,
    0.05
   ],
   [
    -1,
    3,
    -0.00012512301513900592,
    -0.00012098947853206555,
    0.05248876167204375,
    0.049885059530303026,
    0.0004800859965739723,
    0.05
   ],
   [
    -1,
    4,
    -0.0010121867312177235,
    -0.0010089295124968718,
    0.0515550202474413,
    0.04904151308473655,
    -0.00023929699637607193,
    0.05
   ],
   [
    -1,
    5,
    -0.0014785919427173777,
    -0.0014759518393791624,
    0.051064072575099706,
    0.048597840078933814,
    -0.0005660421766113147,
    0.05
   ],
   [
    -1,
    6,
    -0.0016362292571946435,
    -0.0016338293381840641,
    0.05089814031124695,
    0.04844785584817766,
    -0.0006623403668065725,
    0.05
   ]
  ],
  "D3_L20": [
   [
    1,
    1,
    -0.006987233947915006,
    -0.006985323949465228,
    -0.05427341903149063,
    -0.057334598124513,
    -0.007059165146926036,
    -0.05
   ],
   [
    1,
    2,
    -0.0026715380854782644,
    -0.0026698653689011765,
    -0.0501632318999502,
    -0.05280336168646036,
    -0.0028078640659945403,
    -0.05
   ],
   [
    1,
    3,
    -0.001102647430188619,
    -0.0011011908533953524,
    -0.048669050106860845,
    -0.05115625165367782,
    -0.00130549299214161,
    -0.05
   ],
   [
    1,
    4,
    -0.00034118889817758014,
    -0.0003399238136699697,
    -0.047943851399807946,
    -0.05035692039256274,
    -0.0006056487536216209,
    -0.05
   ],
   [
    1,
    5,
    0.00010698444475054636,
    0.00010808121295844482,
    -0.047517019582757,
    -0.049886514602959696,
    -0.00021404018511527164,
    -0.05
   ],
   [
    1,
    6,
    0.00040166262347329784,
    0.000402612879727573,
    -0.047236373657571606,
    -0.049577256016482864,
    2.7875616521606296e-05,
    -0.05
   ],
   [
    1,
    7,
    0.0006060803010978777,
    0.0006069059168031513,
    -0.047041690126847224,
    -0.049362748094241064,
    0.00018293681252272458,
    -0.05
   ],
   [
    1,
    8,
    0.0007471679280208754,
    0.0007478937904387047,
    -0.046907320938868624,
    -0.04921471066590884,
    0.0002801427919159292,
    -0.05
   ],
   [
    1,
    9,
    0.0008332998045172957,
    0.0008339584600801242,
    -0.04682529056840823,
    -0.04912434266449536,
    0.0003339241187014234,
    -0.05
   ],
   [
    1,
    10,
    0.000862796117653053,
    0.0008634305948358645,
    -0.04679719883753941,
    -0.049093396889343256,
    0.00035116682257593774,
    -0.05
   ],
   [
    -1,
    1,
    0.007718266724922929,
    0.007720619356311528,
    0.06075321708602277,
    0.057334598124513,
    0.007802233291833128,
    0.05
   ],
   [
    -1,
    2,
    0.002948852394213608,
    0.0029509031213165574,
    0.05573279445206736,
    0.05280336168646036,
    0.003103427960484019,
    0.05
   ],
   [
    -1,
    3,
    0.001215323027309756,
    0.0012171053882716841,
    0.053908031605569784,
    0.05115625165367782,
    0.0014429129611704264,
    0.05
   ],
   [
    -1,
    4,
    0.00037415879064317537,
    0.00037570517766645863,
    0.053022597948560124,
    0.05035692039256274,
    0.0006694010935162164,
    0.05
   ],
   [
    -1,
    5,
    -0.00012079776688900257,
    -0.00011945815410385948,
    0.05250159244726449,
    0.049886514602959696,
    0.00023657067419986737,
    0.05
   ],
   [
    -1,
    6,
    -0.00044615298242778323,
    -0.00044499307617121975,
    0.052159114194190836,
    0.049577256016482864,
    -3.080988455836859e-05,
    0.05
   ],
   [
    -1,
    7,
    -0.0006717977526689012,
    -0.0006707905893310968,
    0.05192159402227814,
    0.049362748094241064,
    -0.00020219327062891737,
    0.05
   ],
   [
    -1,
    8,
    -0.0008275042353627082,
    -0.0008266192544110019,
    0.05175769290239929,
    0.04921471066590884,
    -0.00030963143262136873,
    0.05
   ],
   [
    -1,
    9,
    -0.0009225460203980962,
    -0.0009217433401563697,
    0.05165764918724314,
    0.04912434266449536,
    -0.0003690739374489352,
    0.05
   ],
   [
    -1,
    10,
    -0.000955090860117247,
    -0.0009543177970767835,
    0.05162339155336704,
    0.049093396889343256,
    -0.0003881316582149012,
    0.05
   ]
  ]
 },
 "R_reduction_full_operator": {
  "D1_L32": {
   "a": 8.881784197001252e-15,
   "full_grip": 3.552713678800501e-15
  },
  "D2_L8": {
   "a": 1.2434497875801753e-14,
   "full_grip": 1.2878587085651816e-14
  },
  "D3_L6": {
   "a": 3.552713678800501e-14,
   "full_grip": 1.6431300764452317e-14
  }
 },
 "S_asymptotics": {
  "2D_L64_ell0_minus_lnL_over_2pi": 0.04875850142855509,
  "2D_L128_ell0_minus_lnL_over_2pi": 0.048763848880206684,
  "2D_L256_ell0_minus_lnL_over_2pi": 0.04876518701309496,
  "2D_L512_ell0_minus_lnL_over_2pi": 0.048765521625585895,
  "2D_L512_ell(j)_minus_(3.10)_form": {
   "2": 0.00428099887944533,
   "4": 0.0009853796100285361,
   "8": 0.00021494543781747089,
   "16": 5.237185890971663e-05,
   "32": 1.489079006741778e-05,
   "64": 3.384913195475514e-05
  },
  "3D_L32_ell0": 0.2456765511865084,
  "3D_L32_ell0_plus_xi3_over_L": 0.2527323299766616,
  "3D_L64_ell0": 0.24920328389380492,
  "3D_L64_ell0_plus_xi3_over_L": 0.2527311732888815,
  "3D_L96_ell0": 0.25037913173616777,
  "3D_L96_ell0_plus_xi3_over_L": 0.2527310579995522,
  "3D_L96_rel_dev_of_ell(j)_from_(3.11)_form": {
   "2": 0.07792726294245632,
   "4": 0.01983248840345087,
   "8": 0.004125206484064414,
   "16": 0.0013984766296878245,
   "24": 0.0035316337626140827
  }
 }
}

Interested in this research?

Join the readers' list and get an e-mail when a new stage or a correction is published. No password, and you can unsubscribe at any time.