This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
Re: little improvement
- From: Brian Gough <bjg at network-theory dot co dot uk>
- To: Fabian <fseoane at wanadoo dot es>
- Cc: gsl-discuss at sources dot redhat dot com
- Date: Fri, 3 Sep 2004 15:23:14 +0100
- Subject: Re: little improvement
- References: <1094044714.3480.12.camel@Eox>
Fabian writes:
> I just thought that if you set s to "1.0 / gsl_complex_abs2 (a)"
this will underflow if |a| > sqrt(DBL_MAX)
See volume 1 of the book "Numerical Computation" by C.Ueberhuber for a
good discussion of IEEE arithmetic.
--
Brian Gough
Network Theory Ltd,
Commercial support for GSL --- http://www.network-theory.co.uk/gsl/