At Wed, 16 Mar 2011 08:03:13 -0400, James Ward wrote: > 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.
I agree this should be improved, looking at the references I'm guessing it comes from Knuth. I've added a note to the TODO file. Any patch and/or test cases would be appreciated. -- Brian Gough _______________________________________________ Bug-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/bug-gsl
