pub fn poly_mul(a: &[f64], b: &[f64]) -> Vec<f64>
Multiply two polynomials (convolution of coefficient arrays).