Update to 1.9:

* What is new in gsl-1.9:

** Fixed the elliptic integrals F,E,P,D so that they have the correct
behavior for phi > pi/2 and phi < 0.  The angular argument is now
valid for all phi.  Also added the complete elliptic integral
gsl_sf_ellint_Pcomp.

** Added a new BFGS minimisation method gsl_multimin_fdfminimizer_vector_bfgs2
based on the algorithm given by R.Fletcher in "Practical Methods of
Optimisation" (Second edition).  This requires substantially fewer
function and gradient evaluations, and supercedes the existing BFGS
minimiser.

** The beta functions gsl_sf_beta_e(a,b) and gsl_sf_lnbeta_e(a,b) now
handle negative arguments a,b.  Added new function gsl_sf_lnbeta_sgn_e
for computing magnitude and sign of negative beta values, analagous to
gsl_sf_lngamma_sgn_e.

** gsl_cheb_eval_mode now uses the same error estimate as
gsl_cheb_eval_err.

** Improved gsl_sf_legendre_sphPlm_e to avoid underflow with large
arguments.

** Added updated Knuth generator, gsl_rng_knuthran2002, from 9th
printing of "The Art of Computer Programming".  Fixes various
weaknesses in the earlier version gsl_rng_knuthran.  See
http://www-cs-faculty.stanford.edu/~knuth/news02.htm

** The functions gsl_multifit_fsolver_set, gsl_multifit_fdfsolver_set
and gsl_multiroot_fsolver_set, gsl_multiroot_fdfsolver_set now have a
const qualifier for the input vector x, reflecting their actual usage.

** gsl_sf_expint_E2(x) now returns the correct value 1 for x==0,
instead of NaN.

** The gsl_ran_gamma function now uses the Marsaglia-Tsang fast gamma
method of gsl_ran_gamma_mt by default.

** The matrix and vector min/max functions now always propagate any
NaNs in their input.

** Prevented NaN occuring for extreme parameters in
gsl_cdf_fdist_{P,Q}inv and gsl_cdf_beta_{P,Q}inv

** Corrected error estimates for the angular reduction functions
gsl_sf_angle_restrict_symm_err and gsl_sf_angle_restrict_pos_err.
Fixed gsl_sf_angle_restrict_pos to avoid possibility of returning
small negative values.  Errors are now reported for out of range
negative arguments as well as positive.  These functions now return
NaN when there would be significant loss of precision.

** Corrected an error in the higher digits of M_PI_4 (this was beyond
the limit of double precision, so double precision results are not
affected).

** gsl_root_test_delta now always returns success if two iterates are
the same, x1==x0.

** A Japanese translation of the reference manual is now available
from the GSL webpage at http://www.gnu.org/software/gsl/ thanks to
Daisuke TOMINAGA.

** Added new functions for basis splines, see the "Basis Splines"
chapter in the GSL Reference Manual for details.

** Added new functions for testing the sign of vectors and matrices,
gsl_vector_ispos, gsl_vector_isneg, gsl_matrix_ispos and
gsl_matrix_isneg.

** Fixed a bug in gsl_sf_lnpoch_e and gsl_sf_lnpoch_sgn_e which caused
the incorrect value 1.0 instead of 0.0 to be returned for x==0.

** Fixed cancellation error in gsl_sf_laguerre_n for n > 1e7 so that
larger arguments can be calculated without loss of precision.

** Improved gsl_sf_zeta_e to return exactly zero for negative even
integers, avoiding less accurate trigonometric reduction.

** Fixed a bug in gsl_sf_zetam1_int_e where 0 was returned instead of
-1 for negative even integer arguments.

** When the differential equation solver gsl_odeiv_apply encounters a
singularity it returns the step-size which caused the error code from
the user-defined function, as opposed to leaving the step-size
unchanged.

** Added support for nonsymmetric eigensystems

** Added Mathieu functions

5 files changed, 10 insertions, 36 deletions

diff --git a/math/gsl/Makefile b/math/gsl/Makefile
index 9a164ac4c91..c5c3b21f6ff 100644
--- a/math/gsl/Makefile
+++ b/math/gsl/Makefile
@@ -1,6 +1,6 @@
-# $NetBSD: Makefile,v 1.25 2006/04/20 18:31:15 adam Exp $
+# $NetBSD: Makefile,v 1.26 2007/02/26 08:05:15 wiz Exp $
 
-DISTNAME=	gsl-1.8
+DISTNAME=	gsl-1.9
 CATEGORIES=	math devel
 MASTER_SITES=	${MASTER_SITE_GNU:=gsl/}
 
@@ -10,7 +10,7 @@ COMMENT=	The GNU Scientific Library
 
 USE_LIBTOOL=		yes
 GNU_CONFIGURE=		yes
