lgamma_r
double lgamma_r(double, int *);
PROTO_NORMAL(lgamma_r);
DEF_NONSTD(lgamma_r);
return lgamma_r(x,&signgam);
return lgamma_r(x,signgamp);
return lgamma_r(x,&signgam);