I have noticed that the builtin gamma function is very accurate and extremely fast. Can someone tell me where to find the source code for the implementation?
gdb skips over the call to the builtin gamma. I assume it is not implemented by the simple Lanczos algorithm in tr1/gamma.tcc because I tried this code and it was about half the speed and less accurate. Thanks, James
