/Users/mourrain/Devel/mmx/algebramix/glue/glue_matrix_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 <basix/int.hpp>
#include <algebramix/vector_unrolled.hpp>
#include <algebramix/vector_simd.hpp>
#include <algebramix/permutation.hpp>
#include <basix/row_tuple.hpp>
#include <algebramix/matrix.hpp>
#include <algebramix/matrix_ring_naive.hpp>
#include <algebramix/matrix_integer.hpp>
#include <algebramix/matrix_quotient.hpp>
#include <algebramix/matrix_complex.hpp>
#include <basix/tuple.hpp>
#include <basix/alias.hpp>
#include <basix/glue.hpp>

#define int_literal(x) as_int (as_string (x))
#define identity_matrix_integer identity_matrix<integer>
#define hilbert_matrix_rational hilbert_matrix<rational>

namespace mmx {
  template<typename C> matrix<C>
  matrix_new (const vector<row_tuple<C> >& t) {
    if (N(t) == 0) return matrix<C> ();
    nat i, j, rows= N(t), cols= N(t[0]);
    C dummy= 0;
    matrix<C> r (dummy, rows, cols);
    for (i=0; i<rows; i++) {
      ASSERT (N(t[i]) == cols, "unequal row lengths");
      for (j=0; j<cols; j++)
        r(i,j)= t[i][j];
    }
    return r;
  }

  template<typename C> matrix<C>
  matrix_new (const vector<C>& t) {
    nat j, rows= 1, cols= N(t);
    C dummy= 0;
    matrix<C> r (dummy, rows, cols);
    for (j=0; j<cols; j++)
      r(0,j)= t[j];
    return r;
  }

  template<typename C> vector<generic>
  wrap_column_reduced_echelon_with_permutation (const matrix<C>& m) {
    permutation permut;
    generic tp=as<generic> (column_reduced_echelon (m, permut));
    return vec (tp, as<generic> (permut));
  }

  template<typename C> vector<generic>
  wrap_column_reduced_echelon_with_transform (const matrix<C>& m) {
    matrix<C> k;
    generic tp=as<generic> (column_reduced_echelon (m, k));
    return vec (tp, as<generic> (k));
  }

  template<typename C> vector<generic>
  wrap_row_reduced_echelon_with_transform (const matrix<C>& m) {
    matrix<C> k;
    generic tp=as<generic> (row_reduced_echelon (m, k));
    return vec (tp, as<generic> (k));
  }
}


namespace mmx {
  static vector<rational>
  GLUE_1 (const vector<int> &arg_1) {
    return as<vector<rational> > (arg_1);
  }
  
  static vector<complex<rational> >
  GLUE_2 (const vector<int> &arg_1) {
    return as<vector<complex<rational> > > (arg_1);
  }
  
  static matrix<rational>
  GLUE_3 (const int &arg_1) {
    return hilbert_matrix_rational (arg_1);
  }
  
  static row_tuple<rational>
  GLUE_4 (const tuple<rational> &arg_1) {
    return row_tuple<rational > (as_vector (arg_1));
  }
  
  static matrix<rational>
  GLUE_5 (const tuple<rational> &arg_1) {
    return matrix_new (as_vector (arg_1));
  }
  
  static matrix<rational>
  GLUE_6 (const tuple<row_tuple<rational> > &arg_1) {
    return matrix_new (as_vector (arg_1));
  }
  
  static matrix<rational>
  GLUE_7 (const tuple<row_tuple<rational> > &arg_1) {
    return matrix_new (as_vector (arg_1));
  }
  
  static int
  GLUE_8 (const matrix<rational> &arg_1) {
    return N (arg_1);
  }
  
  static int
  GLUE_9 (const matrix<rational> &arg_1) {
    return rows (arg_1);
  }
  
  static int
  GLUE_10 (const matrix<rational> &arg_1) {
    return cols (arg_1);
  }
  
  static rational
  GLUE_11 (const matrix<rational> &arg_1, const int &arg_2, const int &arg_3) {
    return arg_1 (arg_2, arg_3);
  }
  
  static alias<rational>
  GLUE_12 (const alias<matrix<rational> > &arg_1, const int &arg_2, const int &arg_3) {
    return alias_access<rational > (arg_1, arg_2, arg_3);
  }
  
  static matrix<rational>
  GLUE_13 (const matrix<rational> &arg_1, const int &arg_2, const int &arg_3, const int &arg_4, const int &arg_5) {
    return range (arg_1, arg_2, arg_3, arg_4, arg_5);
  }
  
