projects
/
libc-test
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
|
inline
| side by side
add regexec REG_NOSUB regression test by Hiltjo Posthuma
[libc-test]
/
src
/
math
/
j0.c
diff --git
a/src/math/j0.c
b/src/math/j0.c
index
ddf4064
..
f937ab4
100644
(file)
--- a/
src/math/j0.c
+++ b/
src/math/j0.c
@@
-1,14
+1,18
@@
+#define _BSD_SOURCE 1
+#define _XOPEN_SOURCE 700
#include <stdint.h>
#include <stdio.h>
#include <stdint.h>
#include <stdio.h>
-#include "
util
.h"
+#include "
mtest
.h"
static struct d_d t[] = {
#include "sanity/j0.h"
static struct d_d t[] = {
#include "sanity/j0.h"
+#include "special/j0.h"
};
int main(void)
{
};
int main(void)
{
+ #pragma STDC FENV_ACCESS ON
double y;
float d;
int e, i, err = 0;
double y;
float d;
int e, i, err = 0;
@@
-16,19
+20,28
@@
int main(void)
for (i = 0; i < sizeof t/sizeof *t; i++) {
p = t + i;
for (i = 0; i < sizeof t/sizeof *t; i++) {
p = t + i;
- setupfenv(p->r);
+
+ if (p->r < 0)
+ continue;
+ fesetround(p->r);
+ feclearexcept(FE_ALL_EXCEPT);
y = j0(p->x);
y = j0(p->x);
- e = getexcept();
+ e = fetestexcept(INEXACT|INVALID|DIVBYZERO|UNDERFLOW|OVERFLOW);
+
if (!checkexcept(e, p->e, p->r)) {
if (!checkexcept(e, p->e, p->r)) {
- printf("%s j0(%a)==%a except: want %s", rstr(p->r), p->x, p->y, estr(p->e));
+ printf("%s:%d: bad fp exception: %s j0(%a)=%a, want %s",
+ p->file, p->line, rstr(p->r), p->x, p->y, estr(p->e));
printf(" got %s\n", estr(e));
err++;
}
d = ulperr(y, p->y, p->dy);
if (!checkulp(d, p->r)) {
printf(" got %s\n", estr(e));
err++;
}
d = ulperr(y, p->y, p->dy);
if (!checkulp(d, p->r)) {
-
printf("
%s j0(%a) want %a got %a ulperr %.3f = %a + %a\n",
-
rstr(p->r), p->x, p->y, y, d, d-p->dy, p->dy);
+
// printf("%s:%d:
%s j0(%a) want %a got %a ulperr %.3f = %a + %a\n",
+
// p->file, p->line,
rstr(p->r), p->x, p->y, y, d, d-p->dy, p->dy);
err++;
err++;
+ // TODO: avoid spamming the output
+ printf(__FILE__ ": known to be broken near zeros\n");
+ break;
}
}
return !!err;
}
}
return !!err;