X-Git-Url: https://git.dogcows.com/gitweb?p=chaz%2Fyoink;a=blobdiff_plain;f=src%2FMoof%2FPlane.hh;h=11cd63ec23f2bf7f6517a8b2d5e937eb6f3832f3;hp=06e4abdce63f59398f961a531942d7b63ec667ea;hb=7e84479de612a4ce287c6f63deb014b447a993ec;hpb=25aefe01ef7dbdb603c51411e04b0d6a6107684f

diff --git a/src/Moof/Plane.hh b/src/Moof/Plane.hh
index 06e4abd..11cd63e 100644
--- a/src/Moof/Plane.hh
+++ b/src/Moof/Plane.hh
@@ -28,15 +28,16 @@
 
 #ifndef _MOOF_PLANE_HH_
 #define _MOOF_PLANE_HH_
-	
+
 #include <Moof/Math.hh>
+#include <Moof/Shape.hh>
 
 
 namespace Mf {
 
 
-class Aabb;
-class Sphere;
+template <int D> class Aabb;
+template <int D> class Sphere;
 
 
 /*
@@ -44,7 +45,7 @@ class Sphere;
  * is normal to the plane.
  */
 
-struct Plane
+struct Plane : public Shape<3>
 {
 	Vector3	normal;
 	Scalar	d;
@@ -65,8 +66,45 @@ struct Plane
 		d(scalar) {}
 
 
+	Scalar intersectRay(const Ray<3>& ray, Ray<3>::Intersection& intersection)
+	{
+		// solve: [(ray.point + t*ray.direction) dot normal] + d = 0
+
+		Scalar denominator = cml::dot(ray.direction, normal);
+
+		// check for parallel condition
+		if (denominator == SCALAR(0.0))
+		{
+			if (isEqual(cml::dot(ray.point, normal), -d))
+			{
+				// the ray lies on the plane
+				intersection.point = ray.point;
+				intersection.normal.set(0.0, 0.0, 0.0);
+				//intersection.normal = normal;
+				return SCALAR(0.0);
+			}
+
+			// no solution
+			return SCALAR(-1.0);
+		}
+
+		Scalar distance = cml::dot(ray.point, normal) + d;
+		Scalar t = -distance / denominator;
+		if (t > SCALAR(0.0))
+		{
+			ray.solve(intersection.point, t);
+
+			if (distance >= 0.0) intersection.normal = normal;
+			else                 intersection.normal = -normal;
+		}
+
+		return t;
+	}
+
+
 	/* Causes the normal of the plane to become normalized.  The scalar may also
-	 * be changed to keep the equation true. */
+	 * be changed to keep the equation true.  Word to the wise: don't normalize
+	 * a plane if the normal is the zero vector. */
 	void normalize()
 	{
 		Scalar mag = normal.length();
@@ -92,8 +130,8 @@ struct Plane
 		else                        return POSITIVE;
 	}
 
-	Halfspace intersects(const Aabb& aabb) const;
-	Halfspace intersects(const Sphere& sphere) const;
+	Halfspace intersects(const Aabb<3>& aabb) const;
+	Halfspace intersects(const Sphere<3>& sphere) const;
 };