changelog

数学公式库--非常不错 The GNU Scientific Library (GSL) is a collection of routines for numerical computing. T

2006-03-26  Brian Gough  <[e-mail]>

	* fermi_dirac.c (fd_neg): initialize s to zero (avoid spurious
	warning from compiler)

2006-02-23  Brian Gough  <[e-mail]>

	* coulomb.c (gsl_sf_coulomb_wave_FG_e): fixed sign of F_lam_min,
	covers case when k_lam_G is nonzero and F_lam_min and F_lam differ
	in sign.

2006-01-21  Brian Gough  <[e-mail]>

	* synchrotron.c (gsl_sf_synchrotron_1_e): added first order
	correction term for the taylor expansion
	(gsl_sf_synchrotron_2_e): added first order correction term for
	the taylor expansion

2006-01-20  Brian Gough  <[e-mail]>

	* bessel_j.c (gsl_sf_bessel_jl_e): limit the use of
	gsl_sf_bessel_Jnu_asympx_e to the range x>100*l*l to satisfy
	the requirement x>>l*l in the asymptotic expansion

	* bessel_In.c (gsl_sf_bessel_In_scaled_e): limit the use of
	gsl_sf_bessel_I_CF1_ser to the region where the continued
	converges with the allowed 20000 terms (|x|<1e7)

2005-12-20  Brian Gough  <[e-mail]>

	* debye.c (gsl_sf_debye_5_e, gsl_sf_debye_6_e): added n=5,6
	(R. J. Mathar)

2005-11-15  Brian Gough  <[e-mail]>

	* dilog.c (dilog_xge0): fix calculation of error estimate

2005-08-04  Brian Gough  <[e-mail]>

	* coulomb.c (gsl_sf_coulomb_wave_sphF_array): fix warning about
	variable shadowing for k

2005-07-29  Brian Gough  <[e-mail]>

	* gamma_inc.c (gsl_sf_gamma_inc_Q_e): use continued fraction close
	to peak of the integrand x > a - sqrt(a)

2005-07-28  Brian Gough  <[e-mail]>

	* elljac.c (gsl_sf_elljac_e): use separate iterations to avoid
	division by zero, new algorithm from Bulirsch avoids inverse trig
	functions.

2005-07-03  Brian Gough  <[e-mail]>

	* legendre_poly.c (gsl_sf_legendre_sphPlm_e): compute exponential
	internally to avoid underflow error when calling gsl_sf_exp_err

2005-07-02  Brian Gough  <[e-mail]>

	* bessel_i.c (gsl_sf_bessel_il_scaled_array): handle x==0 as a
	special case

	* bessel_k.c (gsl_sf_bessel_kl_scaled_array): added lmax==0 as a
	special case

	* bessel_y.c (gsl_sf_bessel_yl_array): added lmax==0 as a special
	case

	* transport.c (gsl_sf_transport_2_e): improve error estimate for
	small x

2005-06-26  Brian Gough  <[e-mail]>

	* test_sf.c (test_jac): added tests over the full period for
	elljac functions

2005-05-23  Brian Gough  <[e-mail]>

	* test_bessel.c (test_bessel): added test for steed(99,1,...)

	* legendre_H3d.c (legendre_H3d_CF1): use hypot
	(legendre_H3d_CF1_ser): use hypot
	(gsl_sf_legendre_H3d_e): use hypot
	(gsl_sf_legendre_H3d_e): use hypot
	(gsl_sf_legendre_H3d_array): use hypot

	* dilog.c (dilogc_unitdisk): use hypot

	* coulomb.c (gsl_sf_coulomb_CL_array): use hypot
	(coulomb_F_recur): use hypot
	(coulomb_G_recur): use hypot
	(coulomb_jwkb): use hypot
	(gsl_sf_coulomb_wave_F_array): use hypot
	(gsl_sf_coulomb_wave_FG_array): use hypot
	(gsl_sf_coulomb_wave_FGp_array): use hypot

	* bessel_j.c (gsl_sf_bessel_jl_steed_array): use hypot

	* bessel.c (gsl_sf_bessel_Inu_scaled_asymp_unif_e): use hypot
	(gsl_sf_bessel_Knu_scaled_asymp_unif_e): use hypot

2005-03-02  Brian Gough  <[e-mail]>

	* coulomb_bound.c (gsl_sf_hydrogenicR_e): don't check for
	underflow when function is known to be zero (e.g. r=0 l>0 or at
	zeroes of the laguerre polynomial).

2004-12-29  Brian Gough  <[e-mail]>

	* 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  <[e-mail]>

	* 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  <[e-mail]>

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

2004-11-12  Brian Gough  <[e-mail]>

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

2004-10-11  Brian Gough  <[e-mail]>

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

2004-06-03  Brian Gough  <[e-mail]>

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

2003-11-30  Brian Gough  <[e-mail]>

	* gsl_sf_log.h: added missing headers for inline functions

2003-08-11  Brian Gough  <[e-mail]>

	* 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  <[e-mail]>

	* 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  <[e-mail]>

	* 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  <[e-mail]>

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

2003-06-02  Brian Gough  <[e-mail]>

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

2003-04-18  Brian Gough  <[e-mail]>

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

2003-04-12  Brian Gough  <[e-mail]>

	* 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  <[e-mail]>

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

2003-01-25  Brian Gough  <[e-mail]>

	* 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  <[e-mail]>

	* 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  <[e-mail]>

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

Sat Jul 13 23:11:37 2002  Brian Gough  <[e-mail]>

	* 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  <[e-mail]>

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

Tue May 28 21:04:29 2002  Brian Gough  <[e-mail]>

	* 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  <[e-mail]>

	* 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  <[e-mail]>

	* 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  <[e-mail]>

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

Tue Jan  8 14:31:04 2002  Brian Gough  <[e-mail]>

	* 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