htswanalysis/MACS/lib/gsl/gsl-1.11/specfunc/ChangeLog

   1 2008-03-17  Brian Gough  <bjg@network-theory.co.uk>
   2
   3         * hyperg_2F1.c (gsl_sf_hyperg_2F1_e): replace duplicate error
   4         check on stat3 by stat4
   5
   6 2008-03-15  Brian Gough  <bjg@network-theory.co.uk>
   7
   8         * bessel.c (gsl_sf_bessel_Jnu_asympx_e): use full asymptotic
   9         series
  10
  11 2008-03-14  Brian Gough  <bjg@network-theory.co.uk>
  12
  13         * bessel_j.c (gsl_sf_bessel_jl_e): increase error estimate by
  14         factor of 2
  15
  16 2008-02-09  Brian Gough  <bjg@network-theory.co.uk>
  17
  18         * exp.c (gsl_sf_exp_e10_e): only use exponent e10 when standard
  19         exp() is out of range
  20         (gsl_sf_exp_mult_e10_e): add full set of error terms
  21
  22 2008-01-10  Brian Gough  <bjg@network-theory.co.uk>
  23
  24         * hyperg_2F1.c (gsl_sf_hyperg_2F1_e): handle case of x==1
  25
  26 2007-10-25  Brian Gough  <bjg@network-theory.co.uk>
  27
  28         * bessel.c (gsl_sf_bessel_J_CF1): handle underflow by rescaling in
  29         recurrence
  30
  31 2007-09-13  Brian Gough  <bjg@network-theory.co.uk>
  32
  33         * ellint.c (gsl_sf_ellint_Pcomp_e): remove test for n <1, no
  34         restriction needed
  35
  36 2007-09-10  Brian Gough  <bjg@network-theory.co.uk>
  37
  38         * expint.c (gsl_sf_expint_En_e): added for support En(x)
  39         (gsl_sf_expint_En_scaled_e): added for support En(x) scaled
  40
  41         * gamma_inc.c (gamma_inc_CF): include finite precision of log term
  42
  43 2007-08-30  Brian Gough  <bjg@network-theory.co.uk>
  44
  45         * psi.c (psi_complex_rhp): accumulate error, to allow for
  46         cancellation effects
  47
  48         * beta.c (gsl_sf_lnbeta_sgn_e): added missing factor of 2 for
  49         error.
  50
  51 2007-08-27  Brian Gough  <bjg@network-theory.co.uk>
  52
  53         * beta_inc.c (gsl_sf_beta_inc_e): handle cases where a<=0 or b<=0
  54
  55 2007-04-27  Brian Gough  <bjg@network-theory.co.uk>
  56
  57         * lambert.c (halley_iteration): increase tolerance to prevent
  58         exceeding max iters due to finite precision
  59
  60 2007-04-04  Brian Gough  <bjg@network-theory.co.uk>
  61
  62         * laguerre.c (gsl_sf_laguerre_n_e): use recursion for a=0 in
  63         addition to a>0
  64
  65 2007-02-17  Brian Gough  <bjg@network-theory.co.uk>
  66
  67         * log.c (gsl_sf_log_e): removed HIDE_INLINE_STATIC
  68
  69         * exp.c (gsl_sf_exp_e): removed HIDE_INLINE_STATIC
  70
  71 2007-02-14  Brian Gough  <bjg@network-theory.co.uk>
  72
  73         * mathieu_charv.c: made solve_cubic static
  74
  75 2007-02-12  Brian Gough  <bjg@network-theory.co.uk>
  76
  77         * mathieu_charv.c (figi): ensure that e[ii] is set when
  78         e2[ii]==0.0 and there is no error condition, as per the original
  79         eispack routine
  80
  81 2007-02-09  Brian Gough  <bjg@network-theory.co.uk>
  82
  83         * ellint.c (gsl_sf_ellint_F_e): do argument reduction for phi>pi/2
  84         (gsl_sf_ellint_E_e): do argument reduction for phi>pi/2
  85         (gsl_sf_ellint_P_e): do argument reduction for phi>pi/2
  86         (gsl_sf_ellint_D_e): do argument reduction for phi>pi/2
  87         (gsl_sf_ellint_Dcomp_e): added complete D integral
  88         (gsl_sf_ellint_Pcomp_e): added complete P integral
  89
  90 2007-01-31  Brian Gough  <bjg@network-theory.co.uk>
  91
  92         * beta.c (gsl_sf_lnbeta_sgn_e): added to support calculations with
  93         negative a,b
  94         (gsl_sf_lnbeta_e): rewritten in terms of gsl_sf_lnbeta_sgn_e
  95         (gsl_sf_beta_e): handle negative a,b
  96
  97         * gamma.c (gsl_sf_lngamma_sgn_e): make error calculations an exact
  98         copy of gsl_sf_lngamma_e (these functions could be merged to avoid
  99         duplication)
 100
 101 2007-01-29  Brian Gough  <bjg@network-theory.co.uk>
 102
 103         * test_legendre.c (test_legendre): added extra test cases for
 104         underflow
 105
 106 2007-01-26  Brian Gough  <bjg@network-theory.co.uk>
 107
 108         * expint.c (expint_E2_impl): handle x==0.0 as a special case
 109         (expint_E2_impl): corrected error term
 110
 111         * gsl_sf_log.h: removed inline version of log
 112
 113         * gsl_sf_exp.h: removed inline version of exp
 114
 115 2007-01-23  Brian Gough  <bjg@network-theory.co.uk>
 116
 117         * hyperg_1F1.c (hyperg_1F1_1_series): increase accuracy by factor
 118         of 4 in sum, tighter convergence condition, increase error
 119         estimate to allow for accumulated roundoff
 120         (hyperg_1F1_1): use series when |x| > |b|
 121         (gsl_sf_hyperg_1F1_e): only use Kummer when |x| < 100 otherwise
 122         exponential takes extreme value
 123
 124         * hyperg.c (gsl_sf_hyperg_1F1_series_e): allow 10000 iterations in
 125         series to extend valid range
 126         (gsl_sf_hyperg_1F1_series_e): increase accuracy by factor of 4 in
 127         sum, tighter convergence condition
 128
 129 2007-01-19  Brian Gough  <bjg@network-theory.co.uk>
 130
 131         * laguerre.c (laguerre_large_n): use the second term in the
 132         asymptotic expansion from Slater, p.73.
 133
 134 2007-01-17  Brian Gough  <bjg@network-theory.co.uk>
 135
 136         * hyperg_1F1.c (hyperg_1F1_largebx): asymptotic expansion for
 137         large b and x, with |x|<|b| from Slater 4.3.7
 138         (hyperg_1F1_1): use new asymptotic expansion for |x|<|b|
 139         (hyperg_1F1_small_a_bgt0): use new asymptotic expansion for |x|<|b|
 140         (hyperg_1F1_renorm_b0): add neglected terms in expansion for AS13.3.7
 141
 142 2007-01-14  Brian Gough  <bjg@network-theory.co.uk>
 143
 144         * legendre_poly.c (gsl_sf_legendre_sphPlm_e): added explicit
 145         computation of error term to allow for case when final term is
 146         zero.
 147
 148 2007-01-12  Brian Gough  <bjg@network-theory.co.uk>
 149
 150         * trig.c (gsl_sf_angle_restrict_symm_err_e): compute edge cases
 151         more reliably, return NaN when total loss of precision
 152         (gsl_sf_angle_restrict_pos_err_e): as above
 153
 154         * legendre_poly.c (gsl_sf_legendre_Pl_e): improve error estimate
 155         for large l by including rounding error at each step of recurrence
 156
 157 2006-10-03  Brian Gough  <bjg@network-theory.co.uk>
 158
 159         * poch.c (gsl_sf_lnpoch_e, gsl_sf_lnpoch_sgn_e): corrected
 160         result->val to 0.0 for x==0, previously returned incorrect value
 161         1.0
 162
 163 2006-09-24  Brian Gough  <bjg@network-theory.co.uk>
 164
 165         * laguerre.c (laguerre_large_n): work with small angles to avoid
 166         cancellation error, computer angular reduction exactly for integer
 167         eta.
 168
 169 2006-09-22  Brian Gough  <bjg@network-theory.co.uk>
 170
 171         * zeta.c (gsl_sf_zeta_e): make sin_term exactly zero for negative
 172         even integers
 173         (gsl_sf_zetam1_int_e): return value is -1 for zetam1_int with
 174         negative even integers
 175
 176 2006-03-26  Brian Gough  <bjg@network-theory.co.uk>
 177
 178         * fermi_dirac.c (fd_neg): initialize s to zero (avoid spurious
 179         warning from compiler)
 180
 181 2006-02-23  Brian Gough  <bjg@network-theory.co.uk>
 182
 183         * coulomb.c (gsl_sf_coulomb_wave_FG_e): fixed sign of F_lam_min,
 184         covers case when k_lam_G is nonzero and F_lam_min and F_lam differ
 185         in sign.
 186
 187 2006-01-21  Brian Gough  <bjg@network-theory.co.uk>
 188
 189         * synchrotron.c (gsl_sf_synchrotron_1_e): added first order
 190         correction term for the taylor expansion
 191         (gsl_sf_synchrotron_2_e): added first order correction term for
 192         the taylor expansion
 193
 194 2006-01-20  Brian Gough  <bjg@network-theory.co.uk>
 195
 196         * bessel_j.c (gsl_sf_bessel_jl_e): limit the use of
 197         gsl_sf_bessel_Jnu_asympx_e to the range x>100*l*l to satisfy
 198         the requirement x>>l*l in the asymptotic expansion
 199
 200         * bessel_In.c (gsl_sf_bessel_In_scaled_e): limit the use of
 201         gsl_sf_bessel_I_CF1_ser to the region where the continued
 202         converges with the allowed 20000 terms (|x|<1e7)
 203
 204 2005-12-20  Brian Gough  <bjg@network-theory.co.uk>
 205
 206         * debye.c (gsl_sf_debye_5_e, gsl_sf_debye_6_e): added n=5,6
 207         (R. J. Mathar)
 208
 209 2005-11-15  Brian Gough  <bjg@network-theory.co.uk>
 210
 211         * dilog.c (dilog_xge0): fix calculation of error estimate
 212
 213 2005-08-04  Brian Gough  <bjg@network-theory.co.uk>
 214
 215         * coulomb.c (gsl_sf_coulomb_wave_sphF_array): fix warning about
 216         variable shadowing for k
 217
 218 2005-07-29  Brian Gough  <bjg@network-theory.co.uk>
 219
 220         * gamma_inc.c (gsl_sf_gamma_inc_Q_e): use continued fraction close
 221         to peak of the integrand x > a - sqrt(a)
 222
 223 2005-07-28  Brian Gough  <bjg@network-theory.co.uk>
 224
 225         * elljac.c (gsl_sf_elljac_e): use separate iterations to avoid
 226         division by zero, new algorithm from Bulirsch avoids inverse trig
 227         functions.
 228
 229 2005-07-03  Brian Gough  <bjg@network-theory.co.uk>
 230
 231         * legendre_poly.c (gsl_sf_legendre_sphPlm_e): compute exponential
 232         internally to avoid underflow error when calling gsl_sf_exp_err
 233
 234 2005-07-02  Brian Gough  <bjg@network-theory.co.uk>
 235
 236         * bessel_i.c (gsl_sf_bessel_il_scaled_array): handle x==0 as a
 237         special case
 238
 239         * bessel_k.c (gsl_sf_bessel_kl_scaled_array): added lmax==0 as a
 240         special case
 241
 242         * bessel_y.c (gsl_sf_bessel_yl_array): added lmax==0 as a special
 243         case
 244
 245         * transport.c (gsl_sf_transport_2_e): improve error estimate for
 246         small x
 247
 248 2005-06-26  Brian Gough  <bjg@network-theory.co.uk>
 249
 250         * test_sf.c (test_jac): added tests over the full period for
 251         elljac functions
 252
 253 2005-05-23  Brian Gough  <bjg@network-theory.co.uk>
 254
 255         * test_bessel.c (test_bessel): added test for steed(99,1,...)
 256
 257         * legendre_H3d.c (legendre_H3d_CF1): use hypot
 258         (legendre_H3d_CF1_ser): use hypot
 259         (gsl_sf_legendre_H3d_e): use hypot
 260         (gsl_sf_legendre_H3d_e): use hypot
 261         (gsl_sf_legendre_H3d_array): use hypot
 262
 263         * dilog.c (dilogc_unitdisk): use hypot
 264
 265         * coulomb.c (gsl_sf_coulomb_CL_array): use hypot
 266         (coulomb_F_recur): use hypot
 267         (coulomb_G_recur): use hypot
 268         (coulomb_jwkb): use hypot
 269         (gsl_sf_coulomb_wave_F_array): use hypot
 270         (gsl_sf_coulomb_wave_FG_array): use hypot
 271         (gsl_sf_coulomb_wave_FGp_array): use hypot
 272
 273         * bessel_j.c (gsl_sf_bessel_jl_steed_array): use hypot
 274
 275         * bessel.c (gsl_sf_bessel_Inu_scaled_asymp_unif_e): use hypot
 276         (gsl_sf_bessel_Knu_scaled_asymp_unif_e): use hypot
 277
 278 2005-03-02  Brian Gough  <bjg@network-theory.co.uk>
 279
 280         * coulomb_bound.c (gsl_sf_hydrogenicR_e): don't check for
 281         underflow when function is known to be zero (e.g. r=0 l>0 or at
 282         zeroes of the laguerre polynomial).
 283
 284 2004-12-29  Brian Gough  <bjg@network-theory.co.uk>
 285
 286         * dilog.c (gsl_sf_complex_dilog_e): use const consistently in
 287         arguments of declaration and definition
 288         (gsl_sf_complex_dilog_xy_e): as above
 289
 290 2004-12-26  Brian Gough  <bjg@network-theory.co.uk>
 291
 292         * gamma_inc.c (gamma_inc_D): improve error estimate for case of
 293         u=x/a to include cancellation errors and only use it when x < 0.5*a
 294         since the cancellation errors are significant for x/a ~ 1
 295
 296 2004-12-23  Brian Gough  <bjg@network-theory.co.uk>
 297
 298         * gsl_sf_coupling.h: fixed declaration to
 299         gsl_sf_coupling_6j_INCORRECT instead of
 300         gsl_sf_coupling_INCORRECT_6j
 301
 302 2004-11-12  Brian Gough  <bjg@network-theory.co.uk>
 303
 304         * psi.c (gsl_sf_psi_1): added missing function definition
 305         gsl_sf_psi_1
 306
 307 2004-10-11  Brian Gough  <bjg@network-theory.co.uk>
 308
 309         * expint.c (gsl_sf_expint_Ei_scaled): fixed call to incorrect
 310         function gsl_sf_expint_Ei_e
 311
 312 2004-06-03  Brian Gough  <bjg@network-theory.co.uk>
 313
 314         * psi.c (psi_n_xg0): added missing return type int
 315
 316 2003-11-30  Brian Gough  <bjg@network-theory.co.uk>
 317
 318         * gsl_sf_log.h: added missing headers for inline functions
 319
 320 2003-08-11  Brian Gough  <bjg@network-theory.co.uk>
 321
 322         * test_sf.c: added preprocessor definitions TEST_FACTOR and
 323         TEST_SIGMA to allow larger tolerances on released versions (to
 324         reduce the number of spurious bug reports).
 325
 326 2003-07-24  Brian Gough  <bjg@network-theory.co.uk>
 327
 328         * gamma.c (gsl_sf_choose_e): declare k as unsigned int instead of
 329         int
 330         (gsl_sf_gamma_e): avoid shadowed declaration of sgn
 331
 332         * bessel_sequence.c (gsl_sf_bessel_sequence_Jnu_e): declare i as
 333         size_t instead of int
 334
 335 2003-06-11  Brian Gough  <bjg@network-theory.co.uk>
 336
 337         * coupling.c (gsl_sf_coupling_6j_INCORRECT_e): fixed typo in test
 338         for two_jf < 0
 339         (gsl_sf_coupling_6j_e): moved the implementation from the
 340         INCORRECT function into the correct function
 341
 342 2003-06-05  Brian Gough  <bjg@network-theory.co.uk>
 343
 344         * test_sf.c (test_coupling): added some regression tests for
 345         coupling functions
 346
 347 2003-06-02  Brian Gough  <bjg@network-theory.co.uk>
 348
 349         * gamma_inc.c (gamma_inc_D): handle x<a, a>10 as a special case
 350
 351 2003-04-18  Brian Gough  <bjg@network-theory.co.uk>
 352
 353         * gamma.c (gsl_sf_gammainv_e): handle any singularities in
 354         gamma(x) up front and return zero directly
 355
 356 2003-04-12  Brian Gough  <bjg@network-theory.co.uk>
 357
 358         * psi.c: changed value of psi(1,1) in table psi_1_table to be
 359         positive ((-)^2 * 1! * zeta(2)), in accordance with Abramowitz &
 360         Stegun 6.4.2.
 361
 362 2003-04-08  G. Jungman
 363
 364     * erfc.c, gsl_sf_erf.h: added gsl_sf_hazard_e() and
 365     gsl_sf_hazard() functions
 366
 367 2003-04-07  G. Jungman
 368
 369     * gamma_inc.c, gsl_sf_gamma.h: added gsl_sf_gamma_inc_e() and
 370     gsl_sf_gamma_inc(), implmentations of the non-normalized incomplete
 371     gamma function
 372
 373 2003-03-09  Brian Gough  <bjg@network-theory.co.uk>
 374
 375         * gsl_sf_legendre.h: added missing const to formal parameter on
 376         gsl_sf_legendre_sphPlm_deriv_array
 377
 378 2003-01-25  Brian Gough  <brian.gough@network-theory.co.uk>
 379
 380         * test_gamma.c (test_gamma): added a test for
 381         gamma_inc_P(10,1e-16) (BUG#12)
 382
 383 Sat Sep  7 15:56:15 2002  Brian Gough  <bjg@network-theory.co.uk>
 384
 385         * test_sf.h (TEST_FACTOR): include an overall factor in the test
 386         tolerances, otherwise there are too many bug reports for minor
 387         failures.
 388
 389 Sun Sep  1 17:34:27 2002  Brian Gough  <bjg@network-theory.co.uk>
 390
 391         * test_legendre.c (test_legendre): increased tolerance on test
 392
 393 Sat Jul 13 23:11:37 2002  Brian Gough  <bjg@network-theory.co.uk>
 394
 395         * ellint.c (gsl_sf_ellint_Kcomp_e): improved error estimate to
 396         take cancellation of y=1-k^2 into account near k=1
 397
 398 Sun Jul  7 21:39:12 2002  Brian Gough  <bjg@network-theory.co.uk>
 399
 400         * test_bessel.c (test_bessel): increased tolerance on test of
 401         bessel_Jn_array
 402
 403 Tue May 28 21:04:29 2002  Brian Gough  <bjg@network-theory.co.uk>
 404
 405         * psi.c (gsl_sf_psi_1piy_e): changed log(y) to log(ay) since
 406         function is even and can be extended to negative y
 407
 408 Mon Jan 28 19:02:42 2002  Brian Gough  <bjg@network-theory.co.uk>
 409
 410         * gamma_inc.c (gamma_inc_Q_CF): express 'small' constant in terms
 411         of GSL_DBL_EPSILON
 412         (gamma_inc_Q_CF): patch for more reliable series from Hans Plesser
 413
 414 Sat Jan 26 17:33:29 2002  Brian Gough  <bjg@network-theory.co.uk>
 415
 416         * test_gamma.c (test_gamma): increased tolerance on a test
 417
 418         * test_hyperg.c (test_hyperg): increased tolerance on a couple of
 419         tests
 420
 421 Fri Jan 18 15:12:30 2002  Brian Gough  <bjg@network-theory.co.uk>
 422
 423         * test_airy.c (test_airy): increased tolerance on test
 424
 425 Tue Jan  8 14:31:04 2002  Brian Gough  <bjg@network-theory.co.uk>
 426
 427         * test_legendre.c (test_legendre): increased tolerance by one
 428         level on _array tests
 429
 430         * hyperg_1F1.c (hyperg_1F1_small_a_bgt0): fix branch for a==1
 431
 432 Fri Dec  7 12:38:40 2001  Brian Gough  <bjg@network-theory.co.uk>
 433
 434         * laguerre.c (laguerre_n_cp): catch internal error, not global
 435         error
 436
 437         * error.h (INTERNAL_OVERFLOW_ERROR): added internal error macros
 438         which do not call the error handler.
 439
 440 Wed Dec  5 19:25:26 2001  Brian Gough  <bjg@network-theory.co.uk>
 441
 442         * gamma.c (gamma_xgthalf): return gamma(x) exactly when x is an
 443         integer
 444
 445         * hyperg_1F1.c (hyperg_1F1_ab_posint):  fix bug in recurrence
 446         initialisation, as done for hyperg_1F1_ab_pos
 447
 448 Thu Oct 18 11:37:25 2001  Brian Gough  <bjg@network-theory.co.uk>
 449
 450         * coulomb.c (gsl_sf_coulomb_CL_array): renamed from
 451         gsl_sf_coulomb_CL_list for consistency
 452
 453 Sat Oct 13 15:55:56 2001  Brian Gough  <bjg@network-theory.co.uk>
 454
 455         * cheb_eval.c (cheb_eval_e): keep track of cancellation error to
 456         prevent underestimates of total error.
 457
 458 Fri Oct 12 16:39:27 2001  Brian Gough  <bjg@network-theory.co.uk>
 459
 460         * test_hyperg.c (test_hyperg): increased tolerance from TOL0 to
 461         TOL1 on one test to allow it to pass under different optimizations
 462
 463 Thu Oct 11 14:21:34 2001  Brian Gough  <bjg@network-theory.co.uk>
 464
 465         * gegenbauer.c (gsl_sf_gegenpoly_n_e): initialize variable gk to
 466         zero to avoid warning
 467
 468         * bessel_i.c (gsl_sf_bessel_il_scaled_e): initialize variable
 469         iellm1 to zero to avoid warning
 470
 471         * bessel_Jnu.c (gsl_sf_bessel_Jnu_e): initialize variable Ynp1 to
 472         zero to avoid warning
 473
 474         * legendre_poly.c (gsl_sf_legendre_sphPlm_e): initialize variables
 475         p_ell, y_ell to zero to avoid warning
 476
 477 Thu Sep  6 10:38:51 2001  Brian Gough  <bjg@network-theory.co.uk>
 478
 479         * clausen.c: added missing initialiser to cheb_series struct
 480
 481         * pow_int.c (gsl_sf_pow_int_e): handle the case x == 0, n < 0
 482
 483         * legendre_poly.c (gsl_sf_legendre_array_size): added missing
 484         static version of this inline function
 485
 486 Wed Aug 15 20:18:43 2001  Brian Gough  <bjg@network-theory.co.uk>
 487
 488         * hyperg_1F1.c (hyperg_1F1_ab_pos): fix bug in recurrence
 489         initialisation
 490         (hyperg_1F1_ab_neg): if all else fails, try the series.
 491
 492 Wed Aug  8 19:55:34 2001  Brian Gough  <bjg@network-theory.co.uk>
 493
 494         * test_sf.c (test_coupling): used analytic expressions for the
 495         exact values to problems with insufficient precision
 496
 497 Thu Jul 12 15:25:04 2001  Brian Gough  <bjg@network-theory.co.uk>
 498
 499         * test_legendre.c (test_legendre): corrected tests to use exact
 500         floating point arguments near singularities.
 501
 502         * legendre_Qn.c (gsl_sf_legendre_Q0_e): removed unnecessary error
 503         terms that I added
 504         (gsl_sf_legendre_Q1_e): removed unnecessary error terms that I
 505         added
 506
 507 Wed Jul 11 16:06:45 2001  Brian Gough  <bjg@network-theory.co.uk>
 508
 509         * error.h: moved domain, overflow and underflow errors into
 510         macros, return Nan for domain error, Inf for overflow.
 511
 512 Tue Jul 10 22:00:55 2001  Brian Gough  <bjg@network-theory.co.uk>
 513
 514         * erfc.c (gsl_sf_log_erfc_e): introduce an additional log(erfc)
 515         branch to prevent loss of accuracy for small x
 516
 517         * test_sf.c (test_zeta): increased tolerance on zeta function test
 518         to take into account differences in rounding
 519
 520 2001-07-10  Brian Gough  <bjg@hppav.local>
 521
 522         * legendre_Qn.c (gsl_sf_legendre_Q0_e): use (log1p(x) - log1p(-x))
 523         instead of log((1+x)/(1-x)), for accuracy
 524         (gsl_sf_legendre_Q1_e): use (log1p(x) - log1p(-x)) instead of
 525         log((1+x)/(1-x)), for accuracy
 526         (gsl_sf_legendre_Q0_e): improve error estimate near singular
 527         points
 528         (gsl_sf_legendre_Q1_e): improve error estimate near singular
 529         points
 530         (gsl_sf_legendre_Q0_e): fixed incorrect factor of 2 in asymptotic
 531         expansion
 532
 533         * check.h (CHECK_UNDERFLOW): provide macro for checking results
 534         for underflow
 535
 536 Fri Jul  6 20:16:19 2001  Brian Gough  <bjg@network-theory.co.uk>
 537
 538         * zeta.c (riemann_zeta1m_slt0): added a function to compute
 539         zeta(1-s) without loss of accuracy for s near zero, as used by the
 540         reflection formula.  This fixes a bug in the accuracy of results
 541         of zeta(-x) for small x, where loss of precision previously
 542         occurred by evaluating 1-x.
 543
 544 2001-07-06  Brian Gough  <bjg@network-theory.co.uk>
 545
 546         * bessel_I0.c (gsl_sf_bessel_I0_scaled_e): added missing scaling
 547         factor of exp(-x) for case x<2sqrt(epsilon)
 548
 549 Thu Jul  5 16:16:13 2001  Brian Gough  <bjg@network-theory.co.uk>
 550
 551         * erfc.c (gsl_sf_erfc_e): corrected case of -10<x<-5 to take into
 552         account that chebyshev fit is defined in terms of positive x.
 553
 554 Tue Jul  3 12:09:33 2001  Brian Gough  <bjg@network-theory.co.uk>
 555
 556         * eval.h (EVAL_DOUBLE): avoid returning the status value as a
 557         numerical result
 558
 559 2001-06-28  Brian Gough  <bjg@network-theory.co.uk>
 560
 561         * elementary.c (gsl_sf_multiply_e): use GSL_COERCE_DBL macro to
 562         deal with extended register problem
 563
 564         * coupling.c: include stdlib.h for prototype of abs()
 565
 566 Wed Jun 27 21:20:22 2001  Brian Gough  <bjg@network-theory.co.uk>
 567
 568         * test_sf.c (test_fermidirac): increased tolerance on test of
 569         gsl_sf_fermi_dirac_int(9,500) by factor of two to allow for MSVC
 570
 571 Tue Jun 26 12:10:17 2001  Brian Gough  <bjg@network-theory.co.uk>
 572
 573         * gsl_sf_gamma.h: added const to prototype, to match function
 574         definition
 575
 576 Sun May  6 13:01:01 2001  Brian Gough  <bjg@network-theory.co.uk>
 577
 578         * test_sf.c: initialize message_buff to null string to prevent
 579         random segmentation faults
 580
 581         * test_sf.h: simplified test macros to reduce stack usage
 582
 583 Mon Apr 30 12:38:15 2001  Brian Gough  <bjg@network-theory.co.uk>
 584
 585         * airy_zero.c bessel_zero.c: zeros are now counted using an
 586         unsigned int
 587
 588         * poly.c: moved poly_eval function into poly/ directory
 589
 590 Wed Apr 25 17:28:48 2001  Brian Gough  <bjg@network-theory.co.uk>
 591
 592         * trig.c (gsl_sf_polar_to_rect): dropped _e from polar/rect
 593         conversion functions
 594
 595 Tue Apr 24 17:05:22 2001  Brian Gough  <bjg@network-theory.co.uk>
 596
 597         * split out chebyshev functions into their own cheb/ directory,
 598         leaving behind an internal cheb struct and cheb evaluation
 599         function.
 600
 601 Mon Apr 23 10:21:06 2001  Brian Gough  <bjg@network-theory.co.uk>
 602
 603         * changed tests for EFAULT into a commented-out macro, since
 604         EFAULT should only apply to invalid non-null pointers, but it
 605         might be useful to keep the tests around for debugging in this
 606         directory.
 607
 608         * unified error handling conventions to _e for error handling
 609         functions and no suffix for plain functions, so _impl functions
 610         are no longer needed.
 611
 612 1999-08-02  Mark Galassi  <rosalia@lanl.gov>
 613
 614         * fermi_dirac.c: took out the use of some "const int" constants
 615         which were being used to size arrays, since this is not portable
 616         (thanks to Bernd Petrovitsch <bernd@ict.tuwien.ac.at> for pointing
 617         this out).
 618
 619 1999-01-02  Mark Galassi  <rosalia@cygnus.com>
 620
 621         * trig.c (gsl_sf_rect_to_polar_impl): introduced an #ifdef
 622         HAVE_HYPOT for systems that do not have the hypot() routine.
 623
 624 Sun Feb 14 20:59:50 1999  Brian Gough  <bjg@netsci.freeserve.co.uk>
 625
 626         * Makefile.am (include_HEADERS): added gsl_sf_result.h
 627
 628 1998-12-05  Mark Galassi  <rosalia@cygnus.com>
 629
 630         * Makefile.am: t-erf.c was commented out, which is a problem
 631         (pointed out by automake-1.3b).  Moved the commented-out t-erf.c
 632         line to the end of the long list of files.
 633
 634         * legendre_poly.c (gsl_sf_legendre_Pl_impl):
 635         * legendre_con.c (gsl_sf_conicalP_xgt1_neg_mu_largetau_impl):
 636         added an ugly fix to the invocation of gsl_sf_bessel_Jn_impl() to
 637         give it two bogus arguments so that it builds.  I hope Jerry fixes
 638         it up soon!
 639
 640         * gsl_sf_airy.h: added an include of gsl_precision.h so that
 641         gsl_prec_t is defined.
 642
 643 Tue Nov 17 17:29:31 1998  Brian Gough  <bjg@vvv.lanl.gov>
 644
 645         * added #include <config.h> to all top-level source files
 646
 647         * chebyshev.c (gsl_sf_cheb_eval_n): changed local gslMIN and
 648         gslMAX to the standard GSL_MIN and GSL_MAX macros
 649
 650 Tue Aug 18 13:36:15 1998  Brian Gough  <bjg@vvv.lanl.gov>
 651
 652         * coulomb.c (C0sq): changed to using gsl_sf_expm1_impl instead of
 653         expm1, since the latter is a GNU extension.
 654
 655 Mon Aug  3 14:43:16 1998  Brian Gough  <bjg@vvv.lanl.gov>
 656
 657         * bessel_amp_phase.h: undefined consts are now declared extern
 658
 659 Mon Jul 13 12:40:27 1998  Brian Gough  <bjg@vvv.lanl.gov>
 660
 661         * gamma.c: changed all the factorial functions to take unsigned
 662         int arguments, since negative values are not allowed.
 663         (gsl_sf_choose): fixed off by one error in call to factorial.
 664
 665 Sun Jul 12 13:21:41 1998  Brian Gough  <bjg@vvv.lanl.gov>
 666
 667         * gsl_sf_legendre.h, gsl_sf_poly.h, gsl_sf_pow_int.h: added
 668         HAVE_INLINE around some inline definitions in the header files
 669
 670         * gamma.c: implemented the natural versions of gsl_sf_lnchoose and
 671         gsl_sf_choose.
 672
 673 Wed Apr 15 21:27:48 1998  Brian Gough  <bjg@vvv.lanl.gov>
 674
 675         * added the return code as the second argument of GSL_WARNING, so
 676         we can filter on the error number or output a standard message for
 677         each one.
 678
 679 1998-02-02  Mark Galassi  <rosalia@cygnus.com>
 680
 681         * Makefile.am (include_HEADERS): added gsl_specfunc.h so it gets
 682         into the distribution.
 683         (INCLUDES): added -I$(top_srcdir) so that bessel.c can find
 684         gsl_math.h