polynomial_bernstein.hpp

Go to the documentation of this file.

#ifndef realroot_ring_bernstein_hpp
#define realroot_ring_bernstein_hpp

#include <realroot/Seq.hpp>
#include <realroot/tensor_bernstein.hpp>
#include <realroot/tensor_bernstein_fcts.hpp>
#include <realroot/monomial.hpp>
#include <realroot/polynomial.hpp>
#include <realroot/ring.hpp>

#define TMPL template<class C,class B>
#define RING ring<C,B>

namespace mmx {
  struct Bernstein {};

  //----------------------------------------------------------------------
  template<class C> 
  struct use<polynomial_of,polynomial<C,with<Bernstein> > >{
    
    typedef tensor::bernstein<C>            rep_t;
    typedef C                               Scalar;
    typedef polynomial<C, with<Bernstein> > Polynomial;
    typedef monom<C>                        Monomial;
    typedef ring<C, Bernstein >             Ring;

  };

  template<class C>
  struct ring<C, Bernstein > {
    typedef tensor::bernstein<C>            rep_t;
    typedef C                               Scalar;
    typedef polynomial<C, with<Bernstein> > Polynomial;
    typedef monom<C>                        Monomial;
    typedef ring<C, Bernstein >             Ring;

    typedef ring<C, Bernstein >             self_t;
    //    typedef self_t                          var_t;

    ring(const char* s) {var = variables(s); }

    static variables var;
    static variables& vars () { return var;}
    static int nbvar() { return var.nbvar(); }

    Polynomial  operator[](int i) const {
      return Polynomial((C)1,1,i);
    }
  };

  //--------------------------------------------------------------------
  template<class C> struct ring_of;
  template<class C> struct ring_of<polynomial<C,with<Bernstein> > > {
    typedef ring<C, Bernstein> Ring;
  } ;
  //--------------------------------------------------------------------
  template<class C> variables ring<C,Bernstein>::var;
  //--------------------------------------------------------------------
  TMPL inline Seq<typename ring<C,Bernstein>::Polynomial> 
  coefficients(const typename ring<C,Bernstein>::Polynomial& pol, int v ) {
      Seq<typename ring<C,Bernstein>::Polynomial> res;
      
      return res;
    }
    
  //====================================================================
} //namespace mmx
#undef TMPL
#undef RING
#undef Polynomial
#endif // realroot_polynom_mdse_hpp