Re-run of release/reproducibility from an extracted copy of release/source.zip (2026-10-18, Python 3.13.5,
NumPy 2.5.2, SciPy 1.18.1, SymPy 1.14.0, mpmath 1.3.0; nice -n 10, at most 4 programs at a time).
Each program was run in its own folder with 'python3 -B <command>', standard output to a new file, compared with
the recorded output (cmp; then diff). Columns: recorded output, command, exit status, wall time, comparison.

check_note_output	check_note.py	exit=0	1s	IDENTICAL
disk_check	disk_check.py	exit=0	201s	DIFFERENT
disk_front	disk_front.py	exit=0	2s	IDENTICAL
disk_nondeg	disk_nondeg.py	exit=0	0s	IDENTICAL
ellipse_check	ellipse_check.py	exit=0	7s	DIFFERENT
ellipse_checks	ellipse_checks.py	exit=0	11s	IDENTICAL
ellipse_freq	ellipse_freq.py	exit=0	1s	IDENTICAL
ellipse_front	ellipse_front.py	exit=0	73s	IDENTICAL
scout_front_copy	scout_front_copy.py	exit=0	8s	IDENTICAL
theoremA_check	theoremA_check.py	exit=0	19s	DIFFERENT
theoremA_geometry	theoremA_geometry.py	exit=0	1s	IDENTICAL
vb_ball_part3	vb_disk_ball.py 3	exit=0	220s	DIFFERENT
vb_ball_part4	vb_disk_ball.py 4	exit=0	395s	DIFFERENT
vb_disk_length	vb_disk_length.py	exit=0	10s	IDENTICAL
vb_disk_part1	vb_disk_ball.py 1	exit=0	6s	DIFFERENT
vb_disk_part2	vb_disk_ball.py 2	exit=0	65s	DIFFERENT
vb_ellipse_freq	vb_ellipse.py freq	exit=0	29s	DIFFERENT
vb_ellipse_front_a1_b0.25	vb_ellipse.py front 1 0.25	exit=0	812s	DIFFERENT
vb_ellipse_front_a1_b0.6	vb_ellipse.py front 1 0.6	exit=0	401s	DIFFERENT
vb_ellipse_front_a1_b0.95	vb_ellipse.py front 1 0.95	exit=0	287s	DIFFERENT
vb_ellipse_front_a2_b1.732	vb_ellipse.py front 2 1.7320508075688772	exit=0	153s	DIFFERENT
vb_ellipse_limits	vb_ellipse.py limits	exit=0	1s	IDENTICAL
vb_freq_wrong_input	vb_freq_wrong_input.py	exit=0	2s	IDENTICAL
vb_theoremA_long	vb_theoremA.py 1.0 100 1000 600	exit=0	151s	DIFFERENT
vb_theoremA_main	vb_theoremA.py 1.0 60 2000 200	exit=0	66s	DIFFERENT
vb_theoremA_pocket	vb_theoremA.py 1.0 40 1500 300 7 pocket	exit=0	57s	DIFFERENT
vb_theoremA_wrong_teeth0.6	vb_theoremA.py 0.6 30 1000 200	exit=0	30s	DIFFERENT

