_coder@spark-0f87:~/Dev/pressure-wave$ uv run simulate_blasts.py --config config.yaml \
> --set blast.cube_extents_m=20 --set blast.cube_segments=9 \
> --set blast.min_dist_m=2.5 --set rays.resolution_arcmin=5 --set rays.max_bounces=5 \
> --set output.file=/tmp/warp_match.csv
2026-07-21 20:59:59,642 - INFO - Config: config.yaml -> output: /tmp/warp_match.csv
2026-07-21 20:59:59,643 - INFO - Cube: center=(0.0, -4.0, 0.0) extents=20.0m segments=9 | rays: 5.0 arcmin, max_bounces=5, fan=sensor_cone
Warp 1.15.0 initialized:
CUDA Toolkit 12.9, Driver 13.0
Devices:
"cpu" : "aarch64"
"cuda:0" : "NVIDIA GB10" (122 GiB, sm_121, mempool enabled)
Kernel cache:
/home/_coder/.cache/warp/1.15.0
2026-07-21 21:00:00,046 - INFO - Warp device: cuda
2026-07-21 21:00:00,080 - INFO - triangulating faces
2026-07-21 21:00:00,114 - INFO - triangulating faces
2026-07-21 21:00:00,121 - INFO - 727 blast points survive the min-distance filter; grid=30x30 rays/explosion (device=cuda).
Module warp_core 2a1bb02 load on device 'cuda:0' took 0.86 ms (cached)
2026-07-21 21:00:00,321 - INFO - Finished config: full_armor
2026-07-21 21:00:00,322 - INFO - Finished config: helmet_only
2026-07-21 21:00:00,323 - INFO - Finished config: vest_only
2026-07-21 21:00:00,323 - INFO - Finished config: no_armor
2026-07-21 21:00:00,324 - INFO - Wrote 727 rows to /tmp/warp_match.csv
_coder@spark-0f87:~/Dev/pressure-wave$ python3 - <<'PY'
> import csv
> def load(p): return {(r['bp_x'],r['bp_y'],r['bp_z']):r for r in csv.DictReader(open(p))}
> w, j = load('/tmp/warp_match.csv'), load('/tmp/bvh.csv')
> cols=['full_armor','helmet_only','vest_only','no_armor']
> import statistics
> for c in cols:
> d=[abs(float(w[k][c])-float(j[k][c])) for k in w if k in j]
> print(f"{c:12s} max|Δ|={max(d):.2e} mean|Δ|={statistics.fmean(d):.2e}")
> PY
full_armor max|Δ|=8.36e-03 mean|Δ|=4.91e-05
helmet_only max|Δ|=8.36e-03 mean|Δ|=5.00e-05
vest_only max|Δ|=1.08e-04 mean|Δ|=3.89e-06
no_armor max|Δ|=1.08e-04 mean|Δ|=3.89e-06