Wilson-Cowan excitatory-inhibitory

PARITYS0 · dim 2

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 →

2D Wilson-Cowan E/I model with separate population sigmoids. Excitatory and inhibitory firing rate dynamics with τ_E=2.5 ms, τ_I=3.75 ms. Canonical minimal neural mass model.

Neuroscience

Problem definition

Wilson & Cowan (1972), Biophys. J.

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 _stable_expit(x):
    """Numerically stable sigmoid 1/(1+exp(-x)), scalar or array."""
    x = np.asarray(x, dtype=float)
    return np.where(
        x >= 0,
        1.0 / (1.0 + np.exp(-x)),
        np.exp(x) / (1.0 + np.exp(x)),
    )

def rhs(t: float, y: np.ndarray) -> np.ndarray:
    r_E, r_I = y[0], y[1]
    input_E = c_EE * G_E * r_E - c_EI * r_I + I_ext * G_E
    input_I = c_IE * r_E
    F_E = float(_stable_expit(sigma_E * (input_E - mu_E)))
    F_I = float(_stable_expit(sigma_I * (input_I - mu_I)))
    dr_E = (-r_E + F_E) / tau_E
    dr_I = (-r_I + F_I) / tau_I
    return np.array([dr_E, dr_I])
Parameters
  • G_E = 1
  • I_ext = 1.5
  • c_EE = 16
  • c_EI = 12
  • c_IE = 15
  • mu_E = 1.3
  • mu_I = 2
  • sigma_E = 4
  • sigma_I = 3.7
  • tau_E = 0.0025
  • tau_I = 0.00375
Initial condition
y(0) = [0.100001, 0.050001]
Horizon
t ∈ [0, 2]

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 LSODA

Noisy best: SciPy BDF

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: Wilson-Cowan excitatory-inhibitory (wilson-cowan-excitatory-inhibitory)

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
1SciPy LSODASciPy
100%
16.03164 ms1.000
2CVODE Adamsexternal
100%
16.023910 ms1.000
3SolvSRK
100%
16.02875 ms1.000
4SciPy RadauSciPy
100%
15.51,20732 ms0.988
5SciPy BDFSciPy
100%
14.147319 ms0.955
6TRBDF2external
100%
13.97274.70 s0.950
7FBDFexternal
100%
12.55534.61 s0.917
8CVODE BDFexternal
100%
11.230811 ms0.886
9Vern9external
100%
10.83,2983.67 s0.876
10Vern7external
100%
9.92,1623.66 s0.855
11SciPy DOP853SciPy
100%
9.01,79023 ms0.832
12Tsit5external
100%
8.02,070863 ms0.809
13SciPy RK45SciPy
100%
7.91,98227 ms0.808
14SciPy RK23SciPy
100%
7.62,25534 ms0.801

At Clean, best balanced arm is SciPy LSODA · 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

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_wilson_cowan_excitatory_inhibitory_2026,
  title        = {Resonix Evidence Portal: Wilson-Cowan excitatory-inhibitory},
  author       = {{Resonix Labs (Canada) Inc.}},
  year         = {2026},
  howpublished = {\url{https://resonix.tech/evidence/problems/wilson-cowan-excitatory-inhibitory}},
  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