numerical evaluation of quotients of gamma functions

>>>>> "Barton" == Barton Willis <willisb at unk.edu> writes:

    Barton> I need to numerically evaluate quotients of products of ?gamma functions. ?In order of?preference,
    Barton> I would like to?

    Barton> (1) avoid overflows, (2) have good accuracy, ?and (3) be fast (much less important). Any
    Barton> suggestions?

Use base 2 logs?  This might be a little more accurate then base e
logs, but not any slower.  Base 2 since floats are base 2, so log(f*2^n)
can be represented in 2 parts very accurately: n + log(f).  Base e
logs don't have this simple separation, and for large values you lose
many bits to hold the integer part.

Ray