/Users/mourrain/Devel/mmx/realroot/glue/glue_polynomial_sparse_rational.cpp

Go to the documentation of this file.

#include <basix/int.hpp>
#include <basix/vector.hpp>
#include <basix/port.hpp>
#include <numerix/integer.hpp>
#include <numerix/rational.hpp>
#include <numerix/floating.hpp>
#include <numerix/kernel.hpp>
#include <realroot/polynomial.hpp>
#include <realroot/polynomial_glue.hpp>
#include <realroot/polynomial_sparse_glue.hpp>
#include <basix/tuple.hpp>
#include <basix/glue.hpp>

#define int_literal(x) as_int (as_string (x))
#define set_of_generic set_of(generic)
#define set_of_double set_of(double)
#define set_of_integer set_of(integer)
#define set_of_rational set_of(rational)
#define set_of_bigfloat set_of(bigfloat)
#define set_of_complex_bigfloat set_of(complex_bigfloat)

namespace mmx {
  static rational
  GLUE_1 (const ring<rational, Sparse, DegRevLex> &arg_1) {
    return sample (arg_1);
  }
  
  static ring<rational, Sparse, DegRevLex>
  GLUE_2 (const ring<rational, Sparse, DegRevLex> &arg_1, const tuple<generic> &arg_2) {
    return ring_sparse_extend_generic (arg_1, as_vector (arg_2));
  }
  
  static ring<rational, Sparse, DegRevLex>
  GLUE_3 (const scalar_set<rational> &arg_1, const tuple<generic> &arg_2) {
    return ring_sparse_generic (arg_1, as_vector (arg_2));
  }
  
  static ring<rational, Sparse, DegRevLex>
  GLUE_4 (const scalar_set<rational> &arg_1, const tuple<generic> &arg_2) {
    return ring_sparse_generic (arg_1, as_vector (arg_2));
  }
  
  static int
  GLUE_5 (const ring<rational, Sparse, DegRevLex>::Polynomial &arg_1) {
    return size (arg_1);
  }
  
  static int
  GLUE_6 (const ring<rational, Sparse, DegRevLex>::Polynomial &arg_1) {
    return degree (arg_1);
  }
  
  static int
  GLUE_7 (const ring<rational, Sparse, DegRevLex>::Polynomial &arg_1) {
    return nbvar (arg_1);
  }
  
  static ring<rational, Sparse, DegRevLex>::Polynomial
  GLUE_8 (const ring<rational, Sparse, DegRevLex>::Polynomial &arg_1) {
    return -arg_1;
  }
  
  static ring<rational, Sparse, DegRevLex>::Polynomial
  GLUE_9 (const ring<rational, Sparse, DegRevLex>::Polynomial &arg_1, const ring<rational, Sparse, DegRevLex>::Polynomial &arg_2) {
    return arg_1 + arg_2;
  }
  
  static ring<rational, Sparse, DegRevLex>::Polynomial
  GLUE_10 (const ring<rational, Sparse, DegRevLex>::Polynomial &arg_1, const ring<rational, Sparse, DegRevLex>::Polynomial &arg_2) {
    return arg_1 - arg_2;
  }
  
  static ring<rational, Sparse, DegRevLex>::Polynomial
  GLUE_11 (const ring<rational, Sparse, DegRevLex>::Polynomial &arg_1, const ring<rational, Sparse, DegRevLex>::Polynomial &arg_2) {
    return arg_1 * arg_2;
  }
  
  static ring<rational, Sparse, DegRevLex>::Polynomial
  GLUE_12 (const ring<rational, Sparse, DegRevLex>::Polynomial &arg_1) {
    return square (arg_1);
  }
  
  static ring<rational, Sparse, DegRevLex>::Polynomial
  GLUE_13 (const ring<rational, Sparse, DegRevLex>::Polynomial &arg_1, const int &arg_2) {
    return pow (arg_1, arg_2);
  }
  
  static ring<rational, Sparse, DegRevLex>::Polynomial
  GLUE_14 (const ring<rational, Sparse, DegRevLex>::Polynomial &arg_1, const int &arg_2) {
    return diff (arg_1, arg_2);
  }
  
  static ring<rational, Sparse, DegRevLex>::Polynomial
  GLUE_15 (const ring<rational, Sparse, DegRevLex>::Polynomial &arg_1, const generic &arg_2) {
    return diff (arg_1, arg_2);
  }
  
