SIDE24 d = 5: exact closed form of the cone moment D_4 (closed forms of c_{d,ref} for d ≤ 5, conditional on #199/#200/#201 and SIDE24)

.github/workflows/side24-d5-closed-form.yml

@@ -0,0 +1,80 @@
+name: SIDE24 d5 closed form
+
+'on':
+  pull_request:
+    paths:
+      - 'coefficients/side24_d5_closed_form_20260930/**'
+      - '.github/workflows/side24-d5-closed-form.yml'
+      - 'imports/lifetime_parent_20260925/UNIFORM_MATRIX_CAP_AND_LIFETIME.md'
+      - 'coefficients/side24_v1/PROOF.md'
+  workflow_dispatch:
+
+permissions:
+  contents: read
+
+jobs:
+  certify-and-reject-mutants:
+    runs-on: ubuntu-latest
+    timeout-minutes: 15
+    steps:
+      - uses: actions/checkout@11d5960a326750d5838078e36cf38b85af677262
+        with:
+          persist-credentials: false
+      - uses: actions/setup-python@a26af69be951a213d495a4c3e4e4022e16d87065
+        with:
+          python-version: '3.11.16'
+      - name: Verify sources and pins, replay both modes, reject mutants
+        env:
+          PYTHONDONTWRITEBYTECODE: '1'
+        run: |
+          set -euo pipefail
+          python -B -S - <<'PY'
+          import concurrent.futures, hashlib, json, pathlib, subprocess, sys
+          root = pathlib.Path('coefficients/side24_d5_closed_form_20260930').resolve()
+          manifest_doc = json.loads((root / 'SOURCE_FILES.json').read_text())
+          manifest = manifest_doc['files']
+          present = sorted(p.name for p in root.iterdir())
+          expected_names = sorted([e['path'] for e in manifest] + ['SOURCE_FILES.json'])
+          if present != expected_names:
+              raise ValueError('packet tree differs from manifest: ' + repr(present))
+          for entry in manifest:
+              path = root / entry['path']
+              if pathlib.PurePosixPath(entry['path']).name != entry['path'] or path.is_symlink():
+                  raise ValueError('flat regular source required: ' + entry['path'])
+              data = path.read_bytes()
+              if len(data) != entry['bytes'] or hashlib.sha256(data).hexdigest() != entry['sha256']:
+                  raise ValueError('source identity mismatch: ' + entry['path'])
+          for pin in manifest_doc['pins_on_main']:
+              if not pin.get('workflow_checked'):
+                  continue
+              p = pathlib.Path(pin['path'])
+              if p.is_symlink() or not p.is_file():
+                  raise ValueError('pinned source missing: ' + pin['path'])
+              data = p.read_bytes()
+              blob = hashlib.sha1(b'blob %d\0' % len(data) + data).hexdigest()
+              if blob != pin['git_blob']:
+                  raise ValueError('pinned source drifted: ' + pin['path'] + ' ' + blob)
+          script = str(root / 'closed_form_d4.py')
+          expected = (root / 'RESULTS.json').read_bytes()
+          mutants = ['shift-variance', 'vandermonde', 'erfc-sign', 'jacobian', 'mehta', 'by-parts', 'one-erf', 'two-erf', 'jk-init', 'gauss-half', 'boundary-term', 'atom-pi', 'cumulative-init', 'cumulative-recurrence']
+          def mutant_rc(args):
+              flags, m = args
+              return m, flags, subprocess.run([sys.executable, *flags, script, '--mutant', m],
+                                              capture_output=True, timeout=600).returncode
+          for flags in (['-B', '-S'], ['-B', '-O', '-S']):
+              out = subprocess.run([sys.executable, *flags, script],
+                                   check=True, capture_output=True, timeout=600).stdout
+              if out != expected:
+                  raise RuntimeError('output differs from RESULTS.json under ' + ' '.join(flags))
+              doc = json.loads(out)
+              if not (doc['passed'] and len(doc['checks']) == 16 and all(doc['checks'].values())):
+                  raise RuntimeError('certification did not pass under ' + ' '.join(flags))
+              if doc['closed_forms']['D_4'] != '6695/54  +  -405/32 sqrt(6)  +  -1375/72 sqrt(21) pi^(-1)  +  -405/16 sqrt(6) pi^(-1) arctan(1/4 sqrt(14))  +  -6695/27 pi^(-1) arctan(1/7 sqrt(21))  +  405/16 sqrt(6) pi^(-1) arctan(1/7 sqrt(14))' or doc['closed_forms']['Z_4'] != '1536 pi^(1)':
+                  raise RuntimeError('closed form differs from the stated one')
+              with concurrent.futures.ThreadPoolExecutor(max_workers=4) as pool:
+                  for m, fl, rc in pool.map(mutant_rc, [(flags, m) for m in mutants]):
+                      if rc != 1:
+                          raise RuntimeError('mutant ' + m + ' not rejected under ' + ' '.join(fl))
+          print(json.dumps({'scientific_effect': 'NONE', 'mathematical_acceptance': False, 'exact_identity_checked': True, 'side24_v1': 'unchanged'}))
+          PY
+          git diff --exit-code