c++ - 在 C++ 代码中提供 qconvex 的输入

Question

对于我的代码，我需要计算一系列点的 Convexhull，由于某些原因，我需要使用qhull库。在这个库中，有一个方法qconvex可以完全满足我的需要。我可以在终端中运行这个命令并得到我想要的。例如让我们假设我有一个像这样的输入points.txt：

2   #dimension
5   #number of points
0 0
1 0
0.5 0.5
1 1
0 1

我可以在终端一中运行这些命令来获得结果：qconvex Fx < points.txt或者cat points.txt | qconvex -Fx输出是：

4 
0
1
3
4

现在我的问题是如何在我的 C++ 代码中迭代地调用这个for命令float **rs_tmp;： qconvex 每次这 10 个点。如何qconvex在我的代码中运行并将其作为输入进行管道rs_tmp传输？宁愿避免将其rs_tmp写入某个临时文件并从中读取，因为我需要我的代码非常快。

float **rs_tmp;
for (int i = 0; i < NUMBER; i++)
{
    for (int j = 0; j < NUMBER; j++)
    {
        rs_tmp = generate_points(label, dect[i], dect[j], fun);
        // HERE I NEED TO CALL QCONVEX SOME HOW
        // THE POINTS ARE STORED IN rs_tmp as 2-Dimensional floating points array
    }
     int size = fun.size();
     for(int i = 0; i < size; ++i)     
     {
         delete[] rs_tmp[i]; 
     }   
     delete[] rs_tmp;         
}

Answer 1

我已经有同样的问题很长一段时间了，只是设法解决了它。Qhull 提供了一个非常简洁的示例，内容丰富。如果您从 git 克隆，您可以在目录 qhull/src/user_eg3/user_eg3.cpp 中看到示例。我花了一点时间来理解他们在做什么，但是一旦你掌握了它的窍门，它实际上很容易。我已经删除了所有额外的选项，但您正在寻找的选项。

/*--------------------------------------------
-user_eg3-  main procedure of user_eg3 application
*/
int main(int argc, char **argv) {
    try{
        return user_eg3(argc, argv);
    }catch(QhullError &e){
        cerr << e.what() << std::endl;
        return e.errorCode();
    }
}//main

int user_eg3(int argc, char **argv)
{
  RboxPoints rbox;
  Qhull qhull;
  line = "";
  std::istringstream is("2 4 1 0 1 1 0 0 0 1"); // To be passed to rbox: produces a unit square where 2 is the dimension 4 is the count and the points follow. This will also accept any valid rbox flags. 
  std::stringstream output;

  rbox.appendPoints(is); // appendPoints accepts either a const char* or an istream object. See libqhullcpp/RboxPoints.h and libqhullcpp/PointCoordinates.h 

  cerr << "@@@@@@@@@@\n" << rbox << "@@@@@@@@@@\n";
  qhull.runQhull(rbox, ""); //you can set any flag you would set for qhull in the terminal inside the ""
  qhull.setOutputStream(&output); //this will take any output stream 
  qhull.outputQhull("m"); //this will take any qhull output option
  cerr << "My Output : " << output.str() << "%%\n";
  qhull.setOutputStream(&cout);
  qhull.outputQhull("n");
  return 0;
}//user_eg3

希望这可以帮助。