/Users/mourrain/Devel/mmx/shape/glue/glue_rational_curve.cpp

Go to the documentation of this file.

#include <basix/double.hpp>
#include <basix/vector.hpp>
#include <basix/port.hpp>
#include <numerix/integer.hpp>
#include <numerix/rational.hpp>
#include <numerix/floating.hpp>
#include <realroot/Interval_glue.hpp>
#include <basix/int.hpp>
#include <numerix/kernel.hpp>
#include <realroot/polynomial.hpp>
#include <realroot/polynomial_glue.hpp>
#include <realroot/polynomial_sparse_glue.hpp>
#include <shape/axel_glue.hpp>
#include <shape/curve_rational_glue.hpp>
#include <basix/alias.hpp>
#include <basix/glue.hpp>

#define double_literal(x) as_double (as_string (x))
#define int_literal(x) as_int (as_string (x))
#define set_of_generic set_of(generic)
#define set_of_double set_of(double)
#define set_of_integer set_of(integer)
#define set_of_rational set_of(rational)
#define set_of_bigfloat set_of(bigfloat)
#define set_of_complex_bigfloat set_of(complex_bigfloat)

namespace mmx {
  static interval<double>
  GLUE_1 (const double &arg_1) {
    return interval<double > (arg_1);
  }
  
  static interval<double>
  GLUE_2 (const double &arg_1, const double &arg_2) {
    return interval_from_pair (arg_1, arg_2);
  }
  
  static interval<double>
  GLUE_3 (const interval<double> &arg_1) {
    return -arg_1;
  }
  
  static interval<double>
  GLUE_4 (const interval<double> &arg_1, const interval<double> &arg_2) {
    return arg_1 + arg_2;
  }
  
  static interval<double>
  GLUE_5 (const interval<double> &arg_1, const interval<double> &arg_2) {
    return arg_1 - arg_2;
  }
  
  static interval<double>
  GLUE_6 (const interval<double> &arg_1, const interval<double> &arg_2) {
    return arg_1 * arg_2;
  }
  
  static interval<double>
  GLUE_7 (const double &arg_1, const interval<double> &arg_2) {
    return arg_1 + arg_2;
  }
  
  static interval<double>
  GLUE_8 (const interval<double> &arg_1, const double &arg_2) {
    return arg_1 + arg_2;
  }
  
  static interval<double>
  GLUE_9 (const double &arg_1, const interval<double> &arg_2) {
    return arg_1 - arg_2;
  }
  
  static interval<double>
  GLUE_10 (const interval<double> &arg_1, const double &arg_2) {
    return arg_1 - arg_2;
  }
  
  static interval<double>
  GLUE_11 (const double &arg_1, const interval<double> &arg_2) {
    return arg_1 * arg_2;
  }
  
  static interval<double>
  GLUE_12 (const interval<double> &arg_1, const double &arg_2) {
    return arg_1 * arg_2;
  }
  
  static interval<double>
  GLUE_13 (const interval<double> &arg_1, const interval<double> &arg_2) {
    return arg_1 / arg_2;
  }
  
  static double
  GLUE_14 (const interval<double> &arg_1) {
    return interval_lower (arg_1);
  }
  
  static double
  GLUE_15 (const interval<double> &arg_1) {
    return interval_upper (arg_1);
  }
  
  static alias<shape_axel>
  GLUE_16 (const alias<shape_axel> &arg_1, const shape_curve_rational &arg_2) {
    return alias_write (arg_1, arg_2);
  }
  
  static shape_curve_rational
  GLUE_17 (const ring<rational, Sparse, DegRevLex>::Polynomial &arg_1, const ring<rational, Sparse, DegRevLex>::Polynomial &arg_2, const ring<rational, Sparse, DegRevLex>::Polynomial &arg_3) {
    return create_curve_rational (arg_1, arg_2, arg_3);
  }
  
  static shape_curve_rational
  GLUE_18 (const interval<double> &arg_1, const ring<rational, Sparse, DegRevLex>::Polynomial &arg_2, const ring<rational, Sparse, DegRevLex>::Polynomial &arg_3, const ring<rational, Sparse, DegRevLex>::Polynomial &arg_4) {
    return create_curve_rational (arg_1, arg_2, arg_3, arg_4);
  }
  
  static shape_curve_rational
  GLUE_19 (const interval<double> &arg_1, const ring<rational, Sparse, DegRevLex>::Polynomial &arg_2, const ring<rational, Sparse, DegRevLex>::Polynomial &arg_3, const ring<rational, Sparse, DegRevLex>::Polynomial &arg_4, const ring<rational, Sparse, DegRevLex>::Polynomial &arg_5) {
    return create_curve_rational (arg_1, arg_2, arg_3, arg_4, arg_5);
  }
  
  static shape_curve_rational
  GLUE_20 (const ring<rational, Sparse, DegRevLex>::Polynomial &arg_1, const ring<rational, Sparse, DegRevLex>::Polynomial &arg_2, const ring<rational, Sparse, DegRevLex>::Polynomial &arg_3, const ring<rational, Sparse, DegRevLex>::Polynomial &arg_4) {
    return create_curve_rational (arg_1, arg_2, arg_3, arg_4);
  }
  
  void
  glue_rational_curve () {
    static bool done = false;
    if (done) return;
    done = true;
    call_glue (string ("glue_double"));
    call_glue (string ("glue_string"));
    call_glue (string ("glue_vector_generic"));
    call_glue (string ("glue_interval"));
    call_glue (string ("glue_polynomial_sparse_rational"));
    call_glue (string ("glue_axel"));
    define_type<interval<double> > (gen (lit ("Interval"), lit ("Double")));
    define ("interval", GLUE_1);
    define ("interval", GLUE_2);
    define ("-", GLUE_3);
    define ("+", GLUE_4);
    define ("-", GLUE_5);
    define ("*", GLUE_6);
    define ("+", GLUE_7);
    define ("+", GLUE_8);
    define ("-", GLUE_9);
    define ("-", GLUE_10);
    define ("*", GLUE_11);
    define ("*", GLUE_12);
    define ("/", GLUE_13);
    define ("lower", GLUE_14);
    define ("upper", GLUE_15);
    define_type<shape_curve_rational > (lit ("RationalCurve"));
    define ("<<", GLUE_16);
    define ("rational_curve", GLUE_17);
    define ("rational_curve", GLUE_18);
    define ("rational_curve", GLUE_19);
    define ("rational_curve", GLUE_20);
  }
}