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 →
Method-of-lines discretisation of 1-D Brusselator PDE, 100 cells (dim=200)
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 _brusselator_1d_rhs(t, y, N=100, A=1.0, B=3.0, alpha=0.02):
"""1D Brusselator reaction-diffusion, method of lines."""
dim = 2 * N
u = y[:N]
v = y[N:]
dx = 1.0 / (N + 1)
dx2 = dx * dx
du = np.zeros(N)
dv = np.zeros(N)
for i in range(N):
u_left = u[i-1] if i > 0 else 1.0
u_right = u[i+1] if i < N-1 else 1.0
v_left = v[i-1] if i > 0 else 3.0
v_right = v[i+1] if i < N-1 else 3.0
du[i] = alpha / dx2 * (u_left - 2*u[i] + u_right) + A + u[i]**2 * v[i] - (B+1)*u[i]
dv[i] = alpha / dx2 * (v_left - 2*v[i] + v_right) + B*u[i] - u[i]**2 * v[i]
return np.concatenate([du, dv])
def rhs(t, y):
return _brusselator_1d_rhs(t, y, N=100)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 BDF
Noisy best: SciPy BDF
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: Brusselator 1-D PDE (N=100) (brusselator-1-d-pde-n-100)
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 BDFSciPy | 100% | 16.0 | 218 | 47 ms | 1.000 |
| 2 | SciPy RadauSciPy | 100% | 16.0 | 234 | 56 ms | 1.000 |
| 3 | SciPy RK45SciPy | 100% | 16.0 | 50 | 10 ms | 1.000 |
| 4 | SciPy LSODASciPy | 100% | 16.0 | 4 | <1 ms | 1.000 |
| 5 | SciPy DOP853SciPy | 100% | 16.0 | 98 | 20 ms | 1.000 |
| 6 | SciPy RK23SciPy | 100% | 16.0 | 26 | 5 ms | 1.000 |
| 7 | CVODE BDFexternal | 100% | 16.0 | 206 | 54 ms | 1.000 |
| 8 | CVODE Adamsexternal | 100% | 16.0 | 206 | 56 ms | 1.000 |
| 9 | Tsit5external | 100% | 16.0 | 24 | 849 ms | 1.000 |
| 10 | SolvSRK | 100% | 16.0 | 214 | 54 ms | 1.000 |
At Clean, best balanced arm is SciPy BDF · SolvSRK survival 100%, SCD 16.0.
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_brusselator_1_d_pde_n_100_2026,
title = {Resonix Evidence Portal: Brusselator 1-D PDE (N=100)},
author = {{Resonix Labs (Canada) Inc.}},
year = {2026},
howpublished = {\url{https://resonix.tech/evidence/problems/brusselator-1-d-pde-n-100}},
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