#include "longdbl.h"
+#include "libc.h"
+
union fshape {
float value;
uint32_t bits;
* ====================================================
*/
+#define _GNU_SOURCE
#include "libm.h"
static const float
-#define _GNU_SOURCE
#include "libm.h"
+double __lgamma_r(double, int *);
+
double lgamma(double x)
{
- return lgamma_r(x, &signgam);
+ return __lgamma_r(x, &signgam);
}
// FIXME
}
-double lgamma_r(double x, int *signgamp)
+double __lgamma_r(double x, int *signgamp)
{
double t,y,z,nadj,p,p1,p2,p3,q,r,w;
int32_t hx;
r = nadj - r;
return r;
}
+
+weak_alias(__lgamma_r, lgamma_r);
-#define _GNU_SOURCE
#include "libm.h"
+float __lgammaf_r(float, int *);
+
float lgammaf(float x)
{
- return lgamma_r(x, &signgam);
+ return __lgammaf_r(x, &signgam);
}
// FIXME
}
-float lgammaf_r(float x, int *signgamp)
+float __lgammaf_r(float x, int *signgamp)
{
float t,y,z,nadj,p,p1,p2,p3,q,r,w;
int32_t hx;
r = nadj - r;
return r;
}
+
+weak_alias(__lgammaf_r, lgammaf_r);
#define _GNU_SOURCE
#include "libm.h"
-long double lgammal(long double x)
-{
- return lgammal_r(x, &signgam);
-}
-
#if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
-long double lgammal_r(long double x, int *sg)
+long double __lgammal_r(long double x, int *sg)
{
- return lgamma_r(x, sg);
+ return __lgamma_r(x, sg);
}
#elif LDBL_MANT_DIG == 64 && LDBL_MAX_EXP == 16384
static const long double
return -y;
}
-long double lgammal_r(long double x, int *sg) {
+long double __lgammal_r(long double x, int *sg) {
long double t, y, z, nadj, p, p1, p2, q, r, w;
int i, ix;
uint32_t se, i0, i1;
return r;
}
#endif
+
+long double lgammal(long double x)
+{
+ return __lgammal_r(x, &signgam);
+}
+
+weak_alias(__lgammal_r, lgammal_r);