solver_mv_bernstein_binary_test.cpp

#include <realroot/solver_bernstein_binary.hpp>

using namespace mmx;
typedef polynomial< double, with<Bernstein> >         Bernstein;
typedef polynomial< double, with<Sparse,DegRevLex> >  Polynomial;
typedef solver<double, MvBernsteinBinaryApproximate>  Solver;
int main()
{
  Seq<Interval<double> > dmn;
  dmn<< Interval<double>(-0.1,1)
     << Interval<double>(-0.1,1);

  Seq<Polynomial> pols; 
  pols<< Polynomial("u^2+v^2-1")
      << Polynomial("u*v-0.4");
  
  std::cout << "Equations:\n  "<<pols<<std::endl;

  std::cout << "Solutions:\n  ";
  std::cout << Solver::solve( pols, dmn ).size() << std::endl;

    
};