pub fn dot_product_f64(a: &[f64], b: &[f64]) -> f64
Compute dot product: sum(a[i] * b[i]) (f64 version)