1

我已经尝试了一段时间以使以下代码正常运行。问题出在 sort(comp) 的比较函数中。它比较 2 个双精度值,但在某些情况下,它会导致程序崩溃并显示消息“调试断言失败!程序:retrace.exe 文件:c:程序文件 (x86)\microsoft visual studio 10.0\vc\include\algorithm Expression:Invalid运算符< "。我 100% 来自 stl 的排序,因为当程序崩溃时,成功执行的最后一行是排序前的那一行。起初我认为这是 double 的一些精度问题,但现在我对此表示怀疑。任何帮助或信息将不胜感激。

Vector startg;
bool comp(const std::pair<Vector, int>& p1, const std::pair<Vector, int>& p2)
{
    return (p1.first-startg).lengthSqr() < (p2.first-startg).lengthSqr();
}

bool Blobs::intersect(const Ray& ray, IntersectionInfo& info) const
{
    startg.set(ray.start.x, ray.start.y, ray.start.z);
    std::vector<std::pair<Vector, int> > spheres_inters;
    //not important stuff
    for(int i=0;i<n;i++)
    {
        Vector H = ray.start - centres[i];
        double A = ray.dir.lengthSqr();
        double B = 2 * dot(H, ray.dir);
        double C = H.lengthSqr() - bounding_radius*bounding_radius;
        double D = B*B - 4*A*C;
        if(D<=0)
            continue;
        double x1 = (-B + sqrt(D)) / (2*A);
        double x2 = (-B - sqrt(D)) / (2*A);
        Vector v1 = ray.start + x1*ray.dir;
        Vector v2 = ray.start + x2*ray.dir;
        spheres_inters.push_back(std::make_pair(v1, i));
        spheres_inters.push_back(std::make_pair(v2, i));
    }
    if(spheres_inters.size()==0)
        return false;
    ///////////////////////////////////////////////////////////////

    std::sort(spheres_inters.begin(), spheres_inters.end(), comp);//THE PROBLEM IS HERE!

    /* the rest of the method.....*/
4

0 回答 0