[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Bug-gsl] Bug in gsl_poly_complex_solve() producing inaccurate roots
From: |
Brian Gough |
Subject: |
Re: [Bug-gsl] Bug in gsl_poly_complex_solve() producing inaccurate roots of polynomials |
Date: |
Tue, 11 Dec 2007 08:40:34 +0000 |
User-agent: |
Wanderlust/2.14.0 (Africa) Emacs/22.1 Mule/5.0 (SAKAKI) |
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