# Re-run of the orthant (plain R^3_+) antinorm search at the parameters of the two final
# certificates, and at weaker rates, to confirm the round-3 observation that it does not close.
# 2026-10-17: uses antinorm_exact.search (same algorithm and parameters as antinorm.search, with the three
# axis points) instead of antinorm.search, whose qhull dual vertices are not reliable for these point sets
# (see antinorm_exact.py and orthant_noaxis.py).  Option --qhull runs the old antinorm.search instead.
import sys, time
if '--qhull' in sys.argv:
    from antinorm import search as _s
    search = lambda *a, **k: _s(*a, **k) + ("(qhull version)",)
else:
    from antinorm_exact import search
for (y, b1, b2) in [(0.345, -0.3622113366699219, -0.35), (0.34, -0.33319140869140623, -0.38),
                    (0.345, -0.45, -0.35), (0.345, -0.60, -0.35), (0.345, -1.0, -0.35)]:
    t0 = time.time()
    ok, Z, r, why = search(y, b1, b2, maxpts=4000, maxround=60, verbose='-v' in sys.argv)
    print(f"y={y} log2 e1={b1} log2 e2={b2}: closed={ok} after {r} rounds, |Z|={len(Z)}; {why} [{time.time() - t0:.0f}s]", flush=True)