  static vector<rational>
  GLUE_14 (const matrix<rational> &arg_1, const int &arg_2) {
    return row (arg_1, arg_2);
  }
  
  static vector<rational>
  GLUE_15 (const matrix<rational> &arg_1, const int &arg_2) {
    return column (arg_1, arg_2);
  }
  
  static matrix<rational>
  GLUE_16 (const matrix<rational> &arg_1) {
    return transpose (arg_1);
  }
  
  static matrix<rational>
  GLUE_17 (const matrix<rational> &arg_1, const matrix<rational> &arg_2) {
    return horizontal_join (arg_1, arg_2);
  }
  
  static matrix<rational>
  GLUE_18 (const matrix<rational> &arg_1, const matrix<rational> &arg_2) {
    return vertical_join (arg_1, arg_2);
  }
  
  static row_tuple<complex<rational> >
  GLUE_19 (const tuple<complex<rational> > &arg_1) {
    return row_tuple<complex<rational> > (as_vector (arg_1));
  }
  
  static matrix<complex<rational> >
  GLUE_20 (const tuple<complex<rational> > &arg_1) {
    return matrix_new (as_vector (arg_1));
  }
  
  static matrix<complex<rational> >
  GLUE_21 (const tuple<row_tuple<complex<rational> > > &arg_1) {
    return matrix_new (as_vector (arg_1));
  }
  
  static matrix<complex<rational> >
  GLUE_22 (const tuple<row_tuple<complex<rational> > > &arg_1) {
    return matrix_new (as_vector (arg_1));
  }
  
  static int
  GLUE_23 (const matrix<complex<rational> > &arg_1) {
    return N (arg_1);
  }
  
  static int
  GLUE_24 (const matrix<complex<rational> > &arg_1) {
    return rows (arg_1);
  }
  
  static int
  GLUE_25 (const matrix<complex<rational> > &arg_1) {
    return cols (arg_1);
  }
  
  static complex<rational>
  GLUE_26 (const matrix<complex<rational> > &arg_1, const int &arg_2, const int &arg_3) {
    return arg_1 (arg_2, arg_3);
  }
  
  static alias<complex<rational> >
  GLUE_27 (const alias<matrix<complex<rational> > > &arg_1, const int &arg_2, const int &arg_3) {
    return alias_access<complex<rational> > (arg_1, arg_2, arg_3);
  }
  
  static matrix<complex<rational> >
  GLUE_28 (const matrix<complex<rational> > &arg_1, const int &arg_2, const int &arg_3, const int &arg_4, const int &arg_5) {
    return range (arg_1, arg_2, arg_3, arg_4, arg_5);
  }
  
  static vector<complex<rational> >
  GLUE_29 (const matrix<complex<rational> > &arg_1, const int &arg_2) {
    return row (arg_1, arg_2);
  }
  
  static vector<complex<rational> >
  GLUE_30 (const matrix<complex<rational> > &arg_1, const int &arg_2) {
    return column (arg_1, arg_2);
  }
  
  static matrix<complex<rational> >
  GLUE_31 (const matrix<complex<rational> > &arg_1) {
    return transpose (arg_1);
  }
  
  static matrix<complex<rational> >
  GLUE_32 (const matrix<complex<rational> > &arg_1, const matrix<complex<rational> > &arg_2) {
    return horizontal_join (arg_1, arg_2);
  }
  
  static matrix<complex<rational> >
  GLUE_33 (const matrix<complex<rational> > &arg_1, const matrix<complex<rational> > &arg_2) {
    return vertical_join (arg_1, arg_2);
  }
  
  static matrix<rational>
  GLUE_34 (const rational &arg_1, const int &arg_2, const int &arg_3) {
    return fill_matrix (arg_1, arg_2, arg_3);
  }
  
  static matrix<rational>
  GLUE_35 (const rational &arg_1, const int &arg_2) {
    return jordan_matrix (arg_1, arg_2);
  }
  
  static matrix<rational>
  GLUE_36 (const vector<rational> &arg_1) {
    return toeplitz_matrix (arg_1);
  }
  
  static matrix<rational>
  GLUE_37 (const vector<rational> &arg_1) {
    return hankel_matrix (arg_1);
  }
  
  static matrix<rational>
  GLUE_38 (const vector<rational> &arg_1, const vector<rational> &arg_2) {
    return tensor_matrix (arg_1, arg_2);
  }
  