  static ring<rational, Sparse, DegRevLex>::Polynomial
  GLUE_16 (const rational &arg_1, const ring<rational, Sparse, DegRevLex>::Polynomial &arg_2) {
    return arg_1 + arg_2;
  }
  
  static ring<rational, Sparse, DegRevLex>::Polynomial
  GLUE_17 (const ring<rational, Sparse, DegRevLex>::Polynomial &arg_1, const rational &arg_2) {
    return arg_1 + arg_2;
  }
  
  static ring<rational, Sparse, DegRevLex>::Polynomial
  GLUE_18 (const rational &arg_1, const ring<rational, Sparse, DegRevLex>::Polynomial &arg_2) {
    return arg_1 - arg_2;
  }
  
  static ring<rational, Sparse, DegRevLex>::Polynomial
  GLUE_19 (const ring<rational, Sparse, DegRevLex>::Polynomial &arg_1, const rational &arg_2) {
    return arg_1 - arg_2;
  }
  
  static ring<rational, Sparse, DegRevLex>::Polynomial
  GLUE_20 (const rational &arg_1, const ring<rational, Sparse, DegRevLex>::Polynomial &arg_2) {
    return arg_1 * arg_2;
  }
  
  static ring<rational, Sparse, DegRevLex>::Polynomial
  GLUE_21 (const ring<rational, Sparse, DegRevLex>::Polynomial &arg_1, const rational &arg_2) {
    return arg_1 * arg_2;
  }
  
  static ring<rational, Sparse, DegRevLex>::Polynomial
  GLUE_22 (const ring<rational, Sparse, DegRevLex>::Polynomial &arg_1, const rational &arg_2) {
    return arg_1 / arg_2;
  }
  
  static ring<rational, Sparse, DegRevLex>::Polynomial
  GLUE_23 (const ring<rational, Sparse, DegRevLex> &arg_1, const rational &arg_2) {
    return polynomial_sparse (arg_1, arg_2);
  }
  
  static ring<rational, Sparse, DegRevLex>::Polynomial
  GLUE_24 (const ring<rational, Sparse, DegRevLex> &arg_1, const rational &arg_2, const int &arg_3, const int &arg_4) {
    return polynomial_sparse (arg_1, arg_2, arg_3, arg_4);
  }
  
  static ring<rational, Sparse, DegRevLex>::Polynomial
  GLUE_25 (const ring<rational, Sparse, DegRevLex> &arg_1, const string &arg_2) {
    return polynomial_sparse (arg_1, arg_2);
  }
  
  static ring<rational, Sparse, DegRevLex>::Polynomial
  GLUE_26 (const ring<rational, Sparse, DegRevLex> &arg_1, const generic &arg_2) {
    return polynomial_sparse (arg_1, arg_2);
  }
  
  static ring<rational, Sparse, DegRevLex>::Polynomial
  GLUE_27 (const ring<rational, Sparse, DegRevLex> &arg_1, const string &arg_2) {
    return polynomial_sparse (arg_1, arg_2);
  }
  
  static ring<rational, Sparse, DegRevLex>::Polynomial
  GLUE_28 (const ring<rational, Sparse, DegRevLex> &arg_1, const generic &arg_2) {
    return polynomial_sparse (arg_1, arg_2);
  }
  
  static vector<generic>
  GLUE_29 (const ring<rational, Sparse, DegRevLex>::Polynomial &arg_1, const int &arg_2) {
    return polynomial_sparse_coefficients (arg_1, arg_2);
  }
  
  static vector<generic>
  GLUE_30 (const ring<rational, Sparse, DegRevLex>::Polynomial &arg_1) {
    return polynomial_sparse_coefficients (arg_1);
  }
  
  static vector<generic>
  GLUE_31 (const ring<rational, Sparse, DegRevLex>::Polynomial &arg_1) {
    return polynomial_sparse_monomials (arg_1);
  }
  
  static ring<rational, Sparse, DegRevLex>::Polynomial
  GLUE_32 (const ring<rational, Sparse, DegRevLex>::Polynomial &arg_1, const ring<rational, Sparse, DegRevLex>::Polynomial &arg_2) {
    return homogenize (arg_1, arg_2);
  }
  
