From 40305f74bd70a575ce73260227ed3b64e0625b13 Mon Sep 17 00:00:00 2001
From: nsz <nsz@port70.net>
Date: Thu, 15 Mar 2012 09:29:53 +0100
Subject: [PATCH] in math.h make lgamma_r and non-double bessel _GNU_SOURCE
 only long double and float bessel functions are no longer xsi extensions

---
 include/math.h     | 38 ++++++++++++++++++--------------------
 src/math/lgamma.c  |  1 +
 src/math/lgammaf.c |  1 +
 src/math/lgammal.c |  1 +
 4 files changed, 21 insertions(+), 20 deletions(-)

diff --git a/include/math.h b/include/math.h
index b9139b02..9c06c967 100644
--- a/include/math.h
+++ b/include/math.h
@@ -349,33 +349,13 @@ double      gamma(double);
 float       gammaf(float);
 long double gammal(long double);
 
-double      lgamma_r(double, int*);
-float       lgammaf_r(float, int*);
-long double lgammal_r(long double, int*);
-
 double      j0(double);
-float       j0f(float);
-long double j0l(long double);
-
 double      j1(double);
-float       j1f(float);
-long double j1l(long double);
-
 double      jn(int, double);
-float       jnf(int, float);
-long double jnl(int, long double);
 
 double      y0(double);
-float       y0f(float);
-long double y0l(long double);
-
 double      y1(double);
-float       y1f(float);
-long double y1l(long double);
-
 double      yn(int, double);
-float       ynf(int, float);
-long double ynl(int, long double);
 #endif
 
 #ifdef _GNU_SOURCE
@@ -386,6 +366,24 @@ long double scalbl(long double, long double);
 void        sincos(double, double*, double*);
 void        sincosf(float, float*, float*);
 void        sincosl(long double, long double*, long double*);
+
+double      lgamma_r(double, int*);
+float       lgammaf_r(float, int*);
+long double lgammal_r(long double, int*);
+
+float       j0f(float);
+long double j0l(long double);
+float       j1f(float);
+long double j1l(long double);
+float       jnf(int, float);
+long double jnl(int, long double);
+
+float       y0f(float);
+long double y0l(long double);
+float       y1f(float);
+long double y1l(long double);
+float       ynf(int, float);
+long double ynl(int, long double);
 #endif
 
 #ifdef __cplusplus
diff --git a/src/math/lgamma.c b/src/math/lgamma.c
index d12462b9..9af7eee4 100644
--- a/src/math/lgamma.c
+++ b/src/math/lgamma.c
@@ -1,3 +1,4 @@
+#define _GNU_SOURCE
 #include "libm.h"
 
 double lgamma(double x)
diff --git a/src/math/lgammaf.c b/src/math/lgammaf.c
index f50f2379..aed98ba4 100644
--- a/src/math/lgammaf.c
+++ b/src/math/lgammaf.c
@@ -1,3 +1,4 @@
+#define _GNU_SOURCE
 #include "libm.h"
 
 float lgammaf(float x)
diff --git a/src/math/lgammal.c b/src/math/lgammal.c
index 603477c9..a33707ad 100644
--- a/src/math/lgammal.c
+++ b/src/math/lgammal.c
@@ -85,6 +85,7 @@
  *
  */
 
+#define _GNU_SOURCE
 #include "libm.h"
 
 long double lgammal(long double x)