Function py_medvedev_map_simulate 

pub(crate) fn py_medvedev_map_simulate<'py>(
    py: Python<'py>,
    x0: f64,
    alpha: f64,
    beta: f64,
    x_threshold: f64,
    n_steps: usize,
    current: f64,
) -> (Bound<'py, PyArray1<f64>>, i64, f64)

Parity contract with sc_neurocore.neurons.models.medvedev_map.MedvedevMapNeuron.simulate: for the same parameters and constant input the returned x trace, upward- crossing spike count, and final x state are bit-identical to the Python reference (the map is exact floating-point arithmetic — a multiply, an add, and a fold into [0, 1); f64::rem_euclid(1.0) equals Python’s x % 1.0 bit-for-bit). This is a one-dimensional map, so there is no y state.