solver_uv_contfrac_extended_test.cpp

#include <string>
#include <fstream>
#include <realroot/GMP.hpp>
#include <realroot/polynomial_tensor.hpp>
#include <realroot/solver_uv_continued_fraction.hpp>
#include <realroot/solver_run.hpp>

int main(int argc, char** argv)
{
    using namespace mmx;

    typedef mmx::GMP::extended_integer                            EZ;
    typedef mmx::GMP::integer                                     Z;
    typedef mmx::GMP::rational                                    Q;
    typedef mmx::GMP::extended_rational                           EQ;
    typedef mmx::polynomial<EZ, mmx::with<mmx::MonomialTensor> >  PEZ;

    // F[x_] := (-4 + 2*Sqrt[7])*x^2 + (3 + 4*Sqrt[7])*x + (5 + 6*Sqrt[7])*1
    // approx. roots
    // {x -> -8.64829}
    // {x -> -1.86892}

    Z root = Z(7);

    PEZ p = EZ(Z(-4), Z(2), root) * PEZ("x^2") +
    EZ(Z(3),  Z(4), root) * PEZ("x") +
    EZ(Z(5),  Z(6), root) * PEZ("1");

    mmx::Seq<mmx::Interval<EQ> > roots;
    mmx::solver<EQ, mmx::ContFrac<mmx::Isolate> >::solve(roots, p);

    std::cout << "num roots: " << roots.size() << std::endl;
    std::cout << roots << std::endl;

    return 0;
}