Univariate polynomials |
Root finding is implemented:
analyziz/solver_aberth.hpp contains an implementation of Aberth's method,
analyziz/solver_floating.hpp uses the latter solver with gradual precision doubling,
analyziz/solver_ball.hpp includes certification with ball coefficients.
Mmx] |
use "analyziz" |
Mmx] |
p: Polynomial Floating == polynomial (-1.0001, 0.0^^4, 1.0) |
:
Mmx] |
bit_precision := 128; v == roots p |
:
Mmx] |
[ eval (p, a) | a in v ] |
:
Mmx] |
p: Polynomial Ball (Floating, Complex Floating) == polynomial (ball (-1.0001), (ball 0.0)^^4, ball 1.0) |
:
Mmx] |
roots p |
: