FMA may return *more precise* results since there is no intermediate rounding when computing (a + b * c). (b * c) is evaluate with infinite precision and only the final result after adding a is rounded to float32. This allows for a FLT_EPSILON difference compared to rounding (b * c) first.