[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Bug-gsl] division by zero in gsl_sf_legendre_sphPlm_e
From: |
Brian Gough |
Subject: |
Re: [Bug-gsl] division by zero in gsl_sf_legendre_sphPlm_e |
Date: |
Wed, 31 Jan 2007 18:20:44 +0000 |
User-agent: |
Wanderlust/2.14.0 (Africa) Emacs/21.3 Mule/5.0 (SAKAKI) |
At Fri, 15 Sep 2006 14:00:06 +0200 (CEST),
address@hidden wrote:
> I attach you a little C program which calls the same function with the
> same parameters. Within the gsl_sf_legendre_sphPlm function there is
> indeed a division by zero, but my C compilers don't complain about it,
> producing an "inf" which does not harm the results.
>
> The division by zero takes place in legendre_poly.c, line 576:
>
> result->err += fabs(y_mm_err/y_mm) * fabs(y_ell);
>
This is now fixed in CVS. Thanks for the bug report.
--
Brian Gough
Network Theory Ltd,
Publishing Free Software Manuals --- http://www.network-theory.co.uk/
[Prev in Thread] |
Current Thread |
[Next in Thread] |
- Re: [Bug-gsl] division by zero in gsl_sf_legendre_sphPlm_e,
Brian Gough <=