param_util.c

   1 #include "param_util.h"
   2
   3 void Param_Vertex_Common_Denominator(Param_Vertices *V)
   4 {
   5     unsigned dim;
   6     Value lcm;
   7     int i;
   8
   9     assert(V->Vertex->NbRows > 0);
  10     dim = V->Vertex->NbColumns-2;
  11
  12     value_init(lcm);
  13
  14     value_assign(lcm, V->Vertex->p[0][dim+1]);
  15     for (i = 1; i < V->Vertex->NbRows; ++i)
  16         value_lcm(V->Vertex->p[i][dim+1], lcm, &lcm);
  17
  18     for (i = 0; i < V->Vertex->NbRows; ++i) {
  19         if (value_eq(V->Vertex->p[i][dim+1], lcm))
  20             continue;
  21         value_division(V->Vertex->p[i][dim+1], lcm, V->Vertex->p[i][dim+1]);
  22         Vector_Scale(V->Vertex->p[i], V->Vertex->p[i],
  23                      V->Vertex->p[i][dim+1], dim+1);
  24         value_assign(V->Vertex->p[i][dim+1], lcm);
  25     }
  26
  27     value_clear(lcm);
  28 }