Function std::num::mul_add
pub fn mul_add<T: Float>(a: T, b: T, c: T) -> T
Fused multiply-add. Computes (a * b) + c
with only one rounding error.
This produces a more accurate result with better performance (on some architectures) than a separate multiplication operation followed by an add.