估计两个伽马函数的比率,当各个函数 return inf 值

Estimate ratio of two gamma functions when individual functions return inf value

我正在估算两个伽马函数的比率。 gamma(x)gamma(y) 的两个估计值都非常大 ( > 10^300) 但两者的比率应该相当小。

from scipy.special import gamma
gamma(x) / gamma(y)

不幸的是,有一点 gamma(x) 太大,scipy returns 是 inf 值。有没有办法阻止这种情况、提高阈值或计算这个比率?

谢谢

假设 x 和 y 是实数且为正,您可以使用以下恒等式:

a / b = e^(ln(a) - ln(b))

我建议使用 gmpy2 进行任意精度计算。它有 gmpy2.lgamma,其中 returns 伽玛函数的对数,您可以使用 gmpy2.exp 将对数形式转换为您想要的比率。