diff options
Diffstat (limited to 'gsl-1.9/specfunc/ChangeLog')
| -rw-r--r-- | gsl-1.9/specfunc/ChangeLog | 620 |
1 files changed, 620 insertions, 0 deletions
diff --git a/gsl-1.9/specfunc/ChangeLog b/gsl-1.9/specfunc/ChangeLog new file mode 100644 index 0000000..68223d6 --- /dev/null +++ b/gsl-1.9/specfunc/ChangeLog @@ -0,0 +1,620 @@ +2007-02-17 Brian Gough <bjg@network-theory.co.uk> + + * log.c (gsl_sf_log_e): removed HIDE_INLINE_STATIC + + * exp.c (gsl_sf_exp_e): removed HIDE_INLINE_STATIC + +2007-02-14 Brian Gough <bjg@network-theory.co.uk> + + * mathieu_charv.c: made solve_cubic static + +2007-02-12 Brian Gough <bjg@network-theory.co.uk> + + * mathieu_charv.c (figi): ensure that e[ii] is set when + e2[ii]==0.0 and there is no error condition, as per the original + eispack routine + +2007-02-09 Brian Gough <bjg@network-theory.co.uk> + + * ellint.c (gsl_sf_ellint_F_e): do argument reduction for phi>pi/2 + (gsl_sf_ellint_E_e): do argument reduction for phi>pi/2 + (gsl_sf_ellint_P_e): do argument reduction for phi>pi/2 + (gsl_sf_ellint_D_e): do argument reduction for phi>pi/2 + (gsl_sf_ellint_Dcomp_e): added complete D integral + (gsl_sf_ellint_Pcomp_e): added complete P integral + +2007-01-31 Brian Gough <bjg@network-theory.co.uk> + + * beta.c (gsl_sf_lnbeta_sgn_e): added to support calculations with + negative a,b + (gsl_sf_lnbeta_e): rewritten in terms of gsl_sf_lnbeta_sgn_e + (gsl_sf_beta_e): handle negative a,b + + * gamma.c (gsl_sf_lngamma_sgn_e): make error calculations an exact + copy of gsl_sf_lngamma_e (these functions could be merged to avoid + duplication) + +2007-01-29 Brian Gough <bjg@network-theory.co.uk> + + * test_legendre.c (test_legendre): added extra test cases for + underflow + +2007-01-26 Brian Gough <bjg@network-theory.co.uk> + + * expint.c (expint_E2_impl): handle x==0.0 as a special case + (expint_E2_impl): corrected error term + + * gsl_sf_log.h: removed inline version of log + + * gsl_sf_exp.h: removed inline version of exp + +2007-01-23 Brian Gough <bjg@network-theory.co.uk> + + * hyperg_1F1.c (hyperg_1F1_1_series): increase accuracy by factor + of 4 in sum, tighter convergence condition, increase error + estimate to allow for accumulated roundoff + (hyperg_1F1_1): use series when |x| > |b| + (gsl_sf_hyperg_1F1_e): only use Kummer when |x| < 100 otherwise + exponential takes extreme value + + * hyperg.c (gsl_sf_hyperg_1F1_series_e): allow 10000 iterations in + series to extend valid range + (gsl_sf_hyperg_1F1_series_e): increase accuracy by factor of 4 in + sum, tighter convergence condition + +2007-01-19 Brian Gough <bjg@network-theory.co.uk> + + * laguerre.c (laguerre_large_n): use the second term in the + asymptotic expansion from Slater, p.73. + +2007-01-17 Brian Gough <bjg@network-theory.co.uk> + + * hyperg_1F1.c (hyperg_1F1_largebx): asymptotic expansion for + large b and x, with |x|<|b| from Slater 4.3.7 + (hyperg_1F1_1): use new asymptotic expansion for |x|<|b| + (hyperg_1F1_small_a_bgt0): use new asymptotic expansion for |x|<|b| + (hyperg_1F1_renorm_b0): add neglected terms in expansion for AS13.3.7 + +2007-01-14 Brian Gough <bjg@network-theory.co.uk> + + * legendre_poly.c (gsl_sf_legendre_sphPlm_e): added explicit + computation of error term to allow for case when final term is + zero. + +2007-01-12 Brian Gough <bjg@network-theory.co.uk> + + * trig.c (gsl_sf_angle_restrict_symm_err_e): compute edge cases + more reliably, return NaN when total loss of precision + (gsl_sf_angle_restrict_pos_err_e): as above + + * legendre_poly.c (gsl_sf_legendre_Pl_e): improve error estimate + for large l by including rounding error at each step of recurrence + +2006-10-03 Brian Gough <bjg@network-theory.co.uk> + + * poch.c (gsl_sf_lnpoch_e, gsl_sf_lnpoch_sgn_e): corrected + result->val to 0.0 for x==0, previously returned incorrect value + 1.0 + +2006-09-24 Brian Gough <bjg@network-theory.co.uk> + + * laguerre.c (laguerre_large_n): work with small angles to avoid + cancellation error, computer angular reduction exactly for integer + eta. + +2006-09-22 Brian Gough <bjg@network-theory.co.uk> + + * zeta.c (gsl_sf_zeta_e): make sin_term exactly zero for negative + even integers + (gsl_sf_zetam1_int_e): return value is -1 for zetam1_int with + negative even integers + +2006-03-26 Brian Gough <bjg@network-theory.co.uk> + + * fermi_dirac.c (fd_neg): initialize s to zero (avoid spurious + warning from compiler) + +2006-02-23 Brian Gough <bjg@network-theory.co.uk> + + * 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 <bjg@network-theory.co.uk> + + * 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 <bjg@network-theory.co.uk> + + * 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 <bjg@network-theory.co.uk> + + * debye.c (gsl_sf_debye_5_e, gsl_sf_debye_6_e): added n=5,6 + (R. J. Mathar) + +2005-11-15 Brian Gough <bjg@network-theory.co.uk> + + * dilog.c (dilog_xge0): fix calculation of error estimate + +2005-08-04 Brian Gough <bjg@network-theory.co.uk> + + * coulomb.c (gsl_sf_coulomb_wave_sphF_array): fix warning about + variable shadowing for k + +2005-07-29 Brian Gough <bjg@network-theory.co.uk> + + * 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 <bjg@network-theory.co.uk> + + * 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 <bjg@network-theory.co.uk> + + * 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 <bjg@network-theory.co.uk> + + * 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 <bjg@network-theory.co.uk> + + * test_sf.c (test_jac): added tests over the full period for + elljac functions + +2005-05-23 Brian Gough <bjg@network-theory.co.uk> + + * 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 <bjg@network-theory.co.uk> + + * 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 <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 |