Pascal Molin on Fri, 24 May 2013 21:45:01 +0200 |
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
Re: Number of real roots of a polynomial |
A solution giving the number of distinct real roots would be equally useful for me. Speed is important in this application, so if there's a good way to do it with the library that would be fine (though I'm just using GP now).Usually I look at polroots(P) and count the number of roots which have an imaginary part which 'looks like' 0. But I'd like a solution that certifies that the zeros are not merely close to being real.Given a univariate polynomial, how can I determine the number of real solutions with gp?
In my case these are cubics in Z[x] so there are analytic criteria available but they seem a bit complicated -- putting the full problem into Mathematica and simplifying gave an _expression_ with LeafCount over 8000 -- that is, the number of operators and function calls was nearly 10^4. No doubt careful storing of partial results can reduce the algorithmic complexity but I don't expect it to be quick in any case (trusting here, somewhat unsteadily, to the reliability of Mathematica).
Charles Greathouse
Analyst/Programmer
Case Western Reserve University