At Fri, 15 Sep 2006 14:00:06 +0200 (CEST), [EMAIL PROTECTED] 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/ _______________________________________________ Bug-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/bug-gsl
