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