Calculating log gamma

Hi

I am trying to calculate the log of the gamma function.
I am using an algorithm that involves a quite big table
of Bernoulli numbers (71), and I can't get any better precision
than 128-129 decimal digits.

I am using a heavily modified cephes math library. I have
incresed precision to 132 bits, but I can't get to that precision
only with the lgamma function. Other functions (sin/cos/log
and many others) achieve 132 correct digits.

Where are good references for this problem? (Papers, programs,
and similar documentation)

I would like to try another approach.

Thanks in advance.


--
jacob navia
jacob at jacob point remcomp point fr
logiciels/informatique
http://www.cs.virginia.edu/~lcc-win32
.

Relevant Pages

  • Re: Calculating log gamma
    ... of Bernoulli numbers, and I can't get any better precision ... -- John Harper, School of Mathematics, Statistics and Computer Science, Victoria University, PO Box 600, Wellington 6140, New Zealand ... I obtain 132 digits precision now. ... jacob navia ...
    (sci.math.num-analysis)
  • Re: Conformance
    ... They give a precision of approx 30 decimal digits ... is called AT COMPILE TIME when the compiler finds an unknown number ... jacob at jacob point remcomp point fr ...
    (comp.std.c)
  • Re: Raising a Number doubt..?
    ... I've to do xtremely large calculations. ... What level of precision to you require? ... The number of decimal digits that can be stored in 1048576000 bits is ... jacob at jacob point remcomp point fr ...
    (comp.lang.c)
  • Re: about c to fortran
    ... > more than 16 significant figures of precision. ... Where d is the number of decimal digits you write and ... no deficiencies and the other way is to make it so complicated ...
    (comp.lang.fortran)
  • Re: Is it time to legitimise REAL*8 etc?
    ... want/need to specify things. ... Decimal digits of precision might seem fine ... but that just isn't a good match the reality of programming. ...
    (comp.lang.fortran)