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

Go to the documentation of this file.

#include <numerix/integer.hpp>
#include <numerix/rational.hpp>
#include <numerix/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/series.hpp>
#include <algebramix/series_elementary.hpp>
#include <algebramix/series_integer.hpp>
#include <algebramix/series_rational.hpp>
#include <algebramix/series_complex.hpp>
#include <algebramix/series_sugar.hpp>
#include <basix/tuple.hpp>
#include <basix/routine.hpp>
#include <basix/glue.hpp>

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 }
#define series_shift_default(s,sh) shift (s, sh, mmx_bit_precision)

namespace mmx {
  static void
  GLUE_1 (const series<rational> &arg_1, const generic &arg_2) {
    set_variable_name (arg_1, arg_2);
  }
  
  static void
  GLUE_2 (const series<rational> &arg_1, const int &arg_2) {
    set_output_order (arg_1, arg_2);
  }
  
  static void
  GLUE_3 (const series<rational> &arg_1, const int &arg_2) {
    set_cancel_order (arg_1, arg_2);
  }
  
  static void
  GLUE_4 (const series<rational> &arg_1, const bool &arg_2) {
    set_formula_output (arg_1, arg_2);
  }
  
  static series<rational>
  GLUE_5 (const tuple<rational> &arg_1) {
    return series<rational > (as_vector (arg_1));
  }
  
  static series<rational>
  GLUE_6 (const rational &arg_1) {
    return series<rational > (arg_1);
  }
  
  static series<rational>
  GLUE_7 (const polynomial<rational> &arg_1) {
    return series<rational > (arg_1);
  }
  
  static iterator<generic>
  GLUE_8 (const series<rational> &arg_1) {
    return as<iterator<generic> > (iterate (arg_1));
  }
  
  static rational
  GLUE_9 (const series<rational> &arg_1, const int &arg_2) {
    return arg_1[arg_2];
  }
  
  static polynomial<rational>
  GLUE_10 (const series<rational> &arg_1, const int &arg_2, const int &arg_3) {
    return range (arg_1, arg_2, arg_3);
  }
  
  static series<rational>
  GLUE_11 (const series<rational> &arg_1) {
    return -arg_1;
  }
  
  static series<rational>
  GLUE_12 (const series<rational> &arg_1) {
    return square (arg_1);
  }
  
  static series<rational>
  GLUE_13 (const series<rational> &arg_1, const series<rational> &arg_2) {
    return arg_1 + arg_2;
  }
  
  static series<rational>
  GLUE_14 (const series<rational> &arg_1, const series<rational> &arg_2) {
    return arg_1 - arg_2;
  }
  
  static series<rational>
  GLUE_15 (const series<rational> &arg_1, const series<rational> &arg_2) {
    return arg_1 * arg_2;
  }
  
  static series<rational>
  GLUE_16 (const rational &arg_1, const series<rational> &arg_2) {
    return arg_1 + arg_2;
  }
  
  static series<rational>
  GLUE_17 (const series<rational> &arg_1, const rational &arg_2) {
    return arg_1 + arg_2;
  }
  
  static series<rational>
  GLUE_18 (const rational &arg_1, const series<rational> &arg_2) {
    return arg_1 - arg_2;
  }
  
  static series<rational>
  GLUE_19 (const series<rational> &arg_1, const rational &arg_2) {
    return arg_1 - arg_2;
  }
  
  static series<rational>
  GLUE_20 (const rational &arg_1, const series<rational> &arg_2) {
    return arg_1 * arg_2;
  }
  
  static series<rational>
  GLUE_21 (const series<rational> &arg_1, const rational &arg_2) {
    return arg_1 * arg_2;
  }
  
  static series<rational>
  GLUE_22 (const series<rational> &arg_1, const int &arg_2) {
    return binpow (arg_1, arg_2);
  }
  
  static series<rational>
  GLUE_23 (const series<rational> &arg_1, const integer &arg_2) {
    return binpow (arg_1, arg_2);
  }
  
