[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/

reply via email to

[Prev in Thread] Current Thread [Next in Thread]