python - 個々の関数が inf 値を返す場合の 2 つのガンマ関数の比率の推定

Question

2 つのガンマ関数の比を推定しています。gamma(x)との両方の推定値gamma(y)は非常に大きい (> 10^300) が、2 つの比率はかなり小さいはずです。

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

残念ながら、 がgamma(x)大きすぎて scipy がinf値を返すポイントがあります。これを停止する方法、しきい値を上げる方法、または代わりにこの比率を計算する方法はありますか?

ありがとう

Answer 1

x と y が実数で正であると仮定すると、次の恒等式を使用できます。

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

gmpy2任意精度の計算に使用することをお勧めします。gmpy2.lgammaガンマ関数の対数を返す があり、対数形式から目的の比率に変換するために使用できますgmpy2.exp。