param_util.c

   1 #include <barvinok/options.h>
   2 #include "param_util.h"
   3 #include "topcom.h"
   4 #include "config.h"
   5
   6 void Param_Vertex_Common_Denominator(Param_Vertices *V)
   7 {
   8     unsigned dim;
   9     Value lcm;
  10     int i;
  11
  12     assert(V->Vertex->NbRows > 0);
  13     dim = V->Vertex->NbColumns-2;
  14
  15     value_init(lcm);
  16
  17     value_assign(lcm, V->Vertex->p[0][dim+1]);
  18     for (i = 1; i < V->Vertex->NbRows; ++i)
  19         value_lcm(V->Vertex->p[i][dim+1], lcm, &lcm);
  20
  21     for (i = 0; i < V->Vertex->NbRows; ++i) {
  22         if (value_eq(V->Vertex->p[i][dim+1], lcm))
  23             continue;
  24         value_division(V->Vertex->p[i][dim+1], lcm, V->Vertex->p[i][dim+1]);
  25         Vector_Scale(V->Vertex->p[i], V->Vertex->p[i],
  26                      V->Vertex->p[i][dim+1], dim+1);
  27         value_assign(V->Vertex->p[i][dim+1], lcm);
  28     }
  29
  30     value_clear(lcm);
  31 }
  32
  33 /* Plug in the parametric vertex Vertex (nvar x (nparam + 2))
  34  * in the constraint constraint (1 + nvar + nparam + 1).
  35  * The result is stored in row (1 + nparam + 1),
  36  * with the denominator in position 0.
  37  */
  38 void Param_Inner_Product(Value *constraint, Matrix *Vertex, Value *row)
  39 {
  40     unsigned nparam = Vertex->NbColumns - 2;
  41     unsigned nvar = Vertex->NbRows;
  42     int j;
  43     Value tmp, tmp2;
  44
  45     value_set_si(row[0], 1);
  46     Vector_Set(row+1, 0, nparam+1);
  47
  48     value_init(tmp);
  49     value_init(tmp2);
  50
  51     for (j = 0 ; j < nvar; ++j) {
  52         value_set_si(tmp, 1);
  53         value_assign(tmp2,  constraint[1+j]);
  54         if (value_ne(row[0], Vertex->p[j][nparam+1])) {
  55             value_assign(tmp, row[0]);
  56             value_lcm(row[0], Vertex->p[j][nparam+1], &row[0]);
  57             value_division(tmp, row[0], tmp);
  58             value_multiply(tmp2, tmp2, row[0]);
  59             value_division(tmp2, tmp2, Vertex->p[j][nparam+1]);
  60         }
  61         Vector_Combine(row+1, Vertex->p[j], row+1, tmp, tmp2, nparam+1);
  62     }
  63     value_set_si(tmp, 1);
  64     Vector_Combine(row+1, constraint+1+nvar, row+1, tmp, row[0], nparam+1);
  65
  66     value_clear(tmp);
  67     value_clear(tmp2);
  68 }
  69
  70 static Param_Polyhedron *PL_P2PP(Polyhedron *Din, Polyhedron *Cin,
  71                                               struct barvinok_options *options)
  72 {
  73     unsigned MaxRays = options->MaxRays;
  74     if (MaxRays & (POL_NO_DUAL | POL_INTEGER))
  75         MaxRays = 0;
  76     return Polyhedron2Param_Domain(Din, Cin, MaxRays);
  77 }
  78
  79
  80 #ifndef POINTS2TRIANGS_PATH
  81 Param_Polyhedron *TC_P2PP(Polyhedron *P, Polyhedron *C,
  82                           struct barvinok_options *options)
  83 {
  84     return NULL;
  85 }
  86 #endif
  87
  88 Param_Polyhedron *Polyhedron2Param_Polyhedron(Polyhedron *P, Polyhedron *C,
  89                                               struct barvinok_options *options)
  90 {
  91     switch(options->chambers) {
  92     case BV_CHAMBERS_POLYLIB:
  93         return PL_P2PP(P, C, options);
  94         break;
  95     case BV_CHAMBERS_TOPCOM:
  96         return TC_P2PP(P, C, options);
  97         break;
  98     default:
  99         assert(0);
 100     }
 101 }