Updated 2022-05-21 0
Viewed 6 times

I am trying to compute the gamma function of an array of complex numbers. The problem is: Gamma is returning NAN for certain elements. I would like to know the reason for that. Also, there is any other alternative library in Python providing a more stable gamma function that could perform the same calculation without returning NAN?

special.gamma(b[dim_l][j][0] + b[dim_l][j][1] * s[:, dim_l])]) 
🔴 No definitive solution yet