Updated gsl to 2.2.

* What is new in gsl-2.2:

** updated gsl_linalg_cholesky_invert to use Level-2 BLAS
   and added function gsl_linalg_pcholesky_invert

** added function gsl_linalg_invtri for inverting
   triangular matrices

** fix GSL_EIGEN_SORT_VAL_{ASC,DESC} for nonsymmetric
   eigensystems (Victor Zverovich)

** added complete orthogonal decomposition routines
   (gsl_linalg_COD)

** bug fix where median calculation wasn't reset in
   gsl_rstat_reset(); added gsl_rstat_quantile_reset() function
   (reported by Pedro Donato)

** added multivariate Gaussian random distribution
   gsl_ran_multivariate_gaussian (Timothée Flutre)

** added functions to estimate the 1-norm reciprocal condition
   number for various matrix factorizations:
     * gsl_linalg_cholesky_rcond
     * gsl_linalg_QRPT_rcond

** added functions gsl_linalg_QRPT_{lssolve,lssolve2} to
   compute least squares solutions with the QRPT decomposition

** added function gsl_permute_matrix()

** added modified Cholesky factorization (gsl_linalg_mcholesky)
   to handle symmetric indefinite matrices

** added pivoted Cholesky factorization (gsl_linalg_pcholesky)
   for ill-conditioned matrices

** rewrote (real) Cholesky decomposition to use
   a Level-2 blas algorithm instead of Level-1. Flop
   count is about the same but the code is much simpler
   and easier to follow

** completely rewritten nonlinear least squares module,
   including support for large problems; the user may
   now control the linear solver used, the trust region
   updating strategy, and the scaling method. In addition,
   support has been added for the geodesic acceleration
   step (Transtrum 2011) which can speed up convergence
   on a wide class of problems.

** added gsl_rstat_rms() for root mean square

** optimized lmniel nonlinear least squares solver
   (bug #46369)

** improved precision in Bessel K0/K1 near x = 2
   (Pavel Holoborodko, bug #47401)

** added support for compressed row storage sparse
   matrices (Alexis Tantet)

** bug fix in convergence check of hypergeometric 2F1
   function (bug #45926)

** added gsl_multilarge_linear_lcurve() to compute
   the L-curve for large linear systems

** updated multilarge normal equations method to use
   new Cholesky scaling for better numerical stability

** added scaling to Cholesky routines to reduce the
   condition number prior to factorization

3 files changed, 25 insertions, 8 deletions

diff --git a/math/gsl/Makefile b/math/gsl/Makefile
index 8323e32d01d..fd44c848bdb 100644
--- a/math/gsl/Makefile
+++ b/math/gsl/Makefile
@@ -1,6 +1,6 @@
-# $NetBSD: Makefile,v 1.36 2015/11/17 10:42:12 wiz Exp $
+# $NetBSD: Makefile,v 1.37 2016/09/01 16:05:22 wiz Exp $
 
-DISTNAME=	gsl-2.1
+DISTNAME=	gsl-2.2
 CATEGORIES=	math devel
 MASTER_SITES=	${MASTER_SITE_GNU:=gsl/}
 
diff --git a/math/gsl/PLIST b/math/gsl/PLIST
index 16ca654f87f..1c83cf74554 100644
--- a/math/gsl/PLIST
+++ b/math/gsl/PLIST
@@ -1,4 +1,4 @@
-@comment $NetBSD: PLIST,v 1.20 2015/11/17 10:42:12 wiz Exp $
+@comment $NetBSD: PLIST,v 1.21 2016/09/01 16:05:22 wiz Exp $
 bin/gsl-config
 bin/gsl-histogram
 bin/gsl-randist
@@ -84,7 +84,9 @@ include/gsl/gsl_monte_plain.h
 include/gsl/gsl_monte_vegas.h
 include/gsl/gsl_multifit.h
 include/gsl/gsl_multifit_nlin.h
+include/gsl/gsl_multifit_nlinear.h
 include/gsl/gsl_multilarge.h
+include/gsl/gsl_multilarge_nlinear.h
 include/gsl/gsl_multimin.h
 include/gsl/gsl_multiroots.h
 include/gsl/gsl_multiset.h
@@ -103,6 +105,21 @@ include/gsl/gsl_permute_float.h
 include/gsl/gsl_permute_int.h
 include/gsl/gsl_permute_long.h
 include/gsl/gsl_permute_long_double.h
+include/gsl/gsl_permute_matrix.h
+include/gsl/gsl_permute_matrix_char.h
+include/gsl/gsl_permute_matrix_complex_double.h
+include/gsl/gsl_permute_matrix_complex_float.h
+include/gsl/gsl_permute_matrix_complex_long_double.h
+include/gsl/gsl_permute_matrix_double.h
+include/gsl/gsl_permute_matrix_float.h
+include/gsl/gsl_permute_matrix_int.h
+include/gsl/gsl_permute_matrix_long.h
+include/gsl/gsl_permute_matrix_long_double.h
+include/gsl/gsl_permute_matrix_short.h
+include/gsl/gsl_permute_matrix_uchar.h
+include/gsl/gsl_permute_matrix_uint.h
+include/gsl/gsl_permute_matrix_ulong.h
+include/gsl/gsl_permute_matrix_ushort.h
 include/gsl/gsl_permute_short.h
 include/gsl/gsl_permute_uchar.h
 include/gsl/gsl_permute_uint.h
diff --git a/math/gsl/distinfo b/math/gsl/distinfo
index 13e518aa946..9393b8fd402 100644
--- a/math/gsl/distinfo
+++ b/math/gsl/distinfo
@@ -1,9 +1,9 @@
-$NetBSD: distinfo,v 1.31 2015/11/17 10:42:12 wiz Exp $
+$NetBSD: distinfo,v 1.32 2016/09/01 16:05:22 wiz Exp $
 
-SHA1 (gsl-2.1.tar.gz) = e20110bf060c4030b26c519a70c6c0e113abc9f7
-RMD160 (gsl-2.1.tar.gz) = 4ccf5425d54c2b0ef76db97f5960563a4bf0e311
-SHA512 (gsl-2.1.tar.gz) = d63660fdfae9c9ec9f4fd1894b0af705297cefb6169973ef10faacd0b8a3d52f537508e9c07f4f9f4253138cc48e5570bb10ef303538226d7febd92e8fd184b0
-Size (gsl-2.1.tar.gz) = 3827265 bytes
+SHA1 (gsl-2.2.tar.gz) = b770a2761839689aaab858bb48c6dbc7f58f0269
+RMD160 (gsl-2.2.tar.gz) = 88d32ac6b2fdaa0e9647330aa3a09f2a29591b67
+SHA512 (gsl-2.2.tar.gz) = 588abc1ef07737f29449d081421fe59f5d57cf84b22ce5d25379f258b7ff7f74145d62cbc5ce273b6f3f80bdd5621766056928b87f58b0b806b0e36180e582ed
+Size (gsl-2.2.tar.gz) = 4473105 bytes
 SHA1 (patch-aa) = 3065e1282ded8b4a1b571f62603c3f49fbde8a09
 SHA1 (patch-ae) = df7c4f8863a1b8c0dc598d6fccd6243406b7ba0d
 SHA1 (patch-af) = 570f98483a5f05862da606c06e59e582048c6fcb