/Users/mourrain/Devel/mmx/numerix/glue/glue_complex_double.cpp

Go to the documentation of this file.

#include <basix/double.hpp>
#include <numerix/integer.hpp>
#include <numerix/rational.hpp>
#include <numerix/complex.hpp>
#include <numerix/complex_double.hpp>
#include <basix/glue.hpp>

#define double_literal(x) as_double (as_string (x))

namespace mmx {
  static complex<double>
  GLUE_1 (const double &arg_1) {
    return complex<double > (arg_1);
  }
  
  static complex<double>
  GLUE_2 (const double &arg_1, const double &arg_2) {
    return complex<double > (arg_1, arg_2);
  }
  
  static complex<double>
  GLUE_3 (const double &arg_1) {
    return complex<double > (arg_1);
  }
  
  static double
  GLUE_4 (const complex<double> &arg_1) {
    return Re (arg_1);
  }
  
  static double
  GLUE_5 (const complex<double> &arg_1) {
    return Im (arg_1);
  }
  
  static complex<double>
  GLUE_6 (const complex<double> &arg_1) {
    return conj (arg_1);
  }
  
  static complex<double>
  GLUE_7 (const complex<double> &arg_1) {
    return times_i (arg_1);
  }
  
  static complex<double>
  GLUE_8 (const complex<double> &arg_1) {
    return over_i (arg_1);
  }
  
  static complex<double>
  GLUE_9 (const complex<double> &arg_1) {
    return -arg_1;
  }
  
  static complex<double>
  GLUE_10 (const complex<double> &arg_1) {
    return square (arg_1);
  }
  
  static complex<double>
  GLUE_11 (const complex<double> &arg_1, const complex<double> &arg_2) {
    return arg_1 + arg_2;
  }
  
  static complex<double>
  GLUE_12 (const complex<double> &arg_1, const complex<double> &arg_2) {
    return arg_1 - arg_2;
  }
  
  static complex<double>
  GLUE_13 (const complex<double> &arg_1, const complex<double> &arg_2) {
    return arg_1 * arg_2;
  }
  
  static complex<double>
  GLUE_14 (const double &arg_1, const complex<double> &arg_2) {
    return arg_1 + arg_2;
  }
  
  static complex<double>
  GLUE_15 (const complex<double> &arg_1, const double &arg_2) {
    return arg_1 + arg_2;
  }
  
  static complex<double>
  GLUE_16 (const double &arg_1, const complex<double> &arg_2) {
    return arg_1 - arg_2;
  }
  
  static complex<double>
  GLUE_17 (const complex<double> &arg_1, const double &arg_2) {
    return arg_1 - arg_2;
  }
  
  static complex<double>
  GLUE_18 (const double &arg_1, const complex<double> &arg_2) {
    return arg_1 * arg_2;
  }
  
  static complex<double>
  GLUE_19 (const complex<double> &arg_1, const double &arg_2) {
    return arg_1 * arg_2;
  }
  
  static complex<double>
  GLUE_20 (const complex<double> &arg_1, const complex<double> &arg_2) {
    return arg_1 / arg_2;
  }
  
  static complex<double>
  GLUE_21 (const double &arg_1, const complex<double> &arg_2) {
    return arg_1 / arg_2;
  }
  
  static complex<double>
  GLUE_22 (const complex<double> &arg_1, const double &arg_2) {
    return arg_1 / arg_2;
  }
  
  static double
  GLUE_23 (const complex<double> &arg_1) {
    return abs (arg_1);
  }
  
  static double
  GLUE_24 (const complex<double> &arg_1) {
    return arg (arg_1);
  }
  
  static complex<double>
  GLUE_25 (const complex<double> &arg_1) {
    return sqrt (arg_1);
  }
  
  static complex<double>
  GLUE_26 (const complex<double> &arg_1) {
    return exp (arg_1);
  }
  
  static complex<double>
  GLUE_27 (const complex<double> &arg_1) {
    return log (arg_1);
  }
  
  static complex<double>
  GLUE_28 (const complex<double> &arg_1, const complex<double> &arg_2) {
    return pow (arg_1, arg_2);
  }
  
  static complex<double>
  GLUE_29 (const complex<double> &arg_1) {
    return cos (arg_1);
  }
  
  static complex<double>
  GLUE_30 (const complex<double> &arg_1) {
    return sin (arg_1);
  }
  