  static series<rational>
  GLUE_24 (const series<rational> &arg_1) {
    return derive (arg_1);
  }
  
  static series<rational>
  GLUE_25 (const series<rational> &arg_1) {
    return xderive (arg_1);
  }
  
  static series<rational>
  GLUE_26 (const series<rational> &arg_1, const int &arg_2) {
    return dilate (arg_1, arg_2);
  }
  
  static series<rational>
  GLUE_27 (const series<rational> &arg_1, const int &arg_2) {
    return lshiftz (arg_1, arg_2);
  }
  
  static series<rational>
  GLUE_28 (const series<rational> &arg_1, const int &arg_2) {
    return rshiftz (arg_1, arg_2);
  }
  
  static series<rational>
  GLUE_29 (const series<rational> &arg_1, const series<rational> &arg_2) {
    return arg_1 / arg_2;
  }
  
  static series<rational>
  GLUE_30 (const rational &arg_1, const series<rational> &arg_2) {
    return arg_1 / arg_2;
  }
  
  static series<rational>
  GLUE_31 (const series<rational> &arg_1, const rational &arg_2) {
    return arg_1 / arg_2;
  }
  
  static series<rational>
  GLUE_32 (const series<rational> &arg_1, const series<rational> &arg_2) {
    return arg_1 / arg_2;
  }
  
  static bool
  GLUE_33 (const series<rational> &arg_1, const series<rational> &arg_2) {
    return divides (arg_1, arg_2);
  }
  
  static series<rational>
  GLUE_34 (const series<rational> &arg_1, const series<rational> &arg_2) {
    return gcd (arg_1, arg_2);
  }
  
  static series<rational>
  GLUE_35 (const series<rational> &arg_1, const series<rational> &arg_2) {
    return lcm (arg_1, arg_2);
  }
  
  static series<rational>
  GLUE_36 (const series<rational> &arg_1) {
    return integrate (arg_1);
  }
  
  static series<rational>
  GLUE_37 (const series<rational> &arg_1, const series<rational> &arg_2) {
    return compose (arg_1, arg_2);
  }
  
  static series<rational>
  GLUE_38 (const series<rational> &arg_1) {
    return reverse (arg_1);
  }
  
  static series<rational>
  GLUE_39 (const series<rational> &arg_1, const rational &arg_2) {
    return q_difference (arg_1, arg_2);
  }
  
  static series<rational>
  GLUE_40 (const series<rational> &arg_1, const rational &arg_2) {
    return series_shift_default (arg_1, arg_2);
  }
  
  static series<rational>
  GLUE_41 (const series<rational> &arg_1, const rational &arg_2, const int &arg_3) {
    return shift (arg_1, arg_2, arg_3);
  }
  
  static bool
  GLUE_42 (const series<rational> &arg_1, const series<rational> &arg_2) {
    return arg_1 <= arg_2;
  }
  
  static bool
  GLUE_43 (const series<rational> &arg_1, const series<rational> &arg_2) {
    return arg_1 >= arg_2;
  }
  
  static bool
  GLUE_44 (const series<rational> &arg_1, const series<rational> &arg_2) {
    return arg_1 < arg_2;
  }
  
  static bool
  GLUE_45 (const series<rational> &arg_1, const series<rational> &arg_2) {
    return arg_1 > arg_2;
  }
  
  static series<rational>
  GLUE_46 (const series<integer> &arg_1) {
    return as<series<rational> > (arg_1);
  }
  
  static series<generic>
  GLUE_47 (const series<rational> &arg_1) {
    return as<series<generic> > (arg_1);
  }
  
  static series<rational>
  GLUE_48 (const series<rational> &arg_1, const series<rational> &arg_2) {
    return pow (arg_1, arg_2);
  }
  
  static series<rational>
  GLUE_49 (const series<rational> &arg_1) {
    return sqrt (arg_1);
  }
  
  static series<rational>
  GLUE_50 (const series<rational> &arg_1) {
    return exp (arg_1);
  }
  
  static series<rational>
  GLUE_51 (const series<rational> &arg_1) {
    return log (arg_1);
  }
  
  static series<rational>
  GLUE_52 (const series<rational> &arg_1) {
    return cos (arg_1);
  }
  
  static series<rational>
  GLUE_53 (const series<rational> &arg_1) {
    return sin (arg_1);
  }
  
  static series<rational>
  GLUE_54 (const series<rational> &arg_1) {
    return tan (arg_1);
  }
  
  static series<rational>
  GLUE_55 (const series<rational> &arg_1) {
    return acos (arg_1);
  }
  
  static series<rational>
  GLUE_56 (const series<rational> &arg_1) {
    return asin (arg_1);
  }
  
  static series<rational>
  GLUE_57 (const series<rational> &arg_1) {
    return atan (arg_1);
  }
  
  static unknown<rational>
  GLUE_58 (const rational &arg_1) {
    return unknown<rational > (arg_1);
  }
  
  static unknown<rational>
  GLUE_59 (const unknown<rational> &arg_1) {
    return -arg_1;
  }
  
  static unknown<rational>
  GLUE_60 (const unknown<rational> &arg_1) {
    return square (arg_1);
  }
  
  static unknown<rational>
  GLUE_61 (const unknown<rational> &arg_1, const unknown<rational> &arg_2) {
    return arg_1 + arg_2;
  }
  
  static unknown<rational>
  GLUE_62 (const unknown<rational> &arg_1, const unknown<rational> &arg_2) {
    return arg_1 - arg_2;
  }
  
  static unknown<rational>
  GLUE_63 (const rational &arg_1, const unknown<rational> &arg_2) {
    return arg_1 * arg_2;
  }
  
  static unknown<rational>
  GLUE_64 (const unknown<rational> &arg_1, const rational &arg_2) {
    return arg_1 * arg_2;
  }
  
  static unknown<rational>
  GLUE_65 (const unknown<rational> &arg_1, const unknown<rational> &arg_2) {
    return arg_1 * arg_2;
  }
  
  static series<rational>
  GLUE_66 (const routine &arg_1, const rational &arg_2) {
    return fixed_point_series (arg_1, arg_2);
  }
  
  static vector<generic>
  GLUE_67 (const routine &arg_1, const vector<rational> &arg_2) {
    return gen_fixed_point_vector_series (arg_1, arg_2);
  }
  
  static series<rational>
  GLUE_68 (const routine &arg_1, const rational &arg_2) {
    return integrate_series (arg_1, arg_2);
  }
  
  static vector<generic>
  GLUE_69 (const routine &arg_1, const vector<rational> &arg_2) {
    return gen_integrate_vector_series (arg_1, arg_2);
  }
  
  static series<rational>
  GLUE_70 (const routine &arg_1, const rational &arg_2) {
    return implicit_series (arg_1, arg_2);
  }
  
  static vector<generic>
  GLUE_71 (const routine &arg_1, const vector<rational> &arg_2) {
    return gen_implicit_vector_series (arg_1, arg_2);
  }
  
  static void
  GLUE_72 (const series<complex<rational> > &arg_1, const generic &arg_2) {
    set_variable_name (arg_1, arg_2);
  }
  
  static void
  GLUE_73 (const series<complex<rational> > &arg_1, const int &arg_2) {
    set_output_order (arg_1, arg_2);
  }
  
  static void
  GLUE_74 (const series<complex<rational> > &arg_1, const int &arg_2) {
    set_cancel_order (arg_1, arg_2);
  }
  
  static void
  GLUE_75 (const series<complex<rational> > &arg_1, const bool &arg_2) {
    set_formula_output (arg_1, arg_2);
  }
  
  static series<complex<rational> >
  GLUE_76 (const tuple<complex<rational> > &arg_1) {
    return series<complex<rational> > (as_vector (arg_1));
  }
  
  static series<complex<rational> >
  GLUE_77 (const complex<rational> &arg_1) {
    return series<complex<rational> > (arg_1);
  }
  
  static series<complex<rational> >
  GLUE_78 (const polynomial<complex<rational> > &arg_1) {
    return series<complex<rational> > (arg_1);
  }
  
  static iterator<generic>
  GLUE_79 (const series<complex<rational> > &arg_1) {
    return as<iterator<generic> > (iterate (arg_1));
  }
  
  static complex<rational>
  GLUE_80 (const series<complex<rational> > &arg_1, const int &arg_2) {
    return arg_1[arg_2];
  }
  
  static polynomial<complex<rational> >
  GLUE_81 (const series<complex<rational> > &arg_1, const int &arg_2, const int &arg_3) {
    return range (arg_1, arg_2, arg_3);
  }
  
  static series<complex<rational> >
  GLUE_82 (const series<complex<rational> > &arg_1) {
    return -arg_1;
  }
  
  static series<complex<rational> >
  GLUE_83 (const series<complex<rational> > &arg_1) {
    return square (arg_1);
  }
  
  static series<complex<rational> >
  GLUE_84 (const series<complex<rational> > &arg_1, const series<complex<rational> > &arg_2) {
    return arg_1 + arg_2;
  }
  
  static series<complex<rational> >
  GLUE_85 (const series<complex<rational> > &arg_1, const series<complex<rational> > &arg_2) {
    return arg_1 - arg_2;
  }
  
  static series<complex<rational> >
  GLUE_86 (const series<complex<rational> > &arg_1, const series<complex<rational> > &arg_2) {
    return arg_1 * arg_2;
  }
  
  static series<complex<rational> >
  GLUE_87 (const complex<rational> &arg_1, const series<complex<rational> > &arg_2) {
    return arg_1 + arg_2;
  }
  
  static series<complex<rational> >
  GLUE_88 (const series<complex<rational> > &arg_1, const complex<rational> &arg_2) {
    return arg_1 + arg_2;
  }
  
  static series<complex<rational> >
  GLUE_89 (const complex<rational> &arg_1, const series<complex<rational> > &arg_2) {
    return arg_1 - arg_2;
  }
  
  static series<complex<rational> >
  GLUE_90 (const series<complex<rational> > &arg_1, const complex<rational> &arg_2) {
    return arg_1 - arg_2;
  }
  
  static series<complex<rational> >
  GLUE_91 (const complex<rational> &arg_1, const series<complex<rational> > &arg_2) {
    return arg_1 * arg_2;
  }
  
  static series<complex<rational> >
  GLUE_92 (const series<complex<rational> > &arg_1, const complex<rational> &arg_2) {
    return arg_1 * arg_2;
  }
  
  static series<complex<rational> >
  GLUE_93 (const series<complex<rational> > &arg_1, const int &arg_2) {
    return binpow (arg_1, arg_2);
  }
  
  static series<complex<rational> >
  GLUE_94 (const series<complex<rational> > &arg_1, const integer &arg_2) {
    return binpow (arg_1, arg_2);
  }
  
  static series<complex<rational> >
  GLUE_95 (const series<complex<rational> > &arg_1) {
    return derive (arg_1);
  }
  
  static series<complex<rational> >
  GLUE_96 (const series<complex<rational> > &arg_1) {
    return xderive (arg_1);
  }
  
  static series<complex<rational> >
  GLUE_97 (const series<complex<rational> > &arg_1, const int &arg_2) {
    return dilate (arg_1, arg_2);
  }
  
  static series<complex<rational> >
  GLUE_98 (const series<complex<rational> > &arg_1, const int &arg_2) {
    return lshiftz (arg_1, arg_2);
  }
  
  static series<complex<rational> >
  GLUE_99 (const series<complex<rational> > &arg_1, const int &arg_2) {
    return rshiftz (arg_1, arg_2);
  }
  
