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

Go to the documentation of this file.

#include <numerix/integer.hpp>
#include <basix/int.hpp>
#include <numerix/modular.hpp>
#include <numerix/modular_integer.hpp>
#include <basix/vector.hpp>
#include <algebramix/vector_unrolled.hpp>
#include <algebramix/vector_simd.hpp>
#include <algebramix/vector_modular.hpp>
#include <algebramix/p_expansion.hpp>
#include <algebramix/p_expansion_modular_integer.hpp>
#include <algebramix/p_adic.hpp>
#include <basix/tuple.hpp>
#include <basix/glue.hpp>

#define int_literal(x) as_int (as_string (x))

#define simple_p_expansion(C) polynomial<C, polynomial_carry_variant_helper<C>::PV>
#define simple_as_p_expansion(C) as_p_expansion<C,Modulus_variant(C),modular_local>


namespace mmx {
  static inline simple_p_expansion(mmx_modular(integer))
  integer_as_p_expansion(const integer& c, const modulus<integer>& p) {
    return simple_as_p_expansion(integer)(c, p); }
}


#define simple_p_adic(C) series<C, series_carry_variant_helper<C>::SV>


namespace mmx {
  static void
  GLUE_1 (const simple_p_adic(mmx_modular(integer) ) &arg_1, const generic &arg_2) {
    set_variable_name (arg_1, arg_2);
  }
  
  static void
  GLUE_2 (const simple_p_adic(mmx_modular(integer) ) &arg_1, const int &arg_2) {
    set_output_order (arg_1, arg_2);
  }
  
  static void
  GLUE_3 (const simple_p_adic(mmx_modular(integer) ) &arg_1, const int &arg_2) {
    set_cancel_order (arg_1, arg_2);
  }
  
  static void
  GLUE_4 (const simple_p_adic(mmx_modular(integer) ) &arg_1, const bool &arg_2) {
    set_formula_output (arg_1, arg_2);
  }
  
  static simple_p_adic(mmx_modular(integer) )
  GLUE_5 (const tuple<mmx_modular(integer) > &arg_1) {
    return (simple_p_adic(mmx_modular(integer) ) (as_vector (arg_1)));
  }
  
  static simple_p_adic(mmx_modular(integer) )
  GLUE_6 (const mmx_modular(integer) &arg_1) {
    return (simple_p_adic(mmx_modular(integer) ) (arg_1));
  }
  
  static iterator<generic>
  GLUE_7 (const simple_p_adic(mmx_modular(integer) ) &arg_1) {
    return as<iterator<generic> > (iterate (arg_1));
  }
  
  static mmx_modular(integer)
  GLUE_8 (const simple_p_adic(mmx_modular(integer) ) &arg_1, const int &arg_2) {
    return arg_1[arg_2];
  }
  
  static simple_p_expansion(mmx_modular(integer) )
  GLUE_9 (const simple_p_adic(mmx_modular(integer) ) &arg_1, const int &arg_2, const int &arg_3) {
    return range (arg_1, arg_2, arg_3);
  }
  
  static simple_p_adic(mmx_modular(integer) )
  GLUE_10 (const simple_p_adic(mmx_modular(integer) ) &arg_1) {
    return -arg_1;
  }
  
  static simple_p_adic(mmx_modular(integer) )
  GLUE_11 (const simple_p_adic(mmx_modular(integer) ) &arg_1) {
    return square (arg_1);
  }
  
  static simple_p_adic(mmx_modular(integer) )
  GLUE_12 (const simple_p_adic(mmx_modular(integer) ) &arg_1, const simple_p_adic(mmx_modular(integer) ) &arg_2) {
    return arg_1 + arg_2;
  }
  
  static simple_p_adic(mmx_modular(integer) )
  GLUE_13 (const simple_p_adic(mmx_modular(integer) ) &arg_1, const simple_p_adic(mmx_modular(integer) ) &arg_2) {
    return arg_1 - arg_2;
  }
  
  static simple_p_adic(mmx_modular(integer) )
  GLUE_14 (const simple_p_adic(mmx_modular(integer) ) &arg_1, const simple_p_adic(mmx_modular(integer) ) &arg_2) {
    return arg_1 * arg_2;
  }
  
  static simple_p_adic(mmx_modular(integer) )
  GLUE_15 (const mmx_modular(integer) &arg_1, const simple_p_adic(mmx_modular(integer) ) &arg_2) {
    return arg_1 + arg_2;
  }
  
