pub fn joint_psth( train_a: &[i32], train_b: &[i32], bin_size: usize, ) -> (Vec<f64>, usize)
Joint PSTH matrix (Aertsen et al. 1989). Returns flattened n×n outer product of mean-subtracted binned counts.