1#define _GNU_SOURCE 1 2#include <stdint.h> 3#include <stdio.h> 4#include "mtest.h" 5 6static struct l_l t[] = { 7#if LDBL_MANT_DIG == 53 8#include "sanity/exp10.h" 9#include "special/exp10.h" 10 11#elif LDBL_MANT_DIG == 64 12#include "sanity/exp10l.h" 13#include "special/exp10l.h" 14 15#endif 16}; 17 18int main(void) 19{ 20 #pragma STDC FENV_ACCESS ON 21 long double y; 22 float d; 23 int e, i, err = 0; 24 struct l_l *p; 25 26 for (i = 0; i < sizeof t/sizeof *t; i++) { 27 p = t + i; 28 29 if (p->r < 0) 30 continue; 31 fesetround(p->r); 32 feclearexcept(FE_ALL_EXCEPT); 33 y = exp10l(p->x); 34 e = fetestexcept(INEXACT|INVALID|DIVBYZERO|UNDERFLOW|OVERFLOW); 35 36 if (!checkexcept(e, p->e, p->r)) { 37 printf("%s:%d: bad fp exception: %s exp10l(%La)=%La, want %s", 38 p->file, p->line, rstr(p->r), p->x, p->y, estr(p->e)); 39 printf(" got %s\n", estr(e)); 40 err++; 41 } 42 d = ulperrl(y, p->y, p->dy); 43 if (!checkulp(d, p->r)) { 44 printf("%s:%d: %s exp10l(%La) want %La got %La ulperr %.3f = %a + %a\n", 45 p->file, p->line, rstr(p->r), p->x, p->y, y, d, d-p->dy, p->dy); 46 err++; 47 } 48 } 49 return !!err; 50} 51