ring_dual.hpp

Go to the documentation of this file.

#ifndef realroot_ring_dual_hpp
#define realroot_ring_dual_hpp

//====================================================================
#include <realroot/polynomial.hpp>
#include <realroot/sparse_dual.hpp>

#define TMPL template<class C, class O> 
#define RING       ring<C,Dual,O>
#define RING_POL   ring<C,Sparse,O>
#define POL        polynomial<C,with<Sparse,O> >
#define DUALPOL    polynomial<C,with<Dual,O> >
//====================================================================
namespace mmx {

  struct Dual {};

  template<class A, class B> struct use;
  template<class C, class V> struct polynomial;
  template<class C, class V> struct with;
 
  template<class C, class O> 
  struct use<polynomial_of,polynomial<C,with<Dual,O> > >{

    typedef sparse::dual<C,O>           rep_t;

    typedef C                           Scalar;
    typedef typename rep_t::Monomial    Monomial;
    typedef polynomial<C,with<Dual,O> > Polynomial;
    typedef ring<C,Dual,O>              Ring;

  };

  template<class C, class O> 
  struct ring<C, Dual, O > {

    typedef RING                        self_t;
    typedef sparse::dual<C,O>           rep_t;

    typedef C                           Scalar;
    typedef typename rep_t::Monomial    Monomial;
    typedef polynomial<C,with<Dual,O> > Polynomial;

    typedef RING                      Ring;
    
    static variables& vars () { return RING_POL::vars();}

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

    int nbvar() const { return  RING_POL::var.nbvar(); }
  };
  

  template<class C, class O> C sample(const RING&rg) { return (C)0;} 

  
  TMPL define_operator_rrr(DUALPOL, POL, DUALPOL, operator*, sparse::mul)
  TMPL define_operator_rrr(POL, DUALPOL, POL, operator*, sparse::mul)

  //----------------------------------------------------------------------
  template<class P, class Q> struct mul_helper;
  template<class C, class O> struct mul_helper<
    polynomial< C, with<Sparse,O> >,
    polynomial< C, with<Dual,  O> > > {
    typedef polynomial< C, with<Dual, O> > TYPE;
  };

  template<class C, class O> struct mul_helper<
    polynomial< C, with<Dual,  O> >,
    polynomial< C, with<Sparse,O> > > {
    typedef polynomial< C, with<Sparse,O> > TYPE;
  };
}
//====================================================================
#undef TMPL
#undef RING
#undef RING_POL
#undef POL
#undef DUALPOL
#undef Polynomial
//====================================================================
#endif //realroot_ring_dual.hpp