2005-05-06  Peter A. Cejchan  <pac7@post.cz>
	native port to Plan9:
	#if --> #ifdef
	#if 0 --> #ifdef NEVER
	no #include in header files
	return GSL_FAILURE adder to funcs that call exits()
	GSL_ERROR --> rewritten with exits()
	header files stripped-off to plain form

2004-12-29  Brian Gough  <bjg@network-theory.co.uk>

	* dilog.c (gsl_sf_complex_dilog_e): use const consistently in
	arguments of declaration and definition
	(gsl_sf_complex_dilog_xy_e): as above

2004-12-26  Brian Gough  <bjg@network-theory.co.uk>

	* gamma_inc.c (gamma_inc_D): improve error estimate for case of
	u=x/a to include cancellation errors and only use it when x < 0.5*a
	since the cancellation errors are significant for x/a ~ 1

2004-12-23  Brian Gough  <bjg@network-theory.co.uk>

	* gsl_sf_coupling.h: fixed declaration to
	gsl_sf_coupling_6j_INCORRECT instead of
	gsl_sf_coupling_INCORRECT_6j

2004-11-12  Brian Gough  <bjg@network-theory.co.uk>

	* psi.c (gsl_sf_psi_1): added missing function definition
	gsl_sf_psi_1

2004-10-11  Brian Gough  <bjg@network-theory.co.uk>

	* expint.c (gsl_sf_expint_Ei_scaled): fixed call to incorrect
	function gsl_sf_expint_Ei_e

2004-06-03  Brian Gough  <bjg@network-theory.co.uk>

	* psi.c (psi_n_xg0): added missing return type int

2003-11-30  Brian Gough  <bjg@network-theory.co.uk>

	* gsl_sf_log.h: added missing headers for inline functions

2003-08-11  Brian Gough  <bjg@network-theory.co.uk>

	* test_sf.c: added preprocessor definitions TEST_FACTOR and
	TEST_SIGMA to allow larger tolerances on released versions (to
	reduce the number of spurious bug reports).

2003-07-24  Brian Gough  <bjg@network-theory.co.uk>

	* gamma.c (gsl_sf_choose_e): declare k as unsigned int instead of
	int
	(gsl_sf_gamma_e): avoid shadowed declaration of sgn

	* bessel_sequence.c (gsl_sf_bessel_sequence_Jnu_e): declare i as
	size_t instead of int

2003-06-11  Brian Gough  <bjg@network-theory.co.uk>

	* coupling.c (gsl_sf_coupling_6j_INCORRECT_e): fixed typo in test
	for two_jf < 0
	(gsl_sf_coupling_6j_e): moved the implementation from the
	INCORRECT function into the correct function

2003-06-05  Brian Gough  <bjg@network-theory.co.uk>

	* test_sf.c (test_coupling): added some regression tests for
	coupling functions

2003-06-02  Brian Gough  <bjg@network-theory.co.uk>

	* gamma_inc.c (gamma_inc_D): handle x<a, a>10 as a special case

2003-04-18  Brian Gough  <bjg@network-theory.co.uk>

	* gamma.c (gsl_sf_gammainv_e): handle any singularities in
	gamma(x) up front and return zero directly

2003-04-12  Brian Gough  <bjg@network-theory.co.uk>

	* psi.c: changed value of psi(1,1) in table psi_1_table to be
	positive ((-)^2 * 1! * zeta(2)), in accordance with Abramowitz &
	Stegun 6.4.2.

2003-04-08  G. Jungman

    * erfc.c, gsl_sf_erf.h: added gsl_sf_hazard_e() and
    gsl_sf_hazard() functions

2003-04-07  G. Jungman

    * gamma_inc.c, gsl_sf_gamma.h: added gsl_sf_gamma_inc_e() and
    gsl_sf_gamma_inc(), implmentations of the non-normalized incomplete
    gamma function

2003-03-09  Brian Gough  <bjg@network-theory.co.uk>

	* gsl_sf_legendre.h: added missing const to formal parameter on
	gsl_sf_legendre_sphPlm_deriv_array

