dpoly.h

   1 #ifndef DPOLY_H
   2 #define DPOLY_H
   3
   4 #include <assert.h>
   5 #include <vector>
   6 #include <gmp.h>
   7 #include <NTL/vec_ZZ.h>
   8 extern "C" {
   9 #include <polylib/polylibgmp.h>
  10 }
  11 #include "conversion.h"
  12
  13 #ifdef NTL_STD_CXX
  14 using namespace NTL;
  15 #endif
  16
  17 class dpoly {
  18 public:
  19     vec_ZZ coeff;
  20     dpoly(int d, ZZ& degree, int offset = 0);
  21     void operator *= (dpoly& f);
  22     void div(dpoly& d, mpq_t count, ZZ& sign);
  23     void div(dpoly& d, mpq_t *count, const mpq_t& factor);
  24 private:
  25     mpq_t *div(dpoly &d) const;
  26     void   clear_div(mpq_t *c) const;
  27 };
  28
  29 struct dpoly_r_term {
  30     int     *powers;
  31     ZZ      coeff;
  32 };
  33
  34 /* len: number of elements in c
  35  * each element in c is the coefficient of a power of t
  36  * in the MacLaurin expansion
  37  */
  38 struct dpoly_r {
  39     std::vector< dpoly_r_term * >       *c;
  40     int len;
  41     int dim;
  42     ZZ denom;
  43
  44     void add_term(int i, int * powers, ZZ& coeff);
  45     dpoly_r(int len, int dim);
  46     dpoly_r(dpoly& num, int dim);
  47     dpoly_r(dpoly& num, dpoly& den, int pos, int dim);
  48     dpoly_r(dpoly_r* num, dpoly& den, int pos, int dim);
  49     ~dpoly_r();
  50     dpoly_r *div(dpoly& d);
  51     void dump(void);
  52 };
  53
  54 #endif