* What is new in gsl-1.7:
** Switched gsl_randist_binomial to use the faster binomial random
variate TPE algorithm by default. The previous binomial variate
algorithm is available as gsl_randist_binomial_knuth. This will
result in a different sequence of binomial variates in programs using
this function.
** Improved the algorithm for gsl_sf_elljac_e to avoid cancellation
errors near quarter periods.
** Fixed the branch selection in gsl_sf_gamma_inc_Q_e to avoid
inaccurate results for large a,x where x~=~a.
** The multilinear fitting functions now have forms which accept a
user-specified tolerance for the SVD cutoff and return the
corresponding effective rank of the design matrix.
** The quadratic solvers in poly/ now handle linear equations
gracefully (i.e. quadratrics with a leading coefficient of zero).
** The output of "make check" now only shows test failures by default,
to reduce the amount of output. Set the environment variable
GSL_TEST_VERBOSE=1 to display all the output. To assist debugging,
the test number of each failure is shown in square brackets at the
line-end [NNNN].
** Fixed bugs in gsl_linalg_SV_decomp_jacobi which caused
incorrect results for some input matrices.
** Bessel, coulomb, dilogarithm and legendre_H3d functions now use
hypot internally to avoid overflow when computing terms like
sqrt(1+x*x).
** The 'Usage' chapter of the reference manual now explains how to
handle deprecated functions using the GSL_DISABLE_DEPRECATED macro.
** The conflicting enum definitions for 'forward' and 'backward' in
gsl_ftt.h and gsl_wavelet.h are deprecated. User code should switch
to the new definitions gsl_fft_forward, gsl_fft_backward,
gsl_wavelet_forward and gsl_wavelet_backward. Selectively define
GSL_DISABLE_DEPRECATED before including the headers to use the new
definitions on either or both modules.
** Fixed an error in the the brent minimisation algorithm. Iterations
should now follow Brent's original description correctly.
** The bound coulomb function gsl_sf_hydrogenicR_e no longer reports
an underflow for exact zeroes of the wavefunction.
** gsl_linalg_SV_decomp_jacobi now reports an error for the
unimplemented case M<N correctly.
** Fixed conformance test for the SYRK and HERK blas functions
gsl_blas_{s,d,c,z}syrk and gsl_blas_{c,z}herk for non-square matrices.
** Configure now checks for presence of ieeefp.h if needed.
** Differential equation solvers now propagate error codes returned
from user-defined functions to the top-level in all cases.
** Sort functions now avoid an infinite loop if Nans are present in
the input vector. The order of nans in the output is undefined,
although other elements will be sorted correctly.
7c48924cc4