Have the type of the epsilon parameter templated.
Refactor the quadratic solver.