Description

   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.