solver_mv_bernstein_proj_test.cpp

#include <realroot/polynomial_bernstein.hpp>
#include <realroot/polynomial_sparse.hpp>
#include <realroot/solver_mv_bernstein.hpp>

using namespace mmx;
typedef polynomial< double, with<Bernstein> > Bernstein;
typedef polynomial< double, with<Sparse,DegRevLex> >  Polynomial;
typedef solver< double, ProjRd<SBD_RDRDL> >  Solver;

int main()
{
    std::vector<double> dmn(4);
    dmn[0] = -0.1; dmn[1] = 1;
    dmn[2] = -0.1; dmn[3] = 1;


    Polynomial a("u^2+v^2-1"), b("u*v-0.4");

    Seq<Polynomial> pols; pols<<a; pols<<b;

    std::cout << "Equations:\n"<<pols<<std::endl;

    std::cout << "Solutions:\n";
    Solver::Solutions sol = solve<double>( pols, ProjRd<SBD_RDRDL>(), dmn );
    std::cout<<sol<< std::endl;

    
};