Problem definition
Jansen & Rit (1995); Ableidinger et al. (2017)
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 _jr_sigmoid(v, v_max=5.0, r=0.56, v0=6.0):
"""Jansen-Rit population sigmoid: S(v) = v_max / (1 + exp(r*(v0 - v)))."""
return v_max * _stable_expit(r * (v - v0))
def rhs(t: float, y: np.ndarray) -> np.ndarray:
x0, x1, x2, x3, x4, x5 = y
S_pyr = _jr_sigmoid(x1 - x2)
S_exc = _jr_sigmoid(C1 * x0)
S_inh = _jr_sigmoid(C3 * x0)
return np.array([
x3,
x4,
x5,
A * a * S_pyr - 2.0 * a * x3 - a2 * x0,
A * a * (p + C2 * S_exc) - 2.0 * a * x4 - a2 * x1,
B * b * C4 * S_inh - 2.0 * b * x5 - b2 * x2,
])- Parameters
- A = 3.25
- B = 22
- C1 = 135
- C2 = 108
- C3 = 33.75
- C4 = 33.75
- a = 100
- a2 = 10000
- b = 50
- b2 = 2500
- p = 220
- v_max = 5
- r = 0.56
- v0 = 6
- Initial condition
- y(0) = [1e-06, 1e-06, 1e-06, 1e-06, 1e-06, 1e-06]
- 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.
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_jansen_rit_normal_alpha_noise_0_0001_sub_clinical_2026,
title = {Resonix Evidence Portal: Jansen-Rit normal alpha — noise σ=0.0001 (sub-clinical)},
author = {{Resonix Labs (Canada) Inc.}},
year = {2026},
howpublished = {\url{https://resonix.tech/evidence/problems/jansen-rit-normal-alpha-noise-0-0001-sub-clinical}},
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