Differences (diff new recorded), for every program not byte-identical:
--- vb_ellipse_front_a2_b1.732
9c9
< elapsed 152.6s
---
> elapsed 139.8s
--- vb_ellipse_front_a1_b0.95
9c9
< elapsed 286.2s
---
> elapsed 266.2s
--- vb_disk_part2
9c9
< elapsed 64.7s
---
> elapsed 58.6s
--- disk_check
9c9
< elapsed 200.6s
---
> elapsed 203.7s
--- vb_ellipse_front_a1_b0.6
9c9
< elapsed 400.1s
---
> elapsed 361.3s
--- vb_theoremA_main
4,5c4,5
<   iterations 197, steps shorter than 1e-7: 0, still active at stop: 0
<   photons 120000, reflections 11329385 (ellipse 5364673, tooth sides 4371520, flat tops 311824, rounded corners 1281368), undefined steps 0
---
>   iterations 175, steps shorter than 1e-7: 0, still active at stop: 0
>   photons 120000, reflections 11499203 (ellipse 5312891, tooth sides 4532376, flat tops 326898, rounded corners 1327038), undefined steps 0
7c7
<   max |x| at axis crossings: 0.865692  (claim: crossings only in I_0=(-7/8,7/8) when teeth at +-1)
---
>   max |x| at axis crossings: 0.865691  (claim: crossings only in I_0=(-7/8,7/8) when teeth at +-1)
9,11c9,11
<   invariant (ii): max |x-intercept| of lines of segments meeting U: 0.999999944; photons with value >=1: 0
<   flat-top reflections: 311824, |x| range [0.875008, 1.000000] (claim: within [7/8,1))
< elapsed 65.6s
---
>   invariant (ii): max |x-intercept| of lines of segments meeting U: 0.999998575; photons with value >=1: 0
>   flat-top reflections: 326898, |x| range [0.875000, 0.999999] (claim: within [7/8,1))
> elapsed 52.7s
--- vb_ball_part3
7c7
< elapsed 219.5s
---
> elapsed 155.1s
--- vb_theoremA_long
4,5c4,5
<   iterations 500, steps shorter than 1e-7: 0, still active at stop: 0
<   photons 100000, reflections 28481818 (ellipse 13409198, tooth sides 10971416, flat tops 828397, rounded corners 3272807), undefined steps 0
---
>   iterations 416, steps shorter than 1e-7: 0, still active at stop: 0
>   photons 100000, reflections 28564691 (ellipse 13382453, tooth sides 11049066, flat tops 833222, rounded corners 3299950), undefined steps 0
9,11c9,11
<   invariant (ii): max |x-intercept| of lines of segments meeting U: 0.999999759; photons with value >=1: 0
<   flat-top reflections: 828397, |x| range [0.875001, 1.000000] (claim: within [7/8,1))
< elapsed 150.9s
---
>   invariant (ii): max |x-intercept| of lines of segments meeting U: 0.999999911; photons with value >=1: 0
>   flat-top reflections: 833222, |x| range [0.875003, 1.000000] (claim: within [7/8,1))
> elapsed 144.5s
--- vb_theoremA_pocket
5,6c5,6
<   iterations 337, steps shorter than 1e-7: 0, still active at stop: 0
<   photons 60000, reflections 8992915 (ellipse 3897604, tooth sides 3631384, flat tops 315383, rounded corners 1148544), undefined steps 0
---
>   iterations 257, steps shorter than 1e-7: 0, still active at stop: 0
>   photons 60000, reflections 8995075 (ellipse 3896811, tooth sides 3632688, flat tops 317226, rounded corners 1148350), undefined steps 0
8c8
<   max |x| at axis crossings: 0.865692  (claim: crossings only in I_0=(-7/8,7/8) when teeth at +-1)
---
>   max |x| at axis crossings: 0.865691  (claim: crossings only in I_0=(-7/8,7/8) when teeth at +-1)
10,12c10,12
<   invariant (ii): max |x-intercept| of lines of segments meeting U: 0.999999830; photons with value >=1: 0
<   flat-top reflections: 315383, |x| range [0.875002, 1.000000] (claim: within [7/8,1))
< elapsed 57.5s
---
>   invariant (ii): max |x-intercept| of lines of segments meeting U: 0.999999276; photons with value >=1: 0
>   flat-top reflections: 317226, |x| range [0.875003, 0.999999] (claim: within [7/8,1))
> elapsed 47.5s
--- vb_theoremA_wrong_teeth0.6
4,11c4,11
<   iterations 329, steps shorter than 1e-7: 0, still active at stop: 0
<   photons 30000, reflections 3822248 (ellipse 1775034, tooth sides 1420372, flat tops 148184, rounded corners 478658), undefined steps 0
<   photons meeting V: 15864; segments leaving D (midpoint test): 0
<   max |x| at axis crossings: 1.999994  (claim: crossings only in I_0=(-7/8,7/8) when teeth at +-1)
<   photons with an axis crossing |x|>=7/8: 23868; with |x|>=1: 20794
<   invariant (ii): max |x-intercept| of lines of segments meeting U: 785566.308236835; photons with value >=1: 20794
<   flat-top reflections: 148184, |x| range [0.475070, 0.724956] (claim: within [7/8,1))
< elapsed 30.6s
---
>   iterations 391, steps shorter than 1e-7: 0, still active at stop: 0
>   photons 30000, reflections 3838717 (ellipse 1781943, tooth sides 1425173, flat tops 147719, rounded corners 483882), undefined steps 0
>   photons meeting V: 16327; segments leaving D (midpoint test): 0
>   max |x| at axis crossings: 2.000000  (claim: crossings only in I_0=(-7/8,7/8) when teeth at +-1)
>   photons with an axis crossing |x|>=7/8: 24278; with |x|>=1: 21219
>   invariant (ii): max |x-intercept| of lines of segments meeting U: 1064551.669675374; photons with value >=1: 21219
>   flat-top reflections: 147719, |x| range [0.475077, 0.724967] (claim: within [7/8,1))
> elapsed 28.1s
--- ellipse_check
0a1,3
> ellipse_check.py:44: IntegrationWarning: Extremely bad integrand behavior occurs at some points of the
>   integration interval.
>   Pn = 1/c*quad(lambda n: 1/np.sqrt(np.sin(n)**2+eps**2), 0, 2*np.pi, limit=200)[0]
29c32
< elapsed 6.3s
---
> elapsed 8.2s
--- vb_disk_part1
7c7
< elapsed 6.3s
---
> elapsed 6.0s
--- theoremA_check
12c12
< elapsed 18.9s
---
> elapsed 23.5s
--- vb_ellipse_freq
41c41
< elapsed 28.4s
---
> elapsed 28.0s
--- vb_ball_part4
7c7
< elapsed 394.9s
---
> elapsed 395.2s
--- vb_ellipse_front_a1_b0.25
9c9
< elapsed 810.9s
---
> elapsed 740.5s

Standard error of ellipse_check.py in the re-run (the recorded output contains the same SciPy warning, with the
path shortened to the file name):
<package>/reproducibility/verification_run_A/ellipse_check.py:44: IntegrationWarning: Extremely bad integrand behavior occurs at some points of the
  integration interval.
  Pn = 1/c*quad(lambda n: 1/np.sqrt(np.sin(n)**2+eps**2), 0, 2*np.pi, limit=200)[0]
