GNU Scientific Library 1.16 released

[Patrick Alken <patrick.alken@Colorado.EDU>]

Version 1.16 of the GNU Scientific Library (GSL) is now available. GSL 
provides a large collection of routines for numerical computing in C.

This release includes significant new features and bug fixes:

   - nonlinear least squares fitting can now compute the Jacobian 
numerically with finite differences if the user doesn't wish to specify 
an analytic Jacobian

   - support for robust linear regression

   - Spearman correlation coefficient

   - numerous bug fixes

The full NEWS file entry is appended below.

The file details for this release are:

ftp://ftp.gnu.org/gnu/gsl/gsl-1.16.tar.gz
ftp://ftp.gnu.org/gnu/gsl/gsl-1.16.tar.gz.sig

The GSL project homepage is http://www.gnu.org/software/gsl/

GSL is free software distributed under the GNU General Public License.

Thanks to everyone who reported bugs and contributed improvements.

Patrick Alken

=====================================
* What is new in gsl-1.16:

** fixed error in gsl_rng_fwrite where uninitialized padding
    bytes were being written (bug #39104)

** fixed error in gsl_block_alloc where padding bytes were not
    properly initialized (bugs #39101,#39102,#39103)

** fixed error in ntuple/test.c where padding bytes were not
    properly initialized (bug #39105)

** fixed triangle selection bug in gsl_sf_coupling_6j_e and
    gsl_sf_coupling_9j_e (bugs #39466 and #29606) (HÃ¥kan Johansson and
    Alexey Illarionov)

** added higher level wrapper routine gsl_multifit_fdfsolver_driver

** converted gsl_multifit_linear_residuals to use dgemv to improve
    efficiency (bug #39153)

** added functions gsl_stats_spearman and gsl_sort_vector2 to compute
    Spearman rank correlation

** added function gsl_poly_dd_hermite_init for Hermite interpolation

** Added support for robust linear least squares

** Added function gsl_linalg_SV_leverage for computing statistical
    leverages from SVD decomposition

** Added support for approximating the Jacobian of nonlinear least
    squares fits using forward finite differences

** Extended gsl_sf_coupling_3j to allow larger range and to handle
    the special case (ja jb jc; 0 0 0)=0 when ja+jb+jc is odd

** Fixed gsl_sf_mathieu_se_array to return zero when the order is zero
    [bug #33679].

** Fixed overflow in gsl_sf_lncosh for large negative x (x<-354).

** Improved gsl_ran_negative_binomial_pdf to avoid underflow/overflow
    for large arguments.

** Multisets now allow k strictly greater than n.

** Fixed gsl_matrix_complex_fwrite/fread failure for noncontiguous
    matrices (Matthias Sitte).