Re: [Bug-gsl] division by zero in gsl_sf_legendre_sphPlm_e

[Brian Gough]

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/