search

VA00 / Generalized-Fermi-Dirac-integrals

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.
MIT License
0 stars 0 forks source link

Theta derivatives for theta=0 #3

Open VA00 opened 2 years ago

VA00 commented 2 years ago

There is some NaN's and discrepancies for theta=0 case. Either Ffermi_derivatives_m_n_quad or reference values from Mma are wrong. Or both...

VA00 commented 2 years ago

Analytical results exist for all derivatives id f=0 and theta=0. This shed some light on the problem. E.g: for k= 0.0 eta= 5.12e+02 theta=0.00e+00 second and third eta derivatives should be equal to 4.4e-223, while code returns -0.0.