gmane.comp.lib.gsl.announce

GNU Scientific Library 1.15 released

Brian Gough — 2011-05-07T11:30:31

-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1

Version 1.15 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 contributions from the following
people:

  - Tuomo Keskitalo: new gsl_odeiv2 differential equation framework
    with improved implicit solvers (this supersedes the existing
    gsl_odeiv framework).

  - Pedro Gonnet: implementation of the CQUAD algorithm for robust
    integration of difficult functions.

  - José Luis García Pallero: error checking for the CBLAS library.

The full NEWS file entry is appended below.

The file details for this release are:

  http://ftpmirror.gnu.org/gsl/gsl-1.15.tar.gz      (3.2 MB)
  http://ftpmirror.gnu.org/gsl/gsl-1.15.tar.gz.sig  (GPG signature)
  d914f84b39a5274b0a589d9b83a66f44cd17ca8e          (sha1 checksum)

These URLs will redirect to a local mirror for faster downloading (the
original files are at ftp://ftp.gnu.org/gnu/gsl).

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.

Brian Gough
(GSL Maintainer)

======================================================================

* What is 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 (José Luis García 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]

-----BEGIN PGP SIGNATURE-----
Version: GnuPG v2.0.17 (GNU/Linux)

iJwEAQECAAYFAk3EgHEACgkQ9U2K2oCCH+qZsgQApinkw/p7BksqruZsbpMyrf8Q
uu6WJWhJszxsg4bGfmDonAc9opL3M9iRbxdJuNUPrxy/T85q9Fx+5n3rSWih2iaV
Na5LWfoXEnx8qi+6vipgNWWa6NVjkrSla5fmYb3s5lGFJ3IqVKl/WH/UMzNqW3DY
7BxK3i0/3aIuPT0xb3o=
=lfI5
-----END PGP SIGNATURE-----

_______________________________________________
Info-gsl mailing list
Info-gsl< at >gnu.org
https://lists.gnu.org/mailman/listinfo/info-gsl

[announce] GSL shell 1.0 released

Francesco Abbate — 2010-12-07T21:01:28

GSL shell 1.0 released

I'm glad to announce, after a lot of work, the final release of GSL
Shell 1.0. This new release brings, among the other things, a lot of
improvements in the graphics functions. Here some of the most
remarkable new features:

    * possibility to put multiple plots on the same windows
    * function to produce animations or interactive plots by using
graphical layers
    * function to save the plot in a file in bitmap format

The graphical rendering code was almost completely rewritten and
optimised for to be efficient both for static plot and animations. You
can find in the download directory also a window binary to try the
program easily but on linux it is easy to compile the program by
yourself. The project home page is http://www.nongnu.org/gsl-shell/ .

We have also a first implementation of a 3D plotting module based on
the JavaScript Pre3d library of Dean Mc Namee. This module give you
the possibility to create simple 3D plots and animations with a simple
interface. Just be aware that, While this module is quite functional
and usable, it could nevertheless replaced in future by a more
efficient implementation. The reason is not the quality of the Pre3d
code, which is excellent in itself, but the usage of JavaScript/Lua
for 3D graphics which is inherently inefficient in term of speed and
memory usage.

From the point of view of core mathematical functions we have also
some interesting new features:

    * new implementation of the interpolation functions with, notably,
the Akima and cubic spline interpolation
    * improvement of the handling of mixed operation between complex
and real matrix

There is also a new function to produce contour plot over a circular
domain. The previous contour functions was only working on rectangular
domains.

There is also a major change also the matrix operations semantic. Now
GSL shell automatically promote a matrix to complex when the other
operands are complex. This improvement includes the arithmetic
operations, the matrix multiplications functions, and the matrix
operations solve() and inv().

You can try the demo examples with the following commands:


or


In the "examples" folder there are a lot of interesting examples.

I hope you will find GSL Shell useful. In any case don't hesitate to
give me a feedback if you find a bug or if you want suggests some
improvements.

You can also contribute to the project by making a donation and by
expressing you preference about the features that should be
implemented in the next releases.

Enjoy
Francesco

release 0.10 of GSL Shell with new contour plotimplementation

Francesco Abbate — 2010-06-26T09:33:57

Hi all,

I'm glad to announce the release 0.10.0 of GSL Shell. This new release
does improve GSL Shell in several ways to make it more reliable and
simple to use.

Here a short list of the more important changes with this release of GSL
Shell:

* A contour plot algorithm is implemented. It is written entirely in
  Lua and it is quite fast and robust. The module name is "contour".

* An “high precision” experimental contour plot is also available. It
  is much more computational expensive but it draws accurate smooths
  curve and it does require the derivatives of the function. The
  module name is "hpcontour".

* To print an expression now you don’t need to write an ‘=’ sign
  before.

* Some serious bugs related to the graphical window system have been
  fixed. Now the plotting system appears to be reliable in all
  situations both on Windows and on X Window.

* More GSL modules implemented:
  - Basis splines
  - Linear Regression
  - Multi variables minimization
  - Eigensystems resolution

The project documentation website address is:

http://www.nongnu.org/gsl-shell/

In the download page you will find the source code, the windows
binaries and the documentation in HTML format.

You can test the modules by using the examples in the "examples"
directory. For example, to test the ODE module you can type:


You can give a look to the examples directory to understand better how
the different modules works. Some examples you may wish to try are:
'plot', 'graphics', 'fractals', 'minimize', 'contour', 'hpcontour'.

All the modules are documented in the online documentation page and
the documentation can also be downloaded in HTML format from the
download page.

The author has made a lot of efforts to ensure the quality of the
software and of the documentation but some bugs can be still present
and the documentation can incomplete or not entirely accurate in some
cases. You can help us by reporting any problem to the author.

We are also looking for volunteers so if you want to help don't hesitate
to contact the author.

Enjoy,
Francesco

PS: We have started a new mailing list about GSL Shell. If you want to
    subscribe here the adress:

    http://lists.nongnu.org/mailman/listinfo/gsl-shell-info

GNU Scientific Library 1.14 released

Brian Gough — 2010-03-11T18:42:51

-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1

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

This is a maintenance release which fixes reported bugs.  Support for
multisets and fixed order Gauss-Legendre integration has been added
and the build system has also been upgraded to the latest automake,
autoconf and libtool.  The full NEWS file entry is appended below.

The file details are:

  http://ftpmirror.gnu.org/gsl/gsl-1.14.tar.gz      (3 MB)
  http://ftpmirror.gnu.org/gsl/gsl-1.14.tar.gz.sig  (GPG signature)
  e1a600e4fe359692e6f0e28b7e12a96681efbe52          (sha1 checksum)

These URLs will redirect to a local mirror for faster downloading (the
original files are at ftp://ftp.gnu.org/gnu/gsl).

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.

Brian Gough
(GSL Maintainer)

======================================================================

* What is new in gsl-1.14:

** Upgraded to latest libtool, autoconf and automake (libtool-2.2.6b,
   autoconf-2.65, automake-1.11.1).  Fixes security hole in 'make
   dist' (see Automake CVE-2009-4029).

** Added support for "multisets", which are similar to permutations
   and combinations.  A multiset is an array of k integers in the
   range 0 to n-1 where each value may occur more than once. Multisets
   can be used to iterate over the indices of a k-th order symmetric
   tensor in n-space. (Rhys Ulerich)

** Added gsl_integrate_glfixed routines for performing fixed order
   Gauss-Legendre integration 
   (Pavel Holoborodko, Konstantin Holoborodko, Rhys Ulerich)

** In the LMDER multi-dimensional fitting routines, the return code
   GSL_ENOPROG is now used instead of GSL_CONTINUE for the case where
   10 or more attempts to find a suitable trial step have been made
   without success. [bug #25383]

** The confluent hypergeometric function gsl_sf_hyperg_U (a,b,x) has
   been extended to support x < 0 for cases not involving
   singularities in 1F1(a,b,x). [bug #27859]

** The F-distribution gsl_ran_fdist_pdf now avoids unnecessary
   underflow and overflow and handles cases where n > 248. [bug
   #28500]

** The SVD routines now use a scaled version of the implicit QR method
   and compute the shift more reliably for values at the limit of
   double precision, avoiding potential NaNs. [bug #28767]

** A compile time error is emitted if the configuration stage prevents
   the functions gsl_isnan and gsl_finite from being defined.

** Added missing dereference in GSL_VECTOR_COMPLEX when not using
   GSL_RANGE_CHECK [bug #28017]

** Improved the range of gsl_sf_hyperg1F1(a,b,x) when a<0,b>0. [bug
   #28718]

** Added macros GSL_MAJOR_VERSION and GSL_MINOR_VERSION in
   <gsl/gsl_version.h>

** Improved gsl_eigen_symmv and gsl_eigen_symm to avoid a potential
   infinite loop when the tridiagonalised matrix contains very small
   values. [bug #28096]

** Added functions gsl_multifit_linear_usvd and
   gsl_multifit_wlinear_usvd for multilinear fitting without
   column-scaling of the fit matrix.

-----BEGIN PGP SIGNATURE-----
Version: GnuPG v1.4.9 (GNU/Linux)

iJwEAQECAAYFAkuZOOYACgkQ9U2K2oCCH+rPdwP/ZzSBLn5RMl10WVK6C60t4Yl5
9+6IURU1J+c1fg3O0s+z1h2wOahl+NSf3KD9h4yZXbVg799fGIgGsPIK68Um7hWo
ge3vR7JwSVZFJls1CTXv5hgZpC1NhK+2zn1O8imkAEwKKoyHXBa5KUrLR+i78NoY
7YqpIGRhUpzqq91CCow=
=d74X
-----END PGP SIGNATURE-----

GSL extension ode-initval2-1.0 released

Tuomo Keskitalo — 2010-01-02T08:19:16

Hello,

FYI: I've released version 1.0 of ode-initval2 (ODE solver) extension 
library to GSL. Comments are welcome. It's available at

http://iki.fi/tuomo.keskitalo/gsl/ode-initval2/

Compared to GSL 1.13 ode-initval, ode-initval2 includes new Adams and 
BDF multistep methods (msadams and msbdf). Implicit steppers rk2imp and 
rk4imp have been modified to use Newton iteration, and implicit Euler 
method (rk1imp) has been added. Also, a new driver level has been added 
to simplify the use of the library.

The main change compared to previous version (ode-initval2-0.9) is in 
the interfaces. Now the stepper, evolve and control objects include a 
pointer to the driver object, through which they can communicate.

PS. I took down my ode-initval-additions git repository, since I have 
not updated it for a while.

Regards,
Tuomo

[ANN] new GSL Shell release with graphics module

Francesco Abbate — 2009-12-27T12:35:29

Hi all,

I'm proud to announce the release 0.9.5 of GSL Shell. This new release
introduce a new powerful and flexible module to produce graphical
plots and diagrams.

The new module is very flexible and let the user produce almost any
kind of graphics to represent data or functions. The most basic usage
is just to plot a function or multiple functions but the most advanced
users will be able to do almost anything by using the powerful
graphical primitives.

The module is based on the excellent Anti-Grain Geometry (AGG) library
by Maxim Shemanarev.

You can find at

https://savannah.nongnu.org/files/?group=gsl-shell

the source code as well as a windows executable. In order to test the
software you can type:


The basic documentation for the module is already available in the
website: http://www.nongnu.org/gsl-shell/graphics.html together with a
couple of examples http://www.nongnu.org/gsl-shell/examples.html .

The compilation with Windows is tricky and there is a bug in a mingw
library that will raise an error during the linkage. If you are
interested in compiling the code under windows please contact me
directly.

I'm sure that there are still a lot of problems and there are a lot of
features that should be added. As usual any suggestion or criticism is
welcome.

I'm also looking for volunteers so if you want to help don't hesitate
to contact the author.

Enjoy,
Francesco

GSL Shell release 0.6

Francesco Abbate — 2009-10-17T07:21:06

Hi all,

I'm glad to announce to the GSL users and developers the release of
GSL Shell 0.6 alpha. This new version features, among other things:

- a complete interface to perform Fast Fourier Transform of real and
complex data with a simple interface
- a complete interface to numerical integration routines, including
weighted functions and indefinite integrals

GSL Shell is current usable for many common tasks without any problem.
A lot of work is still needed to implement some GSL modules but the
GSL modules already implemented are fully functional and documented.

With GSL Shell you can very easily use the GSL routines with the full
power of Lua scripting language. By using Lua you are free to define
even complex functions to perform elaborate computations on structured
data and use, at the same time, the GSL routines for performing
numerical computations.

The documentation system used by GSL Shell is reStructured Text and
Sphinx, both are Python projects. With Sphinx you can produce high
quality HTML and Latex documentation by using a *very* simple textual
format.

For more information on GSL Shell you can visit the GSL shell website,
http://www.nongnu.org/gsl-shell/ or the GSL shell project page at
Savannah, https://savannah.nongnu.org/projects/gsl-shell/ .

You can also checkout the Subversion repository:

svn co svn://svn.sv.gnu.org/gsl-shell

I'm also looking for volunteers to further develop GSL Shell. A lot of
work has to be done to create the interface to many modules like
Random numbers, Eigenvalues, Simulated annealing, Special Functions
and many others. I need also help to improve and complete the
documentation and to develop a module to plot data.

Any feedback will be appreciated, thank you very much in any case for
your attention.

Best regards,
Francesco

GNU Scientific Library 1.13 released

Brian Gough — 2009-09-08T16:07:12

-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1

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

This is a maintenance release, which fixes reported bugs and upgrades
the build system to the latest autoconf and automake.  The full NEWS
file entry is appended below.

The file details are:

  ftp://ftp.gnu.org/gnu/gsl/gsl-1.13.tar.gz      (2.9 MB)
  ftp://ftp.gnu.org/gnu/gsl/gsl-1.13.tar.gz.sig  (GPG signature)
  02db78b9583bc7b2a577da6f45f5dd9f23ef737e       (sha1 checksum)

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.

Brian Gough
(GSL Maintainer)

======================================================================

* What is new in gsl-1.13:

** Upgraded to latest autoconf and automake (autoconf-2.64,
   automake-1.11)

** Fixed the rk4 and bspline allocators to avoid invalid free()
   calls under out of memory conditions. [bug #27194, #27236]

** Fixed a bug in gsl_multimin_fminimizer_nmsimplex2 where the center
   and size of the simplex were not updated on contract-by-best steps,
   causing failures in convergence. [bug #27180]

** Added new functions to set MISER and VEGAS Monte Carlo integration
   parameters, and to examine VEGAS chi-squared value and intermediate
   results.

** Added the function gsl_bspline_greville_abscissa to compute
   Greville abscissae for B-splines.

** The cumulative distribution functions gsl_cdf_gumbel1_{P,Q} should
   now handle a larger range of parameters without underflow and
   overflow.

** The header file gsl_const_cgs.h no longer defines values for
   electromagnetic units.  Applications should use gsl_const_cgsm.h
   instead to obtain the values in the CGS-Magnetic system.  The
   previous values for these units given in gsl_const_cgs.h were
   ill-defined as the type of CGS electromagnetic system was
   unspecified (the values were a mixture of CGS units with the Ampere
   of the MSKA system).  The affected constants are
   GSL_CONST_CGS_BOHR_MAGNETON, GSL_CONST_CGS_ELECTRON_CHARGE,
   GSL_CONST_CGS_ELECTRON_MAGNETIC_MOMENT, GSL_CONST_CGS_FARADAY,
   GSL_CONST_CGS_GAUSS, GSL_CONST_CGS_NUCLEAR_MAGNETON,
   GSL_CONST_CGS_PROTON_MAGNETIC_MOMENT, and GSL_CONST_CGS_ROENTGEN.

** The Pochhammer functions gsl_sf_poch(a,x) and gsl_sf_lnpoch(a,x) now
   handle the special cases where a and a+x are zero or negative
   integers.

** The confluent hypergeometric function gsl_sf_hyperg_U (a,b,x) now
   handles some cases where x=0. The case where 1+a-b is a negative
   integer no longer returns an error [bug #22859] and the incorrect
   termination of the series in certain cases is fixed [bug #26706].

** Added a new function gsl_poly_eval_derivs to evaluate a polynomial
   and its derivatives simultaneously.

** Added a new univariate minimisation algorithm
   gsl_min_fminimizer_quad_golden which is a variant of Brent's
   algorithm with safeguarded step-length adjustment.

** Added a new Nelder-Mead minimiser gsl_multimin_fminimizer_nmsimplex2rand
   which uses a randomly oriented simplex rather than one fixed on
   the coordinate axes [bug #25077]

** The texinfo file now uses the dircategory "Software libraries" from
   the Free Software Directory, as recommended in the Texinfo manual.

** The function gsl_ran_exponential now includes zero in its output
   range. [bug #25039]

** All functions for freeing allocated memory now accept a NULL
   pointer, following the standard C convention for free(). [bug
   #25319]

** The function gsl_sum_levin_u_accel now handles the special case
   c_0 + 0 + 0 + 0 + .... that occurs when summing power series
   c_n*x^n with x=0. [bug #26807]

** The functions gsl_linalg_LU_solve, gsl_linalg_LU_svx,
   gsl_linalg_LU_refine, gsl_linalg_LU_invert and their complex
   equivalents now return an error for singular matrices.

** The multifit LMDER hybrid solvers now check the return code of the
   user-supplied function in the gsl_multifit_fdfsolver_set
   method. [bug #26871]

** Improved the implementation of gsl_ran_discrete_preproc to avoid
   internal errors due to inconsistencies from excess precision on
   some platforms. [bug #26502]

** Corrected gsl_sf_hyperg_2F1(a,b,c,x) to not give a domain error in
   the case where c is a negative integer and the series terminates
   with a finite result.

** The C99 inline keyword is now supported, in addition to the
   previously supported GNU-style inline.

** Modified gsl_poly_complex_solve_cubic and gsl_poly_solve_cubic to
   avoid returning NaNs in cases where excess precision causes a
   change in the number of roots. 

** Fixed incorrect length check in gsl_blas_drotm. [bug #26503]

** Fixed gsl_odeiv_step_gear2 to restore y on step failure

** gsl_odeiv_evolve_apply now restores the correct value of t on step
   failures [bug #26255].

** Using make install prefix=DIR now puts correct paths in package
   config files gsl-config and gsl.pc

** Modified gsl_monte_vegas to work around pow() function inaccuracies
   on MinGW [bug #25413].

** Increased the number of terms in gsl_sf_mathieu_a and
   gsl_sf_mathieu_b to improve convergence in difficult regions [bug
   #25075]

-----BEGIN PGP SIGNATURE-----
Version: GnuPG v1.4.9 (GNU/Linux)

iJwEAQECAAYFAkqmebEACgkQ9U2K2oCCH+pjsQQAzr6z0BY8z/dM3tVHTczbEXo5
+q5Ayil/v7KpeThc58Sy7XeWoLgf0BcgccoZ97rt7YWr4sB1e9Um5N7vIZO2bNo3
4lb9wBr4kL1V/Ldd9QoTC6aIWeXYmVBQ3tGXtsOKl3LufW/QKQLjrBnMihRsF9uk
UH0a97c5JxhQjuSoiso=
=klpC
-----END PGP SIGNATURE-----