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

Go to the documentation of this file.

#include <numerix/integer.hpp>
#include <numerix/rational.hpp>
#include <numerix/floating.hpp>
#include <basix/double.hpp>
#include <numerix/complex.hpp>
#include <numerix/complex_double.hpp>
#include <numerix/ball.hpp>
#include <numerix/ball_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/algebraic.hpp>
#include <algebramix/algebraic_number.hpp>
#include <basix/glue.hpp>

#define double_literal(x) as_double (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 algebraic_real
  GLUE_1 (const polynomial<rational> &arg_1, const mmx_ball(mmx_floating, mmx_floating) &arg_2) {
    return algebraic_real (arg_1, arg_2);
  }
  
  static algebraic_real
  GLUE_2 (const polynomial<rational> &arg_1, const polynomial<rational> &arg_2, const mmx_ball(mmx_floating, mmx_floating) &arg_3) {
    return algebraic_real (arg_1, arg_2, arg_3);
  }
  
  static algebraic_real
  GLUE_3 (const rational &arg_1) {
    return algebraic_real (arg_1);
  }
  
  static mmx_floating
  GLUE_4 (const algebraic_real &arg_1) {
    return as_floating (arg_1);
  }
  
  static mmx_ball(mmx_floating, mmx_floating)
  GLUE_5 (const algebraic_real &arg_1) {
    return as_ball (arg_1);
  }
  
  static polynomial<rational>
  GLUE_6 (const algebraic_real &arg_1) {
    return annihilator (arg_1);
  }
  
  static algebraic_real
  GLUE_7 (const algebraic_real &arg_1) {
    return normalize (arg_1);
  }
  
  static algebraic_real
  GLUE_8 (const algebraic_real &arg_1) {
    return -arg_1;
  }
  
  static algebraic_real
  GLUE_9 (const algebraic_real &arg_1) {
    return square (arg_1);
  }
  
  static algebraic_real
  GLUE_10 (const algebraic_real &arg_1, const algebraic_real &arg_2) {
    return arg_1 + arg_2;
  }
  
  static algebraic_real
  GLUE_11 (const algebraic_real &arg_1, const algebraic_real &arg_2) {
    return arg_1 - arg_2;
  }
  
  static algebraic_real
  GLUE_12 (const algebraic_real &arg_1, const algebraic_real &arg_2) {
    return arg_1 * arg_2;
  }
  
  static algebraic_real
  GLUE_13 (const algebraic_real &arg_1, const algebraic_real &arg_2) {
    return arg_1 / arg_2;
  }
  
  static algebraic_real
  GLUE_14 (const rational &arg_1, const algebraic_real &arg_2) {
    return arg_1 + arg_2;
  }
  
  static algebraic_real
  GLUE_15 (const algebraic_real &arg_1, const rational &arg_2) {
    return arg_1 + arg_2;
  }
  
  static algebraic_real
  GLUE_16 (const rational &arg_1, const algebraic_real &arg_2) {
    return arg_1 - arg_2;
  }
  
  static algebraic_real
  GLUE_17 (const algebraic_real &arg_1, const rational &arg_2) {
    return arg_1 - arg_2;
  }
  
  static algebraic_real
  GLUE_18 (const rational &arg_1, const algebraic_real &arg_2) {
    return arg_1 * arg_2;
  }
  
  static algebraic_real
  GLUE_19 (const algebraic_real &arg_1, const rational &arg_2) {
    return arg_1 * arg_2;
  }
  
  static algebraic_real
  GLUE_20 (const rational &arg_1, const algebraic_real &arg_2) {
    return arg_1 / arg_2;
  }
  
  static algebraic_real
  GLUE_21 (const algebraic_real &arg_1, const rational &arg_2) {
    return arg_1 / arg_2;
  }
  
  static algebraic_real
  GLUE_22 (const algebraic_real &arg_1) {
    return sqrt (arg_1);
  }
  
  static algebraic_real
  GLUE_23 (const algebraic_real &arg_1, const int &arg_2) {
    return root (arg_1, arg_2);
  }
  
  static int
  GLUE_24 (const algebraic_real &arg_1) {
    return sign (arg_1);
  }
  
  static bool
  GLUE_25 (const algebraic_real &arg_1, const algebraic_real &arg_2) {
    return arg_1 < arg_2;
  }
  
  static bool
  GLUE_26 (const algebraic_real &arg_1, const algebraic_real &arg_2) {
    return arg_1 <= arg_2;
  }
  
  static bool
  GLUE_27 (const algebraic_real &arg_1, const algebraic_real &arg_2) {
    return arg_1 > arg_2;
  }
  
  static bool
  GLUE_28 (const algebraic_real &arg_1, const algebraic_real &arg_2) {
    return arg_1 >= arg_2;
  }
  
  static algebraic_number
  GLUE_29 (const polynomial<rational> &arg_1, const mmx_ball(mmx_floating, complex<mmx_floating> ) &arg_2) {
    return algebraic_number (arg_1, arg_2);
  }
  
  static algebraic_number
  GLUE_30 (const polynomial<rational> &arg_1, const polynomial<rational> &arg_2, const mmx_ball(mmx_floating, complex<mmx_floating> ) &arg_3) {
    return algebraic_number (arg_1, arg_2, arg_3);
  }
  
  static algebraic_number
  GLUE_31 (const rational &arg_1) {
    return algebraic_number (arg_1);
  }
  
  static algebraic_number
  GLUE_32 (const algebraic_real &arg_1) {
    return algebraic_number (arg_1);
  }
  
  static complex<mmx_floating>
  GLUE_33 (const algebraic_number &arg_1) {
    return as_floating (arg_1);
  }
  
  static mmx_ball(mmx_floating, complex<mmx_floating> )
  GLUE_34 (const algebraic_number &arg_1) {
    return as_ball (arg_1);
  }
  
  static polynomial<rational>
  GLUE_35 (const algebraic_number &arg_1) {
    return annihilator (arg_1);
  }
  
  static algebraic_number
  GLUE_36 (const algebraic_number &arg_1) {
    return normalize (arg_1);
  }
  
  static algebraic_number
  GLUE_37 (const algebraic_number &arg_1) {
    return -arg_1;
  }
  
  static algebraic_number
  GLUE_38 (const algebraic_number &arg_1) {
    return square (arg_1);
  }
  
  static algebraic_number
  GLUE_39 (const algebraic_number &arg_1, const algebraic_number &arg_2) {
    return arg_1 + arg_2;
  }
  
  static algebraic_number
  GLUE_40 (const algebraic_number &arg_1, const algebraic_number &arg_2) {
    return arg_1 - arg_2;
  }
  
  static algebraic_number
  GLUE_41 (const algebraic_number &arg_1, const algebraic_number &arg_2) {
    return arg_1 * arg_2;
  }
  
  static algebraic_number
  GLUE_42 (const algebraic_number &arg_1, const algebraic_number &arg_2) {
    return arg_1 / arg_2;
  }
  
  static algebraic_number
  GLUE_43 (const rational &arg_1, const algebraic_number &arg_2) {
    return arg_1 + arg_2;
  }
  
  static algebraic_number
  GLUE_44 (const algebraic_number &arg_1, const rational &arg_2) {
    return arg_1 + arg_2;
  }
  
  static algebraic_number
  GLUE_45 (const rational &arg_1, const algebraic_number &arg_2) {
    return arg_1 - arg_2;
  }
  
  static algebraic_number
  GLUE_46 (const algebraic_number &arg_1, const rational &arg_2) {
    return arg_1 - arg_2;
  }
  
  static algebraic_number
  GLUE_47 (const rational &arg_1, const algebraic_number &arg_2) {
    return arg_1 * arg_2;
  }
  
  static algebraic_number
  GLUE_48 (const algebraic_number &arg_1, const rational &arg_2) {
    return arg_1 * arg_2;
  }
  
  static algebraic_number
  GLUE_49 (const rational &arg_1, const algebraic_number &arg_2) {
    return arg_1 / arg_2;
  }
  
  static algebraic_number
  GLUE_50 (const algebraic_number &arg_1, const rational &arg_2) {
    return arg_1 / arg_2;
  }
  
  static algebraic_number
  GLUE_51 (const algebraic_number &arg_1) {
    return sqrt (arg_1);
  }
  
  static algebraic_number
  GLUE_52 (const algebraic_number &arg_1, const int &arg_2) {
    return root (arg_1, arg_2);
  }
  
  static algebraic_number
  GLUE_53 (const algebraic_real &arg_1, const algebraic_real &arg_2) {
    return gaussian (arg_1, arg_2);
  }
  
  static algebraic_real
  GLUE_54 (const algebraic_number &arg_1) {
    return abs (arg_1);
  }
  
  static algebraic_real
  GLUE_55 (const algebraic_number &arg_1) {
    return Re (arg_1);
  }
  
  static algebraic_real
  GLUE_56 (const algebraic_number &arg_1) {
    return Im (arg_1);
  }
  
  static algebraic_number
  GLUE_57 (const algebraic_number &arg_1) {
    return conj (arg_1);
  }
  
  static algebraic_number
  GLUE_58 (const algebraic_number &arg_1) {
    return times_i (arg_1);
  }
  
  static algebraic_number
  GLUE_59 (const algebraic_number &arg_1) {
    return over_i (arg_1);
  }
  
  void
  glue_algebraic_number () {
    static bool done = false;
    if (done) return;
    done = true;
    call_glue (string ("glue_ball_floating"));
    call_glue (string ("glue_polynomial_rational"));
    call_glue (string ("glue_algebraic_generic"));
    define_type<algebraic_real > (lit ("Algebraic_real"));
    define ("algebraic", GLUE_1);
    define ("algebraic", GLUE_2);
    define_converter ("upgrade", GLUE_3, PENALTY_INCLUSION);
    define_converter ("upgrade", GLUE_4, PENALTY_INCLUSION);
    define_converter ("upgrade", GLUE_5, PENALTY_INCLUSION);
    define ("annihilator", GLUE_6);
    define ("normalize", GLUE_7);
    define ("-", GLUE_8);
    define ("square", 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);
    define ("/", GLUE_21);
    define ("sqrt", GLUE_22);
    define ("root", GLUE_23);
    define ("sign", GLUE_24);
    define ("<", GLUE_25);
    define ("<=", GLUE_26);
    define (">", GLUE_27);
    define (">=", GLUE_28);
    define_type<algebraic_number > (lit ("Algebraic_number"));
    define ("algebraic", GLUE_29);
    define ("algebraic", GLUE_30);
    define_converter ("upgrade", GLUE_31, PENALTY_INCLUSION);
    define_converter ("upgrade", GLUE_32, PENALTY_INCLUSION);
    define_converter ("upgrade", GLUE_33, PENALTY_INCLUSION);
    define_converter ("upgrade", GLUE_34, PENALTY_INCLUSION);
    define ("annihilator", GLUE_35);
    define ("normalize", GLUE_36);
    define ("-", GLUE_37);
    define ("square", GLUE_38);
    define ("+", GLUE_39);
    define ("-", GLUE_40);
    define ("*", GLUE_41);
    define ("/", GLUE_42);
    define ("+", GLUE_43);
    define ("+", GLUE_44);
    define ("-", GLUE_45);
    define ("-", GLUE_46);
    define ("*", GLUE_47);
    define ("*", GLUE_48);
    define ("/", GLUE_49);
    define ("/", GLUE_50);
    define ("sqrt", GLUE_51);
    define ("root", GLUE_52);
    define ("complex", GLUE_53);
    define ("abs", GLUE_54);
    define ("Re", GLUE_55);
    define ("Im", GLUE_56);
    define ("conj", GLUE_57);
    define ("times_i", GLUE_58);
    define ("over_i", GLUE_59);
  }
}