Hi, I was looking at the code for gsl_sum. I have not completely understood it yet, but I noticed that whereas levin_utrunc.c explicitly tests for terms that are zero, this is not the case for levin_u.c. As one divides by such terms, it seems that one may potentially divide by zero this way. Am I wrong?
The documentation does not seem to mention that terms should be different from zero. Regards, Michel _______________________________________________ Help-gsl mailing list Help-gsl@gnu.org http://lists.gnu.org/mailman/listinfo/help-gsl
