Dunno if this is a bug, but I was comparing the Complex Arithmetic routines in Section 5.5 of Numerical Recipes, and gsl_complex_div() uses simple complex division while 5.5.5 has a more sophisticated one that avoids underflow/overflow.
gsl_complex_abs() should match 5.5.4 if hypot() is implemented correctly on your machine, and gsl_complex_sqrt() matches 5.5.6 and 5.5.7. Best wishes, Jim Ward _______________________________________________ Bug-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/bug-gsl
