RpolyPlusPlus
RpolyPlusPlus copied to clipboard
A modern c++ root finding algorithm based on the original Jenkins-Traub RPOLY software.
I have been testing the library on some polynomials with integer coefficients. On many of them incorrect roots are found. One of the examples: Coefficients: [1,0,-1,0,1,1,1,1,1,1,1,0,-1,0,1] Roots found: ``` 1.1910...
Another solution is `static constexpr` but I don't think it is compatible with previous C++ versions. `const double` should be enough, without hurting memory allocation.
Under some situations, I have found that the solver is outputting all zeros, below are some coefficients that caused it for me. I have tried putting some of these into...
Fix the build using g++ (GCC) 6.2.1 20160916 (Red Hat 6.2.1-2)