  static complex<double>
  GLUE_31 (const complex<double> &arg_1) {
    return tan (arg_1);
  }
  
  static complex<double>
  GLUE_32 (const complex<double> &arg_1) {
    return acos (arg_1);
  }
  
  static complex<double>
  GLUE_33 (const complex<double> &arg_1) {
    return asin (arg_1);
  }
  
  static complex<double>
  GLUE_34 (const complex<double> &arg_1) {
    return atan (arg_1);
  }
  
  static bool
  GLUE_35 (const complex<double> &arg_1) {
    return is_finite (arg_1);
  }
  
  static bool
  GLUE_36 (const complex<double> &arg_1) {
    return is_infinite (arg_1);
  }
  
  static bool
  GLUE_37 (const complex<double> &arg_1) {
    return is_nan (arg_1);
  }
  
  static complex<double>
  GLUE_38 (const complex<double> &arg_1) {
    return times_infinity (arg_1);
  }
  
  static int
  GLUE_39 (const complex<double> &arg_1) {
    return precision (arg_1);
  }
  
  static int
  GLUE_40 (const complex<double> &arg_1) {
    return exponent (arg_1);
  }
  
  static double
  GLUE_41 (const complex<double> &arg_1) {
    return magnitude (arg_1);
  }
  
  static complex<double>
  GLUE_42 (const complex<double> &arg_1, const int &arg_2) {
    return incexp2 (arg_1, arg_2);
  }
  
  static complex<double>
  GLUE_43 (const complex<double> &arg_1, const int &arg_2) {
    return decexp2 (arg_1, arg_2);
  }
  
  static double
  GLUE_44 (const complex<double> &arg_1) {
    return rounding_error (arg_1);
  }
  
  static double
  GLUE_45 (const complex<double> &arg_1) {
    return additive_error (arg_1);
  }
  
  static double
  GLUE_46 (const complex<double> &arg_1) {
    return multiplicative_error (arg_1);
  }
  
  static double
  GLUE_47 (const complex<double> &arg_1) {
    return elementary_error (arg_1);
  }
  
  static complex<generic>
  GLUE_48 (const complex<double> &arg_1) {
    return as<complex<generic> > (arg_1);
  }
  
  static complex<double>
  GLUE_49 (const complex<rational> &arg_1) {
    return as<complex<double> > (arg_1);
  }
  
  void
  glue_complex_double () {
    static bool done = false;
    if (done) return;
    done = true;
    call_glue (string ("glue_double"));
    call_glue (string ("glue_complex_rational"));
    define_type<complex<double> > (gen (lit ("Complex"), lit ("Double")));
    define ("complex", GLUE_1);
    define ("complex", GLUE_2);
    define_converter ("upgrade", GLUE_3, PENALTY_HOMOMORPHISM);
    define ("Re", GLUE_4);
    define ("Im", GLUE_5);
    define ("conj", GLUE_6);
    define ("times_i", GLUE_7);
    define ("over_i", GLUE_8);
    define ("-", GLUE_9);
    define ("square", 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 ("/", GLUE_22);
    define ("abs", GLUE_23);
    define ("arg", GLUE_24);
    define ("sqrt", GLUE_25);
    define ("exp", GLUE_26);
    define ("log", GLUE_27);
    define ("^", GLUE_28);
    define ("cos", GLUE_29);
    define ("sin", GLUE_30);
    define ("tan", GLUE_31);
    define ("arccos", GLUE_32);
    define ("arcsin", GLUE_33);
    define ("arctan", GLUE_34);
    define ("finite?", GLUE_35);
    define ("infinite?", GLUE_36);
    define ("nan?", GLUE_37);
    define ("times_infinity", GLUE_38);
    define ("precision", GLUE_39);
    define ("exponent", GLUE_40);
    define ("magnitude", GLUE_41);
    define ("increase_exponent", GLUE_42);
    define ("decrease_exponent", GLUE_43);
    define ("rounding_error", GLUE_44);
    define ("additive_error", GLUE_45);
    define ("multiplicative_error", GLUE_46);
    define ("elementary_error", GLUE_47);
    define_converter (":>", GLUE_48, PENALTY_PROMOTE_GENERIC);
    define_converter (":>", GLUE_49, PENALTY_INCLUSION);
  }
}