  static ring<rational, Sparse, DegRevLex>::Polynomial
  GLUE_33 (const ring<rational, Sparse, DegRevLex>::Polynomial &arg_1, const int &arg_2, const ring<rational, Sparse, DegRevLex>::Polynomial &arg_3) {
    return homogenize (arg_1, arg_2, arg_3);
  }
  
  static generic
  GLUE_34 (const ring<rational, Sparse, DegRevLex>::Polynomial &arg_1, const vector<generic> &arg_2) {
    return polynomial_sparse_eval_generic (arg_1, arg_2);
  }
  
  static ring<rational, Sparse, DegRevLex>::Polynomial
  GLUE_35 (const ring<rational, Sparse, DegRevLex> &arg_1, const int &arg_2) {
    return arg_1[arg_2];
  }
  
  static vector<generic>
  GLUE_36 (const ring<rational, Sparse, DegRevLex>::Polynomial &arg_1, const ring<rational, Sparse, DegRevLex>::Polynomial &arg_2, const int &arg_3) {
    return polynomial_sturm_sequence (arg_1, arg_2, arg_3);
  }
  
  static ring<rational, Sparse, DegRevLex>::Polynomial
  GLUE_37 (const ring<rational, Sparse, DegRevLex>::Polynomial &arg_1, const ring<rational, Sparse, DegRevLex>::Polynomial &arg_2, const int &arg_3) {
    return polynomial_resultant (arg_1, arg_2, arg_3);
  }
  
  static ring<rational, Sparse, DegRevLex>::Polynomial
  GLUE_38 (const ring<rational, Sparse, DegRevLex>::Polynomial &arg_1, const ring<rational, Sparse, DegRevLex>::Polynomial &arg_2) {
    return arg_1 / arg_2;
  }
  
  static ring<rational, Sparse, DegRevLex>::Polynomial
  GLUE_39 (const ring<rational, Sparse, DegRevLex>::Polynomial &arg_1, const ring<rational, Sparse, DegRevLex>::Polynomial &arg_2) {
    return arg_1 / arg_2;
  }
  
  static ring<rational, Sparse, DegRevLex>::Polynomial
  GLUE_40 (const ring<rational, Sparse, DegRevLex>::Polynomial &arg_1, const ring<rational, Sparse, DegRevLex>::Polynomial &arg_2) {
    return arg_1 % arg_2;
  }
  
  void
  glue_polynomial_sparse_rational () {
    static bool done = false;
    if (done) return;
    done = true;
    call_glue (string ("glue_int"));
    call_glue (string ("glue_vector_generic"));
    call_glue (string ("glue_string"));
    call_glue (string ("glue_kernel"));
    define_type<ring<rational, Sparse, DegRevLex> > (gen (lit ("MonomialSparseRing"), lit ("Rational")));
    define ("sample", GLUE_1);
    define (".[]", GLUE_2);
    define (".[]", GLUE_3);
    define ("ring", GLUE_4);
    define_type<ring<rational, Sparse, DegRevLex>::Polynomial > (gen (lit ("Polynomial"), gen (lit ("MonomialSparseRing"), lit ("Rational"))));
    define ("#", GLUE_5);
    define ("degree", GLUE_6);
    define ("nbvar", GLUE_7);
    define ("-", GLUE_8);
    define ("+", GLUE_9);
    define ("-", GLUE_10);
    define ("*", GLUE_11);
    define ("square", GLUE_12);
    define ("^", GLUE_13);
    define ("diff", GLUE_14);
    define ("diff", GLUE_15);
    define ("+", GLUE_16);
    define ("+", GLUE_17);
    define ("-", GLUE_18);
    define ("-", GLUE_19);
    define ("*", GLUE_20);
    define ("*", GLUE_21);
    define ("/", GLUE_22);
    define ("polynomial", GLUE_23);
    define ("polynomial", GLUE_24);
    define ("polynomial", GLUE_25);
    define ("polynomial", GLUE_26);
    define ("<<", GLUE_27);
    define ("<<", GLUE_28);
    define ("coefficients", GLUE_29);
    define ("coefficients", GLUE_30);
    define ("monomials", GLUE_31);
    define ("homogenize", GLUE_32);
    define ("homogenize", GLUE_33);
    define ("eval", GLUE_34);
    define (".[]", GLUE_35);
    define ("sturm_seq", GLUE_36);
    define ("resultant", GLUE_37);
    define ("/", GLUE_38);
    define ("quo", GLUE_39);
    define ("rem", GLUE_40);
  }
}