Unexecuted reproducibility supplement

Copy this listing to 02_e_v2.py. Dependencies are Python, NumPy, and SciPy; optional SDPs require CVXPY and SCS. Installed versions are recorded when run.

python en02_e_v2.py –mode synthetic python en02_e_v2.py –mode pilot python en02_e_v2.py –mode pilot –sdp –tex En02-E.tex –compile

The synthetic command requires no radial calculation. The pilot runs synthetic checks first, then the field layer and linked operational layer. The source filename is consistently 02-E.tex.

Array artifacts have companion basis metadata. Raw solver artifacts are saved in the solve wrapper's block, before status or feasibility acceptance. No passing status is written if a requested check fails.

"""En02-E v2. Unexecuted implementation; no results are prepopulated.""" import argparse import datetime import hashlib import itertools as it import json import math from pathlib import Path import shutil import subprocess import sys import traceback import numpy as np import scipy from scipy.integrate import solve_ivp from scipy.special import lambertw

I = np.eye(2,dtype=complex) X = np.array([[0,1],[1,0]],complex) Y = np.array([[0,-1j],[1j,0]],complex) Z = np.diag([1,-1]).astype(complex) H = np.array([[1,1],[1,-1]],complex)/np.sqrt(2) SD = np.diag([1,-1j]).astype(complex) CZ = np.diag([1,1,1,-1]).astype(complex) CX = np.array([[1,0,0,0],[0,1,0,0], [0,0,0,1],[0,0,1,0]],complex) ROOT = None TOL, STOL, GAP = 1e-9, 3e-4, 5e-4 DATA = {"checks":{}, "residuals":{}, "solves":[], "field":{}, "operational":[], "occupancy":[], "certified_numerics":False}

def need(ok,msg): if not ok: raise AssertionError(msg)

def clean(x): if isinstance(x,np.ndarray): return clean(x.tolist()) if isinstance(x,np.generic): return clean(x.item()) if isinstance(x,complex): return "real":clean(x.real),"imag":clean(x.imag) if isinstance(x,float) and not np.isfinite(x): return str(x) if isinstance(x,dict): return str(k):clean(v) for k,v in x.items() if isinstance(x,(list,tuple)): return [clean(v) for v in x] if x is None or isinstance(x,(str,int,float,bool)): return x return repr(x)

def jwrite(name,obj): (ROOT/(name+".json")).write_text( json.dumps(clean(obj),indent=2,allow_nan=False))

def dump(): jwrite("records",DATA)

def save(name,**arrays): np.savez_compressed(ROOT/(name+".npz"),**arrays)

def herm(A): return (A+A.conj().T)/2

def check(name,A,B,tol=TOL): need(name not in DATA["checks"],"duplicate check "+name) A,B = np.asarray(A),np.asarray(B) need(A.shape==B.shape,"shape "+name) err = float(np.max(np.abs(A-B))) DATA["checks"][name] = err need(np.isfinite(err) and err<=tol,name+" residual "+str(err))