  static void
  GLUE_100 (const series<unknown<rational> > &arg_1, const generic &arg_2) {
    set_variable_name (arg_1, arg_2);
  }
  
  static void
  GLUE_101 (const series<unknown<rational> > &arg_1, const int &arg_2) {
    set_output_order (arg_1, arg_2);
  }
  
  static void
  GLUE_102 (const series<unknown<rational> > &arg_1, const int &arg_2) {
    set_cancel_order (arg_1, arg_2);
  }
  
  static void
  GLUE_103 (const series<unknown<rational> > &arg_1, const bool &arg_2) {
    set_formula_output (arg_1, arg_2);
  }
  
  static series<unknown<rational> >
  GLUE_104 (const tuple<unknown<rational> > &arg_1) {
    return series<unknown<rational> > (as_vector (arg_1));
  }
  
  static series<unknown<rational> >
  GLUE_105 (const unknown<rational> &arg_1) {
    return series<unknown<rational> > (arg_1);
  }
  
  static iterator<generic>
  GLUE_106 (const series<unknown<rational> > &arg_1) {
    return as<iterator<generic> > (iterate (arg_1));
  }
  
  static unknown<rational>
  GLUE_107 (const series<unknown<rational> > &arg_1, const int &arg_2) {
    return arg_1[arg_2];
  }
  
  static series<unknown<rational> >
  GLUE_108 (const series<unknown<rational> > &arg_1) {
    return -arg_1;
  }
  
  static series<unknown<rational> >
  GLUE_109 (const series<unknown<rational> > &arg_1) {
    return square (arg_1);
  }
  
  static series<unknown<rational> >
  GLUE_110 (const series<unknown<rational> > &arg_1, const series<unknown<rational> > &arg_2) {
    return arg_1 + arg_2;
  }
  
  static series<unknown<rational> >
  GLUE_111 (const series<unknown<rational> > &arg_1, const series<unknown<rational> > &arg_2) {
    return arg_1 - arg_2;
  }
  
  static series<unknown<rational> >
  GLUE_112 (const series<unknown<rational> > &arg_1, const series<unknown<rational> > &arg_2) {
    return arg_1 * arg_2;
  }
  
  static series<unknown<rational> >
  GLUE_113 (const unknown<rational> &arg_1, const series<unknown<rational> > &arg_2) {
    return arg_1 + arg_2;
  }
  
  static series<unknown<rational> >
  GLUE_114 (const series<unknown<rational> > &arg_1, const unknown<rational> &arg_2) {
    return arg_1 + arg_2;
  }
  
  static series<unknown<rational> >
  GLUE_115 (const unknown<rational> &arg_1, const series<unknown<rational> > &arg_2) {
    return arg_1 - arg_2;
  }
  
  static series<unknown<rational> >
  GLUE_116 (const series<unknown<rational> > &arg_1, const unknown<rational> &arg_2) {
    return arg_1 - arg_2;
  }
  
  static series<unknown<rational> >
  GLUE_117 (const unknown<rational> &arg_1, const series<unknown<rational> > &arg_2) {
    return arg_1 * arg_2;
  }
  
  static series<unknown<rational> >
  GLUE_118 (const series<unknown<rational> > &arg_1, const unknown<rational> &arg_2) {
    return arg_1 * arg_2;
  }
  
  static series<unknown<rational> >
  GLUE_119 (const series<unknown<rational> > &arg_1, const int &arg_2) {
    return binpow (arg_1, arg_2);
  }
  
  static series<unknown<rational> >
  GLUE_120 (const series<unknown<rational> > &arg_1, const integer &arg_2) {
    return binpow (arg_1, arg_2);
  }
  
  static series<unknown<rational> >
  GLUE_121 (const series<unknown<rational> > &arg_1) {
    return derive (arg_1);
  }
  
  static series<unknown<rational> >
  GLUE_122 (const series<unknown<rational> > &arg_1) {
    return xderive (arg_1);
  }
  
