Function py_izhikevich2007_simulate 

pub(crate) fn py_izhikevich2007_simulate<'py>(
    py: Python<'py>,
    v0: f64,
    u0: f64,
    cap: f64,
    k: f64,
    vr: f64,
    vt: f64,
    vpeak: f64,
    a: f64,
    b: f64,
    c: f64,
    d: f64,
    dt: f64,
    n_steps: usize,
    current: f64,
) -> (Bound<'py, PyArray1<f64>>, i64, f64, f64)

Parity contract with sc_neurocore.neurons.models.izhikevich2007.Izhikevich2007Neuron.simulate: for the same parameters and constant input the returned v trace, spike count, and final (v, u) state are bit-identical to the Python RK4 reference (the NeuroML right-hand side k (v-vr)(v-vt)/C is exact arithmetic — products, a sum and a division, no transcendental functions).