/Users/mourrain/Devel/mmx/realroot/include/realroot/polynomial_operators.hpp

Go to the documentation of this file.

#ifndef realroot_polynom_operators_hpp
#define realroot_polynom_operators_hpp
//====================================================================
#include <realroot/texp_sup.hpp>
#include <realroot/texp_expression.hpp>
#include <realroot/texp_operators.hpp>

# define TMPL       template <class C, class Rep, class Ord>
# define TMPLX      template <class C, class Rep, class Ord, class X>
# define VARIANT    with<Rep,Ord>
# define POLYNOMIAL polynomial<C,VARIANT>
# define Scalar     C
//====================================================================
namespace mmx {

TMPL inline void 
add (POLYNOMIAL &r, const POLYNOMIAL &a ) { 
  add (r.rep (), a.rep () ); 
}
TMPL inline void 
add (POLYNOMIAL &r, const POLYNOMIAL &a, const POLYNOMIAL &b) {
  add (r.rep (), a.rep (), b.rep ()); 
}
TMPL inline void
add (POLYNOMIAL &r, const POLYNOMIAL &a, const C & b) {
  add (r.rep (), a.rep (), b); 
}
TMPL inline void
add (POLYNOMIAL &r, const C & a, const POLYNOMIAL &b) {
  add (r.rep(), b.rep(), a); 
}
TMPL inline void
sub (POLYNOMIAL &r, const POLYNOMIAL &a ) {
  sub (r.rep (), a.rep () ); 
}
TMPL inline void
sub (POLYNOMIAL &r, const POLYNOMIAL &a, const POLYNOMIAL &b) {
  sub (r.rep (), a.rep (), b.rep ()); 
}
TMPL inline void
sub (POLYNOMIAL &r, const C & a, const POLYNOMIAL &b) {
  sub (r.rep (), a, b.rep ()); 
}
TMPL inline void
sub (POLYNOMIAL &r, const POLYNOMIAL &a, const C & b) {
  sub (r.rep (), a.rep (), b); 
}
TMPL inline void
mul (POLYNOMIAL &r, const POLYNOMIAL &a ) {
  mul (r.rep (), a.rep () ); 
}
TMPL inline void
mul (POLYNOMIAL &r, const POLYNOMIAL &a, const POLYNOMIAL &b) {
  mul (r.rep (), a.rep (), b.rep ()); 
}
TMPL inline void
mul (POLYNOMIAL &r, const POLYNOMIAL &a, const C & b) {
  mul (r.rep (), a.rep (), b); }
TMPL inline void
mul (POLYNOMIAL &r, const C & a, const POLYNOMIAL &b) {
  mul (r.rep (), b.rep (), a); 
}
TMPL inline void
div (POLYNOMIAL &r, const POLYNOMIAL &a, const POLYNOMIAL &b) {
  div (r.rep (), a.rep (), b.rep ()); 
}
TMPL inline void
div (POLYNOMIAL &r, const POLYNOMIAL &a, const C & b) {
  div (r.rep (), a.rep (), b); 
}
TMPL inline void
rem (POLYNOMIAL &r, const POLYNOMIAL &a, const POLYNOMIAL &b) {
  rem (r.rep (), b.rep (), a.rep); 
}
//====================================================================
struct operators_of {};
TMPL struct use<operators_of,POLYNOMIAL> {
    static inline void 
    add (POLYNOMIAL &r, const POLYNOMIAL &a ) { 
      add (r.rep (), a.rep () ); 
    }
    static inline void 
    add (POLYNOMIAL &r, const POLYNOMIAL &a, const POLYNOMIAL &b) {
      add (r.rep (), a.rep (), b.rep ()); 
    }
    static inline void
    add (POLYNOMIAL &r, const POLYNOMIAL &a, const C & b) {
      add (r.rep (), a.rep (), b); 
    }
    static inline void
    add (POLYNOMIAL &r, const C & a, const POLYNOMIAL &b) {
      add (r.rep(), b.rep(), a); 
    }
    static inline void
    sub (POLYNOMIAL &r, const POLYNOMIAL &a ) {
      sub (r.rep (), a.rep () ); 
    }
    static inline void
    sub (POLYNOMIAL &r, const POLYNOMIAL &a, const POLYNOMIAL &b) {
      sub (r.rep (), a.rep (), b.rep ()); 
    }
    static inline void
    sub (POLYNOMIAL &r, const C & a, const POLYNOMIAL &b) {
      sub (r.rep (), a, b.rep ()); 
    }
    static inline void
    sub (POLYNOMIAL &r, const POLYNOMIAL &a, const C & b) {
      sub (r.rep (), a.rep (), b); 
    }
    static inline void
    mul (POLYNOMIAL &r, const POLYNOMIAL &a ) {
      mul (r.rep (), a.rep () ); 
    }
    static inline void
    mul (POLYNOMIAL &r, const POLYNOMIAL &a, const POLYNOMIAL &b) {
      mul (r.rep (), a.rep (), b.rep ()); 
}
    static inline void
    mul (POLYNOMIAL &r, const POLYNOMIAL &a, const C & b) {
      mul (r.rep (), a.rep (), b); }
    static inline void
    mul (POLYNOMIAL &r, const C & a, const POLYNOMIAL &b) {
      mul (r.rep (), b.rep (), a); 
    }
    static inline void
    div (POLYNOMIAL &r, const POLYNOMIAL &a, const POLYNOMIAL &b) {
      div (r.rep (), a.rep (), b.rep ()); 
    }
    static inline void
    div (POLYNOMIAL &r, const POLYNOMIAL &a, const C & b) {
      div (r.rep (), a.rep (), b); 
    }
    static inline void
    rem (POLYNOMIAL &r, const POLYNOMIAL &a, const POLYNOMIAL &b) {
      rem (r.rep (), b.rep (), a.rep); 
    }
};
//====================================================================
TMPL inline POLYNOMIAL
operator-(const POLYNOMIAL& p) { 
  POLYNOMIAL r(p); mul(r.rep(),(typename POLYNOMIAL::value_type)(-1)); return r;
}

//----------------------------------------------------------------------
# define OP use<operators_of,typename POLYNOMIAL::rep_t>
TMPL  define_operator_rrr(POLYNOMIAL,POLYNOMIAL,POLYNOMIAL,operator+,OP::add)
TMPLX define_operator_r_r(POLYNOMIAL,X         ,POLYNOMIAL,operator+,OP::add)
TMPLX define_operator_rr_(POLYNOMIAL,POLYNOMIAL,X         ,operator+,OP::add)

TMPL  define_operator_rrr(POLYNOMIAL,POLYNOMIAL,POLYNOMIAL,operator-,OP::sub)
TMPLX define_operator_rr_(POLYNOMIAL,POLYNOMIAL,X         ,operator-,OP::sub)
TMPLX define_operator_r_r(POLYNOMIAL,X         ,POLYNOMIAL,operator-,OP::sub)

TMPL  define_operator_rrr(POLYNOMIAL,POLYNOMIAL,POLYNOMIAL,operator*,OP::mul)
TMPLX define_operator_rr_(POLYNOMIAL,POLYNOMIAL,X         ,operator*,OP::mul)
TMPLX define_operator_r_r(POLYNOMIAL,X         ,POLYNOMIAL,operator*,OP::mul)

TMPL  define_operator_rrr(POLYNOMIAL,POLYNOMIAL,POLYNOMIAL,operator/,OP::div)
TMPLX define_operator_rr_(POLYNOMIAL,POLYNOMIAL,X,operator/,OP::div)

TMPL  define_operator_rrr(POLYNOMIAL,POLYNOMIAL,POLYNOMIAL,operator%,OP::rem)
# undef OP
//----------------------------------------------------------------------
TMPL POLYNOMIAL 
operator^(const POLYNOMIAL &a, unsigned n) {
  return pow (a, n); }


TMPL bool 
operator== (const POLYNOMIAL &mp, const C & c) {
  return mp.rep() == c; 
}

TMPL bool 
operator== (const POLYNOMIAL &mp, unsigned n) {
  return mp.rep() == (C)n; 
}


  template<class C, class V, class X> 
bool operator !=( const polynomial<C,V>& p, const X& q ) {return !(p==q);}

//====================================================================== 
} //namespace mmx
# undef TMPL
# undef TMPLX
# undef VARIANT
# undef POLYNOMIAL
# undef Scalar
//====================================================================
#endif