  static matrix<rational>
  GLUE_39 (const vector<rational> &arg_1) {
    return vandermonde (arg_1);
  }
  
  static matrix<rational>
  GLUE_40 (const matrix<rational> &arg_1) {
    return -arg_1;
  }
  
  static matrix<rational>
  GLUE_41 (const matrix<rational> &arg_1) {
    return square (arg_1);
  }
  
  static matrix<rational>
  GLUE_42 (const matrix<rational> &arg_1, const matrix<rational> &arg_2) {
    return arg_1 + arg_2;
  }
  
  static matrix<rational>
  GLUE_43 (const matrix<rational> &arg_1, const matrix<rational> &arg_2) {
    return arg_1 - arg_2;
  }
  
  static matrix<rational>
  GLUE_44 (const matrix<rational> &arg_1, const matrix<rational> &arg_2) {
    return arg_1 * arg_2;
  }
  
  static matrix<rational>
  GLUE_45 (const rational &arg_1, const matrix<rational> &arg_2) {
    return arg_1 + arg_2;
  }
  
  static matrix<rational>
  GLUE_46 (const matrix<rational> &arg_1, const rational &arg_2) {
    return arg_1 + arg_2;
  }
  
  static matrix<rational>
  GLUE_47 (const rational &arg_1, const matrix<rational> &arg_2) {
    return arg_1 - arg_2;
  }
  
  static matrix<rational>
  GLUE_48 (const matrix<rational> &arg_1, const rational &arg_2) {
    return arg_1 - arg_2;
  }
  
  static matrix<rational>
  GLUE_49 (const rational &arg_1, const matrix<rational> &arg_2) {
    return arg_1 * arg_2;
  }
  
  static matrix<rational>
  GLUE_50 (const matrix<rational> &arg_1, const rational &arg_2) {
    return arg_1 * arg_2;
  }
  
  static vector<rational>
  GLUE_51 (const matrix<rational> &arg_1, const vector<rational> &arg_2) {
    return arg_1 * arg_2;
  }
  
  static vector<rational>
  GLUE_52 (const vector<rational> &arg_1, const matrix<rational> &arg_2) {
    return arg_1 * arg_2;
  }
  
  static matrix<rational>
  GLUE_53 (const rational &arg_1, const matrix<rational> &arg_2) {
    return arg_1 / arg_2;
  }
  
  static matrix<rational>
  GLUE_54 (const matrix<rational> &arg_1, const rational &arg_2) {
    return arg_1 / arg_2;
  }
  
  static matrix<rational>
  GLUE_55 (const matrix<rational> &arg_1, const matrix<rational> &arg_2) {
    return krylov (arg_1, arg_2);
  }
  
  static rational
  GLUE_56 (const matrix<rational> &arg_1) {
    return det (arg_1);
  }
  
  static matrix<rational>
  GLUE_57 (const matrix<rational> &arg_1) {
    return row_echelon (arg_1);
  }
  
  static matrix<rational>
  GLUE_58 (const matrix<rational> &arg_1) {
    return column_echelon (arg_1);
  }
  
  static matrix<rational>
  GLUE_59 (const matrix<rational> &arg_1) {
    return row_reduced_echelon (arg_1);
  }
  
  static matrix<rational>
  GLUE_60 (const matrix<rational> &arg_1) {
    return column_reduced_echelon (arg_1);
  }
  
  static vector<generic>
  GLUE_61 (const matrix<rational> &arg_1) {
    return wrap_row_reduced_echelon_with_transform (arg_1);
  }
  
  static vector<generic>
  GLUE_62 (const matrix<rational> &arg_1) {
    return wrap_column_reduced_echelon_with_transform (arg_1);
  }
  
  static vector<generic>
  GLUE_63 (const matrix<rational> &arg_1) {
    return wrap_column_reduced_echelon_with_permutation (arg_1);
  }
  
  static matrix<rational>
  GLUE_64 (const matrix<rational> &arg_1) {
    return kernel (arg_1);
  }
  
  static matrix<rational>
  GLUE_65 (const matrix<rational> &arg_1) {
    return image (arg_1);
  }
  
  static int
  GLUE_66 (const matrix<rational> &arg_1) {
    return rank (arg_1);
  }
  
  static matrix<rational>
  GLUE_67 (const matrix<rational> &arg_1) {
    return invert (arg_1);
  }
  
