|
Packit |
67cb25 |
/* specfunc/gsl_sf_fermi_dirac.h
|
|
Packit |
67cb25 |
*
|
|
Packit |
67cb25 |
* Copyright (C) 1996, 1997, 1998, 1999, 2000 Gerard Jungman
|
|
Packit |
67cb25 |
*
|
|
Packit |
67cb25 |
* This program is free software; you can redistribute it and/or modify
|
|
Packit |
67cb25 |
* it under the terms of the GNU General Public License as published by
|
|
Packit |
67cb25 |
* the Free Software Foundation; either version 3 of the License, or (at
|
|
Packit |
67cb25 |
* your option) any later version.
|
|
Packit |
67cb25 |
*
|
|
Packit |
67cb25 |
* This program is distributed in the hope that it will be useful, but
|
|
Packit |
67cb25 |
* WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
Packit |
67cb25 |
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
|
|
Packit |
67cb25 |
* General Public License for more details.
|
|
Packit |
67cb25 |
*
|
|
Packit |
67cb25 |
* You should have received a copy of the GNU General Public License
|
|
Packit |
67cb25 |
* along with this program; if not, write to the Free Software
|
|
Packit |
67cb25 |
* Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
|
|
Packit |
67cb25 |
*/
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
/* Author: G. Jungman */
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
#ifndef __GSL_SF_FERMI_DIRAC_H__
|
|
Packit |
67cb25 |
#define __GSL_SF_FERMI_DIRAC_H__
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
#include <gsl/gsl_sf_result.h>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
#undef __BEGIN_DECLS
|
|
Packit |
67cb25 |
#undef __END_DECLS
|
|
Packit |
67cb25 |
#ifdef __cplusplus
|
|
Packit |
67cb25 |
# define __BEGIN_DECLS extern "C" {
|
|
Packit |
67cb25 |
# define __END_DECLS }
|
|
Packit |
67cb25 |
#else
|
|
Packit |
67cb25 |
# define __BEGIN_DECLS /* empty */
|
|
Packit |
67cb25 |
# define __END_DECLS /* empty */
|
|
Packit |
67cb25 |
#endif
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
__BEGIN_DECLS
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
/* Complete Fermi-Dirac Integrals:
|
|
Packit |
67cb25 |
*
|
|
Packit |
67cb25 |
* F_j(x) := 1/Gamma[j+1] Integral[ t^j /(Exp[t-x] + 1), {t,0,Infinity}]
|
|
Packit |
67cb25 |
*
|
|
Packit |
67cb25 |
*
|
|
Packit |
67cb25 |
* Incomplete Fermi-Dirac Integrals:
|
|
Packit |
67cb25 |
*
|
|
Packit |
67cb25 |
* F_j(x,b) := 1/Gamma[j+1] Integral[ t^j /(Exp[t-x] + 1), {t,b,Infinity}]
|
|
Packit |
67cb25 |
*/
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
/* Complete integral F_{-1}(x) = e^x / (1 + e^x)
|
|
Packit |
67cb25 |
*
|
|
Packit |
67cb25 |
* exceptions: GSL_EUNDRFLW
|
|
Packit |
67cb25 |
*/
|
|
Packit |
67cb25 |
int gsl_sf_fermi_dirac_m1_e(const double x, gsl_sf_result * result);
|
|
Packit |
67cb25 |
double gsl_sf_fermi_dirac_m1(const double x);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
/* Complete integral F_0(x) = ln(1 + e^x)
|
|
Packit |
67cb25 |
*
|
|
Packit |
67cb25 |
* exceptions: GSL_EUNDRFLW
|
|
Packit |
67cb25 |
*/
|
|
Packit |
67cb25 |
int gsl_sf_fermi_dirac_0_e(const double x, gsl_sf_result * result);
|
|
Packit |
67cb25 |
double gsl_sf_fermi_dirac_0(const double x);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
/* Complete integral F_1(x)
|
|
Packit |
67cb25 |
*
|
|
Packit |
67cb25 |
* exceptions: GSL_EUNDRFLW, GSL_EOVRFLW
|
|
Packit |
67cb25 |
*/
|
|
Packit |
67cb25 |
int gsl_sf_fermi_dirac_1_e(const double x, gsl_sf_result * result);
|
|
Packit |
67cb25 |
double gsl_sf_fermi_dirac_1(const double x);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
/* Complete integral F_2(x)
|
|
Packit |
67cb25 |
*
|
|
Packit |
67cb25 |
* exceptions: GSL_EUNDRFLW, GSL_EOVRFLW
|
|
Packit |
67cb25 |
*/
|
|
Packit |
67cb25 |
int gsl_sf_fermi_dirac_2_e(const double x, gsl_sf_result * result);
|
|
Packit |
67cb25 |
double gsl_sf_fermi_dirac_2(const double x);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
/* Complete integral F_j(x)
|
|
Packit |
67cb25 |
* for integer j
|
|
Packit |
67cb25 |
*
|
|
Packit |
67cb25 |
* exceptions: GSL_EUNDRFLW, GSL_EOVRFLW
|
|
Packit |
67cb25 |
*/
|
|
Packit |
67cb25 |
int gsl_sf_fermi_dirac_int_e(const int j, const double x, gsl_sf_result * result);
|
|
Packit |
67cb25 |
double gsl_sf_fermi_dirac_int(const int j, const double x);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
/* Complete integral F_{-1/2}(x)
|
|
Packit |
67cb25 |
*
|
|
Packit |
67cb25 |
* exceptions: GSL_EUNDRFLW, GSL_EOVRFLW
|
|
Packit |
67cb25 |
*/
|
|
Packit |
67cb25 |
int gsl_sf_fermi_dirac_mhalf_e(const double x, gsl_sf_result * result);
|
|
Packit |
67cb25 |
double gsl_sf_fermi_dirac_mhalf(const double x);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
/* Complete integral F_{1/2}(x)
|
|
Packit |
67cb25 |
*
|
|
Packit |
67cb25 |
* exceptions: GSL_EUNDRFLW, GSL_EOVRFLW
|
|
Packit |
67cb25 |
*/
|
|
Packit |
67cb25 |
int gsl_sf_fermi_dirac_half_e(const double x, gsl_sf_result * result);
|
|
Packit |
67cb25 |
double gsl_sf_fermi_dirac_half(const double x);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
/* Complete integral F_{3/2}(x)
|
|
Packit |
67cb25 |
*
|
|
Packit |
67cb25 |
* exceptions: GSL_EUNDRFLW, GSL_EOVRFLW
|
|
Packit |
67cb25 |
*/
|
|
Packit |
67cb25 |
int gsl_sf_fermi_dirac_3half_e(const double x, gsl_sf_result * result);
|
|
Packit |
67cb25 |
double gsl_sf_fermi_dirac_3half(const double x);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
/* Incomplete integral F_0(x,b) = ln(1 + e^(b-x)) - (b-x)
|
|
Packit |
67cb25 |
*
|
|
Packit |
67cb25 |
* exceptions: GSL_EUNDRFLW, GSL_EDOM
|
|
Packit |
67cb25 |
*/
|
|
Packit |
67cb25 |
int gsl_sf_fermi_dirac_inc_0_e(const double x, const double b, gsl_sf_result * result);
|
|
Packit |
67cb25 |
double gsl_sf_fermi_dirac_inc_0(const double x, const double b);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
__END_DECLS
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
#endif /* __GSL_SF_FERMI_DIRAC_H__ */
|