/* 1 1 log gamma function
*/
double nf_lgamma(const Node *n, const Cell *c)
{
  return lgamma( eval_tree( Right(n), c) );
}