polyhedron_sample.c

   1 #include <barvinok/options.h>
   2 #include <barvinok/sample.h>
   3 #include "config.h"
   4 #include "argp.h"
   5 #include "progname.h"
   6
   7 int main(int argc, char **argv)
   8 {
   9     Matrix *M;
  10     Polyhedron *P;
  11     Vector *sample;
  12     struct barvinok_options *bv_options = barvinok_options_new_with_defaults();
  13
  14     set_program_name(argv[0]);
  15     argp_parse(&barvinok_argp, argc, argv, 0, 0, bv_options);
  16
  17     M = Matrix_Read();
  18     P = Constraints2Polyhedron(M, bv_options->MaxRays);
  19     Matrix_Free(M);
  20
  21     sample = Polyhedron_Sample(P, bv_options);
  22     if (sample) {
  23         assert(in_domain(P, sample->p));
  24         Vector_Print(stdout, P_VALUE_FMT, sample);
  25         Vector_Free(sample);
  26     }
  27
  28     Polyhedron_Free(P);
  29     barvinok_options_free(bv_options);
  30 }