Re: Bug in gsl_poly_complex_solve() producing inaccurate roots of polynomials

[Brian Gough <bjg@network-theory.co.uk>]

At Mon, 10 Dec 2007 16:41:31 -0600 (CST),
salman wrote:

> I am using gsl_poly_complex_solve() to determine the roots of some basic 
> polynomials. The results I am getting are fairly off from the correct 
> roots. Here is the relevant section of the code (deg_L is the degree of 
> the polynomial I am considering):
> 
>    gsl_poly_complex_workspace *u=gsl_poly_complex_workspace_alloc(deg_L+1);
>    gsl_poly_complex_solve(coeffs, deg_L+1, u, roots);
>    gsl_poly_complex_workspace_free(u);
>    for(i=0; i<deg_L; i++)
>      printf("\t %.5f, %.5f\n", roots[i], roots[2*i+1]);
> 

Thanks for your email.  Please can you send a complete example program
so we can reproduce the problem.  Also, details of which version of
GSL you are using and operating system etc.

-- 
Brian Gough