polyroots-fortran
polyroots-fortran copied to clipboard
Modern Fortran library for finding the roots of real and complex polynomial equations
Digging through Netlib I've uncovered one more polynomial root finder: https://netlib.org/port/prop.upd/drpoly.f According to the comments it uses an adaptation of the Jenkins-Traub method. It is part of the [PORT Subroutine...
Look into the [eiscor](https://github.com/eiscor/eiscor) algorithm. It's a somewhat large library, so not sure if it makes sense to absorb all of that code into this library. Or maybe the code...
Should refactor all the routines so that: * dummy variable names are the same * all accept the coefficients in the same order (currently, some use increasing powers, other decreasing).
For cubics and quartics, a Fortran solver was published recently in ACM TOMS: > Flocke, N. (2015). Algorithm 954: An accurate and efficient cubic and quartic equation solver for physical...
I have extracted the polynomials from an [early issue](http://algol60.org/acm/003.pdf) of ACM TOMS (in case anyone is interested in the Algol code, a copy is available [here](https://github.com/JvanKatwijk/algol-60-compiler/blob/master/bairstow-e.alg)): > Grau, A. A....