More data concerning the precision of bigfloat calculation of gamma
Subject: More data concerning the precision of bigfloat calculation of gamma
From: Richard Fateman
Date: Sat, 11 Oct 2008 09:23:50 -0700
If you feel that the accuracy is not full precision, you know that you can
increase fpprec during the calculation to get more correct digits.
This is what (for example) Mathematica does.
If time is not important, you can compute the result in precision n+5 and
n+10. If they agree to n digits,
then they are (very probably) correct digits.
Testing at random inputs is usually not very effective in determining the
"worst case" errors.
Knowing what the algorithm does, and testing at critical points, helps.
RJF