polyhedron_sample.c

   1 #include <polylib/polylibgmp.h>
   2 #include <barvinok/sample.h>
   3 #include "config.h"
   4
   5 #ifdef HAVE_GROWING_CHERNIKOVA
   6 #define MAXRAYS    (POL_NO_DUAL | POL_INTEGER)
   7 #else
   8 #define MAXRAYS  600
   9 #endif
  10
  11 int main(int argc, char **argv)
  12 {
  13     Matrix *M;
  14     Polyhedron *P;
  15     Vector *sample;
  16
  17     M = Matrix_Read();
  18     P = Constraints2Polyhedron(M, MAXRAYS);
  19     Matrix_Free(M);
  20
  21     sample = Polyhedron_Sample(P, MAXRAYS);
  22     if (sample) {
  23         Vector_Print(stdout, P_VALUE_FMT, sample);
  24         Vector_Free(sample);
  25     }
  26
  27     Polyhedron_Free(P);
  28 }