pub fn gpfa( trains: &[&[i32]], n_latents: usize, bin_ms: f64, dt: f64, max_iter: usize, tol: f64, seed: u64, ) -> GpfaResult
Extract smooth latent trajectories from parallel spike trains via EM.