No 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 →
7-body gravitational system; non-stiff, high-dimensional, long integration
Canonical benchmark implementation
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.
def _pleiades_rhs(t, y):
n = 7
x = y[:n]
z = y[n:2*n]
vx = y[2*n:3*n]
vz = y[3*n:4*n]
m = np.arange(1, n + 1, dtype=float)
ax = np.zeros(n)
az = np.zeros(n)
for i in range(n):
for j in range(n):
if i != j:
rij3 = ((x[i] - x[j])**2 + (z[i] - z[j])**2) ** 1.5
ax[i] += m[j] * (x[j] - x[i]) / rij3
az[i] += m[j] * (z[j] - z[i]) / rij3
return np.concatenate([vx, vz, ax, az])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: high
Default noise: low
Recommendation snapshot
Clean best: SciPy Radau
Noisy best: SciPy RK23
Coverage
14 solver arms · clean + 5 noise levels
Ranked on survival, precision, and speed
20 seeds/cell default · 14 arms · TRL 4–5 · simulation-lab validated · this page: Pleiades 7-body problem (pleiades-7-body-problem)
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
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 | SciPy RadauSciPy | 100% | 8.3 | 11,193 | 957 ms | 0.817 |
| 2 | Vern7external | 100% | 7.5 | 2,772 | 3.96 s | 0.797 |
| 3 | SolvSRK | 100% | 7.5 | 10,133 | 397 ms | 0.797 |
| 4 | Vern9external | 100% | 7.2 | 2,402 | 4.12 s | 0.790 |
| 5 | Tsit5external | 100% | 6.6 | 2,964 | 1.14 s | 0.776 |
| 6 | SciPy DOP853SciPy | 100% | 6.5 | 2,222 | 157 ms | 0.775 |
| 7 | SciPy RK45SciPy | 100% | 5.8 | 2,474 | 176 ms | 0.756 |
| 8 | SciPy RK23SciPy | 100% | 5.5 | 14,348 | 1.06 s | 0.751 |
| 9 | SciPy LSODASciPy | 100% | 5.4 | 1,847 | 127 ms | 0.747 |
| 10 | CVODE Adamsexternal | 100% | 4.5 | 1,702 | 126 ms | 0.727 |
| 11 | FBDFexternal | 100% | 4.4 | 2,558 | 5.25 s | 0.724 |
| 12 | SciPy BDFSciPy | 100% | 4.3 | 2,995 | 302 ms | 0.722 |
| 13 | CVODE BDFexternal | 100% | 4.1 | 2,350 | 172 ms | 0.717 |
| 14 | TRBDF2external | 100% | 1.4 | 8,914 | 5.97 s | 0.653 |
At Clean, best balanced arm is SciPy Radau · SolvSRK survival 100%, SCD 7.5.
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 →
SolvScout · free
This page shows one published benchmark cell. SolvScout fingerprints your ODE, compares it to the full corpus, and recommends a solver with the same survival / precision / speed ranking you see here — including when a SciPy arm wins.
SolvSRK · 30-day trial
SolvSRK is the stiffness-adaptive integrator behind the SolvSRK column in these tables. Create an account, activate a machine, and take a 30-day trial — same binary you'd ship after purchase.
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_pleiades_7_body_problem_2026,
title = {Resonix Evidence Portal: Pleiades 7-body problem},
author = {{Resonix Labs (Canada) Inc.}},
year = {2026},
howpublished = {\url{https://resonix.tech/evidence/problems/pleiades-7-body-problem}},
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.}
}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