  static matrix<rational>
  GLUE_68 (const matrix<rational> &arg_1) {
    return abs (arg_1);
  }
  
  static row_tuple<rational>
  GLUE_69 (const row_tuple<integer> &arg_1) {
    return as<row_tuple<rational> > (arg_1);
  }
  
  static matrix<rational>
  GLUE_70 (const matrix<integer> &arg_1) {
    return as<matrix<rational> > (arg_1);
  }
  
  static row_tuple<complex<rational> >
  GLUE_71 (const row_tuple<integer> &arg_1) {
    return as<row_tuple<complex<rational> > > (arg_1);
  }
  
  static matrix<complex<rational> >
  GLUE_72 (const matrix<integer> &arg_1) {
    return as<matrix<complex<rational> > > (arg_1);
  }
  
  static row_tuple<generic>
  GLUE_73 (const row_tuple<rational> &arg_1) {
    return as<row_tuple<generic> > (arg_1);
  }
  
  static matrix<generic>
  GLUE_74 (const matrix<rational> &arg_1) {
    return as<matrix<generic> > (arg_1);
  }
  
  static row_tuple<complex<rational> >
  GLUE_75 (const row_tuple<rational> &arg_1) {
    return as<row_tuple<complex<rational> > > (arg_1);
  }
  
  static matrix<complex<rational> >
  GLUE_76 (const matrix<rational> &arg_1) {
    return as<matrix<complex<rational> > > (arg_1);
  }
  
  static row_tuple<generic>
  GLUE_77 (const row_tuple<complex<rational> > &arg_1) {
    return as<row_tuple<generic> > (arg_1);
  }
  
  static matrix<generic>
  GLUE_78 (const matrix<complex<rational> > &arg_1) {
    return as<matrix<generic> > (arg_1);
  }
  
  static matrix<complex<rational> >
  GLUE_79 (const complex<rational> &arg_1, const int &arg_2, const int &arg_3) {
    return fill_matrix (arg_1, arg_2, arg_3);
  }
  
  static matrix<complex<rational> >
  GLUE_80 (const complex<rational> &arg_1, const int &arg_2) {
    return jordan_matrix (arg_1, arg_2);
  }
  
  static matrix<complex<rational> >
  GLUE_81 (const vector<complex<rational> > &arg_1) {
    return toeplitz_matrix (arg_1);
  }
  
  static matrix<complex<rational> >
  GLUE_82 (const vector<complex<rational> > &arg_1) {
    return hankel_matrix (arg_1);
  }
  
  static matrix<complex<rational> >
  GLUE_83 (const vector<complex<rational> > &arg_1, const vector<complex<rational> > &arg_2) {
    return tensor_matrix (arg_1, arg_2);
  }
  
  static matrix<complex<rational> >
  GLUE_84 (const vector<complex<rational> > &arg_1) {
    return vandermonde (arg_1);
  }
  
  static matrix<complex<rational> >
  GLUE_85 (const matrix<complex<rational> > &arg_1) {
    return -arg_1;
  }
  
  static matrix<complex<rational> >
  GLUE_86 (const matrix<complex<rational> > &arg_1) {
    return square (arg_1);
  }
  
  static matrix<complex<rational> >
  GLUE_87 (const matrix<complex<rational> > &arg_1, const matrix<complex<rational> > &arg_2) {
    return arg_1 + arg_2;
  }
  
  static matrix<complex<rational> >
  GLUE_88 (const matrix<complex<rational> > &arg_1, const matrix<complex<rational> > &arg_2) {
    return arg_1 - arg_2;
  }
  
  static matrix<complex<rational> >
  GLUE_89 (const matrix<complex<rational> > &arg_1, const matrix<complex<rational> > &arg_2) {
    return arg_1 * arg_2;
  }
  
  static matrix<complex<rational> >
  GLUE_90 (const complex<rational> &arg_1, const matrix<complex<rational> > &arg_2) {
    return arg_1 + arg_2;
  }
  
  static matrix<complex<rational> >
  GLUE_91 (const matrix<complex<rational> > &arg_1, const complex<rational> &arg_2) {
    return arg_1 + arg_2;
  }
  
  static matrix<complex<rational> >
  GLUE_92 (const complex<rational> &arg_1, const matrix<complex<rational> > &arg_2) {
    return arg_1 - arg_2;
  }
  
