gsl: update to 2.6.

* What is new in gsl-2.6:

** add BLAS calls for the following functions:
     - gsl_vector_memcpy
     - gsl_vector_scale
     - gsl_matrix_memcpy
     - gsl_matrix_transpose_memcpy
     - gsl_matrix_tricpy
     - gsl_matrix_transpose_tricpy

** deprecated functions gsl_linalg_complex_householder_hm and
   gsl_linalg_complex_householder_mh

** add unit tests for gsl_linalg_symmtd and gsl_linalg_hermtd

** multilarge TSQR algorithm has been converted to use the new Level 3 QR decomposition

** nonlinear least squares Cholesky solver now uses the new Level 3 BLAS
   method; the old modified Cholesky solver is still available under
   gsl_multifit_nlinear_solver_mcholesky and gsl_multilarge_nlinear_solver_mcholesky

** implemented Level 3 BLAS versions of several linear algebra routines:
     - Triangular matrix inversion
     - Cholesky decomposition and inversion (real and complex)
     - LU decomposition and inversion (real and complex)
     - QR decomposition (courtesy of Julien Langou)
     - Generalized symmetric/hermitian eigensystem reduction to standard form

** removed deprecated function gsl_linalg_hessenberg()

** renamed gsl_interp2d_eval_e_extrap() to gsl_interp2d_eval_extrap_e()
   to match documentation (reported by D. Lebrun-Grandie)

** renamed some of the gsl_sf_hermite functions to be more consistent
   with rest of the library, and deprecated old function names

** updated gsl_sf_hermite_func() to use a newer algorithm
   due to B. Bunck which is more stable for large x; also added
   gsl_sf_hermite_func_fast() which uses the faster Cauchy integral
   algorithm in the same paper by Bunck

** add gsl_vector_axpby()

** add un-pivoted LDLT decomposition and its banded
   variant (gsl_linalg_ldlt_* and gsl_linalg_ldlt_band_*)

** add binary search tree module (gsl_bst); based on GNU libavl

** remove -u flag to gsl-histogram

** updated spmatrix module
   - added routines and data structures for all types (float,uint,char,...)
   - added gsl_spmatrix_scale_columns() and gsl_spmatrix_scale_rows()
   - added gsl_spmatrix_add_to_dense()
   - more efficient reallocation of COO/triplet matrices (no longer rebuilds binary tree)
   - enhanced test suite
   - added gsl_spmatrix_min_index()

** add routines for banded Cholesky decomposition (gsl_linalg_cholesky_band_*)

** documented gsl_linalg_LQ routines and added gsl_linalg_LQ_lssolve()

4 files changed, 28 insertions, 10 deletions

diff --git a/math/gsl/Makefile b/math/gsl/Makefile
index cccdbfeabc7..d30d827b45b 100644
--- a/math/gsl/Makefile
+++ b/math/gsl/Makefile
@@ -1,6 +1,6 @@
-# $NetBSD: Makefile,v 1.42 2020/01/18 23:32:51 rillig Exp $
+# $NetBSD: Makefile,v 1.43 2021/05/31 10:35:17 wiz Exp $
 
-DISTNAME=	gsl-2.5
+DISTNAME=	gsl-2.6
 CATEGORIES=	math devel
 MASTER_SITES=	${MASTER_SITE_GNU:=gsl/}
 
diff --git a/math/gsl/PLIST b/math/gsl/PLIST
index 49c8fd8411d..2c251b20045 100644
--- a/math/gsl/PLIST
+++ b/math/gsl/PLIST
@@ -1,4 +1,4 @@
-@comment $NetBSD: PLIST,v 1.23 2018/10/02 12:10:29 adam Exp $
+@comment $NetBSD: PLIST,v 1.24 2021/05/31 10:35:17 wiz Exp $
 bin/gsl-config
 bin/gsl-histogram
 bin/gsl-randist
@@ -20,6 +20,10 @@ include/gsl/gsl_block_uint.h
 include/gsl/gsl_block_ulong.h
 include/gsl/gsl_block_ushort.h
 include/gsl/gsl_bspline.h
+include/gsl/gsl_bst.h
+include/gsl/gsl_bst_avl.h
+include/gsl/gsl_bst_rb.h
+include/gsl/gsl_bst_types.h
 include/gsl/gsl_cblas.h
 include/gsl/gsl_cdf.h
 include/gsl/gsl_chebyshev.h
