El vie, 03-09-2004 a las 16:23, Brian Gough escribió: > 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. thanks for the reply, i'll have to study more :)