From da54e2f2c27f68c08f4db7551c57b2c1136dc778 Mon Sep 17 00:00:00 2001 From: Alberto Ruiz Date: Fri, 22 Feb 2008 13:08:36 +0000 Subject: rm modified gsl_sf_exp,log --- lib/Numeric/GSL/Special/autoall.sh | 4 +- lib/Numeric/GSL/Special/gsl_sf_exp.h | 146 ----------------------------------- lib/Numeric/GSL/Special/gsl_sf_log.h | 90 --------------------- 3 files changed, 2 insertions(+), 238 deletions(-) delete mode 100644 lib/Numeric/GSL/Special/gsl_sf_exp.h delete mode 100644 lib/Numeric/GSL/Special/gsl_sf_log.h (limited to 'lib/Numeric/GSL') diff --git a/lib/Numeric/GSL/Special/autoall.sh b/lib/Numeric/GSL/Special/autoall.sh index 1581e38..da349a4 100644 --- a/lib/Numeric/GSL/Special/autoall.sh +++ b/lib/Numeric/GSL/Special/autoall.sh @@ -19,7 +19,7 @@ rep ', coupling_6j_INCORRECT_e\n, coupling_6j_INCORRECT\n' '' Coupling.hs ./auto.hs elementary ./auto.hs ellint ./auto.hs erf -./auto.hs exp gsl_sf_exp.h +./auto.hs exp rep ', exp\n' ', Numeric.GSL.Special.Exp.exp\n' Exp.hs ./auto.hs expint ./auto.hs fermi_dirac @@ -29,7 +29,7 @@ rep ', exp\n' ', Numeric.GSL.Special.Exp.exp\n' Exp.hs ./auto.hs laguerre ./auto.hs lambert ./auto.hs legendre gsl_sf_legendre.h -./auto.hs log gsl_sf_log.h +./auto.hs log rep ', log\n' ', Numeric.GSL.Special.Log.log\n' Log.hs ./auto.hs pow_int ./auto.hs psi diff --git a/lib/Numeric/GSL/Special/gsl_sf_exp.h b/lib/Numeric/GSL/Special/gsl_sf_exp.h deleted file mode 100644 index b1f0d89..0000000 --- a/lib/Numeric/GSL/Special/gsl_sf_exp.h +++ /dev/null @@ -1,146 +0,0 @@ -/* specfunc/gsl_sf_exp.h - * - * Copyright (C) 1996, 1997, 1998, 1999, 2000, 2004 Gerard Jungman - * - * This program is free software; you can redistribute it and/or modify - * it under the terms of the GNU General Public License as published by - * the Free Software Foundation; either version 2 of the License, or (at - * your option) any later version. - * - * This program is distributed in the hope that it will be useful, but - * WITHOUT ANY WARRANTY; without even the implied warranty of - * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU - * General Public License for more details. - * - * You should have received a copy of the GNU General Public License - * along with this program; if not, write to the Free Software - * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA. - */ - -/* Author: G. Jungman */ - -#ifndef __GSL_SF_EXP_H__ -#define __GSL_SF_EXP_H__ - -#include -#include - -#undef __BEGIN_DECLS -#undef __END_DECLS -#ifdef __cplusplus -# define __BEGIN_DECLS extern "C" { -# define __END_DECLS } -#else -# define __BEGIN_DECLS /* empty */ -# define __END_DECLS /* empty */ -#endif - -__BEGIN_DECLS - -/* Provide an exp() function with GSL semantics, - * i.e. with proper error checking, etc. - * - * exceptions: GSL_EOVRFLW, GSL_EUNDRFLW - */ -int gsl_sf_exp_e(const double x, gsl_sf_result * result); -double gsl_sf_exp(const double x); - - -/* Exp(x) - * - * exceptions: GSL_EOVRFLW, GSL_EUNDRFLW - */ -int gsl_sf_exp_e10_e(const double x, gsl_sf_result_e10 * result); - - -/* Exponentiate and multiply by a given factor: y * Exp(x) - * - * exceptions: GSL_EOVRFLW, GSL_EUNDRFLW - */ -int gsl_sf_exp_mult_e(const double x, const double y, gsl_sf_result * result); -double gsl_sf_exp_mult(const double x, const double y); - - -/* Exponentiate and multiply by a given factor: y * Exp(x) - * - * exceptions: GSL_EOVRFLW, GSL_EUNDRFLW - */ -int gsl_sf_exp_mult_e10_e(const double x, const double y, gsl_sf_result_e10 * result); - - -/* exp(x)-1 - * - * exceptions: GSL_EOVRFLW - */ -int gsl_sf_expm1_e(const double x, gsl_sf_result * result); -double gsl_sf_expm1(const double x); - - -/* (exp(x)-1)/x = 1 + x/2 + x^2/(2*3) + x^3/(2*3*4) + ... - * - * exceptions: GSL_EOVRFLW - */ -int gsl_sf_exprel_e(const double x, gsl_sf_result * result); -double gsl_sf_exprel(const double x); - - -/* 2(exp(x)-1-x)/x^2 = 1 + x/3 + x^2/(3*4) + x^3/(3*4*5) + ... - * - * exceptions: GSL_EOVRFLW - */ -int gsl_sf_exprel_2_e(double x, gsl_sf_result * result); -double gsl_sf_exprel_2(const double x); - - -/* Similarly for the N-th generalization of - * the above. The so-called N-relative exponential - * - * exprel_N(x) = N!/x^N (exp(x) - Sum[x^k/k!, {k,0,N-1}]) - * = 1 + x/(N+1) + x^2/((N+1)(N+2)) + ... - * = 1F1(1,1+N,x) - */ -int gsl_sf_exprel_n_e(const int n, const double x, gsl_sf_result * result); -double gsl_sf_exprel_n(const int n, const double x); - - -/* Exponentiate a quantity with an associated error. - */ -int gsl_sf_exp_err_e(const double x, const double dx, gsl_sf_result * result); - -/* Exponentiate a quantity with an associated error. - */ -int gsl_sf_exp_err_e10_e(const double x, const double dx, gsl_sf_result_e10 * result); - - -/* Exponentiate and multiply by a given factor: y * Exp(x), - * for quantities with associated errors. - * - * exceptions: GSL_EOVRFLW, GSL_EUNDRFLW - */ -int gsl_sf_exp_mult_err_e(const double x, const double dx, const double y, const double dy, gsl_sf_result * result); - - -/* Exponentiate and multiply by a given factor: y * Exp(x), - * for quantities with associated errors. - * - * exceptions: GSL_EOVRFLW, GSL_EUNDRFLW - */ -int gsl_sf_exp_mult_err_e10_e(const double x, const double dx, const double y, const double dy, gsl_sf_result_e10 * result); - -__END_DECLS - - -#ifdef HAVE_INLINE -#include -#include - -__BEGIN_DECLS - - - -__END_DECLS - -#endif /* HAVE_INLINE */ - - -#endif /* __GSL_SF_EXP_H__ */ diff --git a/lib/Numeric/GSL/Special/gsl_sf_log.h b/lib/Numeric/GSL/Special/gsl_sf_log.h deleted file mode 100644 index 5225d05..0000000 --- a/lib/Numeric/GSL/Special/gsl_sf_log.h +++ /dev/null @@ -1,90 +0,0 @@ -/* specfunc/gsl_sf_log.h - * - * Copyright (C) 1996, 1997, 1998, 1999, 2000, 2004 Gerard Jungman - * - * This program is free software; you can redistribute it and/or modify - * it under the terms of the GNU General Public License as published by - * the Free Software Foundation; either version 2 of the License, or (at - * your option) any later version. - * - * This program is distributed in the hope that it will be useful, but - * WITHOUT ANY WARRANTY; without even the implied warranty of - * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU - * General Public License for more details. - * - * You should have received a copy of the GNU General Public License - * along with this program; if not, write to the Free Software - * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA. - */ - -/* Author: G. Jungman */ - -#ifndef __GSL_SF_LOG_H__ -#define __GSL_SF_LOG_H__ - -#include - -#undef __BEGIN_DECLS -#undef __END_DECLS -#ifdef __cplusplus -# define __BEGIN_DECLS extern "C" { -# define __END_DECLS } -#else -# define __BEGIN_DECLS /* empty */ -# define __END_DECLS /* empty */ -#endif - -__BEGIN_DECLS - - -/* Provide a logarithm function with GSL semantics. - * - * exceptions: GSL_EDOM - */ -int gsl_sf_log_e(const double x, gsl_sf_result * result); -double gsl_sf_log(const double x); - - -/* Log(|x|) - * - * exceptions: GSL_EDOM - */ -int gsl_sf_log_abs_e(const double x, gsl_sf_result * result); -double gsl_sf_log_abs(const double x); - - -/* Complex Logarithm - * exp(lnr + I theta) = zr + I zi - * Returns argument in [-pi,pi]. - * - * exceptions: GSL_EDOM - */ -int gsl_sf_complex_log_e(const double zr, const double zi, gsl_sf_result * lnr, gsl_sf_result * theta); - - -/* Log(1 + x) - * - * exceptions: GSL_EDOM - */ -int gsl_sf_log_1plusx_e(const double x, gsl_sf_result * result); -double gsl_sf_log_1plusx(const double x); - - -/* Log(1 + x) - x - * - * exceptions: GSL_EDOM - */ -int gsl_sf_log_1plusx_mx_e(const double x, gsl_sf_result * result); -double gsl_sf_log_1plusx_mx(const double x); - - -#ifdef HAVE_INLINE -#include -#include - -#endif /* HAVE_INLINE */ - - -__END_DECLS - -#endif /* __GSL_SF_LOG_H__ */ -- cgit v1.2.3