pub fn mad<T: Float, V: FloatOps<T>>(m: V, a: V, b: V) -> V
returns fused multiply add m * a + b
m * a + b