lgammal
y = expl(lgammal(x) / 2);
__weak_reference(lgamma, lgammal);
long double lgammal(long double);