Function py_cazelles_map_simulate 

pub(crate) fn py_cazelles_map_simulate<'py>(
    py: Python<'py>,
    x0: f64,
    y0: f64,
    a: f64,
    epsilon: f64,
    sigma: f64,
    x_threshold: f64,
    n_steps: usize,
    current: f64,
) -> (Bound<'py, PyArray1<f64>>, i64, f64, f64)

N-step Cazelles-Courbage-Rabinovich (2001) bursting-map simulation.

Parity contract with sc_neurocore.neurons.models.cazelles_map.CazellesMapNeuron.simulate: for the same parameters and constant input the returned x trace, spike count, and final (x, y) state are bit-identical to the Python reference (the map is exact floating-point arithmetic, no transcendental functions).