Соотношение оценки двух гамма-функций, когда отдельные функции возвращают значение inf

Question

Я оцениваю соотношение двух гамма-функций. Обе оценки для gamma(x) и gamma(y) очень велики (> 10^300), но соотношение двух должно быть довольно небольшим.

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

К сожалению, есть точка, где gamma(x) слишком велико, а SciPy возвращает inf значение. Есть ли способ остановить это, увеличить порог или, наоборот, рассчитать это соотношение?

Благодаря

Answer 1

Предполагая, что х и у вещественны и положительны вы можете использовать следующее тождество:

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

Я могу предложить использовать gmpy2 для вычисления произвольной точности. Он имеет gmpy2.lgamma, который возвращает логарифм функции гамма, и вы можете использовать gmpy2.exp для преобразования из формы логарифма в желаемое соотношение.