big float gamma



In bffac.mac, there is a big float gamma function bffac (author Bill
Gosper).
I was testing this function with complex inputs -- the output isn't in
rectangular form (it's a big mess). So I was wondering: Is bffac
algorithmically
correct for (all) complex inputs? If so, maybe we could insert some calls
to rectform
into the code? Example:

(%i49) float2bf : true$
(%i50) load(bffac)$

(%i51) bffac(-5.0 + %i,25);

< long mess deleted >

(%i52) rectform(%);
(%o52) 1.262948529867786b-2-7.163743828758051b-4*%i

BW