X-Git-Url: http://nsz.repo.hu/git/?a=blobdiff_plain;f=src%2Fmath%2Flgammal_r.c;h=848f48d464954303a24923ca6db72fbd3a968a24;hb=32fca61b208425e16262e4d007dabad44bcd872b;hp=2ebe612aa22f41cf4745b716875aab13e313fe55;hpb=d22f4f43182fe1cc70b343b6a52cb8e3a54d1f44;p=libc-test

diff --git a/src/math/lgammal_r.c b/src/math/lgammal_r.c
index 2ebe612..848f48d 100644
--- a/src/math/lgammal_r.c
+++ b/src/math/lgammal_r.c
@@ -1,13 +1,16 @@
+#define _GNU_SOURCE 1
 #include <stdint.h>
 #include <stdio.h>
-#include "util.h"
+#include "mtest.h"
 
 static struct l_li t[] = {
 #if LDBL_MANT_DIG == 53
-#include "sanity/lgammal_r.h"
+#include "sanity/lgamma_r.h"
+#include "special/lgamma_r.h"
 
 #elif LDBL_MANT_DIG == 64
 #include "sanity/lgammal_r.h"
+#include "special/lgammal_r.h"
 
 #endif
 };
@@ -38,7 +41,8 @@ int main(void)
 			err++;
 		}
 		d = ulperrl(y, p->y, p->dy);
-		if (!checkulp(d, p->r) || yi != p->i) {
+		// TODO: 2 ulp errors allowed
+		if ((p->r==RN && fabs(d)>2) || (!isnan(p->x) && p->x!=-inf && !(p->e&DIVBYZERO) && yi != p->i)) {
 			printf("%s:%d: %s lgammal_r(%La) want %La,%lld got %La,%d ulperr %.3f = %a + %a\n",
 				p->file, p->line, rstr(p->r), p->x, p->y, p->i, y, yi, d, d-p->dy, p->dy);
 			err++;