2003-01-25  Brian Gough  <brian.gough@network-theory.co.uk>

	* test_gamma.c (test_gamma): added a test for
	gamma_inc_P(10,1e-16) (BUG#12)

Sat Sep  7 15:56:15 2002  Brian Gough  <bjg@network-theory.co.uk>

	* test_sf.h (TEST_FACTOR): include an overall factor in the test
 	tolerances, otherwise there are too many bug reports for minor
 	failures.

Sun Sep  1 17:34:27 2002  Brian Gough  <bjg@network-theory.co.uk>

	* test_legendre.c (test_legendre): increased tolerance on test

Sat Jul 13 23:11:37 2002  Brian Gough  <bjg@network-theory.co.uk>

	* ellint.c (gsl_sf_ellint_Kcomp_e): improved error estimate to
 	take cancellation of y=1-k^2 into account near k=1 

Sun Jul  7 21:39:12 2002  Brian Gough  <bjg@network-theory.co.uk>

	* test_bessel.c (test_bessel): increased tolerance on test of
 	bessel_Jn_array

Tue May 28 21:04:29 2002  Brian Gough  <bjg@network-theory.co.uk>

	* psi.c (gsl_sf_psi_1piy_e): changed log(y) to log(ay) since
 	function is even and can be extended to negative y

Mon Jan 28 19:02:42 2002  Brian Gough  <bjg@network-theory.co.uk>

	* gamma_inc.c (gamma_inc_Q_CF): express 'small' constant in terms
 	of GSL_DBL_EPSILON
	(gamma_inc_Q_CF): patch for more reliable series from Hans Plesser

Sat Jan 26 17:33:29 2002  Brian Gough  <bjg@network-theory.co.uk>

	* test_gamma.c (test_gamma): increased tolerance on a test

	* test_hyperg.c (test_hyperg): increased tolerance on a couple of
 	tests

Fri Jan 18 15:12:30 2002  Brian Gough  <bjg@network-theory.co.uk>

	* test_airy.c (test_airy): increased tolerance on test

Tue Jan  8 14:31:04 2002  Brian Gough  <bjg@network-theory.co.uk>

	* test_legendre.c (test_legendre): increased tolerance by one
 	level on _array tests

	* hyperg_1F1.c (hyperg_1F1_small_a_bgt0): fix branch for a==1

Fri Dec  7 12:38:40 2001  Brian Gough  <bjg@network-theory.co.uk>

	* laguerre.c (laguerre_n_cp): catch internal error, not global
 	error

	* error.h (INTERNAL_OVERFLOW_ERROR): added internal error macros
 	which do not call the error handler.

Wed Dec  5 19:25:26 2001  Brian Gough  <bjg@network-theory.co.uk>

	* gamma.c (gamma_xgthalf): return gamma(x) exactly when x is an
 	integer

	* hyperg_1F1.c (hyperg_1F1_ab_posint):  fix bug in recurrence
 	initialisation, as done for hyperg_1F1_ab_pos

Thu Oct 18 11:37:25 2001  Brian Gough  <bjg@network-theory.co.uk>

	* coulomb.c (gsl_sf_coulomb_CL_array): renamed from
 	gsl_sf_coulomb_CL_list for consistency

Sat Oct 13 15:55:56 2001  Brian Gough  <bjg@network-theory.co.uk>

	* cheb_eval.c (cheb_eval_e): keep track of cancellation error to
 	prevent underestimates of total error.

Fri Oct 12 16:39:27 2001  Brian Gough  <bjg@network-theory.co.uk>

	* test_hyperg.c (test_hyperg): increased tolerance from TOL0 to
 	TOL1 on one test to allow it to pass under different optimizations

Thu Oct 11 14:21:34 2001  Brian Gough  <bjg@network-theory.co.uk>

	* gegenbauer.c (gsl_sf_gegenpoly_n_e): initialize variable gk to
 	zero to avoid warning

	* bessel_i.c (gsl_sf_bessel_il_scaled_e): initialize variable
 	iellm1 to zero to avoid warning

	* bessel_Jnu.c (gsl_sf_bessel_Jnu_e): initialize variable Ynp1 to
 	zero to avoid warning

	* legendre_poly.c (gsl_sf_legendre_sphPlm_e): initialize variables
 	p_ell, y_ell to zero to avoid warning

Thu Sep  6 10:38:51 2001  Brian Gough  <bjg@network-theory.co.uk>

	* clausen.c: added missing initialiser to cheb_series struct

	* pow_int.c (gsl_sf_pow_int_e): handle the case x == 0, n < 0

	* legendre_poly.c (gsl_sf_legendre_array_size): added missing
 	static version of this inline function

Wed Aug 15 20:18:43 2001  Brian Gough  <bjg@network-theory.co.uk>

	* hyperg_1F1.c (hyperg_1F1_ab_pos): fix bug in recurrence
 	initialisation
	(hyperg_1F1_ab_neg): if all else fails, try the series.

Wed Aug  8 19:55:34 2001  Brian Gough  <bjg@network-theory.co.uk>

	* test_sf.c (test_coupling): used analytic expressions for the
 	exact values to problems with insufficient precision

Thu Jul 12 15:25:04 2001  Brian Gough  <bjg@network-theory.co.uk>

	* test_legendre.c (test_legendre): corrected tests to use exact
	floating point arguments near singularities.

	* legendre_Qn.c (gsl_sf_legendre_Q0_e): removed unnecessary error
 	terms that I added
	(gsl_sf_legendre_Q1_e): removed unnecessary error terms that I
 	added

Wed Jul 11 16:06:45 2001  Brian Gough  <bjg@network-theory.co.uk>

	* error.h: moved domain, overflow and underflow errors into
 	macros, return Nan for domain error, Inf for overflow.

Tue Jul 10 22:00:55 2001  Brian Gough  <bjg@network-theory.co.uk>

	* erfc.c (gsl_sf_log_erfc_e): introduce an additional log(erfc)
 	branch to prevent loss of accuracy for small x

	* test_sf.c (test_zeta): increased tolerance on zeta function test
 	to take into account differences in rounding

2001-07-10  Brian Gough  <bjg@hppav.local>

	* legendre_Qn.c (gsl_sf_legendre_Q0_e): use (log1p(x) - log1p(-x))
	instead of log((1+x)/(1-x)), for accuracy
	(gsl_sf_legendre_Q1_e): use (log1p(x) - log1p(-x)) instead of
	log((1+x)/(1-x)), for accuracy
	(gsl_sf_legendre_Q0_e): improve error estimate near singular
	points
	(gsl_sf_legendre_Q1_e): improve error estimate near singular
	points
	(gsl_sf_legendre_Q0_e): fixed incorrect factor of 2 in asymptotic
	expansion

	* check.h (CHECK_UNDERFLOW): provide macro for checking results
	for underflow

Fri Jul  6 20:16:19 2001  Brian Gough  <bjg@network-theory.co.uk>

	* zeta.c (riemann_zeta1m_slt0): added a function to compute
 	zeta(1-s) without loss of accuracy for s near zero, as used by the
 	reflection formula.  This fixes a bug in the accuracy of results
 	of zeta(-x) for small x, where loss of precision previously
 	occurred by evaluating 1-x.

2001-07-06  Brian Gough  <bjg@network-theory.co.uk>

	* bessel_I0.c (gsl_sf_bessel_I0_scaled_e): added missing scaling
	factor of exp(-x) for case x<2sqrt(epsilon)

Thu Jul  5 16:16:13 2001  Brian Gough  <bjg@network-theory.co.uk>

	* erfc.c (gsl_sf_erfc_e): corrected case of -10<x<-5 to take into
 	account that chebyshev fit is defined in terms of positive x.

Tue Jul  3 12:09:33 2001  Brian Gough  <bjg@network-theory.co.uk>

	* eval.h (EVAL_DOUBLE): avoid returning the status value as a
 	numerical result

2001-06-28  Brian Gough  <bjg@network-theory.co.uk>

	* elementary.c (gsl_sf_multiply_e): use GSL_COERCE_DBL macro to
	deal with extended register problem

	* coupling.c: include stdlib.h for prototype of abs()

Wed Jun 27 21:20:22 2001  Brian Gough  <bjg@network-theory.co.uk>

	* test_sf.c (test_fermidirac): increased tolerance on test of
 	gsl_sf_fermi_dirac_int(9,500) by factor of two to allow for MSVC

Tue Jun 26 12:10:17 2001  Brian Gough  <bjg@network-theory.co.uk>

	* gsl_sf_gamma.h: added const to prototype, to match function
 	definition

Sun May  6 13:01:01 2001  Brian Gough  <bjg@network-theory.co.uk>

	* test_sf.c: initialize message_buff to null string to prevent
 	random segmentation faults

	* test_sf.h: simplified test macros to reduce stack usage

Mon Apr 30 12:38:15 2001  Brian Gough  <bjg@network-theory.co.uk>

	* airy_zero.c bessel_zero.c: zeros are now counted using an
 	unsigned int

	* poly.c: moved poly_eval function into poly/ directory

Wed Apr 25 17:28:48 2001  Brian Gough  <bjg@network-theory.co.uk>

	* trig.c (gsl_sf_polar_to_rect): dropped _e from polar/rect
 	conversion functions

Tue Apr 24 17:05:22 2001  Brian Gough  <bjg@network-theory.co.uk>

	* split out chebyshev functions into their own cheb/ directory,
 	leaving behind an internal cheb struct and cheb evaluation
 	function.

Mon Apr 23 10:21:06 2001  Brian Gough  <bjg@network-theory.co.uk>

	* changed tests for EFAULT into a commented-out macro, since
 	EFAULT should only apply to invalid non-null pointers, but it
 	might be useful to keep the tests around for debugging in this
 	directory.

	* unified error handling conventions to _e for error handling
 	functions and no suffix for plain functions, so _impl functions
 	are no longer needed.

1999-08-02  Mark Galassi  <rosalia@lanl.gov>

	* fermi_dirac.c: took out the use of some "const int" constants
	which were being used to size arrays, since this is not portable
	(thanks to Bernd Petrovitsch <bernd@ict.tuwien.ac.at> for pointing
	this out).

1999-01-02  Mark Galassi  <rosalia@cygnus.com>

	* trig.c (gsl_sf_rect_to_polar_impl): introduced an #ifdef
	HAVE_HYPOT for systems that do not have the hypot() routine.

Sun Feb 14 20:59:50 1999  Brian Gough  <bjg@netsci.freeserve.co.uk>

	* Makefile.am (include_HEADERS): added gsl_sf_result.h

1998-12-05  Mark Galassi  <rosalia@cygnus.com>

	* Makefile.am: t-erf.c was commented out, which is a problem
	(pointed out by automake-1.3b).  Moved the commented-out t-erf.c
	line to the end of the long list of files.

	* legendre_poly.c (gsl_sf_legendre_Pl_impl):
	* legendre_con.c (gsl_sf_conicalP_xgt1_neg_mu_largetau_impl):
	added an ugly fix to the invocation of gsl_sf_bessel_Jn_impl() to
	give it two bogus arguments so that it builds.  I hope Jerry fixes
	it up soon!

	* gsl_sf_airy.h: added an include of gsl_precision.h so that
	gsl_prec_t is defined.

Tue Nov 17 17:29:31 1998  Brian Gough  <bjg@vvv.lanl.gov>

	* added #include <config.h> to all top-level source files

	* chebyshev.c (gsl_sf_cheb_eval_n): changed local gslMIN and
 	gslMAX to the standard GSL_MIN and GSL_MAX macros

Tue Aug 18 13:36:15 1998  Brian Gough  <bjg@vvv.lanl.gov>

	* coulomb.c (C0sq): changed to using gsl_sf_expm1_impl instead of
 	expm1, since the latter is a GNU extension.

Mon Aug  3 14:43:16 1998  Brian Gough  <bjg@vvv.lanl.gov>

	* bessel_amp_phase.h: undefined consts are now declared extern

Mon Jul 13 12:40:27 1998  Brian Gough  <bjg@vvv.lanl.gov>

	* gamma.c: changed all the factorial functions to take unsigned
 	int arguments, since negative values are not allowed.
	(gsl_sf_choose): fixed off by one error in call to factorial.

Sun Jul 12 13:21:41 1998  Brian Gough  <bjg@vvv.lanl.gov>

	* gsl_sf_legendre.h, gsl_sf_poly.h, gsl_sf_pow_int.h: added
 	HAVE_INLINE around some inline definitions in the header files

	* gamma.c: implemented the natural versions of gsl_sf_lnchoose and
	gsl_sf_choose.

Wed Apr 15 21:27:48 1998  Brian Gough  <bjg@vvv.lanl.gov>

	* added the return code as the second argument of GSL_WARNING, so
 	we can filter on the error number or output a standard message for
 	each one.

1998-02-02  Mark Galassi  <rosalia@cygnus.com>

	* Makefile.am (include_HEADERS): added gsl_specfunc.h so it gets
	into the distribution.
	(INCLUDES): added -I$(top_srcdir) so that bessel.c can find
	gsl_math.h