  static matrix<complex<rational> >
  GLUE_93 (const matrix<complex<rational> > &arg_1, const complex<rational> &arg_2) {
    return arg_1 - arg_2;
  }
  
  static matrix<complex<rational> >
  GLUE_94 (const complex<rational> &arg_1, const matrix<complex<rational> > &arg_2) {
    return arg_1 * arg_2;
  }
  
  static matrix<complex<rational> >
  GLUE_95 (const matrix<complex<rational> > &arg_1, const complex<rational> &arg_2) {
    return arg_1 * arg_2;
  }
  
  static vector<complex<rational> >
  GLUE_96 (const matrix<complex<rational> > &arg_1, const vector<complex<rational> > &arg_2) {
    return arg_1 * arg_2;
  }
  
  static vector<complex<rational> >
  GLUE_97 (const vector<complex<rational> > &arg_1, const matrix<complex<rational> > &arg_2) {
    return arg_1 * arg_2;
  }
  
  static matrix<complex<rational> >
  GLUE_98 (const complex<rational> &arg_1, const matrix<complex<rational> > &arg_2) {
    return arg_1 / arg_2;
  }
  
  static matrix<complex<rational> >
  GLUE_99 (const matrix<complex<rational> > &arg_1, const complex<rational> &arg_2) {
    return arg_1 / arg_2;
  }
  
  static matrix<complex<rational> >
  GLUE_100 (const matrix<complex<rational> > &arg_1, const matrix<complex<rational> > &arg_2) {
    return krylov (arg_1, arg_2);
  }
  
  static complex<rational>
  GLUE_101 (const matrix<complex<rational> > &arg_1) {
    return det (arg_1);
  }
  
  static matrix<complex<rational> >
  GLUE_102 (const matrix<complex<rational> > &arg_1) {
    return row_echelon (arg_1);
  }
  
  static matrix<complex<rational> >
  GLUE_103 (const matrix<complex<rational> > &arg_1) {
    return column_echelon (arg_1);
  }
  
  static matrix<complex<rational> >
  GLUE_104 (const matrix<complex<rational> > &arg_1) {
    return row_reduced_echelon (arg_1);
  }
  
  static matrix<complex<rational> >
  GLUE_105 (const matrix<complex<rational> > &arg_1) {
    return column_reduced_echelon (arg_1);
  }
  
  static vector<generic>
  GLUE_106 (const matrix<complex<rational> > &arg_1) {
    return wrap_row_reduced_echelon_with_transform (arg_1);
  }
  
  static vector<generic>
  GLUE_107 (const matrix<complex<rational> > &arg_1) {
    return wrap_column_reduced_echelon_with_transform (arg_1);
  }
  
  static vector<generic>
  GLUE_108 (const matrix<complex<rational> > &arg_1) {
    return wrap_column_reduced_echelon_with_permutation (arg_1);
  }
  
  static matrix<complex<rational> >
  GLUE_109 (const matrix<complex<rational> > &arg_1) {
    return kernel (arg_1);
  }
  
  static matrix<complex<rational> >
  GLUE_110 (const matrix<complex<rational> > &arg_1) {
    return image (arg_1);
  }
  
  static int
  GLUE_111 (const matrix<complex<rational> > &arg_1) {
    return rank (arg_1);
  }
  
  static matrix<complex<rational> >
  GLUE_112 (const matrix<complex<rational> > &arg_1) {
    return invert (arg_1);
  }
  
