diff options
author | schnetter <schnetter@0f96c2a6-8d4f-0410-8182-94f2e8281398> | 2008-10-10 16:22:36 +0000 |
---|---|---|
committer | schnetter <schnetter@0f96c2a6-8d4f-0410-8182-94f2e8281398> | 2008-10-10 16:22:36 +0000 |
commit | 5eeb2401d99c2458960db1d6a25527f9c5f07a78 (patch) | |
tree | dd5712ef14e29b48e09b92da2f9924e990c1e9d3 /src/gsl_sf_erf.c | |
parent | dd9e8d6928f79c2b3abad885a9d68c1ea436c744 (diff) |
git-svn-id: http://svn.aei.mpg.de/numrel/AEIThorns/GSLFortran/trunk@2 0f96c2a6-8d4f-0410-8182-94f2e8281398
Diffstat (limited to 'src/gsl_sf_erf.c')
-rw-r--r-- | src/gsl_sf_erf.c | 118 |
1 files changed, 118 insertions, 0 deletions
diff --git a/src/gsl_sf_erf.c b/src/gsl_sf_erf.c new file mode 100644 index 0000000..1b38a76 --- /dev/null +++ b/src/gsl_sf_erf.c @@ -0,0 +1,118 @@ +#include <gsl/gsl_sf_erf.h> + +#include "cctk.h" + + + +/* Complementary Error Function + * erfc(x) := 2/Sqrt[Pi] Integrate[Exp[-t^2], {t,x,Infinity}] + * + * exceptions: none + */ +int CCTK_FCALL +CCTK_FNAME(gsl_sf_erfc_e) (double const * restrict const x, + gsl_sf_result * restrict const result) +{ + return gsl_sf_erfc_e (* x, result); +} + +double CCTK_FCALL +CCTK_FNAME(gsl_sf_erfc) (double const * restrict const x) +{ + return gsl_sf_erfc (* x); +} + + +/* Log Complementary Error Function + * + * exceptions: none + */ +int CCTK_FCALL +CCTK_FNAME(gsl_sf_log_erfc_e) (double const * restrict const x, + gsl_sf_result * restrict const result) +{ + return gsl_sf_log_erfc_e (* x, result); +} + +double CCTK_FCALL +CCTK_FNAME(gsl_sf_log_erfc) (double const * restrict const x) +{ + return gsl_sf_log_erfc (* x); +} + + +/* Error Function + * erf(x) := 2/Sqrt[Pi] Integrate[Exp[-t^2], {t,0,x}] + * + * exceptions: none + */ +int CCTK_FCALL +CCTK_FNAME(gsl_sf_erf_e) (double const * restrict const x, + gsl_sf_result * restrict const result) +{ + return gsl_sf_erf_e (* x, result); +} + +double CCTK_FCALL +CCTK_FNAME(gsl_sf_erf) (double const * restrict const x) +{ + return gsl_sf_erf (* x); +} + + +/* Probability functions: + * Z(x) : Abramowitz+Stegun 26.2.1 + * Q(x) : Abramowitz+Stegun 26.2.3 + * + * exceptions: none + */ +int CCTK_FCALL +CCTK_FNAME(gsl_sf_erf_Z_e) (double const * restrict const x, + gsl_sf_result * restrict const result) +{ + return gsl_sf_erf_Z_e (* x, result); +} + +int CCTK_FCALL +CCTK_FNAME(gsl_sf_erf_Q_e) (double const * restrict const x, + gsl_sf_result * restrict const result) +{ + return gsl_sf_erf_Q_e (* x, result); +} + +double CCTK_FCALL +CCTK_FNAME(gsl_sf_erf_Z) (double const * restrict const x) +{ + return gsl_sf_erf_Z (* x); +} + +double CCTK_FCALL +CCTK_FNAME(gsl_sf_erf_Q) (double const * restrict const x) +{ + return gsl_sf_erf_Q (* x); +} + + +#if 0 +/* Does not exist in older versions of GSL*/ + +/* Hazard function, also known as the inverse Mill's ratio. + * + * H(x) := Z(x)/Q(x) + * = Sqrt[2/Pi] Exp[-x^2 / 2] / Erfc[x/Sqrt[2]] + * + * exceptions: GSL_EUNDRFLW + */ +int CCTK_FCALL +CCTK_FNAME(gsl_sf_hazard_e) (double const * restrict const x, + gsl_sf_result * restrict const result) +{ + return gsl_sf_hazard_e (* x, result); +} + +double CCTK_FCALL +CCTK_FNAME(gsl_sf_hazard) (double const * restrict const x) +{ + return gsl_sf_hazard (* x); +} +#endif |