double lgamma_r(double x, List<int> signp) { double result = lgamma(x); signp[0] = signgam; return result; }