Provide reliable, accurate down to several ULP's and fast numerical framework to compute generalized Fermi-Dirac integrals aiming at full floating-point parameter coverage, generality, precision control and speed.
For very small theta<10^-100 or very large theta>10^100 NaN's or huge errors produced. This is not serious problem in practice, but annoying, and spoils the test results.
For very small theta<10^-100 or very large theta>10^100 NaN's or huge errors produced. This is not serious problem in practice, but annoying, and spoils the test results.