Logarithm of the factorial function.
Avoids complications with factorials of large numbers overflowing 64-bit numbers.
n - non-negative integer number
lf - logarithm of the factorial of n
Double precision overflow is a persistent problem with Clebsch-Gordan coefficients and other objects that involve factorials.
Version 2.1, authors: Ilya Kuprov