This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
Re: Possible bug in QUADPACK routine
- From: Brian Gough <bjg at network-theory dot co dot uk>
- To: Paul Barclay <pbarclay at caltech dot edu>
- Cc: gsl-discuss at sources dot redhat dot com
- Date: Sun, 10 Feb 2002 12:40:54 +0000 (GMT)
- Subject: Re: Possible bug in QUADPACK routine
- References: <Pine.GSO.4.42.0202092214430.20453-100000@inky>
Paul Barclay writes:
> I've just begun using GSL, and so far I've only tried out some of
> the quadrature routines (qgsl_integration_qags and
> gsl_integration_qagiu). Numerically they seem fine, but they have
> a nasty habit of causing a segmentation fault + core dump rather
> then return an error code whenever they are not able to do the
> integral.
Hi,
Thanks for your comments. The default error handling behavior is to
call abort() from gsl_error() but it is configurable -- see chapter 3
in the manual. If there is a core dump for any other reason please
send a bug report.
> One addition which would be very valuable I think are Bessel
> functions capable of taking complex arguments. Currently the
> fortran AMOS implementation available on Netlib is the only free
> code I've been able to find to do this.
There is a developers guide on the website at http://sources.redhat.com/gsl/
regards
--
Brian Gough
----------------------------------------------------------------------
Network Theory Ltd Phone: 0117 3179309 (+44 117 3179309)
15 Royal Park WWW: http://www.network-theory.co.uk/
Bristol BS8 3AL Email: bjg@network-theory.co.uk
United Kingdom
----------------------------------------------------------------------