The fma functions compute (x * y) + z, rounded as one ternary
    operation: the value is computed as if to infinite precision and
    rounded once to the result format, according to the rounding mode
    characterized by the value of FLT_ROUNDS.