test/math/signgam.c

   1 #define _XOPEN_SOURCE 600
   2 #include <math.h>
   3 #include <float.h>
   4 #include <stdio.h>
   5
   6 double zero = 0.0;
   7 double mzero;
   8
   9 int main(void)
  10 {
  11         double d;
  12         int errors = 0;
  13
  14         mzero = copysign(zero, -1.0);
  15
  16         d = lgamma(zero);
  17         printf("%g %d\n", d, signgam);
  18         errors += !(d == HUGE_VAL);
  19         errors += !(signgam == 1);
  20
  21         d = lgamma(mzero);
  22         printf("%g %d\n", d, signgam);
  23         errors += !(d == HUGE_VAL);
  24         errors += !(signgam == -1);
  25
  26         printf("Errors: %d\n", errors);
  27         return errors;
  28 }