Parma Polyhedra库:顶点枚举和浮点运算

时间:2018-04-02 08:44:43

标签: c++ discrete-mathematics polyhedra

我正在使用Parma Polyhedra库(PPL)来执行给定多面体的顶点枚举。这已经讨论过here。但是,我无法弄清楚如何在计算中使用有理数而不是整数。

下面的代码生成一个线段[0.3,3.7](一维凸面多面体),PPL返回两个整数{0,3},但我想要有理数{0.3,3.7}。我如何建议PPL使用有理数(浮点运算)?

    #include <cstddef>
    #include <stdio.h>
    #include "ppl.hh"

    using namespace Parma_Polyhedra_Library;

    int main() {
      Variable x(0);
      C_Polyhedron ph(1);
      ph.refine_with_constraint( x <= 3.7);
      ph.refine_with_constraint( x >= 0.3);

      Generator_System gs = ph.generators();
      for(Generator_System::const_iterator it = gs.begin(); it != gs.end(); it++) {
        const Generator& g = *it;
        std::cout << g.coefficient(x) << std::endl;
      }
      return 0;
    }

1 个答案:

答案 0 :(得分:0)

ppl期望在其约束的定义中使用整数。你必须用公分母来衡量你的不平等。