  static simple_p_adic(mmx_modular(integer) )
  GLUE_16 (const simple_p_adic(mmx_modular(integer) ) &arg_1, const mmx_modular(integer) &arg_2) {
    return arg_1 + arg_2;
  }
  
  static simple_p_adic(mmx_modular(integer) )
  GLUE_17 (const mmx_modular(integer) &arg_1, const simple_p_adic(mmx_modular(integer) ) &arg_2) {
    return arg_1 - arg_2;
  }
  
  static simple_p_adic(mmx_modular(integer) )
  GLUE_18 (const simple_p_adic(mmx_modular(integer) ) &arg_1, const mmx_modular(integer) &arg_2) {
    return arg_1 - arg_2;
  }
  
  static simple_p_adic(mmx_modular(integer) )
  GLUE_19 (const mmx_modular(integer) &arg_1, const simple_p_adic(mmx_modular(integer) ) &arg_2) {
    return arg_1 * arg_2;
  }
  
  static simple_p_adic(mmx_modular(integer) )
  GLUE_20 (const simple_p_adic(mmx_modular(integer) ) &arg_1, const mmx_modular(integer) &arg_2) {
    return arg_1 * arg_2;
  }
  
  static simple_p_adic(mmx_modular(integer) )
  GLUE_21 (const simple_p_adic(mmx_modular(integer) ) &arg_1, const int &arg_2) {
    return binpow (arg_1, arg_2);
  }
  
  static simple_p_adic(mmx_modular(integer) )
  GLUE_22 (const simple_p_adic(mmx_modular(integer) ) &arg_1, const simple_p_adic(mmx_modular(integer) ) &arg_2) {
    return arg_1 / arg_2;
  }
  
  static simple_p_adic(mmx_modular(integer) )
  GLUE_23 (const mmx_modular(integer) &arg_1, const simple_p_adic(mmx_modular(integer) ) &arg_2) {
    return arg_1 / arg_2;
  }
  
  static simple_p_adic(mmx_modular(integer) )
  GLUE_24 (const simple_p_adic(mmx_modular(integer) ) &arg_1, const mmx_modular(integer) &arg_2) {
    return arg_1 / arg_2;
  }
  
  static simple_p_adic(mmx_modular(integer) )
  GLUE_25 (const simple_p_adic(mmx_modular(integer) ) &arg_1, const simple_p_adic(mmx_modular(integer) ) &arg_2) {
    return gcd (arg_1, arg_2);
  }
  
  static simple_p_adic(mmx_modular(integer) )
  GLUE_26 (const simple_p_adic(mmx_modular(integer) ) &arg_1, const simple_p_adic(mmx_modular(integer) ) &arg_2) {
    return lcm (arg_1, arg_2);
  }
  
  static simple_p_adic(mmx_modular(integer) )
  GLUE_27 (const simple_p_adic(mmx_modular(integer) ) &arg_1, const int &arg_2) {
    return separable_root (arg_1, arg_2);
  }
  
  static simple_p_adic(mmx_modular(integer) )
  GLUE_28 (const simple_p_adic(mmx_modular(integer) ) &arg_1) {
    return pth_root (arg_1);
  }
  
  void
  glue_p_adic_modular_integer () {
    static bool done = false;
    if (done) return;
    done = true;
    call_glue (string ("glue_integer"));
    call_glue (string ("glue_modular_integer"));
    call_glue (string ("glue_p_expansion_modular_integer"));
    define_type<simple_p_adic(mmx_modular(integer) ) > (gen (lit ("P_adic"), gen (lit ("Modular"), lit ("Integer"))));
    define ("set_variable_name", GLUE_1);
    define ("set_output_order", GLUE_2);
    define ("set_cancel_order", GLUE_3);
    define ("set_formula_output", GLUE_4);
    define ("p_adic", GLUE_5);
    define_converter ("upgrade", GLUE_6, PENALTY_INCLUSION);
    define_converter (":>", GLUE_7, PENALTY_CAST);
    define (".[]", GLUE_8);
    define (".[]", GLUE_9);
    define ("-", GLUE_10);
    define ("square", 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 ("/", GLUE_23);
    define ("/", GLUE_24);
    define ("gcd", GLUE_25);
    define ("lcm", GLUE_26);
    define ("separable_root", GLUE_27);
    define ("pth_root", GLUE_28);
  }
}