In order to get rid of the overflows due to factorials in wigner d-function (wigd), we have represented the factorial as a Gamma function which relates to the factorials as follows: G(n)=(n-1)! The Gamma function is defined in BSD library. The multiplication in factorials sections \'93pre1\'94 and \'93deform\'94 can be replaced by the summation if we take the logarithm of these values. Then, to get back to \'93pre2=pre1/denorm\'94 used in calculation of wigd, we have to take the exponent of log(pre1/denorm). \
0 commit comments