/Users/mourrain/Devel/mmx/algebramix/glue/glue_algebraic_generic.cpp

Go to the documentation of this file.

#include <basix/vector.hpp>
#include <algebramix/polynomial.hpp>
#include <algebramix/polynomial_polynomial.hpp>
#include <algebramix/polynomial_schonhage.hpp>
#include <algebramix/algebraic.hpp>
#include <basix/glue.hpp>

namespace mmx {
    template<typename C> polynomial<C>
    polynomial_reverse (const vector<C>& v) {
      return polynomial<C> (reverse (v)); }

    template<typename C> polynomial<modular<modulus<C>, modular_local> >
    as_polynomial_modular (const polynomial<C>& f, const modulus<C>& p) {
      modular<modulus<C>, modular_local>::set_modulus (p);
      return as<polynomial<modular<modulus<C>, modular_local> > > (f); }

    template<typename C> vector<generic>
    wrap_subresultants (const polynomial<C>& f, const polynomial<C>& g) {
      return as<vector<generic> > (subresultants (f, g)); }

  }
namespace mmx { POLYNOMIAL_GENERIC_USES_SCHONHAGE }

namespace mmx {
  static algebraic<generic>
  GLUE_1 (const polynomial<generic> &arg_1) {
    return algebraic<generic > (arg_1);
  }
  
  static algebraic<generic>
  GLUE_2 (const polynomial<generic> &arg_1, const polynomial<generic> &arg_2) {
    return algebraic<generic > (arg_1, arg_2);
  }
  
  static polynomial<generic>
  GLUE_3 (const algebraic<generic> &arg_1) {
    return annihilator (arg_1);
  }
  
  static algebraic<generic>
  GLUE_4 (const algebraic<generic> &arg_1) {
    return normalize (arg_1);
  }
  
  static algebraic<generic>
  GLUE_5 (const algebraic<generic> &arg_1) {
    return -arg_1;
  }
  
  static algebraic<generic>
  GLUE_6 (const algebraic<generic> &arg_1) {
    return square (arg_1);
  }
  
  static algebraic<generic>
  GLUE_7 (const algebraic<generic> &arg_1, const algebraic<generic> &arg_2) {
    return arg_1 + arg_2;
  }
  
  static algebraic<generic>
  GLUE_8 (const algebraic<generic> &arg_1, const algebraic<generic> &arg_2) {
    return arg_1 - arg_2;
  }
  
  static algebraic<generic>
  GLUE_9 (const algebraic<generic> &arg_1, const algebraic<generic> &arg_2) {
    return arg_1 * arg_2;
  }
  
  static algebraic<generic>
  GLUE_10 (const algebraic<generic> &arg_1, const algebraic<generic> &arg_2) {
    return arg_1 / arg_2;
  }
  
  void
  glue_algebraic_generic () {
    static bool done = false;
    if (done) return;
    done = true;
    call_glue (string ("glue_polynomial_generic"));
    define_type<algebraic<generic> > (gen (lit ("Algebraic"), lit ("Generic")));
    define ("algebraic", GLUE_1);
    define ("algebraic", GLUE_2);
    define ("annihilator", GLUE_3);
    define ("normalize", GLUE_4);
    define ("-", GLUE_5);
    define ("square", GLUE_6);
    define ("+", GLUE_7);
    define ("-", GLUE_8);
    define ("*", GLUE_9);
    define ("/", GLUE_10);
  }
}