2D Proportional Navigation (Pure PN, N=3)
PARITYS0 · dim 4No clear winner. The survival gap is under 10 percentage points and the balanced-score gap is under 0.05, so neither SolvSRK nor the best baseline clears the win threshold. Either works - choose on cost, licensing, or integration effort. All verdicts →
2D proportional navigation engagement. Interceptor steers via pure PN (N=3) against a constant-velocity target. No actuator dynamics. Baseline engagement geometry.
Problem definition
Zarchan, 'Tactical and Strategic Missile Guidance' (6th ed., 2012), Ch. 2
Canonical RHS excerpt from the registered callable used for this benchmark cell. Expand it to verify the state equations; it is not a standalone runnable fixture.
Show canonical RHS excerpt
def _clamp(x: float, lo: float, hi: float) -> float:
return max(lo, min(hi, x))
def rhs(t, y):
rx, ry = y[0], y[1]
vx_i, vy_i = y[2], y[3]
d = np.empty(4)
R = np.sqrt(rx * rx + ry * ry)
if R < _RANGE_GUARD:
d[:] = 0.0
return d
lam = np.arctan2(ry, rx)
vx_rel = vx_i - v_tgt_x
vy_rel = vy_i - v_tgt_y
R_dot = (rx * vx_rel + ry * vy_rel) / R
lam_dot = (rx * vy_rel - ry * vx_rel) / (R * R)
lam_dot = _clamp(lam_dot, -_LOS_RATE_CLAMP, _LOS_RATE_CLAMP)
V_c = -R_dot
a_cmd = N * V_c * lam_dot
perp_x = -np.sin(lam)
perp_y = np.cos(lam)
d[0] = vx_i - v_tgt_x
d[1] = vy_i - v_tgt_y
d[2] = a_cmd * perp_x
d[3] = a_cmd * perp_y
return d- Parameters
- N = 3
- _LOS_RATE_CLAMP = 0.5
- _RANGE_GUARD = 0.1
- v_tgt_x = -200
- v_tgt_y = 0
- Initial condition
- y(0) = [5000, 2000, -250, 0]
- Horizon
- t ∈ [0, 20]
Canonical RHS excerpt captured from the same registered callable used for the published benchmark. Frozen closure values are summarized below; helper imports and solver settings are intentionally omitted.
Fingerprint
Spread: extreme
Default noise: low
Recommendation snapshot
Clean best: Vern9
Noisy best: SciPy LSODA
Coverage
14 solver arms · clean + 5 noise levels
Ranked on survival, precision, and speed
Versions & freeze
Methodology →- Freeze
- 2026-08-13
- libsolvsrk
- 2.3.0
- SciPy
- 1.14
- SUNDIALS
- CVODE (bundled backend)
20 seeds/cell default · 14 arms · TRL 4–5 · simulation-lab validated · this page: 2D Proportional Navigation (Pure PN, N=3) (2d-proportional-navigation-pure-pn-n-3)
Governed SolvTune benchmark freeze; per-arm medians only. RHS definitions and raw trial rows are not published.
Self-reported by Resonix Labs · not independently verified
Results matrix
Pick an objective and a noise level to rank all arms on survival, median SCD, median nfev, and median wall time. Medians across seeds.
Objective
Best overall trade-off of survival, precision, and speed.
Noise level
| # | Solver | Survival | SCD | nfev | Wall | Score |
|---|---|---|---|---|---|---|
| 1 | Vern9external | 100% | 13.9 | 114 | 4.10 s | 0.949 |
| 2 | Vern7external | 100% | 13.0 | 92 | 4.05 s | 0.928 |
| 3 | SolvSRK | 100% | 11.9 | 213 | 30 ms | 0.901 |
| 4 | SciPy RadauSciPy | 100% | 11.7 | 245 | 6 ms | 0.897 |
| 5 | Tsit5external | 100% | 10.8 | 66 | 803 ms | 0.875 |
| 6 | SciPy DOP853SciPy | 100% | 10.7 | 74 | <1 ms | 0.873 |
| 7 | SciPy RK45SciPy | 100% | 9.9 | 74 | <1 ms | 0.856 |
| 8 | SciPy LSODASciPy | 100% | 9.2 | 81 | <1 ms | 0.837 |
| 9 | FBDFexternal | 100% | 9.1 | 110 | 5.20 s | 0.835 |
| 10 | SciPy RK23SciPy | 100% | 8.8 | 386 | 4 ms | 0.829 |
| 11 | SciPy BDFSciPy | 100% | 8.5 | 108 | 4 ms | 0.820 |
| 12 | CVODE Adamsexternal | 100% | 8.2 | 54 | 7 ms | 0.815 |
| 13 | CVODE BDFexternal | 100% | 8.1 | 65 | 7 ms | 0.811 |
| 14 | TRBDF2external | 100% | 6.2 | 453 | 5.52 s | 0.767 |
At Clean, best balanced arm is Vern9 · SolvSRK survival 100%, SCD 11.9.
Values are medians across seeds, measured by Resonix Labs on Resonix hardware and not independently verified; nfev and wall are on reference lab hardware (indicative). Under injected noise only SolvSRK and the SciPy arms are run. How we measure accuracy → · Verification status →
Cite this page
Replace the access date. Pin the freeze ID and library versions when comparing against a later export. Cite it as what it is - a self-reported vendor benchmark, not an independently verified result. The note field says so; please keep it.
@misc{resonix_evidence_2d_proportional_navigation_pure_pn_n_3_2026,
title = {Resonix Evidence Portal: 2D Proportional Navigation (Pure PN, N=3)},
author = {{Resonix Labs (Canada) Inc.}},
year = {2026},
howpublished = {\url{https://resonixusa.com/evidence/problems/2d-proportional-navigation-pure-pn-n-3}},
note = {Self-reported vendor benchmark; internally generated by Resonix Labs and not independently verified. Accessed YYYY-MM-DD. Freeze 2026-08-13; libsolvsrk 2.3.0; SciPy 1.14.}
}Related
TRL 4–5 · simulation-lab validated · 398 problems · 14 solver arms · clean + 5 noise levels
Freeze: 2026-08-13 · scipy 1.14 · libsolvsrk 2.3.0 · Methodology
Self-reported by Resonix Labs · not independently verified · Verification status