fn fit_linear_per_dim(trajectory: &[(i64, &[f32])]) -> (Vec<f64>, Vec<f64>)
Fit a linear model y_d = slope_d * t + intercept_d for each dimension.