[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Bug-gsl] gsl_complex_div
|
From: |
Brian Gough |
|
Subject: |
Re: [Bug-gsl] gsl_complex_div |
|
Date: |
Tue, 29 Mar 2011 23:45:18 +0100 |
|
User-agent: |
Wanderlust/2.15.6 (Almost Unreal) Emacs/23.2 Mule/6.0 (HANACHIRUSATO) |
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