Linear eigsweep 28D κ=3e4

PARITYS3 · dim 28

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 →

Diagonal linear ODE ẏ=diag(λ)·y with 28 eigenvalues log-spaced from −1 to −3e+04. Stiffness ratio = κ. Boundary mapping: stiffness × dimension sweep.

Synthetic & linear-stiff benchmarks

Problem definition

Enright & Pryce (1987) ACM TOMS; Hairer & Wanner (1996) Solving ODEs II, §IV.10

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 rhs(t: float, y: np.ndarray) -> np.ndarray:
    return lam * y
Parameters
  • lam = [-1, -1.46493820619, -2.14604394796, -3.14382177153, -4.60550462657, -6.74677968626, …] [shape=(28,), min=-30000, max=-1]
Initial condition
y(0) = repeat(1, 28) [shape=(28,)]
Horizon
t ∈ [0, 50]

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: none

Recommendation snapshot

Clean best: SolvSRK

Noisy best: SolvSRK

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: Linear eigsweep 28D κ=3e4 (linear-eigsweep-28d-3e4)

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

#SolverSurvivalSCDnfevWallScore
1SolvSRK
100%
3.5153,583972 ms0.702
2SciPy BDFSciPy
100%
0.01,11137 ms0.619
3SciPy RadauSciPy
100%
0.03,33364 ms0.619
4SciPy RK45SciPy
100%
0.03,171,65619.10 s0.619
5SciPy LSODASciPy
100%
0.02,6547 ms0.619
6SciPy DOP853SciPy
100%
0.02,815,57415.15 s0.619
7SciPy RK23SciPy
100%
0.01,792,10013.53 s0.619
8CVODE BDFexternal
100%
0.094512 ms0.619
9CVODE Adamsexternal
100%
0.06,42045 ms0.619
Tsit5external
0%

At Clean, best balanced arm is SolvSRK.

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

Profile your problem for 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

Run the winner on your machine

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.

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_linear_eigsweep_28d_3e4_2026,
  title        = {Resonix Evidence Portal: Linear eigsweep 28D κ=3e4},
  author       = {{Resonix Labs (Canada) Inc.}},
  year         = {2026},
  howpublished = {\url{https://resonix.tech/evidence/problems/linear-eigsweep-28d-3e4}},
  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