Brusselator
PARITYS1 · dim 2No 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 chemical oscillator; mild stiffness, widely used benchmark
Problem definition
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.
Show canonical RHS excerpt
def _brusselator_rhs(t, y, A=1.0, B=3.0):
x, z = y
return np.array([A + x**2 * z - (B + 1) * x, B * x - x**2 * z])
def rhs(t, y):
return _brusselator_rhs(t, y)- Parameters
- A = 1
- B = 3
- Initial condition
- y(0) = [1, 1]
- 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: low
Default noise: none
Recommendation snapshot
Clean best: SciPy Radau
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: Brusselator (brusselator)
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 | SciPy RadauSciPy | 100% | 10.6 | 6,025 | 115 ms | 0.872 |
| 2 | Vern9external | 100% | 10.4 | 1,442 | 4.01 s | 0.866 |
| 3 | SolvSRK | 100% | 9.6 | 3,782 | 9 ms | 0.848 |
| 4 | SciPy DOP853SciPy | 100% | 8.7 | 1,298 | 8 ms | 0.827 |
| 5 | Vern7external | 100% | 8.4 | 1,572 | 3.93 s | 0.820 |
| 6 | SciPy RK45SciPy | 100% | 8.1 | 1,718 | 12 ms | 0.811 |
| 7 | Tsit5external | 100% | 7.9 | 1,668 | 849 ms | 0.808 |
| 8 | CVODE Adamsexternal | 100% | 7.4 | 784 | 12 ms | 0.794 |
| 9 | SciPy RK23SciPy | 100% | 7.1 | 7,499 | 64 ms | 0.788 |
| 10 | SciPy LSODASciPy | 100% | 6.9 | 957 | 4 ms | 0.784 |
| 11 | SciPy BDFSciPy | 100% | 6.8 | 2,465 | 73 ms | 0.781 |
| 12 | FBDFexternal | 100% | 6.5 | 1,972 | 5.05 s | 0.774 |
| 13 | CVODE BDFexternal | 100% | 6.5 | 900 | 13 ms | 0.773 |
| 14 | TRBDF2external | 100% | 4.4 | 6,491 | 5.17 s | 0.725 |
At Clean, best balanced arm is SciPy Radau · SolvSRK survival 100%, SCD 9.6.
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_brusselator_2026,
title = {Resonix Evidence Portal: Brusselator},
author = {{Resonix Labs (Canada) Inc.}},
year = {2026},
howpublished = {\url{https://resonixusa.com/evidence/problems/brusselator}},
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