glue/glue_algebraic_number.cpp

Go to the documentation of this file.

#include <basix/int.hpp>
#include <basix/vector.hpp>
#include <basix/port.hpp>
#include <basix/literal.hpp>
#include <numerix/integer.hpp>
#include <numerix/modular.hpp>
#include <numerix/modular_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 <algebramix/vector_unrolled.hpp>
#include <algebramix/vector_simd.hpp>
#include <algebramix/vector_modular.hpp>
#include <basix/compound.hpp>
#include <basix/mmx_syntax.hpp>
#include <basix/lisp_syntax.hpp>
#include <basix/cpp_syntax.hpp>
#include <basix/syntactic.hpp>
#include <algebramix/polynomial.hpp>
#include <algebramix/polynomial_polynomial.hpp>
#include <algebramix/polynomial_integer.hpp>
#include <algebramix/polynomial_rational.hpp>
#include <algebramix/polynomial_modular.hpp>
#include <algebramix/polynomial_modular_integer.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 int_literal(x) as_int (as_string (x))
#define is_generic_literal is<literal>
#define gen_literal_apply(f,v) gen (as<generic> (f), v)
#define gen_literal_access(f,v) access (as<generic> (f), v)
#define double_literal(x) as_double (as_string (x))
#define is_generic_compound is<compound>
#define compound_arguments(x) cdr (as_vector (x))
#define gen_compound_apply(f,v) gen (as<generic> (f), v)
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 rational &arg_1) {
    return algebraic_real (arg_1);
  }
  
  static mmx_floating
  GLUE_2 (const algebraic_real &arg_1) {
    return as_floating (arg_1);
  }
  
  static algebraic_real
  GLUE_3 (const algebraic_real &arg_1, const algebraic_real &arg_2) {
    return lcommon (arg_1, arg_2);
  }
  
  static algebraic_real
  GLUE_4 (const algebraic_real &arg_1, const algebraic_real &arg_2) {
    return rcommon (arg_1, arg_2);
  }
  
  static polynomial<rational>
  GLUE_5 (const algebraic_real &arg_1) {
    return annihilator (arg_1);
  }
  
  static algebraic_real
  GLUE_6 (const algebraic_real &arg_1) {
    return normalize (arg_1);
  }
  
  static algebraic_real
  GLUE_7 (const algebraic_real &arg_1) {
    return -arg_1;
  }
  
  static algebraic_real
  GLUE_8 (const algebraic_real &arg_1) {
    return square (arg_1);
  }
  
  static algebraic_real
  GLUE_9 (const algebraic_real &arg_1, const algebraic_real &arg_2) {
    return arg_1 + arg_2;
  }
  
  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 rational &arg_1, const algebraic_real &arg_2) {
    return arg_1 + arg_2;
  }
  
  static algebraic_real
  GLUE_14 (const algebraic_real &arg_1, const rational &arg_2) {
    return arg_1 + arg_2;
  }
  
  static algebraic_real
  GLUE_15 (const rational &arg_1, const algebraic_real &arg_2) {
    return arg_1 - arg_2;
  }
  
  static algebraic_real
  GLUE_16 (const algebraic_real &arg_1, const rational &arg_2) {
    return arg_1 - arg_2;
  }
  
  static algebraic_real
  GLUE_17 (const rational &arg_1, const algebraic_real &arg_2) {
    return arg_1 * arg_2;
  }
  
  static algebraic_real
  GLUE_18 (const algebraic_real &arg_1, const rational &arg_2) {
    return arg_1 * arg_2;
  }
  
  static algebraic_real
  GLUE_19 (const rational &arg_1, const algebraic_real &arg_2) {
    return arg_1 / arg_2;
  }
  
  static algebraic_real
  GLUE_20 (const algebraic_real &arg_1, const rational &arg_2) {
    return arg_1 / arg_2;
  }
  
  static algebraic_real
  GLUE_21 (const algebraic_real &arg_1) {
    return sqrt (arg_1);
  }
  
  static algebraic_real
  GLUE_22 (const algebraic_real &arg_1, const int &arg_2) {
    return root (arg_1, arg_2);
  }
  
  static int
  GLUE_23 (const algebraic_real &arg_1) {
    return sign (arg_1);
  }
  
  static bool
  GLUE_24 (const algebraic_real &arg_1, const algebraic_real &arg_2) {
    return arg_1 == arg_2;
  }
  
  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 bool
  GLUE_29 (const algebraic_real &arg_1, const algebraic_real &arg_2) {
    return arg_1 >= arg_2;
  }
  
  static bool
  GLUE_30 (const algebraic_real &arg_1, const rational &arg_2) {
    return arg_1 == arg_2;
  }
  
  static bool
  GLUE_31 (const algebraic_real &arg_1, const rational &arg_2) {
    return arg_1 != arg_2;
  }
  
  static bool
  GLUE_32 (const algebraic_real &arg_1, const rational &arg_2) {
    return arg_1 < arg_2;
  }
  
  static bool
  GLUE_33 (const algebraic_real &arg_1, const rational &arg_2) {
    return arg_1 <= arg_2;
  }
  
  static bool
  GLUE_34 (const algebraic_real &arg_1, const rational &arg_2) {
    return arg_1 > arg_2;
  }
  
  static bool
  GLUE_35 (const algebraic_real &arg_1, const rational &arg_2) {
    return arg_1 >= arg_2;
  }
  
  static bool
  GLUE_36 (const rational &arg_1, const algebraic_real &arg_2) {
    return arg_1 == arg_2;
  }
  
  static bool
  GLUE_37 (const rational &arg_1, const algebraic_real &arg_2) {
    return arg_1 != arg_2;
  }
  
  static bool
  GLUE_38 (const rational &arg_1, const algebraic_real &arg_2) {
    return arg_1 < arg_2;
  }
  
  static bool
  GLUE_39 (const rational &arg_1, const algebraic_real &arg_2) {
    return arg_1 <= arg_2;
  }
  
  static bool
  GLUE_40 (const rational &arg_1, const algebraic_real &arg_2) {
    return arg_1 > arg_2;
  }
  
  static bool
  GLUE_41 (const rational &arg_1, const algebraic_real &arg_2) {
    return arg_1 >= arg_2;
  }
  
  static algebraic_number
  GLUE_42 (const rational &arg_1) {
    return algebraic_number (arg_1);
  }
  
  static algebraic_number
  GLUE_43 (const algebraic_real &arg_1) {
    return algebraic_number (arg_1);
  }
  
  static complex<mmx_floating>
  GLUE_44 (const algebraic_number &arg_1) {
    return as_floating (arg_1);
  }
  
  static algebraic_number
  GLUE_45 (const algebraic_number &arg_1, const algebraic_number &arg_2) {
    return lcommon (arg_1, arg_2);
  }
  
  static algebraic_number
  GLUE_46 (const algebraic_number &arg_1, const algebraic_number &arg_2) {
    return rcommon (arg_1, arg_2);
  }
  
  static polynomial<rational>
  GLUE_47 (const algebraic_number &arg_1) {
    return annihilator (arg_1);
  }
  
  static algebraic_number
  GLUE_48 (const algebraic_number &arg_1) {
    return normalize (arg_1);
  }
  
  static algebraic_number
  GLUE_49 (const algebraic_number &arg_1) {
    return -arg_1;
  }
  
  static algebraic_number
  GLUE_50 (const algebraic_number &arg_1) {
    return square (arg_1);
  }
  
  static algebraic_number
  GLUE_51 (const algebraic_number &arg_1, const algebraic_number &arg_2) {
    return arg_1 + arg_2;
  }
  
  static algebraic_number
  GLUE_52 (const algebraic_number &arg_1, const algebraic_number &arg_2) {
    return arg_1 - arg_2;
  }
  
  static algebraic_number
  GLUE_53 (const algebraic_number &arg_1, const algebraic_number &arg_2) {
    return arg_1 * arg_2;
  }
  
  static algebraic_number
  GLUE_54 (const algebraic_number &arg_1, const algebraic_number &arg_2) {
    return arg_1 / arg_2;
  }
  
  static algebraic_number
  GLUE_55 (const rational &arg_1, const algebraic_number &arg_2) {
    return arg_1 + arg_2;
  }
  
  static algebraic_number
  GLUE_56 (const algebraic_number &arg_1, const rational &arg_2) {
    return arg_1 + arg_2;
  }
  
  static algebraic_number
  GLUE_57 (const rational &arg_1, const algebraic_number &arg_2) {
    return arg_1 - arg_2;
  }
  
  static algebraic_number
  GLUE_58 (const algebraic_number &arg_1, const rational &arg_2) {
    return arg_1 - arg_2;
  }
  
  static algebraic_number
  GLUE_59 (const rational &arg_1, const algebraic_number &arg_2) {
    return arg_1 * arg_2;
  }
  
  static algebraic_number
  GLUE_60 (const algebraic_number &arg_1, const rational &arg_2) {
    return arg_1 * arg_2;
  }
  
  static algebraic_number
  GLUE_61 (const rational &arg_1, const algebraic_number &arg_2) {
    return arg_1 / arg_2;
  }
  
  static algebraic_number
  GLUE_62 (const algebraic_number &arg_1, const rational &arg_2) {
    return arg_1 / arg_2;
  }
  
  static algebraic_number
  GLUE_63 (const algebraic_number &arg_1) {
    return sqrt (arg_1);
  }
  
  static algebraic_number
  GLUE_64 (const algebraic_number &arg_1, const int &arg_2) {
    return root (arg_1, arg_2);
  }
  
  static bool
  GLUE_65 (const algebraic_number &arg_1, const algebraic_number &arg_2) {
    return arg_1 == arg_2;
  }
  
  static bool
  GLUE_66 (const algebraic_number &arg_1, const algebraic_number &arg_2) {
    return arg_1 != arg_2;
  }
  
  static bool
  GLUE_67 (const algebraic_number &arg_1, const rational &arg_2) {
    return arg_1 == arg_2;
  }
  
  static bool
  GLUE_68 (const algebraic_number &arg_1, const rational &arg_2) {
    return arg_1 != arg_2;
  }
  
  static bool
  GLUE_69 (const rational &arg_1, const algebraic_number &arg_2) {
    return arg_1 == arg_2;
  }
  
  static bool
  GLUE_70 (const rational &arg_1, const algebraic_number &arg_2) {
    return arg_1 != arg_2;
  }
  
  static algebraic_number
  GLUE_71 (const algebraic_real &arg_1, const algebraic_real &arg_2) {
    return gaussian (arg_1, arg_2);
  }
  
  static algebraic_real
  GLUE_72 (const algebraic_number &arg_1) {
    return abs (arg_1);
  }
  
  static algebraic_real
  GLUE_73 (const algebraic_number &arg_1) {
    return Re (arg_1);
  }
  
  static algebraic_real
  GLUE_74 (const algebraic_number &arg_1) {
    return Im (arg_1);
  }
  
  static algebraic_number
  GLUE_75 (const algebraic_number &arg_1) {
    return conj (arg_1);
  }
  
  static algebraic_number
  GLUE_76 (const algebraic_number &arg_1) {
    return times_i (arg_1);
  }
  
  static algebraic_number
  GLUE_77 (const algebraic_number &arg_1) {
    return over_i (arg_1);
  }
  
  static algebraic_real
  GLUE_78 (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_79 (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 mmx_ball(mmx_floating, mmx_floating)
  GLUE_80 (const algebraic_real &arg_1) {
    return as_ball (arg_1);
  }
  
  static algebraic_number
  GLUE_81 (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_82 (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 mmx_ball(mmx_floating, complex<mmx_floating> )
  GLUE_83 (const algebraic_number &arg_1) {
    return as_ball (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_converter ("upgrade", GLUE_1, PENALTY_INCLUSION);
    define_converter ("upgrade", GLUE_2, PENALTY_INCLUSION);
    define ("lcommon", GLUE_3);
    define ("rcommon", GLUE_4);
    define ("annihilator", GLUE_5);
    define ("normalize", 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);
    define ("sqrt", GLUE_21);
    define ("root", GLUE_22);
    define ("sign", GLUE_23);
    define ("=", GLUE_24);
    define ("!=", GLUE_25);
    define ("<", GLUE_26);
    define ("<=", GLUE_27);
    define (">", GLUE_28);
    define (">=", GLUE_29);
    define ("=", GLUE_30);
    define ("!=", GLUE_31);
    define ("<", GLUE_32);
    define ("<=", GLUE_33);
    define (">", GLUE_34);
    define (">=", GLUE_35);
    define ("=", GLUE_36);
    define ("!=", GLUE_37);
    define ("<", GLUE_38);
    define ("<=", GLUE_39);
    define (">", GLUE_40);
    define (">=", GLUE_41);
    define_type<algebraic_number > (lit ("Algebraic_number"));
    define_converter ("upgrade", GLUE_42, PENALTY_INCLUSION);
    define_converter ("upgrade", GLUE_43, PENALTY_INCLUSION);
    define_converter ("upgrade", GLUE_44, PENALTY_INCLUSION);
    define ("lcommon", GLUE_45);
    define ("rcommon", GLUE_46);
    define ("annihilator", GLUE_47);
    define ("normalize", GLUE_48);
    define ("-", GLUE_49);
    define ("square", GLUE_50);
    define ("+", GLUE_51);
    define ("-", GLUE_52);
    define ("*", GLUE_53);
    define ("/", GLUE_54);
    define ("+", GLUE_55);
    define ("+", GLUE_56);
    define ("-", GLUE_57);
    define ("-", GLUE_58);
    define ("*", GLUE_59);
    define ("*", GLUE_60);
    define ("/", GLUE_61);
    define ("/", GLUE_62);
    define ("sqrt", GLUE_63);
    define ("root", GLUE_64);
    define ("=", GLUE_65);
    define ("!=", GLUE_66);
    define ("=", GLUE_67);
    define ("!=", GLUE_68);
    define ("=", GLUE_69);
    define ("!=", GLUE_70);
    define ("complex", GLUE_71);
    define ("abs", GLUE_72);
    define ("Re", GLUE_73);
    define ("Im", GLUE_74);
    define ("conj", GLUE_75);
    define ("times_i", GLUE_76);
    define ("over_i", GLUE_77);
    define ("algebraic", GLUE_78);
    define ("algebraic", GLUE_79);
    define_converter ("upgrade", GLUE_80, PENALTY_INCLUSION);
    define ("algebraic", GLUE_81);
    define ("algebraic", GLUE_82);
    define_converter ("upgrade", GLUE_83, PENALTY_INCLUSION);
  }
}

---