From 06696f40afe58a231e2531c8bf6d9f0dadb92e51 Mon Sep 17 00:00:00 2001 From: Jakob Kaivo Date: Sun, 3 Mar 2019 21:25:50 -0500 Subject: outline details from C18 annex F --- src/math/fma.c | 24 +++++++++++++++++++++++- 1 file changed, 23 insertions(+), 1 deletion(-) (limited to 'src/math/fma.c') diff --git a/src/math/fma.c b/src/math/fma.c index 18b059d1..845afea4 100644 --- a/src/math/fma.c +++ b/src/math/fma.c @@ -1,10 +1,32 @@ # define TGSOURCE "fma.c" #include "_tgmath.h" #include +#include "fenv.h" TYPE TGFN(fma)(TYPE x, TYPE y, TYPE z) { - return x - y - z; + int classx = fpclassify(x); + int classy = fpclassify(y); + /* int classz = fpclassify(z); */ + + if (classx == FP_INFINITE && classy == FP_ZERO) { + feraiseexcept(FE_INVALID); + return NAN; + } + + if (classx == FP_ZERO && classy == FP_INFINITE) { + feraiseexcept(FE_INVALID); + return NAN; + } + + /* + if (x * y == -z && fpclassify(x*y) == classz == FP_INFINITE) { + feraiseexcept(FE_INVALID); + return NAN; + } + */ + + return x * y + z; } /* -- cgit v1.2.1