check_int ("gamma (0) sets signgam to 1", signgam, 1);
check_int ("gamma (1) sets signgam to 1", signgam, 1);
gamma_test (void){  int save_lib_version = _LIB_VERSION;  errno = 0;  FUNC(gamma) (1);  if (errno == ENOSYS)    /* Function not implemented.  */    return;  feclearexcept (FE_ALL_EXCEPT);  _LIB_VERSION = _SVID_;  check_isinfp ("gamma (+inf) == +inf", FUNC(gamma) (plus_infty));  check_isinfp_exc ("gamma (0) == +inf plus divide by zero exception",		    FUNC(gamma) (0), DIVIDE_BY_ZERO_EXCEPTION);  check_isinfp_exc ("gamma (x) == +inf plus divide by zero exception for integer x <= 0",		    FUNC(gamma) (-3), DIVIDE_BY_ZERO_EXCEPTION);  check_isnan_exc ("gamma (-inf) == NaN plus invalid exception",                   FUNC(gamma) (minus_infty), INVALID_EXCEPTION);  signgam = 0;  check ("gamma (1) == 0", FUNC(gamma) (1), 0);  check_int ("gamma (0) sets signgam to 1", signgam, 1);  signgam = 0;  check ("gamma (3) == M_LN2", FUNC(gamma) (3), M_LN2l);  check_int ("gamma (3) sets signgam to 1", signgam, 1);  signgam = 0;  check_eps ("gamma (0.5) == log(sqrt(pi))", FUNC(gamma) (0.5),             FUNC(log) (FUNC(sqrt) (M_PIl)), CHOOSE (0, 1e-15, 1e-7));  check_int ("gamma (0.5) sets signgam to 1", signgam, 1);  signgam = 0;  check_eps ("gamma (-0.5) == log(2*sqrt(pi))", FUNC(gamma) (-0.5),             FUNC(log) (2*FUNC(sqrt) (M_PIl)), CHOOSE (0, 1e-15, 0));  check_int ("gamma (-0.5) sets signgam to -1", signgam, -1);  _LIB_VERSION = _IEEE_;  check_isinfp ("gamma (+inf) == +inf", FUNC(gamma) (plus_infty));  check_isnan_exc ("gamma (0) == NaN plus invalid exception",                    FUNC(gamma) (0), INVALID_EXCEPTION);  check_isnan_exc_ext ("gamma (x) == NaN plus invalid exception for integer x <= 0",                        FUNC(gamma) (-2), INVALID_EXCEPTION, -2);  check_isnan_exc ("gamma (-inf) == NaN plus invalid exception",                   FUNC(gamma) (minus_infty), INVALID_EXCEPTION);#ifdef TODO  check_eps ("gamma (0.5) == sqrt(pi)", FUNC(gamma) (0.5), FUNC(sqrt) (M_PIl),	     CHOOSE (0, 5e-16, 2e-7));#endif  check_eps ("gamma (-0.5) == -2*sqrt(pi)", FUNC(gamma) (-0.5),	     -2*FUNC(sqrt) (M_PIl), CHOOSE (0, 5e-16, 3e-7));  check ("gamma (1) == 1", FUNC(gamma) (1), 1);  check ("gamma (4) == 6", FUNC(gamma) (4), 6);  check_eps ("gamma (0.7) == 1.29805...", FUNC(gamma) (0.7),	     1.29805533264755778568L, CHOOSE(0, 3e-16, 2e-7));  check ("gamma (1.2) == 0.91816...", FUNC(gamma) (1.2), 0.91816874239976061064L);  _LIB_VERSION = save_lib_version;}