与 Boost::geometry 的多边形相交严重的性能下降答案

【问题标题】：Polygon intersection with Boost::geometry severe performance deterioration与 Boost::geometry 的多边形相交严重的性能下降
【发布时间】：2013-04-17 21:38:32
【问题描述】：

我有一个粒子系统，我正在使用boost::geometry 将我的椭圆粒子近似为多边形，然后使用库的交集函数找到重叠区域。我正在计算一个“内部”和“外部”椭圆（多边形）区域来为每个粒子-粒子相互作用分配一个“势”。

我的潜在功能是这样的：

    double Potential(Cell* current, Cell* next)
{
    double areaRep, areaAtt;
    double distance = Distance(current,next);
    double A1 = current->getLength();
    double B1 = A1/2.0;
    double theta1 = current->getTheta(); //*180.0/M_PI
    double x1 = current->getCurrX();
    double y1 = current->getCurrY();
    double A2 = next->getLength();
    double B2 = A2/2.0;
    double theta2 = next->getTheta();
    double x2 = next->getCurrX();
    double y2 = next->getCurrY();
    polygon_2d poly1, poly2, poly3, poly4;
    double lamda1, lamda2;
    lamda1 = 0.0005; lamda2 = 0.00001;
    if(distance < 2.0*1.5*A1) {
        ellipse2poly(theta1, A1, B1, x1, y1, &poly1);
        ellipse2poly(theta2, A2, B2, x2, y2, &poly2);
        areaRep = getOverlapingAreaPoly(poly1,poly2);
        ellipse2poly(theta1, 1.5*A1, 1.5*B1, x1, y1, &poly3);
        ellipse2poly(theta2, 1.5*A2, 1.5*B2, x2, y2, &poly4);
        areaAtt = getOverlapingAreaPoly(poly3, poly4);
        return (lamda1*areaRep - lamda2*areaAtt); 
    }
    else
        return 0.0;
}

“多边形化”功能是：

int ellipse2poly(double theta, double A1, double B1, double H1, double K1, polygon_2d *po)
{
    using namespace boost::geometry;
    polygon_2d  poly;
    const int n = 20;
    double angle = theta; // cell orientation
    double a = A1; // Long semi-axis length
    double b = B1; // short semi-axis length
    double xc = H1; // current X position
    double yc = K1; // current Y position

    if(!n)
    {
        std::cout << "error ellipse(): n should be >0\n" <<std::endl;
        return 0;
    }
    double t = 0;
    int i = 0;
    double coor[2*n+1][2];

    double x, y;
    double step = M_PI/(double)n;
    double sinphi = sin(angle);
    double cosphi = cos(angle);
    for(i=0; i<2*n+1; i++)
    {   
        x = xc + a*cos(t)*cosphi - b*sin(t)*sinphi;
        y = yc + a*cos(t)*sinphi + b*sin(t)*cosphi;
        coor[i][0] = x;
        coor[i][1] = y;
        t += step;
    }
    assign_points(poly, coor);
    correct(poly);
    *po = poly;
    return 1;
}

而返回区域为：

double getOverlapingAreaPoly(polygon_2d poly, polygon_2d poly2)
{
    point_2d cent; //centre of overlaping area
    double overAreaPoly = 0.0;
    typedef std::vector<polygon_2d > polygon_list;
    polygon_list v;

    intersection(poly,poly2,v);
    for (polygon_list::const_iterator it = v.begin(); it != v.end(); ++it)
    {
        centroid(*it, cent);
        overAreaPoly = area(*it);
    }

    return overAreaPoly;
}

只要不是针对同一个细胞（粒子），就会为每个细胞（粒子）调用该函数。以前，使用另一种方法，我的算法的一次迭代对于 100 个粒子的一次迭代大约需要 43 毫秒。现在大约需要 1 分钟（！！！），所以我想我做错了什么！

我只在 Win7 64bit 下的 MSVC2012 中对此进行了测试。我将使用 Qt 4.7.4 报告 Linux Mint。

编辑：我已经用 Qt 4.7.4 在 Linux Mint 上进行了测试，它运行得非常合理；每次迭代可能需要 90-100 毫秒，这很好。不知道win7出了什么问题...

【问题讨论】：

你的 Potential() 函数是固定的还是你愿意修改它？您能否举例说明您打算如何利用潜力？
其实我修好了。我刚刚在 Visual Studio 中启动了一个新项目，复制了源文件和头文件，重新编译，它再次非常合理地工作了！我不知道为什么实际上是这样......为了回答你的问题，我使用粒子之间距离的归一化向量来确定相互作用的方向，然后乘以势能。
感谢您的澄清。很高兴听到它现在工作正常:-)

标签： c++ performance boost boost-geometry

【解决方案1】：

我实际上已经修复了它。我在 Visual Studio 中启动了一个新项目并复制了所有源文件和头文件，重新编译，现在一切运行顺利。我想从根本上改变代码和添加/减去东西一定会产生一些影响......

【讨论】：