  static series<unknown<rational> >
  GLUE_123 (const series<unknown<rational> > &arg_1, const int &arg_2) {
    return dilate (arg_1, arg_2);
  }
  
  static series<unknown<rational> >
  GLUE_124 (const series<unknown<rational> > &arg_1, const int &arg_2) {
    return lshiftz (arg_1, arg_2);
  }
  
  static series<unknown<rational> >
  GLUE_125 (const series<unknown<rational> > &arg_1, const int &arg_2) {
    return rshiftz (arg_1, arg_2);
  }
  
  static series<complex<rational> >
  GLUE_126 (const series<complex<rational> > &arg_1, const series<complex<rational> > &arg_2) {
    return arg_1 / arg_2;
  }
  
  static series<complex<rational> >
  GLUE_127 (const complex<rational> &arg_1, const series<complex<rational> > &arg_2) {
    return arg_1 / arg_2;
  }
  
  static series<complex<rational> >
  GLUE_128 (const series<complex<rational> > &arg_1, const complex<rational> &arg_2) {
    return arg_1 / arg_2;
  }
  
  static series<complex<rational> >
  GLUE_129 (const series<complex<rational> > &arg_1, const series<complex<rational> > &arg_2) {
    return arg_1 / arg_2;
  }
  
  static bool
  GLUE_130 (const series<complex<rational> > &arg_1, const series<complex<rational> > &arg_2) {
    return divides (arg_1, arg_2);
  }
  
  static series<complex<rational> >
  GLUE_131 (const series<complex<rational> > &arg_1, const series<complex<rational> > &arg_2) {
    return gcd (arg_1, arg_2);
  }
  
  static series<complex<rational> >
  GLUE_132 (const series<complex<rational> > &arg_1, const series<complex<rational> > &arg_2) {
    return lcm (arg_1, arg_2);
  }
  
  static series<complex<rational> >
  GLUE_133 (const series<complex<rational> > &arg_1) {
    return integrate (arg_1);
  }
  
  static series<complex<rational> >
  GLUE_134 (const series<complex<rational> > &arg_1, const series<complex<rational> > &arg_2) {
    return compose (arg_1, arg_2);
  }
  
  static series<complex<rational> >
  GLUE_135 (const series<complex<rational> > &arg_1) {
    return reverse (arg_1);
  }
  
  static series<complex<rational> >
  GLUE_136 (const series<complex<rational> > &arg_1, const complex<rational> &arg_2) {
    return q_difference (arg_1, arg_2);
  }
  
  static series<complex<rational> >
  GLUE_137 (const series<complex<rational> > &arg_1, const complex<rational> &arg_2) {
    return series_shift_default (arg_1, arg_2);
  }
  
  static series<complex<rational> >
  GLUE_138 (const series<complex<rational> > &arg_1, const complex<rational> &arg_2, const int &arg_3) {
    return shift (arg_1, arg_2, arg_3);
  }
  
  static series<complex<rational> >
  GLUE_139 (const series<integer> &arg_1) {
    return as<series<complex<rational> > > (arg_1);
  }
  
  static series<unknown<rational> >
  GLUE_140 (const series<integer> &arg_1) {
    return as<series<unknown<rational> > > (arg_1);
  }
  
  static series<complex<rational> >
  GLUE_141 (const series<rational> &arg_1) {
    return as<series<complex<rational> > > (arg_1);
  }
  
  static series<unknown<rational> >
  GLUE_142 (const series<rational> &arg_1) {
    return as<series<unknown<rational> > > (arg_1);
  }
  
  static series<generic>
  GLUE_143 (const series<complex<rational> > &arg_1) {
    return as<series<generic> > (arg_1);
  }
  
  static series<generic>
  GLUE_144 (const series<unknown<rational> > &arg_1) {
    return as<series<generic> > (arg_1);
  }
  
