fn log_linear_fit(signal: &[f64], echo_times: &[f64]) -> (f64, f64)
Log-linear R2* fit: log(S) = log(S0) - R2* * TE Solves via normal equations.