/*	$OpenBSD: lgamma.c,v 1.3 2016/10/23 18:46:03 otto Exp $	*/

/*
 * Written by Martynas Venckus.  Public domain
 */

#define __POSIX_VISIBLE 201403

#include <assert.h>
#include <math.h>
#include <float.h>


int
main(int argc, char *argv[])
{
	assert(isnan(lgamma(NAN)));
	assert(isnan(lgammaf(NAN)));

	signgam = 0;
	assert(lgamma(-HUGE_VAL) == HUGE_VAL && signgam == 1);
	signgam = 0;
	assert(lgammaf(-HUGE_VALF) == HUGE_VALF && signgam == 1);
	signgam = 0;
	assert(lgammal(-HUGE_VALL) == HUGE_VALL && signgam == 1);

	signgam = 0;
	assert(lgamma(HUGE_VAL) == HUGE_VAL && signgam == 1);
	signgam = 0;
	assert(lgammaf(HUGE_VALF) == HUGE_VALF && signgam == 1);
	signgam = 0;
	assert(lgammal(HUGE_VALL) == HUGE_VALL && signgam == 1);

	signgam = 0;
	assert(lgamma(-0.0) == HUGE_VAL && signgam == -1);
	signgam = 0;
	assert(lgammaf(-0.0F) == HUGE_VALF && signgam == -1);
	signgam = 0;
	assert(lgammal(-0.0L) == HUGE_VALL && signgam == -1);

	signgam = 0;
	assert(lgamma(0.0) == HUGE_VAL && signgam == 1);
	signgam = 0;
	assert(lgammaf(0.0F) == HUGE_VALF && signgam == 1);
	signgam = 0;
	assert(lgammal(0.0L) == HUGE_VALL && signgam == 1);

	signgam = 0;
	assert(lgamma(1.0) == 0.0 && signgam == 1);
	signgam = 0;
	assert(lgammaf(1.0F) == 0.0F && signgam == 1);
	signgam = 0;
	assert(lgammal(1.0L) == 0.0L && signgam == 1);

	signgam = 0;
	assert(fabs(lgamma(3.0) - M_LN2) < DBL_EPSILON && signgam == 1);
	signgam = 0;
	assert(fabsf(lgammaf(3.0F) - (float)M_LN2) < FLT_EPSILON && signgam == 1);
	signgam = 0;
	assert(fabsl(lgammal(3.0L) - M_LN2l) < LDBL_EPSILON && signgam == 1);

	return (0);
}