Concerning Gamma, the main reason I haven't implemented it yet is that

it requires determining some magic constants that were only given to me

for a lower precision than what I need (10 digits instead of 16).

For the Gamma function the book

„Numerical Recipes: The art of Scientific Computing“ (3rd edition);

William H. Press, Saul A. Teukolsky, William T. Vetterling, Brian P. Flannery;

Cambridge University Press

might be of interest as it describes the implementation of a Gamma function with error smaller than < 1e-15.

Ben