[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [Bug-gsl] Bug in polynomial root finder

From: Munagala Ramanath
Subject: Re: [Bug-gsl] Bug in polynomial root finder
Date: Sun, 6 Nov 2011 08:38:22 -0800 (PST)

Thanks for the response.

I have run it on more than 1.6 million polynomials of order 15 and this is
the only one on which it has problems.

I've also run it on some polynomials of order 18 and again it had no problems.

So it seems the order is not the problem, just some corner case with this
particular polynomial.


From: Brian Gough <address@hidden>
To: Munagala Ramanath <address@hidden>
Cc: "address@hidden" <address@hidden>
Sent: Friday, November 4, 2011 3:37 AM
Subject: Re: [Bug-gsl] Bug in polynomial root finder

At Sat, 29 Oct 2011 08:22:53 -0700 (PDT),
Munagala Ramanath wrote:
> # shows bug in GSL polynomial root finder
> p = [1, -4, 2, 10, -17, 10, 6, -16, 12, -16, 16, -8, 28, -8, -48, 32]
> roots = GSL::Poly[ p.reverse ].solve

Thanks for the email.  I think that is just a limitation of the solver
algorithm rather than a bug. The algorithm does not work for
arbitrarily high orders and 15th order is probably too high.

Brian Gough

reply via email to

[Prev in Thread] Current Thread [Next in Thread]