The lgamma functions return the logarithm of the absolute value
    of gamma of x, or ln(|G(x)|), where G is the gamma function.

    The sign of gamma of x is returned in the external integer
    variable signgam. The x argument cannot be 0, a negative integer,
    or Infinity.