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

Go to the documentation of this file.

#include <basix/vector.hpp>
#include <basix/int.hpp>
#include <numerix/integer.hpp>
#include <numerix/modular.hpp>
#include <numerix/modular_integer.hpp>
#include <algebramix/vector_unrolled.hpp>
#include <algebramix/vector_simd.hpp>
#include <algebramix/vector_modular.hpp>
#include <algebramix/polynomial.hpp>
#include <algebramix/polynomial_polynomial.hpp>
#include <algebramix/polynomial_integer.hpp>
#include <algebramix/polynomial_modular.hpp>
#include <algebramix/polynomial_modular_integer.hpp>
#include <algebramix/polynomial_schonhage.hpp>
#include <basix/tuple.hpp>
#include <basix/glue.hpp>

#define int_literal(x) as_int (as_string (x))
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 polynomial<mmx_modular(integer) >
  GLUE_1 (const tuple<mmx_modular(integer) > &arg_1) {
    return polynomial_reverse (as_vector (arg_1));
  }
  
  static polynomial<mmx_modular(integer) >
  GLUE_2 (const tuple<mmx_modular(integer) > &arg_1) {
    return polynomial<mmx_modular(integer) > (as_vector (arg_1));
  }
  
  static void
  GLUE_3 (const polynomial<mmx_modular(integer) > &arg_1, const generic &arg_2) {
    set_variable_name (arg_1, arg_2);
  }
  
  static polynomial<mmx_modular(integer) >
  GLUE_4 (const mmx_modular(integer) &arg_1) {
    return polynomial<mmx_modular(integer) > (arg_1);
  }
  
  static iterator<generic>
  GLUE_5 (const polynomial<mmx_modular(integer) > &arg_1) {
    return as<iterator<generic> > (iterate (arg_1));
  }
  
  static int
  GLUE_6 (const polynomial<mmx_modular(integer) > &arg_1) {
    return N (arg_1);
  }
  
  static int
  GLUE_7 (const polynomial<mmx_modular(integer) > &arg_1) {
    return deg (arg_1);
  }
  
  static mmx_modular(integer)
  GLUE_8 (const polynomial<mmx_modular(integer) > &arg_1, const int &arg_2) {
    return arg_1[arg_2];
  }
  
  static polynomial<mmx_modular(integer) >
  GLUE_9 (const polynomial<mmx_modular(integer) > &arg_1) {
    return -arg_1;
  }
  
  static polynomial<mmx_modular(integer) >
  GLUE_10 (const polynomial<mmx_modular(integer) > &arg_1) {
    return square (arg_1);
  }
  
  static polynomial<mmx_modular(integer) >
  GLUE_11 (const polynomial<mmx_modular(integer) > &arg_1, const polynomial<mmx_modular(integer) > &arg_2) {
    return arg_1 + arg_2;
  }
  
  static polynomial<mmx_modular(integer) >
  GLUE_12 (const polynomial<mmx_modular(integer) > &arg_1, const polynomial<mmx_modular(integer) > &arg_2) {
    return arg_1 - arg_2;
  }
  
  static polynomial<mmx_modular(integer) >
  GLUE_13 (const polynomial<mmx_modular(integer) > &arg_1, const polynomial<mmx_modular(integer) > &arg_2) {
    return arg_1 * arg_2;
  }
  
  static polynomial<mmx_modular(integer) >
  GLUE_14 (const mmx_modular(integer) &arg_1, const polynomial<mmx_modular(integer) > &arg_2) {
    return arg_1 + arg_2;
  }
  
  static polynomial<mmx_modular(integer) >
  GLUE_15 (const polynomial<mmx_modular(integer) > &arg_1, const mmx_modular(integer) &arg_2) {
    return arg_1 + arg_2;
  }
  
  static polynomial<mmx_modular(integer) >
  GLUE_16 (const mmx_modular(integer) &arg_1, const polynomial<mmx_modular(integer) > &arg_2) {
    return arg_1 - arg_2;
  }
  
  static polynomial<mmx_modular(integer) >
  GLUE_17 (const polynomial<mmx_modular(integer) > &arg_1, const mmx_modular(integer) &arg_2) {
    return arg_1 - arg_2;
  }
  
  static polynomial<mmx_modular(integer) >
  GLUE_18 (const mmx_modular(integer) &arg_1, const polynomial<mmx_modular(integer) > &arg_2) {
    return arg_1 * arg_2;
  }
  
  static polynomial<mmx_modular(integer) >
  GLUE_19 (const polynomial<mmx_modular(integer) > &arg_1, const mmx_modular(integer) &arg_2) {
    return arg_1 * arg_2;
  }
  
  static polynomial<mmx_modular(integer) >
  GLUE_20 (const polynomial<mmx_modular(integer) > &arg_1, const int &arg_2) {
    return binpow (arg_1, arg_2);
  }
  
  static polynomial<mmx_modular(integer) >
  GLUE_21 (const polynomial<mmx_modular(integer) > &arg_1, const int &arg_2) {
    return lshiftz (arg_1, arg_2);
  }
  
  static polynomial<mmx_modular(integer) >
  GLUE_22 (const polynomial<mmx_modular(integer) > &arg_1, const int &arg_2) {
    return rshiftz (arg_1, arg_2);
  }
  
  static polynomial<mmx_modular(integer) >
  GLUE_23 (const polynomial<mmx_modular(integer) > &arg_1) {
    return derive (arg_1);
  }
  
  static polynomial<mmx_modular(integer) >
  GLUE_24 (const polynomial<mmx_modular(integer) > &arg_1) {
    return xderive (arg_1);
  }
  
  static mmx_modular(integer)
  GLUE_25 (const polynomial<mmx_modular(integer) > &arg_1, const mmx_modular(integer) &arg_2) {
    return evaluate (arg_1, arg_2);
  }
  
  static vector<mmx_modular(integer) >
  GLUE_26 (const polynomial<mmx_modular(integer) > &arg_1, const vector<mmx_modular(integer) > &arg_2) {
    return evaluate (arg_1, arg_2);
  }
  
  static mmx_modular(integer)
  GLUE_27 (const polynomial<mmx_modular(integer) > &arg_1, const mmx_modular(integer) &arg_2) {
    return evaluate (arg_1, arg_2);
  }
  
  static vector<mmx_modular(integer) >
  GLUE_28 (const polynomial<mmx_modular(integer) > &arg_1, const vector<mmx_modular(integer) > &arg_2) {
    return evaluate (arg_1, arg_2);
  }
  
  static polynomial<mmx_modular(integer) >
  GLUE_29 (const polynomial<integer> &arg_1, const modulus<integer> &arg_2) {
    return as_polynomial_modular (arg_1, arg_2);
  }
  
  static polynomial<mmx_modular(integer) >
  GLUE_30 (const polynomial<mmx_modular(integer) > &arg_1, const mmx_modular(integer) &arg_2) {
    return arg_1 / arg_2;
  }
  
  static polynomial<mmx_modular(integer) >
  GLUE_31 (const polynomial<mmx_modular(integer) > &arg_1, const polynomial<mmx_modular(integer) > &arg_2) {
    return arg_1 / arg_2;
  }
  
  static polynomial<mmx_modular(integer) >
  GLUE_32 (const polynomial<mmx_modular(integer) > &arg_1, const polynomial<mmx_modular(integer) > &arg_2) {
    return quo (arg_1, arg_2);
  }
  
  static polynomial<mmx_modular(integer) >
  GLUE_33 (const polynomial<mmx_modular(integer) > &arg_1, const polynomial<mmx_modular(integer) > &arg_2) {
    return rem (arg_1, arg_2);
  }
  
  static bool
  GLUE_34 (const polynomial<mmx_modular(integer) > &arg_1, const polynomial<mmx_modular(integer) > &arg_2) {
    return divides (arg_1, arg_2);
  }
  
  static polynomial<mmx_modular(integer) >
  GLUE_35 (const polynomial<mmx_modular(integer) > &arg_1, const polynomial<mmx_modular(integer) > &arg_2, const int &arg_3) {
    return subresultant (arg_1, arg_2, arg_3);
  }
  
  static vector<generic>
  GLUE_36 (const polynomial<mmx_modular(integer) > &arg_1, const polynomial<mmx_modular(integer) > &arg_2) {
    return wrap_subresultants (arg_1, arg_2);
  }
  
  static mmx_modular(integer)
  GLUE_37 (const polynomial<mmx_modular(integer) > &arg_1, const polynomial<mmx_modular(integer) > &arg_2) {
    return resultant (arg_1, arg_2);
  }
  
  static mmx_modular(integer)
  GLUE_38 (const polynomial<mmx_modular(integer) > &arg_1) {
    return discriminant (arg_1);
  }
  
  static polynomial<mmx_modular(integer) >
  GLUE_39 (const polynomial<mmx_modular(integer) > &arg_1) {
    return integrate (arg_1);
  }
  
  static polynomial<mmx_modular(integer) >
  GLUE_40 (const polynomial<mmx_modular(integer) > &arg_1, const polynomial<mmx_modular(integer) > &arg_2) {
    return compose (arg_1, arg_2);
  }
  
  static polynomial<mmx_modular(integer) >
  GLUE_41 (const polynomial<mmx_modular(integer) > &arg_1, const mmx_modular(integer) &arg_2) {
    return q_difference (arg_1, arg_2);
  }
  
  static polynomial<mmx_modular(integer) >
  GLUE_42 (const polynomial<mmx_modular(integer) > &arg_1, const int &arg_2) {
    return dilate (arg_1, arg_2);
  }
  
  static polynomial<mmx_modular(integer) >
  GLUE_43 (const vector<mmx_modular(integer) > &arg_1) {
    return annulator (arg_1);
  }
  
  static polynomial<mmx_modular(integer) >
  GLUE_44 (const vector<mmx_modular(integer) > &arg_1, const vector<mmx_modular(integer) > &arg_2) {
    return interpolate (arg_1, arg_2);
  }
  
  static polynomial<mmx_modular(integer) >
  GLUE_45 (const polynomial<mmx_modular(integer) > &arg_1, const mmx_modular(integer) &arg_2) {
    return shift (arg_1, arg_2);
  }
  
  static polynomial<mmx_modular(integer) >
  GLUE_46 (const polynomial<mmx_modular(integer) > &arg_1) {
    return graeffe (arg_1);
  }
  
  static polynomial<generic>
  GLUE_47 (const polynomial<mmx_modular(integer) > &arg_1) {
    return as<polynomial<generic> > (arg_1);
  }
  
  void
  glue_polynomial_modular_integer () {
    static bool done = false;
    if (done) return;
    done = true;
    call_glue (string ("glue_vector_generic"));
    call_glue (string ("glue_vector_modular_integer"));
    call_glue (string ("glue_polynomial_integer"));
    define_type<polynomial<mmx_modular(integer) > > (gen (lit ("Polynomial"), gen (lit ("Modular"), lit ("Integer"))));
    define ("poly", GLUE_1);
    define ("polynomial", GLUE_2);
    define ("set_variable_name", GLUE_3);
    define_converter ("upgrade", GLUE_4, PENALTY_INCLUSION);
    define_converter (":>", GLUE_5, PENALTY_CAST);
    define ("#", GLUE_6);
    define ("deg", GLUE_7);
    define (".[]", GLUE_8);
    define ("-", GLUE_9);
    define ("square", GLUE_10);
    define ("+", GLUE_11);
    define ("-", GLUE_12);
    define ("*", GLUE_13);
    define ("+", GLUE_14);
    define ("+", GLUE_15);
    define ("-", GLUE_16);
    define ("-", GLUE_17);
    define ("*", GLUE_18);
    define ("*", GLUE_19);
    define ("^", GLUE_20);
    define ("<<", GLUE_21);
    define (">>", GLUE_22);
    define ("derive", GLUE_23);
    define ("xderive", GLUE_24);
    define ("eval", GLUE_25);
    define ("eval", GLUE_26);
    define ("evaluate", GLUE_27);
    define ("evaluate", GLUE_28);
    define ("mod", GLUE_29);
    define ("/", GLUE_30);
    define ("div", GLUE_31);
    define ("quo", GLUE_32);
    define ("rem", GLUE_33);
    define ("divides?", GLUE_34);
    define ("subresultant", GLUE_35);
    define ("subresultants", GLUE_36);
    define ("resultant", GLUE_37);
    define ("discriminant", GLUE_38);
    define ("integrate", GLUE_39);
    define ("@", GLUE_40);
    define ("q_difference", GLUE_41);
    define ("dilate", GLUE_42);
    define ("annulator", GLUE_43);
    define ("interpolate", GLUE_44);
    define ("shift", GLUE_45);
    define ("graeffe", GLUE_46);
    define_converter (":>", GLUE_47, PENALTY_PROMOTE_GENERIC);
  }
}