Update GSL to version 1.15

New in gsl-1.15:

** Added Tuomo Keskitalo's new ode branch ode-initval2 with a
   gsl_odeiv2 prefix. This provides proper support for implicit
   solvers.  It is intended to be the new default for differential
   equations.  The existing gsl_odeiv routines will be retained for
   binary compatibility but their interface will be deprecated.

** Added new gsl_integrate_cquad routines for robust integration of
   difficult functions using the doubly-adaptive CQUAD algorithm
   (Pedro Gonnet).

** Added error checking to CBLAS functions (Jose Luis Garcia Pallero)

** Added a new function gsl_integration_glfixed_point to return
   ordered Gauss-Legendre points and weights contained within a
   gsl_integration_glfixed_table [bug #32237].

** Added a new function gsl_interp_type_min_size to return the size of
   an interpolation type.

** Added a function gsl_pow_uint(x,n) to compute x^n for unsigned
   exponents (needed when n exceeds the range of signed integers).

** Added new routine gsl_linalg_complex_cholesky_invert to handle the
   matrix inversion for complex Cholesky decompositions (Huan Wu).

** Added the functions gsl_vector_equal(x,y) and gsl_matrix_equal(x,y)
   for testing equality of two vectors or matrices.

** Added function gsl_eigen_nonsymmv_params to control the balancing
   transformation for eigenvector calculations. Balancing is now
   turned off by default for gsl_eigen_nonsymmv.

** It is now possible to choose an alternative cblas library via
   pkg-config using the GSL_CBLAS_LIB environment variable or
   the pkg-config --define-variable option.

** The jacobi method gsl_eigen_jacobi now uses the norm of the
   off-diagonal elements for its convergence criterion, as in
   algorithm 8.4.3 of Golub and van Loan.

** The newton multiroot solvers now return an error when a singular
   jacobian is detected.

** The interpolation functions now return NaN and when x is out of range,
   instead of extrapolating.

** The gsl_multimin_fdfsolver multidimensional minimisers now return
   GSL_ENOPROG immediately if the generated trial point does not
   differ from the initial point (to machine precision), avoiding
   unnecessary further iterations.

** Extended the range of gsl_sf_bessel_lnKnu_e by rescaling
   intermediate results to avoid internal overflows [bug #31528].

** Improved the result of gsl_sf_atanint_e for large arguments by
   adding the first order 1/x correction term. [bug #29562]

** Fixed the gsl_rng_ranlxs generators to enforce a maximum seed value
   of 2^31-1.  Larger seed values caused out of range values to be
   returned.

** Fixed gsl_ran_chisq_pdf(x,nu) to return correct result of 1/2
   instead of 0 when x=0 and nu=2, and +inf when x=0 and nu<2.

** Fixed gsl_pow_int(x,n) to avoid an infinite loop when n=INT_MIN due
   to wrapping of signed integers.

** Fixed gsl_sf_hyperg_2F1(a,b,c,x) to avoid returning NaN for
   arguments |a|>10. [bug #24812]

** Avoid spurious underflow return code in gsl_sf_beta_inc_e when
   intermediate underflow does not affect the result. [bug #30933]

** Avoid segfault in Chebyshev series derivatives gsl_cheb_calc_deriv
   for n=1. [bug #29139]

4 files changed, 16 insertions, 15 deletions

diff --git a/math/gsl/Makefile b/math/gsl/Makefile
index ca1bd641a41..fb86ab2fd8d 100644
--- a/math/gsl/Makefile
+++ b/math/gsl/Makefile
@@ -1,6 +1,6 @@
-# $NetBSD: Makefile,v 1.31 2010/03/15 10:45:28 wiz Exp $
+# $NetBSD: Makefile,v 1.32 2011/05/08 20:35:25 asau Exp $
 
-DISTNAME=	gsl-1.14
+DISTNAME=	gsl-1.15
 CATEGORIES=	math devel
 MASTER_SITES=	${MASTER_SITE_GNU:=gsl/}
 
diff --git a/math/gsl/PLIST b/math/gsl/PLIST
index e6274f0ebfa..9391c7220ec 100644
--- a/math/gsl/PLIST
+++ b/math/gsl/PLIST
@@ -1,4 +1,4 @@
-@comment $NetBSD: PLIST,v 1.17 2010/03/15 10:45:28 wiz Exp $
+@comment $NetBSD: PLIST,v 1.18 2011/05/08 20:35:25 asau Exp $
 bin/gsl-config
 bin/gsl-histogram
 bin/gsl-randist
@@ -89,6 +89,7 @@ include/gsl/gsl_multiset.h
 include/gsl/gsl_nan.h
 include/gsl/gsl_ntuple.h
 include/gsl/gsl_odeiv.h
+include/gsl/gsl_odeiv2.h
 include/gsl/gsl_permutation.h
 include/gsl/gsl_permute.h
 include/gsl/gsl_permute_char.h
diff --git a/math/gsl/distinfo b/math/gsl/distinfo
index 6c21011202f..779fad631fa 100644
--- a/math/gsl/distinfo
+++ b/math/gsl/distinfo
@@ -1,8 +1,8 @@
-$NetBSD: distinfo,v 1.24 2010/03/15 10:45:28 wiz Exp $
+$NetBSD: distinfo,v 1.25 2011/05/08 20:35:25 asau Exp $
 
-SHA1 (gsl-1.14.tar.gz) = e1a600e4fe359692e6f0e28b7e12a96681efbe52
-RMD160 (gsl-1.14.tar.gz) = c89a30a0b38a0c83d73ded19180d52acd56e6a52
-Size (gsl-1.14.tar.gz) = 3128269 bytes
-SHA1 (patch-aa) = 7166f1ca7ec25214e8985f0b012f9acaf1b59675
+SHA1 (gsl-1.15.tar.gz) = d914f84b39a5274b0a589d9b83a66f44cd17ca8e
+RMD160 (gsl-1.15.tar.gz) = f34aa2a6514228f7df42fc94638286bcce5cd94c
+Size (gsl-1.15.tar.gz) = 3281134 bytes
+SHA1 (patch-aa) = 5db4acc481ec04f9a129cc2da4acb4dd669f9a49
 SHA1 (patch-ae) = ae00f8c0e18377c66672bcd0099cc49f151b1a07
 SHA1 (patch-af) = 6c0fe8679cf511f9455630bbd42bfba1f293c35d
diff --git a/math/gsl/patches/patch-aa b/math/gsl/patches/patch-aa
index 4e977a57864..72931222e8f 100644
--- a/math/gsl/patches/patch-aa
+++ b/math/gsl/patches/patch-aa
@@ -1,13 +1,13 @@
-$NetBSD: patch-aa,v 1.9 2008/05/02 08:49:50 adam Exp $
+$NetBSD: patch-aa,v 1.10 2011/05/08 20:35:25 asau Exp $
 
---- configure.orig	2008-03-29 17:23:19.000000000 +0100
-+++ configure
-@@ -10715,7 +10715,7 @@ fi
+--- configure.orig	2011-05-05 01:10:35.000000000 +0400
++++ configure	2011-05-09 00:00:51.000000000 +0400
+@@ -11598,7 +11598,7 @@
  
  
  GSL_CFLAGS="-I$includedir"
 -GSL_LIBS="-L$libdir -lgsl"
 +GSL_LIBS="-Wl,-R$libdir -L$libdir -lgsl"
- 
- 
- 
+ LIBM=
+ case $host in
+ *-*-beos* | *-*-cegcc* | *-*-cygwin* | *-*-haiku* | *-*-pw32* | *-*-darwin*)