code: def lgamma(x, context=None):
      
      return _apply_function_in_current_context(
          BigFloat,
          lambda rop, op, rnd: mpfr.mpfr_lgamma(rop, op, rnd)[0],
          (BigFloat._implicit_convert(x),),
          context,
      )
docs: Return the logarithm of the absolute value of the Gamma function at x.