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

Go to the documentation of this file.

#include <numerix/integer.hpp>
#include <numerix/rational.hpp>
#include <numerix/complex.hpp>
#include <basix/vector.hpp>
#include <algebramix/polynomial.hpp>
#include <algebramix/polynomial_polynomial.hpp>
#include <algebramix/polynomial_integer.hpp>
#include <algebramix/polynomial_rational.hpp>
#include <algebramix/polynomial_complex.hpp>
#include <algebramix/polynomial_schonhage.hpp>
#include <algebramix/quotient.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 }
#define simple_quotient(C) quotient<C,C >

namespace mmx {
  static simple_quotient(polynomial<rational> )
  GLUE_1 (const polynomial<rational> &arg_1, const polynomial<rational> &arg_2) {
    return (simple_quotient(polynomial<rational> ) (arg_1, arg_2));
  }
  
  static simple_quotient(polynomial<rational> )
  GLUE_2 (const polynomial<rational> &arg_1) {
    return (simple_quotient(polynomial<rational> ) (arg_1));
  }
  
  static simple_quotient(polynomial<rational> )
  GLUE_3 (const polynomial<rational> &arg_1, const polynomial<rational> &arg_2) {
    return (simple_quotient(polynomial<rational> ) (arg_1, arg_2));
  }
  
  static simple_quotient(polynomial<rational> )
  GLUE_4 (const polynomial<rational> &arg_1) {
    return (simple_quotient(polynomial<rational> ) (arg_1));
  }
  
  static polynomial<rational>
  GLUE_5 (const simple_quotient(polynomial<rational> ) &arg_1) {
    return numerator (arg_1);
  }
  
  static polynomial<rational>
  GLUE_6 (const simple_quotient(polynomial<rational> ) &arg_1) {
    return denominator (arg_1);
  }
  
  static simple_quotient(polynomial<rational> )
  GLUE_7 (const simple_quotient(polynomial<rational> ) &arg_1) {
    return -arg_1;
  }
  
  static simple_quotient(polynomial<rational> )
  GLUE_8 (const simple_quotient(polynomial<rational> ) &arg_1) {
    return square (arg_1);
  }
  
  static simple_quotient(polynomial<rational> )
  GLUE_9 (const simple_quotient(polynomial<rational> ) &arg_1, const simple_quotient(polynomial<rational> ) &arg_2) {
    return arg_1 + arg_2;
  }
  
  static simple_quotient(polynomial<rational> )
  GLUE_10 (const simple_quotient(polynomial<rational> ) &arg_1, const simple_quotient(polynomial<rational> ) &arg_2) {
    return arg_1 - arg_2;
  }
  
  static simple_quotient(polynomial<rational> )
  GLUE_11 (const simple_quotient(polynomial<rational> ) &arg_1, const simple_quotient(polynomial<rational> ) &arg_2) {
    return arg_1 * arg_2;
  }
  
  static simple_quotient(polynomial<rational> )
  GLUE_12 (const simple_quotient(polynomial<rational> ) &arg_1, const simple_quotient(polynomial<rational> ) &arg_2) {
    return arg_1 / arg_2;
  }
  
  static simple_quotient(polynomial<rational> )
  GLUE_13 (const polynomial<rational> &arg_1, const simple_quotient(polynomial<rational> ) &arg_2) {
    return arg_1 + arg_2;
  }
  
  static simple_quotient(polynomial<rational> )
  GLUE_14 (const simple_quotient(polynomial<rational> ) &arg_1, const polynomial<rational> &arg_2) {
    return arg_1 + arg_2;
  }
  
  static simple_quotient(polynomial<rational> )
  GLUE_15 (const polynomial<rational> &arg_1, const simple_quotient(polynomial<rational> ) &arg_2) {
    return arg_1 - arg_2;
  }
  
  static simple_quotient(polynomial<rational> )
  GLUE_16 (const simple_quotient(polynomial<rational> ) &arg_1, const polynomial<rational> &arg_2) {
    return arg_1 - arg_2;
  }
  
  static simple_quotient(polynomial<rational> )
  GLUE_17 (const polynomial<rational> &arg_1, const simple_quotient(polynomial<rational> ) &arg_2) {
    return arg_1 * arg_2;
  }
  
  static simple_quotient(polynomial<rational> )
  GLUE_18 (const simple_quotient(polynomial<rational> ) &arg_1, const polynomial<rational> &arg_2) {
    return arg_1 * arg_2;
  }
  
  static simple_quotient(polynomial<rational> )
  GLUE_19 (const polynomial<rational> &arg_1, const simple_quotient(polynomial<rational> ) &arg_2) {
    return arg_1 / arg_2;
  }
  
  static simple_quotient(polynomial<rational> )
  GLUE_20 (const simple_quotient(polynomial<rational> ) &arg_1, const polynomial<rational> &arg_2) {
    return arg_1 / arg_2;
  }
  
  void
  glue_quotient_polynomial_rational () {
    static bool done = false;
    if (done) return;
    done = true;
    call_glue (string ("glue_polynomial_rational"));
    define_type<simple_quotient(polynomial<rational> ) > (gen (lit ("Quotient"), gen (lit ("Polynomial"), lit ("Rational"))));
    define ("quotient", GLUE_1);
    define ("quotient", GLUE_2);
    define ("/", GLUE_3);
    define_converter ("upgrade", GLUE_4, PENALTY_INCLUSION);
    define ("numerator", GLUE_5);
    define ("denominator", GLUE_6);
    define ("-", GLUE_7);
    define ("square", GLUE_8);
    define ("+", GLUE_9);
    define ("-", 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);
  }
}