aboutsummaryrefslogtreecommitdiff
path: root/sysdeps/x86/fpu/s_ffma.c
blob: bcd9edc74f91f5940645e3292c2bf932769d4f8c (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
/* Fused multiply-add of double value, narrowing the result to float.
   x86 version.
   Copyright (C) 2021-2022 Free Software Foundation, Inc.
   This file is part of the GNU C Library.

   The GNU C Library is free software; you can redistribute it and/or
   modify it under the terms of the GNU Lesser General Public
   License as published by the Free Software Foundation; either
   version 2.1 of the License, or (at your option) any later version.

   The GNU C Library is distributed in the hope that it will be useful,
   but WITHOUT ANY WARRANTY; without even the implied warranty of
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
   Lesser General Public License for more details.

   You should have received a copy of the GNU Lesser General Public
   License along with the GNU C Library; if not, see
   <https://www.gnu.org/licenses/>.  */

#define f32fmaf64 __hide_f32fmaf64
#define f32fmaf32x __hide_f32fmaf32x
#define ffmal __hide_ffmal
#include <math.h>
#undef f32fmaf64
#undef f32fmaf32x
#undef ffmal

#include <math-narrow.h>

#ifndef __FP_FAST_FMA
/* Depending on the details of the glibc configuration, fma might use
   either SSE or 387 arithmetic; ensure that both parts of the
   floating-point state are handled in the round-to-odd code.  If
   __FP_FAST_FMA is defined, that implies that the compiler is using
   SSE floating point and that the fma call will be inlined, so the
   x86 macros will work with only the SSE state and that is
   sufficient.  */
# undef libc_feholdexcept_setround
# define libc_feholdexcept_setround	default_libc_feholdexcept_setround
# undef libc_feupdateenv_test
# define libc_feupdateenv_test		default_libc_feupdateenv_test
#endif

float
__ffma (double x, double y, double z)
{
  NARROW_FMA_ROUND_TO_ODD (x, y, z, float, union ieee754_double, , mantissa1,
			   false);
}
libm_alias_float_double (fma)