@@ -214,6 +218,20 @@ include/gsl/gsl_splinalg.h
 include/gsl/gsl_spline.h
 include/gsl/gsl_spline2d.h
 include/gsl/gsl_spmatrix.h
+include/gsl/gsl_spmatrix_char.h
+include/gsl/gsl_spmatrix_complex_double.h
+include/gsl/gsl_spmatrix_complex_float.h
+include/gsl/gsl_spmatrix_complex_long_double.h
+include/gsl/gsl_spmatrix_double.h
+include/gsl/gsl_spmatrix_float.h
+include/gsl/gsl_spmatrix_int.h
+include/gsl/gsl_spmatrix_long.h
+include/gsl/gsl_spmatrix_long_double.h
+include/gsl/gsl_spmatrix_short.h
+include/gsl/gsl_spmatrix_uchar.h
+include/gsl/gsl_spmatrix_uint.h
+include/gsl/gsl_spmatrix_ulong.h
+include/gsl/gsl_spmatrix_ushort.h
 include/gsl/gsl_statistics.h
 include/gsl/gsl_statistics_char.h
 include/gsl/gsl_statistics_double.h
diff --git a/math/gsl/distinfo b/math/gsl/distinfo
index 968e101dd5f..6b7585c051c 100644
--- a/math/gsl/distinfo
+++ b/math/gsl/distinfo
@@ -1,9 +1,9 @@
-$NetBSD: distinfo,v 1.36 2018/10/02 12:10:29 adam Exp $
+$NetBSD: distinfo,v 1.37 2021/05/31 10:35:17 wiz Exp $
 
-SHA1 (gsl-2.5.tar.gz) = 9b3a32ba4f4e79fc6164fdc62e56003ae6e991e5
-RMD160 (gsl-2.5.tar.gz) = 2ec06d425adf5992191ff1c0dd992f3bae3af064
-SHA512 (gsl-2.5.tar.gz) = 5b4c5c023f9029ce220f4e09aa4d0234fed94d42ac224d58fda095fe0532d54237a7c33278f8b5d0ba051f6004486edb38d0adb4fcb49337a8c1d8a18cf4a24a
-Size (gsl-2.5.tar.gz) = 7303881 bytes
+SHA1 (gsl-2.6.tar.gz) = 9273164b6bdf60d0577518a1c1310eff6659e3dd
+RMD160 (gsl-2.6.tar.gz) = ac5a8879e03ae93d60342df4a0d6082fab02170c
+SHA512 (gsl-2.6.tar.gz) = 0be8240715f0b86aba2c63d9f12da4dba4719d4e350e9308d279e0dd3b2f0519ea26fd2e38a17f3e8cf43aacbaa2455207a7ca0d6c305f3b8725e8ece2250a74
+Size (gsl-2.6.tar.gz) = 7398548 bytes
 SHA1 (patch-aa) = 3065e1282ded8b4a1b571f62603c3f49fbde8a09
 SHA1 (patch-ae) = df7c4f8863a1b8c0dc598d6fccd6243406b7ba0d
 SHA1 (patch-af) = 570f98483a5f05862da606c06e59e582048c6fcb
diff --git a/math/gsl/hacks.mk b/math/gsl/hacks.mk
index 74492607b68..1fbfcc1e6f4 100644
--- a/math/gsl/hacks.mk
+++ b/math/gsl/hacks.mk
@@ -1,4 +1,4 @@
-# $NetBSD: hacks.mk,v 1.1 2006/03/31 18:57:28 jlam Exp $
+# $NetBSD: hacks.mk,v 1.2 2021/05/31 10:35:17 wiz Exp $
 
 .include "../../mk/compiler.mk"
 
@@ -6,7 +6,7 @@
 ### At least gcc-3.3.5 failed to compile gsl-1.6/statistics/lag1.c
 ### with -mcpu=v9 on a solaris/sparc system.
 ###
-.if (${OPSYS} == "SunOS") && !empty(MACHINE_ARCH:Msparc) && \
+.if (${OPSYS} == "SunOS") && ${MACHINE_ARCH} == sparc && \
     !empty(CC_VERSION:M*gcc*)
 PKG_HACKS+=		solaris-sparc-ice
 BUILDLINK_TRANSFORM+=	rm:-mcpu=v9