glue/glue_floating.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/alias.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)

namespace mmx {
  static mmx_floating
  GLUE_1 (const literal &arg_1) {
    return make_literal_floating (arg_1);
  }
  
  static mmx_floating
  GLUE_2 (const int &arg_1) {
    return mmx_floating (arg_1);
  }
  
  static mmx_floating
  GLUE_3 (const integer &arg_1) {
    return mmx_floating (arg_1);
  }
  
  static mmx_floating
  GLUE_4 (const rational &arg_1) {
    return mmx_floating (arg_1);
  }
  
  static mmx_floating
  GLUE_5 (const int &arg_1) {
    return mmx_floating (arg_1);
  }
  
  static mmx_floating
  GLUE_6 (const double &arg_1) {
    return mmx_floating (arg_1);
  }
  
  static mmx_floating
  GLUE_7 (const integer &arg_1) {
    return mmx_floating (arg_1);
  }
  
  static mmx_floating
  GLUE_8 (const rational &arg_1) {
    return mmx_floating (arg_1);
  }
  
  static mmx_floating
  GLUE_9 (const string &arg_1) {
    return mmx_floating (arg_1);
  }
  
  static int
  GLUE_10 (const mmx_floating &arg_1) {
    return as_int (arg_1);
  }
  
  static double
  GLUE_11 (const mmx_floating &arg_1) {
    return as_double (arg_1);
  }
  
  static integer
  GLUE_12 (const mmx_floating &arg_1) {
    return as_integer (arg_1);
  }
  
  static string
  GLUE_13 (const mmx_floating &arg_1) {
    return as_string (arg_1);
  }
  
  static mmx_floating
  GLUE_14 (const mmx_floating &arg_1, const mmx_floating &arg_2) {
    return uniform_deviate (arg_1, arg_2);
  }
  
  static mmx_floating
  GLUE_15 (const mmx_floating &arg_1) {
    return -arg_1;
  }
  
  static mmx_floating
  GLUE_16 (const mmx_floating &arg_1) {
    return square (arg_1);
  }
  
  static mmx_floating
  GLUE_17 (const mmx_floating &arg_1, const mmx_floating &arg_2) {
    return arg_1 + arg_2;
  }
  
  static mmx_floating
  GLUE_18 (const mmx_floating &arg_1, const mmx_floating &arg_2) {
    return arg_1 - arg_2;
  }
  
  static mmx_floating
  GLUE_19 (const mmx_floating &arg_1, const mmx_floating &arg_2) {
    return arg_1 * arg_2;
  }
  
  static mmx_floating
  GLUE_20 (const mmx_floating &arg_1, const mmx_floating &arg_2) {
    return arg_1 / arg_2;
  }
  
  static mmx_floating
  GLUE_21 (const mmx_floating &arg_1) {
    return sqrt (arg_1);
  }
  
  static mmx_floating
  GLUE_22 (const mmx_floating &arg_1, const mmx_floating &arg_2) {
    return pow (arg_1, arg_2);
  }
  
  static mmx_floating
  GLUE_23 (const mmx_floating &arg_1) {
    return exp (arg_1);
  }
  
  static mmx_floating
  GLUE_24 (const mmx_floating &arg_1) {
    return log (arg_1);
  }
  
  static mmx_floating
  GLUE_25 (const mmx_floating &arg_1) {
    return cos (arg_1);
  }
  
  static mmx_floating
  GLUE_26 (const mmx_floating &arg_1) {
    return sin (arg_1);
  }
  
  static mmx_floating
  GLUE_27 (const mmx_floating &arg_1) {
    return tan (arg_1);
  }
  
  static mmx_floating
  GLUE_28 (const mmx_floating &arg_1) {
    return acos (arg_1);
  }
  
  static mmx_floating
  GLUE_29 (const mmx_floating &arg_1) {
    return asin (arg_1);
  }
  
  static mmx_floating
  GLUE_30 (const mmx_floating &arg_1) {
    return atan (arg_1);
  }
  
  static bool
  GLUE_31 (const mmx_floating &arg_1, const mmx_floating &arg_2) {
    return arg_1 < arg_2;
  }
  
  static bool
  GLUE_32 (const mmx_floating &arg_1, const mmx_floating &arg_2) {
    return arg_1 <= arg_2;
  }
  
  static bool
  GLUE_33 (const mmx_floating &arg_1, const mmx_floating &arg_2) {
    return arg_1 > arg_2;
  }
  
  static bool
  GLUE_34 (const mmx_floating &arg_1, const mmx_floating &arg_2) {
    return arg_1 >= arg_2;
  }
  
  static int
  GLUE_35 (const mmx_floating &arg_1) {
    return sign (arg_1);
  }
  
  static mmx_floating
  GLUE_36 (const mmx_floating &arg_1) {
    return abs (arg_1);
  }
  
  static mmx_floating
  GLUE_37 (const mmx_floating &arg_1, const mmx_floating &arg_2) {
    return min (arg_1, arg_2);
  }
  
  static mmx_floating
  GLUE_38 (const mmx_floating &arg_1, const mmx_floating &arg_2) {
    return max (arg_1, arg_2);
  }
  
  static mmx_floating
  GLUE_39 (const mmx_floating &arg_1, const mmx_floating &arg_2) {
    return min (arg_1, arg_2);
  }
  
  static mmx_floating
  GLUE_40 (const mmx_floating &arg_1, const mmx_floating &arg_2) {
    return max (arg_1, arg_2);
  }
  
  static mmx_floating
  GLUE_41 (const mmx_floating &arg_1) {
    return floor (arg_1);
  }
  
  static mmx_floating
  GLUE_42 (const mmx_floating &arg_1) {
    return ceil (arg_1);
  }
  
