#ifndef SGIntersect_HXX
#define SGIntersect_HXX
+#include <algorithm>
+
template<typename T>
inline bool
intersects(const SGSphere<T>& s1, const SGSphere<T>& s2)
{
if (sphere.empty())
return false;
- // Is more or less trivially included in the next tests
- // if (box.empty())
- // return false;
-
- if (sphere.getCenter().x() < box.getMin().x() - sphere.getRadius())
- return false;
- if (sphere.getCenter().y() < box.getMin().y() - sphere.getRadius())
- return false;
- if (sphere.getCenter().z() < box.getMin().z() - sphere.getRadius())
+ if (box.empty())
return false;
- if (box.getMax().x() + sphere.getRadius() < sphere.getCenter().x())
- return false;
- if (box.getMax().y() + sphere.getRadius() < sphere.getCenter().y())
- return false;
- if (box.getMax().z() + sphere.getRadius() < sphere.getCenter().z())
- return false;
-
- return true;
+ SGVec3<T1> closest = box.getClosestPoint(sphere.getCenter());
+ return distSqr(closest, SGVec3<T1>(sphere.getCenter())) <= sphere.getRadius2();
}
// make it symmetric
template<typename T1, typename T2>
// The negative numerator for the \alpha expression
T num = plane.getPositiveDist();
- num -= dot(plane.getNormal(), lineSegment.getOrigin());
+ num -= dot(plane.getNormal(), lineSegment.getStart());
// If the numerator is zero, we have the lines origin included in the plane
if (fabs(num) <= SGLimits<T>::min()) {
- dst = lineSegment.getOrigin();
+ dst = lineSegment.getStart();
return true;
}
if (1 < alpha)
return false;
- dst = lineSegment.getOrigin() + alpha*lineSegment.getDirection();
+ dst = lineSegment.getStart() + alpha*lineSegment.getDirection();
return true;
}
// make it symmetric
T u = b*e - c*d;
T v = b*d - a*e;
+/*
// Regions
// \2|
// \|
// 3 |0\ 1
//----------
// 4 | 5 \ 6
+*/
if (u + v <= det) {
if (u < 0) {
T u = b*e - c*d;
T v = b*d - a*e;
+/*
// Regions
// \2|
// \|
// 3 |0\ 1
//----------
// 4 | 5 \ 6
+*/
if (u + v <= det) {
if (u < 0) {