pub fn expansion_product(a: &[f64], b: &[f64]) -> Vec<f64>
Multiply two expansions exactly.
Each component of b scales a exactly; the partial expansions are then accumulated without rounding away any component.
b
a