  static mmx_floating
  GLUE_43 (const mmx_floating &arg_1) {
    return trunc (arg_1);
  }
  
  static mmx_floating
  GLUE_44 (const mmx_floating &arg_1) {
    return round (arg_1);
  }
  
  static bool
  GLUE_45 (const mmx_floating &arg_1) {
    return is_finite (arg_1);
  }
  
  static bool
  GLUE_46 (const mmx_floating &arg_1) {
    return is_infinite (arg_1);
  }
  
  static bool
  GLUE_47 (const mmx_floating &arg_1) {
    return is_nan (arg_1);
  }
  
  static mmx_floating
  GLUE_48 (const mmx_floating &arg_1) {
    return times_infinity (arg_1);
  }
  
  static int
  GLUE_49 (const mmx_floating &arg_1) {
    return precision (arg_1);
  }
  
  static mmx_floating
  GLUE_50 (const mmx_floating &arg_1, const int &arg_2) {
    return change_precision (arg_1, arg_2);
  }
  
  static mmx_floating
  GLUE_51 (const mmx_floating &arg_1) {
    return next_above (arg_1);
  }
  
  static mmx_floating
  GLUE_52 (const mmx_floating &arg_1) {
    return next_below (arg_1);
  }
  
  static int
  GLUE_53 (const mmx_floating &arg_1) {
    return exponent (arg_1);
  }
  
  static double
  GLUE_54 (const mmx_floating &arg_1) {
    return magnitude (arg_1);
  }
  
  static mmx_floating
  GLUE_55 (const mmx_floating &arg_1, const int &arg_2) {
    return incexp2 (arg_1, arg_2);
  }
  
  static mmx_floating
  GLUE_56 (const mmx_floating &arg_1, const int &arg_2) {
    return decexp2 (arg_1, arg_2);
  }
  
  static mmx_floating
  GLUE_57 (const mmx_floating &arg_1) {
    return rounding_error (arg_1);
  }
  
  static mmx_floating
  GLUE_58 (const mmx_floating &arg_1) {
    return additive_error (arg_1);
  }
  
  static mmx_floating
  GLUE_59 (const mmx_floating &arg_1) {
    return multiplicative_error (arg_1);
  }
  
  static mmx_floating
  GLUE_60 (const mmx_floating &arg_1) {
    return elementary_error (arg_1);
  }
  
  void
  glue_floating () {
    static bool done = false;
    if (done) return;
    done = true;
    call_glue (string ("glue_rational"));
    static alias<int> mmx_significant_digits_alias = global_alias (((int&) mmx_significant_digits));
    define_constant<alias<int> > ("significant_digits", mmx_significant_digits_alias);
    static alias<int> mmx_bit_precision_alias = global_alias (((int&) mmx_bit_precision));
    define_constant<alias<int> > ("bit_precision", mmx_bit_precision_alias);
    static alias<int> mmx_discrepancy_alias = global_alias (((int&) mmx_discrepancy));
    define_constant<alias<int> > ("discrepancy", mmx_discrepancy_alias);
    static alias<bool> mmx_pretty_exponents_alias = global_alias (((bool&) mmx_pretty_exponents));
    define_constant<alias<bool> > ("pretty_exponents", mmx_pretty_exponents_alias);
    define_type<mmx_floating > (lit ("Floating"));
    define ("literal_floating", GLUE_1);
    define_converter ("upgrade", GLUE_2, PENALTY_INCLUSION);
    define_converter ("upgrade", GLUE_3, PENALTY_INCLUSION);
    define_converter ("upgrade", GLUE_4, PENALTY_INCLUSION);
    define ("as_floating", GLUE_5);
    define ("as_floating", GLUE_6);
    define ("as_floating", GLUE_7);
    define ("as_floating", GLUE_8);
    define ("as_floating", GLUE_9);
    define ("as_int", GLUE_10);
    define ("as_double", GLUE_11);
    define ("as_integer", GLUE_12);
    define ("as_string", GLUE_13);
    define ("uniform_deviate", GLUE_14);
    define ("-", GLUE_15);
    define ("square", GLUE_16);
    define ("+", GLUE_17);
    define ("-", GLUE_18);
    define ("*", GLUE_19);
    define ("/", GLUE_20);
    define ("sqrt", GLUE_21);
    define ("^", GLUE_22);
    define ("exp", GLUE_23);
    define ("log", GLUE_24);
    define ("cos", GLUE_25);
    define ("sin", GLUE_26);
    define ("tan", GLUE_27);
    define ("arccos", GLUE_28);
    define ("arcsin", GLUE_29);
    define ("arctan", GLUE_30);
    define ("<", GLUE_31);
    define ("<=", GLUE_32);
    define (">", GLUE_33);
    define (">=", GLUE_34);
    define ("sign", GLUE_35);
    define ("abs", GLUE_36);
    define ("min", GLUE_37);
    define ("max", GLUE_38);
    define ("inf", GLUE_39);
    define ("sup", GLUE_40);
    define ("floor", GLUE_41);
    define ("ceil", GLUE_42);
    define ("trunc", GLUE_43);
    define ("round", GLUE_44);
    define ("finite?", GLUE_45);
    define ("infinite?", GLUE_46);
    define ("nan?", GLUE_47);
    define ("times_infinity", GLUE_48);
    define ("precision", GLUE_49);
    define ("change_precision", GLUE_50);
    define ("next_above", GLUE_51);
    define ("next_below", GLUE_52);
    define ("exponent", GLUE_53);
    define ("magnitude", GLUE_54);
    define ("increase_exponent", GLUE_55);
    define ("decrease_exponent", GLUE_56);
    define ("rounding_error", GLUE_57);
    define ("additive_error", GLUE_58);
    define ("multiplicative_error", GLUE_59);
    define ("elementary_error", GLUE_60);
  }
}

---