A problem with tgamma function

[Lev Bishop]

On 10/12/07, Angelo Graziosi wrote:
>
> Lev Bishop wrote:
>
> > I'm sure glibc and newlib would both appreciate a good algorithm for
> > tgamma(), if you felt like submitting one...
>
> It seems that a good algorithm has yet been coded in the GAMMA
> implementation of CERNLIB
> (http://wwwasdoc.web.cern.ch/wwwasdoc/shortwrupsdir/c302/top.html,
> http://cernlib.web.cern.ch/cernlib/download/2006_source/src/mathlib/gen/c/gamma64.F),
> because with that implementation the test I posted passes...

One immediately obvious problem with that implementation is that it
doesn't handle negative parameter.

Lev

--
Unsubscribe info:      http://cygwin.com/ml/#unsubscribe-simple
Problem reports:       http://cygwin.com/problems.html
Documentation:         http://cygwin.com/docs.html
FAQ:                   http://cygwin.com/faq/