-INFO_FILES=		# PLIST
+INFO_FILES=		yes
 
 TEST_TARGET=		check
 
diff --git a/math/gsl/PLIST b/math/gsl/PLIST
index 70e4c285905..877a376ec5c 100644
--- a/math/gsl/PLIST
+++ b/math/gsl/PLIST
@@ -1,4 +1,4 @@
-@comment $NetBSD: PLIST,v 1.13 2006/03/31 23:56:28 jlam Exp $
+@comment $NetBSD: PLIST,v 1.14 2007/02/26 08:05:15 wiz Exp $
 bin/gsl-config
 bin/gsl-histogram
 bin/gsl-randist
@@ -19,6 +19,7 @@ include/gsl/gsl_block_uchar.h
 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_cblas.h
 include/gsl/gsl_cdf.h
 include/gsl/gsl_chebyshev.h
@@ -146,6 +147,7 @@ include/gsl/gsl_sf_laguerre.h
 include/gsl/gsl_sf_lambert.h
 include/gsl/gsl_sf_legendre.h
 include/gsl/gsl_sf_log.h
+include/gsl/gsl_sf_mathieu.h
 include/gsl/gsl_sf_pow_int.h
 include/gsl/gsl_sf_psi.h
 include/gsl/gsl_sf_result.h
diff --git a/math/gsl/distinfo b/math/gsl/distinfo
index bbdaa86cc0e..34daf5a762d 100644
--- a/math/gsl/distinfo
+++ b/math/gsl/distinfo
@@ -1,11 +1,9 @@
-$NetBSD: distinfo,v 1.19 2007/01/01 14:32:29 rillig Exp $
+$NetBSD: distinfo,v 1.20 2007/02/26 08:05:16 wiz Exp $
 
-SHA1 (gsl-1.8.tar.gz) = 0dabdf2d024ae3b29207e485c5413676f3fa136f
-RMD160 (gsl-1.8.tar.gz) = f3e7f4622be56635f84c6d889bf8316dbac1edf4
-Size (gsl-1.8.tar.gz) = 2419825 bytes
+SHA1 (gsl-1.9.tar.gz) = 7541504ed9864b1a278f9392830ecac9e0146258
+RMD160 (gsl-1.9.tar.gz) = faca66d02876135c31c2f3585baa1931488e1593
+Size (gsl-1.9.tar.gz) = 2574939 bytes
 SHA1 (patch-aa) = 9803d533a5fab6667b3e1f21dbe21a1871ea0b3a
-SHA1 (patch-ab) = 90e974419437aaf3c4c0d50809ccec65f0a12a28
 SHA1 (patch-ac) = 2c11b0d3f89bb5764518124e1a56ecdf5dd29468
-SHA1 (patch-ad) = 1f0e7455c120a31ddbe00c3a05ca2d3d451566ea
 SHA1 (patch-ae) = 016b1b8911acf8fd20bd09d615890053a5807ec4
 SHA1 (patch-af) = e25d461cdc9400f5df017fb65aa395e1a447f9ca
diff --git a/math/gsl/patches/patch-ab b/math/gsl/patches/patch-ab
deleted file mode 100644
index a780bc3b53b..00000000000
--- a/math/gsl/patches/patch-ab
+++ /dev/null
@@ -1,13 +0,0 @@
-$NetBSD: patch-ab,v 1.5 2006/10/22 19:51:05 schwarz Exp $
-
---- ./randist/gausszig.c.orig	2006-10-22 10:54:48.000000000 +0200
-+++ ./randist/gausszig.c	2006-10-22 10:55:05.000000000 +0200
-@@ -164,7 +164,7 @@
- 
- 
- double
--gsl_ran_gaussian_ziggurat (const gsl_rng * r, double sigma)
-+gsl_ran_gaussian_ziggurat (const gsl_rng * r, const double sigma)
- {
-   unsigned long int i, j;
-   int sign;
diff --git a/math/gsl/patches/patch-ad b/math/gsl/patches/patch-ad
deleted file mode 100644
index 7d5e12e5ae4..00000000000
--- a/math/gsl/patches/patch-ad
+++ /dev/null
@@ -1,13 +0,0 @@
-$NetBSD: patch-ad,v 1.3 2006/10/22 19:51:05 schwarz Exp $
-
---- cdf/betainv.c.orig	2006-10-22 15:30:57.000000000 +0200
-+++ cdf/betainv.c	2006-10-22 15:31:37.000000000 +0200
-@@ -152,7 +152,7 @@
- }
- 
- double
--gsl_cdf_beta_Qinv (double Q, double a, double b)
-+gsl_cdf_beta_Qinv (const double Q, const double a, const double b)
- {
- 
-   if (Q < 0.0 || Q > 1.0)