Closed emsr closed 5 years ago
Is the quadratic stepper used in Jenkins-Traub or others? Is this a factorization or could it be made into one? It would be cool/useful to make these solvers composable.
I put the Laguerre solver together. It can be stepped as well as run in one go.
I started thinking about quadratic_t{Tp b; Tp c}; OTOH I will want to hand this into the little quadratic solver that takes std::array.
Done.
The testing of these two components is successful. Quadratic needs a complex quadratic solver. Both need a full CRBond style driver.