I am running GSL version 1.6 on a Macintosh Powerbook G4 (15 inch),
running OS 10.3. I am using gcc version 3.3.
Here's a short program that I called test.c that contains the bug:
I compiled it with:
gcc test.c -lgsl -lgslcblas -lm
***
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <string.h>
#include <gsl/gsl_roots.h>
#include <gsl/gsl_errno.h>
#include <gsl/gsl_math.h>
int main (int argc, char *argv[])
{
printf("Testing gsl_root_test_delta: %d\n", gsl_root_test_delta(0.,
0., 0., 1e-3));
}
***
Note that previously, I wrote:
Dear bug-gsl,
I think that I have found a bug in the function gsl_root_test_delta.
When x1 and x0 are both 0.0, gsl_root_test_delta returns -2, indicating
failure to converge.
For instance, the following line of code:
printf("Testing gsl_root_test_delta: %d\n", gsl_root_test_delta(0.,
0., 0., 1e-3));
results in the following output:
Testing gsl_root_test_delta: -2
Best,
Dave Spiegel
_______________________________________________
Bug-gsl mailing list
[email protected]
http://lists.gnu.org/mailman/listinfo/bug-gsl