  void
  glue_matrix_rational () {
    static bool done = false;
    if (done) return;
    done = true;
    call_glue (string ("glue_vector_rational"));
    call_glue (string ("glue_matrix_integer"));
    define_converter (":>", GLUE_1, PENALTY_CAST);
    define_converter (":>", GLUE_2, PENALTY_CAST);
    define_type<row_tuple<rational> > (gen (lit ("Row"), lit ("Rational")));
    define_type<row_tuple<complex<rational> > > (gen (lit ("Row"), gen (lit ("Complex"), lit ("Rational"))));
    define_type<matrix<rational> > (gen (lit ("Matrix"), lit ("Rational")));
    define_type<matrix<complex<rational> > > (gen (lit ("Matrix"), gen (lit ("Complex"), lit ("Rational"))));
    define ("hilbert_matrix", GLUE_3);
    define ("(.)", GLUE_4);
    define ("matrix", GLUE_5);
    define ("matrix", GLUE_6);
    define ("[]", GLUE_7);
    define ("#", GLUE_8);
    define ("rows", GLUE_9);
    define ("columns", GLUE_10);
    define (".[]", GLUE_11);
    define (".[]", GLUE_12);
    define (".[]", GLUE_13);
    define ("row", GLUE_14);
    define ("column", GLUE_15);
    define ("transpose", GLUE_16);
    define ("horizontal_join", GLUE_17);
    define ("vertical_join", GLUE_18);
    define ("(.)", GLUE_19);
    define ("matrix", GLUE_20);
    define ("matrix", GLUE_21);
    define ("[]", GLUE_22);
    define ("#", GLUE_23);
    define ("rows", GLUE_24);
    define ("columns", GLUE_25);
    define (".[]", GLUE_26);
    define (".[]", GLUE_27);
    define (".[]", GLUE_28);
    define ("row", GLUE_29);
    define ("column", GLUE_30);
    define ("transpose", GLUE_31);
    define ("horizontal_join", GLUE_32);
    define ("vertical_join", GLUE_33);
    define ("fill_matrix", GLUE_34);
    define ("jordan_matrix", GLUE_35);
    define ("toeplitz_matrix", GLUE_36);
    define ("hankel_matrix", GLUE_37);
    define ("tensor_matrix", GLUE_38);
    define ("vandermonde", GLUE_39);
    define ("-", GLUE_40);
    define ("square", GLUE_41);
    define ("+", GLUE_42);
    define ("-", GLUE_43);
    define ("*", GLUE_44);
    define ("+", GLUE_45);
    define ("+", GLUE_46);
    define ("-", GLUE_47);
    define ("-", GLUE_48);
    define ("*", GLUE_49);
    define ("*", GLUE_50);
    define ("*", GLUE_51);
    define ("*", GLUE_52);
    define ("/", GLUE_53);
    define ("/", GLUE_54);
    define ("krylov", GLUE_55);
    define ("det", GLUE_56);
    define ("row_echelon", GLUE_57);
    define ("column_echelon", GLUE_58);
    define ("row_reduced_echelon", GLUE_59);
    define ("column_reduced_echelon", GLUE_60);
    define ("row_reduced_echelon_with_transform", GLUE_61);
    define ("column_reduced_echelon_with_transform", GLUE_62);
    define ("column_reduced_echelon_with_permutation", GLUE_63);
    define ("ker", GLUE_64);
    define ("im", GLUE_65);
    define ("rank", GLUE_66);
    define ("invert", GLUE_67);
    define ("abs", GLUE_68);
    define_converter (":>", GLUE_69, PENALTY_INCLUSION);
    define_converter (":>", GLUE_70, PENALTY_INCLUSION);
    define_converter (":>", GLUE_71, PENALTY_INCLUSION);
    define_converter (":>", GLUE_72, PENALTY_INCLUSION);
    define_converter (":>", GLUE_73, PENALTY_PROMOTE_GENERIC);
    define_converter (":>", GLUE_74, PENALTY_PROMOTE_GENERIC);
    define_converter (":>", GLUE_75, PENALTY_INCLUSION);
    define_converter (":>", GLUE_76, PENALTY_INCLUSION);
    define_converter (":>", GLUE_77, PENALTY_PROMOTE_GENERIC);
    define_converter (":>", GLUE_78, PENALTY_PROMOTE_GENERIC);
    define ("fill_matrix", GLUE_79);
    define ("jordan_matrix", GLUE_80);
    define ("toeplitz_matrix", GLUE_81);
    define ("hankel_matrix", GLUE_82);
    define ("tensor_matrix", GLUE_83);
    define ("vandermonde", GLUE_84);
    define ("-", GLUE_85);
    define ("square", 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 ("*", GLUE_95);
    define ("*", GLUE_96);
    define ("*", GLUE_97);
    define ("/", GLUE_98);
    define ("/", GLUE_99);
    define ("krylov", GLUE_100);
    define ("det", GLUE_101);
    define ("row_echelon", GLUE_102);
    define ("column_echelon", GLUE_103);
    define ("row_reduced_echelon", GLUE_104);
    define ("column_reduced_echelon", GLUE_105);
    define ("row_reduced_echelon_with_transform", GLUE_106);
    define ("column_reduced_echelon_with_transform", GLUE_107);
    define ("column_reduced_echelon_with_permutation", GLUE_108);
    define ("ker", GLUE_109);
    define ("im", GLUE_110);
    define ("rank", GLUE_111);
    define ("invert", GLUE_112);
  }
}