TEST(math, lgammaf_r_17471883) {
int sign;
sign = 0;
ASSERT_FLOAT_EQ(HUGE_VALF, lgammaf_r(0.0f, &sign));
ASSERT_EQ(1, sign);
sign = 0;
ASSERT_FLOAT_EQ(HUGE_VALF, lgammaf_r(-0.0f, &sign));
ASSERT_EQ(-1, sign);
}
int main(void)
{
#pragma STDC FENV_ACCESS ON
int yi;
double y;
float d;
int e, i, err = 0;
struct f_fi *p;
for (i = 0; i < sizeof t/sizeof *t; i++) {
p = t + i;
if (p->r < 0)
continue;
fesetround(p->r);
feclearexcept(FE_ALL_EXCEPT);
y = lgammaf_r(p->x, &yi);
e = fetestexcept(INEXACT|INVALID|DIVBYZERO|UNDERFLOW|OVERFLOW);
if (!checkexcept(e, p->e, p->r)) {
printf("%s:%d: bad fp exception: %s lgammaf_r(%a)=%a,%lld, want %s",
p->file, p->line, rstr(p->r), p->x, p->y, p->i, estr(p->e));
printf(" got %s\n", estr(e));
err++;
}
d = ulperrf(y, p->y, p->dy);
if (!checkulp(d, p->r) || (!isnan(p->x) && p->x!=-inf && !(p->e&DIVBYZERO) && yi != p->i)) {
printf("%s:%d: %s lgammaf_r(%a) want %a,%lld got %a,%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++;
}
}
return !!err;
}
static TACommandVerdict lgammaf_r_cmd(TAThread thread,TAInputStream stream)
{
float x, res;
int signp;
x = readFloat(&stream);
START_TARGET_OPERATION(thread);
errno = 0;
res = lgammaf_r(x, &signp);
END_TARGET_OPERATION(thread);
writeInt(thread, errno);
writeFloat(thread, res);
writeInt(thread, signp);
sendResponse(thread);
return taDefaultVerdict;
}
float
lgammaf(float x)
{
return lgammaf_r(x,&signgam);
}
extern float _swift_Darwin_lgammaf_r(float x, int *psigngam) {
return lgammaf_r(x, psigngam);
}
TEST(math, lgammaf_r) {
int sign;
ASSERT_FLOAT_EQ(logf(24.0f), lgammaf_r(5.0f, &sign));
ASSERT_EQ(1, sign);
}
float
gammaf_r(float x, int *signgamp) /* wrapper lgammaf_r */
{
return lgammaf_r(x,signgamp);
}
