polyhedron_sample.c

   1 #include <barvinok/options.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     struct barvinok_options *bv_options = barvinok_options_new_with_defaults();
  17
  18     M = Matrix_Read();
  19     P = Constraints2Polyhedron(M, bv_options->MaxRays);
  20     Matrix_Free(M);
  21
  22     sample = Polyhedron_Sample(P, bv_options);
  23     if (sample) {
  24         assert(in_domain(P, sample->p));
  25         Vector_Print(stdout, P_VALUE_FMT, sample);
  26         Vector_Free(sample);
  27     }
  28
  29     Polyhedron_Free(P);
  30     free(bv_options);
  31 }