math: add an expm1l special case

[libc-test] / src / math / fmal.c

diff --git a/src/math/fmal.c b/src/math/fmal.c
index 801a7a3..791d5d3 100644 (file)
--- a/src/math/fmal.c
+++ b/src/math/fmal.c
@@ -1,12 +1,14 @@
 #include <stdint.h>
 #include <stdio.h>
-#include "util.h"
+#include "mtest.h"
 
 static struct lll_l t[] = {
 #if LDBL_MANT_DIG == 53
 #include "sanity/fma.h"
+#include "special/fma.h"
 #elif LDBL_MANT_DIG == 64
 #include "sanity/fmal.h"
+#include "special/fmal.h"
 #endif
 };
 
@@ -28,14 +30,19 @@ int main(void)
                y = fmal(p->x, p->x2, p->x3);
                e = fetestexcept(INEXACT|INVALID|DIVBYZERO|UNDERFLOW|OVERFLOW);
 
-               if (!checkexcept(e, p->e, p->r)) {
+               /* do not check inexact by default */
+#if defined CHECK_INEXACT || defined CHECK_INEXACT_OMISSION
+               if (!checkexceptall(e, p->e, p->r)) {
+#else
+               if (!checkexceptall(e|INEXACT, p->e|INEXACT, p->r)) {
+#endif
                        printf("%s:%d: bad fp exception: %s fmal(%La,%La,%La)=%La, want %s",
                                p->file, p->line, rstr(p->r), p->x, p->x2, p->x3, p->y, estr(p->e));
                        printf(" got %s\n", estr(e));
                        err++;
                }
                d = ulperrl(y, p->y, p->dy);
-               if (!checkulp(d, p->r)) {
+               if (!checkcr(y, p->y, p->r)) {
                        printf("%s:%d: %s fmal(%La,%La,%La) want %La got %La ulperr %.3f = %a + %a\n",
                                p->file, p->line, rstr(p->r), p->x, p->x2, p->x3, p->y, y, d, d-p->dy, p->dy);
                        err++;