-----Urspr?ngliche Nachricht-----
Von: raymond.toy at ericsson.com [mailto:raymond.toy at ericsson.com]
> As best as I can tell, you've done a very straightforward translation of
> the original Maxima implementation. Does this really make much
> difference? I would think we could make it quite a bit faster if we
> used the "native" bfloat operations instead of using mul/add which call
> the simplifier. (But maybe the major cost is calling $bern?. I didn't
> measure it.)
> Also, do you know the algorithm used? It looks like the asymptotic
> series for log(gamma(z)) (see A&S 6.1.40).
> Can't we arrange it to autoload bffac.mac when needed?
You are right. It is a straightforward translation. The alternative is to call
cbffac. The point is that the Gamma function must support the Complex Bigfloat
to allow further work on other special functions.
By the way, I have studied the Complex Bigfloat routines suggested by David
Billinghurst in a post he sent on January 2006 to the mailing list. Perhaps this
routines could be a starting point for more "native" bigfloat implementation of
Complex Bigfloat algorithm.
Dieter Kaiser