verify_series.cc

   1 #include <barvinok/util.h>
   2 #include "skewed_genfun.h"
   3 #include "verify_series.h"
   4
   5 struct check_poly_gf_data {
   6     struct check_poly_data   cp;
   7     Polyhedron              *S;
   8     const skewed_gen_fun    *gf;
   9     int                      exist;
  10 };
  11
  12 static int cp_gf(const struct check_poly_data *data, int nparam, Value *z,
  13                  const struct verify_options *options)
  14 {
  15     Value c, tmp, one;
  16     struct check_poly_gf_data* gf_data = (struct check_poly_gf_data*) data;
  17     const skewed_gen_fun *gf = gf_data->gf;
  18     int exist = gf_data->exist;
  19     Polyhedron *S = gf_data->S;
  20
  21     value_init(c);
  22     value_init(tmp);
  23     value_init(one);
  24     value_set_si(one, 1);
  25
  26     /* Computes the coefficient */
  27     gf->coefficient(z, &c, options->barvinok);
  28
  29     /* Manually count the number of points */
  30     count_points_e(1, S, exist, nparam, data->z, &tmp);
  31
  32     check_poly_print(value_eq(tmp, c), nparam, z, tmp, one, c, one,
  33                      "EP", "count", "EP eval", options);
  34
  35     value_clear(c);
  36     value_clear(tmp);
  37     value_clear(one);
  38     return 1;
  39 }
  40
  41 int check_poly_gf(Polyhedron *S, Polyhedron *CS, skewed_gen_fun *gf,
  42                     int exist, int nparam, int pos, Value *z,
  43                     const struct verify_options *options)
  44 {
  45     struct check_poly_gf_data data;
  46     data.cp.z = z;
  47     data.cp.check = cp_gf;
  48     data.S = S;
  49     data.gf = gf;
  50     data.exist = exist;
  51     return check_poly(CS, &data.cp, nparam, pos, z+S->Dimension-nparam+1,
  52                         options);
  53 }