polyhedron_sample.c

   1 #include <assert.h>
   2 #include <barvinok/options.h>
   3 #include <barvinok/sample.h>
   4 #include "config.h"
   5 #include "argp.h"
   6 #include "progname.h"
   7
   8 int main(int argc, char **argv)
   9 {
  10     Matrix *M;
  11     Polyhedron *P;
  12     Vector *sample;
  13     struct barvinok_options *bv_options = barvinok_options_new_with_defaults();
  14
  15     set_program_name(argv[0]);
  16     argp_parse(&barvinok_argp, argc, argv, 0, 0, bv_options);
  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     barvinok_options_free(bv_options);
  31 }