decomposer.h

   1 #ifndef DECOMPOSER_H
   2 #define DECOMPOSER_H
   3
   4 #include <gmp.h>
   5 #include <barvinok/polylib.h>
   6 #include <barvinok/options.h>
   7
   8 struct signed_cone {
   9     signed_cone(const mat_ZZ& rays, int sign, unsigned long det,
  10                 int *closed = NULL) :
  11                         C(NULL), rays(rays), sign(sign), det(det), closed(closed) {}
  12     signed_cone(Polyhedron *C, const mat_ZZ& rays, int sign, unsigned long det = 0,
  13                 int *closed = NULL) :
  14                         C(C), rays(rays), sign(sign), det(det), closed(closed) {}
  15     Polyhedron *C;
  16     const mat_ZZ& rays;
  17     int sign;
  18     unsigned long det;
  19     /* facet not containing ray is closed */
  20     int *closed;
  21 };
  22
  23 struct signed_cone_consumer {
  24     virtual void handle(const signed_cone& sc, barvinok_options *options) = 0;
  25 };
  26
  27 struct vertex_decomposer {
  28     Polyhedron *P;
  29     unsigned nbV;       // number of vertices
  30     Param_Vertices *V;  // current vertex
  31     int vert;           // current vertex index
  32     signed_cone_consumer& scc;
  33
  34     vertex_decomposer(Polyhedron *P, unsigned nbV, signed_cone_consumer& scc) :
  35                         P(P), nbV(nbV), scc(scc) {}
  36     void decompose_at_vertex(Param_Vertices *V, int _i, barvinok_options *options);
  37 };
  38
  39 void barvinok_decompose(Polyhedron *C, signed_cone_consumer& scc,
  40                         barvinok_options *options);
  41
  42 #endif