lgamma_r
double lgamma_r(double, int *);
extern double lgamma_r(double, int *);
return (lgamma_r(x, signgamp));