  static series<complex<rational> >
  GLUE_145 (const routine &arg_1, const complex<rational> &arg_2) {
    return fixed_point_series (arg_1, arg_2);
  }
  
  static vector<generic>
  GLUE_146 (const routine &arg_1, const vector<complex<rational> > &arg_2) {
    return gen_fixed_point_vector_series (arg_1, arg_2);
  }
  
  static series<complex<rational> >
  GLUE_147 (const routine &arg_1, const complex<rational> &arg_2) {
    return integrate_series (arg_1, arg_2);
  }
  
  static vector<generic>
  GLUE_148 (const routine &arg_1, const vector<complex<rational> > &arg_2) {
    return gen_integrate_vector_series (arg_1, arg_2);
  }
  
  static series<complex<rational> >
  GLUE_149 (const routine &arg_1, const complex<rational> &arg_2) {
    return implicit_series (arg_1, arg_2);
  }
  
  static vector<generic>
  GLUE_150 (const routine &arg_1, const vector<complex<rational> > &arg_2) {
    return gen_implicit_vector_series (arg_1, arg_2);
  }
  
  void
  glue_series_rational () {
    static bool done = false;
    if (done) return;
    done = true;
    call_glue (string ("glue_polynomial_rational"));
    call_glue (string ("glue_series_integer"));
    define_type<series<rational> > (gen (lit ("Series"), lit ("Rational")));
    define ("set_variable_name", GLUE_1);
    define ("set_output_order", GLUE_2);
    define ("set_cancel_order", GLUE_3);
    define ("set_formula_output", GLUE_4);
    define ("series", GLUE_5);
    define_converter ("upgrade", GLUE_6, PENALTY_INCLUSION);
    define_converter ("upgrade", GLUE_7, PENALTY_INCLUSION);
    define_converter (":>", GLUE_8, PENALTY_CAST);
    define (".[]", GLUE_9);
    define (".[]", GLUE_10);
    define ("-", GLUE_11);
    define ("square", 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 ("^", GLUE_23);
    define ("derive", GLUE_24);
    define ("xderive", GLUE_25);
    define ("dilate", GLUE_26);
    define ("<<", GLUE_27);
    define (">>", GLUE_28);
    define ("/", GLUE_29);
    define ("/", GLUE_30);
    define ("/", GLUE_31);
    define ("div", GLUE_32);
    define ("divides?", GLUE_33);
    define ("gcd", GLUE_34);
    define ("lcm", GLUE_35);
    define ("integrate", GLUE_36);
    define ("@", GLUE_37);
    define ("reverse", GLUE_38);
    define ("q_difference", GLUE_39);
    define ("shift", GLUE_40);
    define ("shift", GLUE_41);
    define ("<=", GLUE_42);
    define (">=", GLUE_43);
    define ("<", GLUE_44);
    define (">", GLUE_45);
    define_converter (":>", GLUE_46, PENALTY_INCLUSION);
    define_converter (":>", GLUE_47, PENALTY_PROMOTE_GENERIC);
    define ("^", GLUE_48);
    define ("sqrt", GLUE_49);
    define ("exp", GLUE_50);
    define ("log", GLUE_51);
    define ("cos", GLUE_52);
    define ("sin", GLUE_53);
    define ("tan", GLUE_54);
    define ("arccos", GLUE_55);
    define ("arcsin", GLUE_56);
    define ("arctan", GLUE_57);
    define_type<unknown<rational> > (gen (lit ("Unknown"), lit ("Rational")));
    define_converter ("upgrade", GLUE_58, PENALTY_INCLUSION);
    define ("-", GLUE_59);
    define ("square", GLUE_60);
    define ("+", GLUE_61);
    define ("-", GLUE_62);
    define ("*", GLUE_63);
    define ("*", GLUE_64);
    define ("*", GLUE_65);
    define ("fixed_point_series", GLUE_66);
    define ("fixed_point_series", GLUE_67);
    define ("integrate_series", GLUE_68);
    define ("integrate_series", GLUE_69);
    define ("implicit_series", GLUE_70);
    define ("implicit_series", GLUE_71);
    define_type<series<complex<rational> > > (gen (lit ("Series"), gen (lit ("Complex"), lit ("Rational"))));
    define_type<series<unknown<rational> > > (gen (lit ("Series"), gen (lit ("Unknown"), lit ("Rational"))));
    define ("set_variable_name", GLUE_72);
    define ("set_output_order", GLUE_73);
    define ("set_cancel_order", GLUE_74);
    define ("set_formula_output", GLUE_75);
    define ("series", GLUE_76);
    define_converter ("upgrade", GLUE_77, PENALTY_INCLUSION);
    define_converter ("upgrade", GLUE_78, PENALTY_INCLUSION);
    define_converter (":>", GLUE_79, PENALTY_CAST);
    define (".[]", GLUE_80);
    define (".[]", GLUE_81);
    define ("-", GLUE_82);
    define ("square", GLUE_83);
    define ("+", GLUE_84);
    define ("-", GLUE_85);
    define ("*", GLUE_86);
    define ("+", GLUE_87);
    define ("+", GLUE_88);
    define ("-", GLUE_89);
    define ("-", GLUE_90);
    define ("*", GLUE_91);
    define ("*", GLUE_92);
    define ("^", GLUE_93);
    define ("^", GLUE_94);
    define ("derive", GLUE_95);
    define ("xderive", GLUE_96);
    define ("dilate", GLUE_97);
    define ("<<", GLUE_98);
    define (">>", GLUE_99);
    define ("set_variable_name", GLUE_100);
    define ("set_output_order", GLUE_101);
    define ("set_cancel_order", GLUE_102);
    define ("set_formula_output", GLUE_103);
    define ("series", GLUE_104);
    define_converter ("upgrade", GLUE_105, PENALTY_INCLUSION);
    define_converter (":>", GLUE_106, PENALTY_CAST);
    define (".[]", GLUE_107);
    define ("-", GLUE_108);
    define ("square", GLUE_109);
    define ("+", GLUE_110);
    define ("-", GLUE_111);
    define ("*", GLUE_112);
    define ("+", GLUE_113);
    define ("+", GLUE_114);
    define ("-", GLUE_115);
    define ("-", GLUE_116);
    define ("*", GLUE_117);
    define ("*", GLUE_118);
    define ("^", GLUE_119);
    define ("^", GLUE_120);
    define ("derive", GLUE_121);
    define ("xderive", GLUE_122);
    define ("dilate", GLUE_123);
    define ("<<", GLUE_124);
    define (">>", GLUE_125);
    define ("/", GLUE_126);
    define ("/", GLUE_127);
    define ("/", GLUE_128);
    define ("div", GLUE_129);
    define ("divides?", GLUE_130);
    define ("gcd", GLUE_131);
    define ("lcm", GLUE_132);
    define ("integrate", GLUE_133);
    define ("@", GLUE_134);
    define ("reverse", GLUE_135);
    define ("q_difference", GLUE_136);
    define ("shift", GLUE_137);
    define ("shift", GLUE_138);
    define_converter (":>", GLUE_139, PENALTY_INCLUSION);
    define_converter (":>", GLUE_140, PENALTY_INCLUSION);
    define_converter (":>", GLUE_141, PENALTY_INCLUSION);
    define_converter (":>", GLUE_142, PENALTY_INCLUSION);
    define_converter (":>", GLUE_143, PENALTY_PROMOTE_GENERIC);
    define_converter (":>", GLUE_144, PENALTY_PROMOTE_GENERIC);
    define ("fixed_point_series", GLUE_145);
    define ("fixed_point_series", GLUE_146);
    define ("integrate_series", GLUE_147);
    define ("integrate_series", GLUE_148);
    define ("implicit_series", GLUE_149);
    define ("implicit_series", GLUE_150);
  }
}