|
Packit |
67cb25 |
2011-10-10 Rhys Ulerich <rhys.ulerich@gmail.com>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* hyperg_U.c (gsl_sf_hyperg_U_int_e,gsl_sf_hyperg_U_e):
|
|
Packit |
67cb25 |
Always initialize the gsl_sf_result_e10 instances.
|
|
Packit |
67cb25 |
Thanks to Victor Zverovich for the bug report and patch.
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2011-10-10 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* coupling.c (gsl_sf_coupling_3j_e): compute 3j factors using logs
|
|
Packit |
67cb25 |
to allow larger range
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2011-09-20 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* coupling.c (triangle_selection_fails): extend triangle selection
|
|
Packit |
67cb25 |
to handle all permutations
|
|
Packit |
67cb25 |
(gsl_sf_coupling_3j_e): special case for (ja jb jc; 0 0 0) = 0
|
|
Packit |
67cb25 |
when ja+jb+jc=odd
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2011-08-10 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* mathieu_radfunc.c (gsl_sf_mathieu_Mc): set odd functions to zero
|
|
Packit |
67cb25 |
for order=0, initialise fn to zero before use
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2011-07-21 Brian Gough <bjg@gnu.org>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* mathieu_angfunc.c (gsl_sf_mathieu_se_array): handle the case
|
|
Packit |
67cb25 |
where q=0
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2011-07-15 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* trig.c (gsl_sf_lncosh_e): handle x symmetrically for middle
|
|
Packit |
67cb25 |
region
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2010-11-11 Brian Gough <bjg@gnu.org>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* ellint.c (gsl_sf_ellint_RC_e, gsl_sf_ellint_RD_e)
|
|
Packit |
67cb25 |
(gsl_sf_ellint_RF_e, gsl_sf_ellint_RJ_e): introduce a limit of
|
|
Packit |
67cb25 |
10000 iterations to avoid infinite loops
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* bessel_Knu.c (gsl_sf_bessel_Knu_scaled_e10_e): alternative
|
|
Packit |
67cb25 |
version of Knu_scaled function to allow greater range for
|
|
Packit |
67cb25 |
gsl_sf_bessel_lnKnu_e.
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2010-10-29 Brian Gough <bjg@gnu.org>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* hyperg_U.c (hyperg_U_small_a_bgt0): corrected result->err for
|
|
Packit |
67cb25 |
case where a==0.0
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2010-08-31 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* beta_inc.c (gsl_sf_beta_inc_e): ignore underflow error when term
|
|
Packit |
67cb25 |
is subtracted from 1
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2010-08-27 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* hyperg_2F1.c (gsl_sf_hyperg_2F1_e): use ap and bp consistently
|
|
Packit |
67cb25 |
in large c and large b cases, previously mixed a,b and ap,bp.
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2010-04-15 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* atanint.c (gsl_sf_atanint_e): added missing term 1/ax for large x
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2010-02-25 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* hyperg_U.c (hyperg_U_negx): handle the case where x<0
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2010-01-23 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* hyperg_1F1.c (gsl_sf_hyperg_1F1_e): use Kummer transformation
|
|
Packit |
67cb25 |
for larger range of x when b>a and a<0
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2009-07-17 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* hyperg_U.c (hyperg_U_series): use a rearrangement of the factor
|
|
Packit |
67cb25 |
in front of the series to avoid incorrect termination when the
|
|
Packit |
67cb25 |
the leading term is zero.
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2009-07-16 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* poch.c (gsl_sf_lnpoch_sgn_e): handle negative and zero values of
|
|
Packit |
67cb25 |
a, e.g. where gamma(a) and/or gamma(a+x) is infinite
|
|
Packit |
67cb25 |
(gsl_sf_poch_e): handle the case where lnpoch_sgn returns
|
|
Packit |
67cb25 |
GSL_NEGINF when the Pochhammer ratio is zero.
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* test_gamma.c (test_gamma): added tests for poch(a,x) with
|
|
Packit |
67cb25 |
negative arguments
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2009-07-12 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* hyperg_U.c (hyperg_U_int_bge1): avoid using the series when
|
|
Packit |
67cb25 |
1+a-b is zero or a negative integer.
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2009-07-11 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* hyperg_U.c (gsl_sf_hyperg_U_int_e10_e, gsl_sf_hyperg_U_e10_e)
|
|
Packit |
67cb25 |
(hyperg_U_origin, hyperg_U_int_origin): added special case for
|
|
Packit |
67cb25 |
U(a,b,z=0).
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2009-07-09 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* mathieu_workspace.c (gsl_sf_mathieu_free): handle NULL argument
|
|
Packit |
67cb25 |
in free
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2009-05-13 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* hyperg_2F1.c (gsl_sf_hyperg_2F1_e): fix condition on a and b
|
|
Packit |
67cb25 |
when c is a negative integer (either a or b must cause
|
|
Packit |
67cb25 |
cancellation of the series)
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2009-01-14 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* mathieu_workspace.c (gsl_sf_mathieu_alloc): increase number of
|
|
Packit |
67cb25 |
terms
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* mathieu_charv.c (gsl_sf_mathieu_a): increase number of terms
|
|
Packit |
67cb25 |
(gsl_sf_mathieu_b): increase number of terms
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2008-12-04 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* gamma_inc.c (gamma_inc_D): propagate cancellation error
|
|
Packit |
67cb25 |
in (x-a)/x for x close to a
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2008-12-03 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* exp.c (exprel_n_CF): changed N to double, to allow non-integer
|
|
Packit |
67cb25 |
usage for gamma_inc, double error factor to allow for two
|
|
Packit |
67cb25 |
parameters in recurrence.
|
|
Packit |
67cb25 |
(gsl_sf_exprel_n_CF_e): exported function to allow calls to
|
|
Packit |
67cb25 |
exprel_n_CF
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* gamma_inc.c (gamma_inc_P_series): improved convergence condition
|
|
Packit |
67cb25 |
using estimate of the remainder of the series, added continued
|
|
Packit |
67cb25 |
fraction as a fallback, increased nmax to 10000
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2008-08-26 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* ellint.c (gsl_sf_ellint_Kcomp_e): corrected taylor expansion for
|
|
Packit |
67cb25 |
k close to 1.
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2008-07-03 Brian Gough <bjg@hp2.network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* Makefile.am (INCLUDES): use top_srcdir instead of top_builddir
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2008-06-27 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* legendre_poly.c (gsl_sf_legendre_array_size): removed inline
|
|
Packit |
67cb25 |
version of this function in gsl_sf_legendre.h
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2008-03-17 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* hyperg_2F1.c (gsl_sf_hyperg_2F1_e): replace duplicate error
|
|
Packit |
67cb25 |
check on stat3 by stat4
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2008-03-15 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* bessel.c (gsl_sf_bessel_Jnu_asympx_e): use full asymptotic
|
|
Packit |
67cb25 |
series
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2008-03-14 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* bessel_j.c (gsl_sf_bessel_jl_e): increase error estimate by
|
|
Packit |
67cb25 |
factor of 2
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2008-02-09 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* exp.c (gsl_sf_exp_e10_e): only use exponent e10 when standard
|
|
Packit |
67cb25 |
exp() is out of range
|
|
Packit |
67cb25 |
(gsl_sf_exp_mult_e10_e): add full set of error terms
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2008-01-10 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* hyperg_2F1.c (gsl_sf_hyperg_2F1_e): handle case of x==1
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2007-10-25 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* bessel.c (gsl_sf_bessel_J_CF1): handle underflow by rescaling in
|
|
Packit |
67cb25 |
recurrence
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2007-09-13 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* ellint.c (gsl_sf_ellint_Pcomp_e): remove test for n <1, no
|
|
Packit |
67cb25 |
restriction needed
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2007-09-10 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* expint.c (gsl_sf_expint_En_e): added for support En(x)
|
|
Packit |
67cb25 |
(gsl_sf_expint_En_scaled_e): added for support En(x) scaled
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* gamma_inc.c (gamma_inc_CF): include finite precision of log term
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2007-08-30 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* psi.c (psi_complex_rhp): accumulate error, to allow for
|
|
Packit |
67cb25 |
cancellation effects
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* beta.c (gsl_sf_lnbeta_sgn_e): added missing factor of 2 for
|
|
Packit |
67cb25 |
error.
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2007-08-27 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* beta_inc.c (gsl_sf_beta_inc_e): handle cases where a<=0 or b<=0
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2007-04-27 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* lambert.c (halley_iteration): increase tolerance to prevent
|
|
Packit |
67cb25 |
exceeding max iters due to finite precision
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2007-04-04 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* laguerre.c (gsl_sf_laguerre_n_e): use recursion for a=0 in
|
|
Packit |
67cb25 |
addition to a>0
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2007-02-17 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* log.c (gsl_sf_log_e): removed HIDE_INLINE_STATIC
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* exp.c (gsl_sf_exp_e): removed HIDE_INLINE_STATIC
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2007-02-14 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* mathieu_charv.c: made solve_cubic static
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2007-02-12 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* mathieu_charv.c (figi): ensure that e[ii] is set when
|
|
Packit |
67cb25 |
e2[ii]==0.0 and there is no error condition, as per the original
|
|
Packit |
67cb25 |
eispack routine
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2007-02-09 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* ellint.c (gsl_sf_ellint_F_e): do argument reduction for phi>pi/2
|
|
Packit |
67cb25 |
(gsl_sf_ellint_E_e): do argument reduction for phi>pi/2
|
|
Packit |
67cb25 |
(gsl_sf_ellint_P_e): do argument reduction for phi>pi/2
|
|
Packit |
67cb25 |
(gsl_sf_ellint_D_e): do argument reduction for phi>pi/2
|
|
Packit |
67cb25 |
(gsl_sf_ellint_Dcomp_e): added complete D integral
|
|
Packit |
67cb25 |
(gsl_sf_ellint_Pcomp_e): added complete P integral
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2007-01-31 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* beta.c (gsl_sf_lnbeta_sgn_e): added to support calculations with
|
|
Packit |
67cb25 |
negative a,b
|
|
Packit |
67cb25 |
(gsl_sf_lnbeta_e): rewritten in terms of gsl_sf_lnbeta_sgn_e
|
|
Packit |
67cb25 |
(gsl_sf_beta_e): handle negative a,b
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* gamma.c (gsl_sf_lngamma_sgn_e): make error calculations an exact
|
|
Packit |
67cb25 |
copy of gsl_sf_lngamma_e (these functions could be merged to avoid
|
|
Packit |
67cb25 |
duplication)
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2007-01-29 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* test_legendre.c (test_legendre): added extra test cases for
|
|
Packit |
67cb25 |
underflow
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2007-01-26 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* expint.c (expint_E2_impl): handle x==0.0 as a special case
|
|
Packit |
67cb25 |
(expint_E2_impl): corrected error term
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* gsl_sf_log.h: removed inline version of log
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* gsl_sf_exp.h: removed inline version of exp
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2007-01-23 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* hyperg_1F1.c (hyperg_1F1_1_series): increase accuracy by factor
|
|
Packit |
67cb25 |
of 4 in sum, tighter convergence condition, increase error
|
|
Packit |
67cb25 |
estimate to allow for accumulated roundoff
|
|
Packit |
67cb25 |
(hyperg_1F1_1): use series when |x| > |b|
|
|
Packit |
67cb25 |
(gsl_sf_hyperg_1F1_e): only use Kummer when |x| < 100 otherwise
|
|
Packit |
67cb25 |
exponential takes extreme value
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* hyperg.c (gsl_sf_hyperg_1F1_series_e): allow 10000 iterations in
|
|
Packit |
67cb25 |
series to extend valid range
|
|
Packit |
67cb25 |
(gsl_sf_hyperg_1F1_series_e): increase accuracy by factor of 4 in
|
|
Packit |
67cb25 |
sum, tighter convergence condition
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2007-01-19 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* laguerre.c (laguerre_large_n): use the second term in the
|
|
Packit |
67cb25 |
asymptotic expansion from Slater, p.73.
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2007-01-17 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* hyperg_1F1.c (hyperg_1F1_largebx): asymptotic expansion for
|
|
Packit |
67cb25 |
large b and x, with |x|<|b| from Slater 4.3.7
|
|
Packit |
67cb25 |
(hyperg_1F1_1): use new asymptotic expansion for |x|<|b|
|
|
Packit |
67cb25 |
(hyperg_1F1_small_a_bgt0): use new asymptotic expansion for |x|<|b|
|
|
Packit |
67cb25 |
(hyperg_1F1_renorm_b0): add neglected terms in expansion for AS13.3.7
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2007-01-14 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* legendre_poly.c (gsl_sf_legendre_sphPlm_e): added explicit
|
|
Packit |
67cb25 |
computation of error term to allow for case when final term is
|
|
Packit |
67cb25 |
zero.
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2007-01-12 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* trig.c (gsl_sf_angle_restrict_symm_err_e): compute edge cases
|
|
Packit |
67cb25 |
more reliably, return NaN when total loss of precision
|
|
Packit |
67cb25 |
(gsl_sf_angle_restrict_pos_err_e): as above
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* legendre_poly.c (gsl_sf_legendre_Pl_e): improve error estimate
|
|
Packit |
67cb25 |
for large l by including rounding error at each step of recurrence
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2006-10-03 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* poch.c (gsl_sf_lnpoch_e, gsl_sf_lnpoch_sgn_e): corrected
|
|
Packit |
67cb25 |
result->val to 0.0 for x==0, previously returned incorrect value
|
|
Packit |
67cb25 |
1.0
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2006-09-24 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* laguerre.c (laguerre_large_n): work with small angles to avoid
|
|
Packit |
67cb25 |
cancellation error, computer angular reduction exactly for integer
|
|
Packit |
67cb25 |
eta.
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2006-09-22 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* zeta.c (gsl_sf_zeta_e): make sin_term exactly zero for negative
|
|
Packit |
67cb25 |
even integers
|
|
Packit |
67cb25 |
(gsl_sf_zetam1_int_e): return value is -1 for zetam1_int with
|
|
Packit |
67cb25 |
negative even integers
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2006-03-26 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* fermi_dirac.c (fd_neg): initialize s to zero (avoid spurious
|
|
Packit |
67cb25 |
warning from compiler)
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2006-02-23 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* coulomb.c (gsl_sf_coulomb_wave_FG_e): fixed sign of F_lam_min,
|
|
Packit |
67cb25 |
covers case when k_lam_G is nonzero and F_lam_min and F_lam differ
|
|
Packit |
67cb25 |
in sign.
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2006-01-21 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* synchrotron.c (gsl_sf_synchrotron_1_e): added first order
|
|
Packit |
67cb25 |
correction term for the taylor expansion
|
|
Packit |
67cb25 |
(gsl_sf_synchrotron_2_e): added first order correction term for
|
|
Packit |
67cb25 |
the taylor expansion
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2006-01-20 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* bessel_j.c (gsl_sf_bessel_jl_e): limit the use of
|
|
Packit |
67cb25 |
gsl_sf_bessel_Jnu_asympx_e to the range x>100*l*l to satisfy
|
|
Packit |
67cb25 |
the requirement x>>l*l in the asymptotic expansion
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* bessel_In.c (gsl_sf_bessel_In_scaled_e): limit the use of
|
|
Packit |
67cb25 |
gsl_sf_bessel_I_CF1_ser to the region where the continued
|
|
Packit |
67cb25 |
converges with the allowed 20000 terms (|x|<1e7)
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2005-12-20 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* debye.c (gsl_sf_debye_5_e, gsl_sf_debye_6_e): added n=5,6
|
|
Packit |
67cb25 |
(R. J. Mathar)
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2005-11-15 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* dilog.c (dilog_xge0): fix calculation of error estimate
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2005-08-04 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* coulomb.c (gsl_sf_coulomb_wave_sphF_array): fix warning about
|
|
Packit |
67cb25 |
variable shadowing for k
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2005-07-29 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* gamma_inc.c (gsl_sf_gamma_inc_Q_e): use continued fraction close
|
|
Packit |
67cb25 |
to peak of the integrand x > a - sqrt(a)
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2005-07-28 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* elljac.c (gsl_sf_elljac_e): use separate iterations to avoid
|
|
Packit |
67cb25 |
division by zero, new algorithm from Bulirsch avoids inverse trig
|
|
Packit |
67cb25 |
functions.
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2005-07-03 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* legendre_poly.c (gsl_sf_legendre_sphPlm_e): compute exponential
|
|
Packit |
67cb25 |
internally to avoid underflow error when calling gsl_sf_exp_err
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2005-07-02 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* bessel_i.c (gsl_sf_bessel_il_scaled_array): handle x==0 as a
|
|
Packit |
67cb25 |
special case
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* bessel_k.c (gsl_sf_bessel_kl_scaled_array): added lmax==0 as a
|
|
Packit |
67cb25 |
special case
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* bessel_y.c (gsl_sf_bessel_yl_array): added lmax==0 as a special
|
|
Packit |
67cb25 |
case
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* transport.c (gsl_sf_transport_2_e): improve error estimate for
|
|
Packit |
67cb25 |
small x
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2005-06-26 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* test_sf.c (test_jac): added tests over the full period for
|
|
Packit |
67cb25 |
elljac functions
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2005-05-23 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* test_bessel.c (test_bessel): added test for steed(99,1,...)
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* legendre_H3d.c (legendre_H3d_CF1): use hypot
|
|
Packit |
67cb25 |
(legendre_H3d_CF1_ser): use hypot
|
|
Packit |
67cb25 |
(gsl_sf_legendre_H3d_e): use hypot
|
|
Packit |
67cb25 |
(gsl_sf_legendre_H3d_e): use hypot
|
|
Packit |
67cb25 |
(gsl_sf_legendre_H3d_array): use hypot
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* dilog.c (dilogc_unitdisk): use hypot
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* coulomb.c (gsl_sf_coulomb_CL_array): use hypot
|
|
Packit |
67cb25 |
(coulomb_F_recur): use hypot
|
|
Packit |
67cb25 |
(coulomb_G_recur): use hypot
|
|
Packit |
67cb25 |
(coulomb_jwkb): use hypot
|
|
Packit |
67cb25 |
(gsl_sf_coulomb_wave_F_array): use hypot
|
|
Packit |
67cb25 |
(gsl_sf_coulomb_wave_FG_array): use hypot
|
|
Packit |
67cb25 |
(gsl_sf_coulomb_wave_FGp_array): use hypot
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* bessel_j.c (gsl_sf_bessel_jl_steed_array): use hypot
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* bessel.c (gsl_sf_bessel_Inu_scaled_asymp_unif_e): use hypot
|
|
Packit |
67cb25 |
(gsl_sf_bessel_Knu_scaled_asymp_unif_e): use hypot
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2005-03-02 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* coulomb_bound.c (gsl_sf_hydrogenicR_e): don't check for
|
|
Packit |
67cb25 |
underflow when function is known to be zero (e.g. r=0 l>0 or at
|
|
Packit |
67cb25 |
zeroes of the laguerre polynomial).
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2004-12-29 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* dilog.c (gsl_sf_complex_dilog_e): use const consistently in
|
|
Packit |
67cb25 |
arguments of declaration and definition
|
|
Packit |
67cb25 |
(gsl_sf_complex_dilog_xy_e): as above
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2004-12-26 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* gamma_inc.c (gamma_inc_D): improve error estimate for case of
|
|
Packit |
67cb25 |
u=x/a to include cancellation errors and only use it when x < 0.5*a
|
|
Packit |
67cb25 |
since the cancellation errors are significant for x/a ~ 1
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2004-12-23 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* gsl_sf_coupling.h: fixed declaration to
|
|
Packit |
67cb25 |
gsl_sf_coupling_6j_INCORRECT instead of
|
|
Packit |
67cb25 |
gsl_sf_coupling_INCORRECT_6j
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2004-11-12 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* psi.c (gsl_sf_psi_1): added missing function definition
|
|
Packit |
67cb25 |
gsl_sf_psi_1
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2004-10-11 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* expint.c (gsl_sf_expint_Ei_scaled): fixed call to incorrect
|
|
Packit |
67cb25 |
function gsl_sf_expint_Ei_e
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2004-06-03 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* psi.c (psi_n_xg0): added missing return type int
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2003-11-30 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* gsl_sf_log.h: added missing headers for inline functions
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2003-08-11 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* test_sf.c: added preprocessor definitions TEST_FACTOR and
|
|
Packit |
67cb25 |
TEST_SIGMA to allow larger tolerances on released versions (to
|
|
Packit |
67cb25 |
reduce the number of spurious bug reports).
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2003-07-24 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* gamma.c (gsl_sf_choose_e): declare k as unsigned int instead of
|
|
Packit |
67cb25 |
int
|
|
Packit |
67cb25 |
(gsl_sf_gamma_e): avoid shadowed declaration of sgn
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* bessel_sequence.c (gsl_sf_bessel_sequence_Jnu_e): declare i as
|
|
Packit |
67cb25 |
size_t instead of int
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2003-06-11 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* coupling.c (gsl_sf_coupling_6j_INCORRECT_e): fixed typo in test
|
|
Packit |
67cb25 |
for two_jf < 0
|
|
Packit |
67cb25 |
(gsl_sf_coupling_6j_e): moved the implementation from the
|
|
Packit |
67cb25 |
INCORRECT function into the correct function
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2003-06-05 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* test_sf.c (test_coupling): added some regression tests for
|
|
Packit |
67cb25 |
coupling functions
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2003-06-02 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* gamma_inc.c (gamma_inc_D): handle x<a, a>10 as a special case
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2003-04-18 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* gamma.c (gsl_sf_gammainv_e): handle any singularities in
|
|
Packit |
67cb25 |
gamma(x) up front and return zero directly
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2003-04-12 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* psi.c: changed value of psi(1,1) in table psi_1_table to be
|
|
Packit |
67cb25 |
positive ((-)^2 * 1! * zeta(2)), in accordance with Abramowitz &
|
|
Packit |
67cb25 |
Stegun 6.4.2.
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2003-04-08 G. Jungman
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* erfc.c, gsl_sf_erf.h: added gsl_sf_hazard_e() and
|
|
Packit |
67cb25 |
gsl_sf_hazard() functions
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2003-04-07 G. Jungman
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* gamma_inc.c, gsl_sf_gamma.h: added gsl_sf_gamma_inc_e() and
|
|
Packit |
67cb25 |
gsl_sf_gamma_inc(), implmentations of the non-normalized incomplete
|
|
Packit |
67cb25 |
gamma function
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2003-03-09 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* gsl_sf_legendre.h: added missing const to formal parameter on
|
|
Packit |
67cb25 |
gsl_sf_legendre_sphPlm_deriv_array
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2003-01-25 Brian Gough <brian.gough@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* test_gamma.c (test_gamma): added a test for
|
|
Packit |
67cb25 |
gamma_inc_P(10,1e-16) (BUG#12)
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Sat Sep 7 15:56:15 2002 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* test_sf.h (TEST_FACTOR): include an overall factor in the test
|
|
Packit |
67cb25 |
tolerances, otherwise there are too many bug reports for minor
|
|
Packit |
67cb25 |
failures.
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Sun Sep 1 17:34:27 2002 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* test_legendre.c (test_legendre): increased tolerance on test
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Sat Jul 13 23:11:37 2002 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* ellint.c (gsl_sf_ellint_Kcomp_e): improved error estimate to
|
|
Packit |
67cb25 |
take cancellation of y=1-k^2 into account near k=1
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Sun Jul 7 21:39:12 2002 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* test_bessel.c (test_bessel): increased tolerance on test of
|
|
Packit |
67cb25 |
bessel_Jn_array
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Tue May 28 21:04:29 2002 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* psi.c (gsl_sf_psi_1piy_e): changed log(y) to log(ay) since
|
|
Packit |
67cb25 |
function is even and can be extended to negative y
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Mon Jan 28 19:02:42 2002 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* gamma_inc.c (gamma_inc_Q_CF): express 'small' constant in terms
|
|
Packit |
67cb25 |
of GSL_DBL_EPSILON
|
|
Packit |
67cb25 |
(gamma_inc_Q_CF): patch for more reliable series from Hans Plesser
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Sat Jan 26 17:33:29 2002 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* test_gamma.c (test_gamma): increased tolerance on a test
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* test_hyperg.c (test_hyperg): increased tolerance on a couple of
|
|
Packit |
67cb25 |
tests
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Fri Jan 18 15:12:30 2002 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* test_airy.c (test_airy): increased tolerance on test
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Tue Jan 8 14:31:04 2002 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* test_legendre.c (test_legendre): increased tolerance by one
|
|
Packit |
67cb25 |
level on _array tests
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* hyperg_1F1.c (hyperg_1F1_small_a_bgt0): fix branch for a==1
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Fri Dec 7 12:38:40 2001 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* laguerre.c (laguerre_n_cp): catch internal error, not global
|
|
Packit |
67cb25 |
error
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* error.h (INTERNAL_OVERFLOW_ERROR): added internal error macros
|
|
Packit |
67cb25 |
which do not call the error handler.
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Wed Dec 5 19:25:26 2001 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* gamma.c (gamma_xgthalf): return gamma(x) exactly when x is an
|
|
Packit |
67cb25 |
integer
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* hyperg_1F1.c (hyperg_1F1_ab_posint): fix bug in recurrence
|
|
Packit |
67cb25 |
initialisation, as done for hyperg_1F1_ab_pos
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Thu Oct 18 11:37:25 2001 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* coulomb.c (gsl_sf_coulomb_CL_array): renamed from
|
|
Packit |
67cb25 |
gsl_sf_coulomb_CL_list for consistency
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Sat Oct 13 15:55:56 2001 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* cheb_eval.c (cheb_eval_e): keep track of cancellation error to
|
|
Packit |
67cb25 |
prevent underestimates of total error.
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Fri Oct 12 16:39:27 2001 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* test_hyperg.c (test_hyperg): increased tolerance from TOL0 to
|
|
Packit |
67cb25 |
TOL1 on one test to allow it to pass under different optimizations
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Thu Oct 11 14:21:34 2001 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* gegenbauer.c (gsl_sf_gegenpoly_n_e): initialize variable gk to
|
|
Packit |
67cb25 |
zero to avoid warning
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* bessel_i.c (gsl_sf_bessel_il_scaled_e): initialize variable
|
|
Packit |
67cb25 |
iellm1 to zero to avoid warning
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* bessel_Jnu.c (gsl_sf_bessel_Jnu_e): initialize variable Ynp1 to
|
|
Packit |
67cb25 |
zero to avoid warning
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* legendre_poly.c (gsl_sf_legendre_sphPlm_e): initialize variables
|
|
Packit |
67cb25 |
p_ell, y_ell to zero to avoid warning
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Thu Sep 6 10:38:51 2001 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* clausen.c: added missing initialiser to cheb_series struct
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* pow_int.c (gsl_sf_pow_int_e): handle the case x == 0, n < 0
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* legendre_poly.c (gsl_sf_legendre_array_size): added missing
|
|
Packit |
67cb25 |
static version of this inline function
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Wed Aug 15 20:18:43 2001 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* hyperg_1F1.c (hyperg_1F1_ab_pos): fix bug in recurrence
|
|
Packit |
67cb25 |
initialisation
|
|
Packit |
67cb25 |
(hyperg_1F1_ab_neg): if all else fails, try the series.
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Wed Aug 8 19:55:34 2001 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* test_sf.c (test_coupling): used analytic expressions for the
|
|
Packit |
67cb25 |
exact values to problems with insufficient precision
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Thu Jul 12 15:25:04 2001 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* test_legendre.c (test_legendre): corrected tests to use exact
|
|
Packit |
67cb25 |
floating point arguments near singularities.
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* legendre_Qn.c (gsl_sf_legendre_Q0_e): removed unnecessary error
|
|
Packit |
67cb25 |
terms that I added
|
|
Packit |
67cb25 |
(gsl_sf_legendre_Q1_e): removed unnecessary error terms that I
|
|
Packit |
67cb25 |
added
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Wed Jul 11 16:06:45 2001 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* error.h: moved domain, overflow and underflow errors into
|
|
Packit |
67cb25 |
macros, return Nan for domain error, Inf for overflow.
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Tue Jul 10 22:00:55 2001 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* erfc.c (gsl_sf_log_erfc_e): introduce an additional log(erfc)
|
|
Packit |
67cb25 |
branch to prevent loss of accuracy for small x
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* test_sf.c (test_zeta): increased tolerance on zeta function test
|
|
Packit |
67cb25 |
to take into account differences in rounding
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2001-07-10 Brian Gough <bjg@hppav.local>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* legendre_Qn.c (gsl_sf_legendre_Q0_e): use (log1p(x) - log1p(-x))
|
|
Packit |
67cb25 |
instead of log((1+x)/(1-x)), for accuracy
|
|
Packit |
67cb25 |
(gsl_sf_legendre_Q1_e): use (log1p(x) - log1p(-x)) instead of
|
|
Packit |
67cb25 |
log((1+x)/(1-x)), for accuracy
|
|
Packit |
67cb25 |
(gsl_sf_legendre_Q0_e): improve error estimate near singular
|
|
Packit |
67cb25 |
points
|
|
Packit |
67cb25 |
(gsl_sf_legendre_Q1_e): improve error estimate near singular
|
|
Packit |
67cb25 |
points
|
|
Packit |
67cb25 |
(gsl_sf_legendre_Q0_e): fixed incorrect factor of 2 in asymptotic
|
|
Packit |
67cb25 |
expansion
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* check.h (CHECK_UNDERFLOW): provide macro for checking results
|
|
Packit |
67cb25 |
for underflow
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Fri Jul 6 20:16:19 2001 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* zeta.c (riemann_zeta1m_slt0): added a function to compute
|
|
Packit |
67cb25 |
zeta(1-s) without loss of accuracy for s near zero, as used by the
|
|
Packit |
67cb25 |
reflection formula. This fixes a bug in the accuracy of results
|
|
Packit |
67cb25 |
of zeta(-x) for small x, where loss of precision previously
|
|
Packit |
67cb25 |
occurred by evaluating 1-x.
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2001-07-06 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* bessel_I0.c (gsl_sf_bessel_I0_scaled_e): added missing scaling
|
|
Packit |
67cb25 |
factor of exp(-x) for case x<2sqrt(epsilon)
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Thu Jul 5 16:16:13 2001 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* erfc.c (gsl_sf_erfc_e): corrected case of -10
|
|
Packit |
67cb25 |
account that chebyshev fit is defined in terms of positive x.
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Tue Jul 3 12:09:33 2001 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* eval.h (EVAL_DOUBLE): avoid returning the status value as a
|
|
Packit |
67cb25 |
numerical result
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
2001-06-28 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* elementary.c (gsl_sf_multiply_e): use GSL_COERCE_DBL macro to
|
|
Packit |
67cb25 |
deal with extended register problem
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* coupling.c: include stdlib.h for prototype of abs()
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Wed Jun 27 21:20:22 2001 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* test_sf.c (test_fermidirac): increased tolerance on test of
|
|
Packit |
67cb25 |
gsl_sf_fermi_dirac_int(9,500) by factor of two to allow for MSVC
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Tue Jun 26 12:10:17 2001 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* gsl_sf_gamma.h: added const to prototype, to match function
|
|
Packit |
67cb25 |
definition
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Sun May 6 13:01:01 2001 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* test_sf.c: initialize message_buff to null string to prevent
|
|
Packit |
67cb25 |
random segmentation faults
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* test_sf.h: simplified test macros to reduce stack usage
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Mon Apr 30 12:38:15 2001 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* airy_zero.c bessel_zero.c: zeros are now counted using an
|
|
Packit |
67cb25 |
unsigned int
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* poly.c: moved poly_eval function into poly/ directory
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Wed Apr 25 17:28:48 2001 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* trig.c (gsl_sf_polar_to_rect): dropped _e from polar/rect
|
|
Packit |
67cb25 |
conversion functions
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Tue Apr 24 17:05:22 2001 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* split out chebyshev functions into their own cheb/ directory,
|
|
Packit |
67cb25 |
leaving behind an internal cheb struct and cheb evaluation
|
|
Packit |
67cb25 |
function.
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Mon Apr 23 10:21:06 2001 Brian Gough <bjg@network-theory.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* changed tests for EFAULT into a commented-out macro, since
|
|
Packit |
67cb25 |
EFAULT should only apply to invalid non-null pointers, but it
|
|
Packit |
67cb25 |
might be useful to keep the tests around for debugging in this
|
|
Packit |
67cb25 |
directory.
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* unified error handling conventions to _e for error handling
|
|
Packit |
67cb25 |
functions and no suffix for plain functions, so _impl functions
|
|
Packit |
67cb25 |
are no longer needed.
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
1999-08-02 Mark Galassi <rosalia@lanl.gov>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* fermi_dirac.c: took out the use of some "const int" constants
|
|
Packit |
67cb25 |
which were being used to size arrays, since this is not portable
|
|
Packit |
67cb25 |
(thanks to Bernd Petrovitsch <bernd@ict.tuwien.ac.at> for pointing
|
|
Packit |
67cb25 |
this out).
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
1999-01-02 Mark Galassi <rosalia@cygnus.com>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* trig.c (gsl_sf_rect_to_polar_impl): introduced an #ifdef
|
|
Packit |
67cb25 |
HAVE_HYPOT for systems that do not have the hypot() routine.
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Sun Feb 14 20:59:50 1999 Brian Gough <bjg@netsci.freeserve.co.uk>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* Makefile.am (include_HEADERS): added gsl_sf_result.h
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
1998-12-05 Mark Galassi <rosalia@cygnus.com>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* Makefile.am: t-erf.c was commented out, which is a problem
|
|
Packit |
67cb25 |
(pointed out by automake-1.3b). Moved the commented-out t-erf.c
|
|
Packit |
67cb25 |
line to the end of the long list of files.
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* legendre_poly.c (gsl_sf_legendre_Pl_impl):
|
|
Packit |
67cb25 |
* legendre_con.c (gsl_sf_conicalP_xgt1_neg_mu_largetau_impl):
|
|
Packit |
67cb25 |
added an ugly fix to the invocation of gsl_sf_bessel_Jn_impl() to
|
|
Packit |
67cb25 |
give it two bogus arguments so that it builds. I hope Jerry fixes
|
|
Packit |
67cb25 |
it up soon!
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* gsl_sf_airy.h: added an include of gsl_precision.h so that
|
|
Packit |
67cb25 |
gsl_prec_t is defined.
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Tue Nov 17 17:29:31 1998 Brian Gough <bjg@vvv.lanl.gov>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* added #include <config.h> to all top-level source files
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* chebyshev.c (gsl_sf_cheb_eval_n): changed local gslMIN and
|
|
Packit |
67cb25 |
gslMAX to the standard GSL_MIN and GSL_MAX macros
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Tue Aug 18 13:36:15 1998 Brian Gough <bjg@vvv.lanl.gov>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* coulomb.c (C0sq): changed to using gsl_sf_expm1_impl instead of
|
|
Packit |
67cb25 |
expm1, since the latter is a GNU extension.
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Mon Aug 3 14:43:16 1998 Brian Gough <bjg@vvv.lanl.gov>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* bessel_amp_phase.h: undefined consts are now declared extern
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Mon Jul 13 12:40:27 1998 Brian Gough <bjg@vvv.lanl.gov>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* gamma.c: changed all the factorial functions to take unsigned
|
|
Packit |
67cb25 |
int arguments, since negative values are not allowed.
|
|
Packit |
67cb25 |
(gsl_sf_choose): fixed off by one error in call to factorial.
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Sun Jul 12 13:21:41 1998 Brian Gough <bjg@vvv.lanl.gov>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* gsl_sf_legendre.h, gsl_sf_poly.h, gsl_sf_pow_int.h: added
|
|
Packit |
67cb25 |
HAVE_INLINE around some inline definitions in the header files
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* gamma.c: implemented the natural versions of gsl_sf_lnchoose and
|
|
Packit |
67cb25 |
gsl_sf_choose.
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
Wed Apr 15 21:27:48 1998 Brian Gough <bjg@vvv.lanl.gov>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* added the return code as the second argument of GSL_WARNING, so
|
|
Packit |
67cb25 |
we can filter on the error number or output a standard message for
|
|
Packit |
67cb25 |
each one.
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
1998-02-02 Mark Galassi <rosalia@cygnus.com>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
* Makefile.am (include_HEADERS): added gsl_specfunc.h so it gets
|
|
Packit |
67cb25 |
into the distribution.
|
|
Packit |
67cb25 |
(INCLUDES): added -I$(top_srcdir) so that bessel.c can find
|
|
Packit |
67cb25 |
gsl_math.h
|