X-Git-Url: https://git.dogcows.com/gitweb?p=chaz%2Fyoink;a=blobdiff_plain;f=src%2FMoof%2FOctree.hh;h=88dae7ac6492e587d1bb133516a27b4343d98735;hp=1c57ff2d3317fe83d10aa6972093683cb634556f;hb=64bd443538f57ad1bdff6c6b35953e72141129b2;hpb=329a48e4c4c2f5f2904b913938fc53154c48b825

diff --git a/src/Moof/Octree.hh b/src/Moof/Octree.hh
index 1c57ff2..88dae7a 100644
--- a/src/Moof/Octree.hh
+++ b/src/Moof/Octree.hh
@@ -29,6 +29,7 @@
 #ifndef _MOOF_OCTREE_HH_
 #define _MOOF_OCTREE_HH_
 
+#include <algorithm>
 #include <list>
 
 #include <boost/shared_ptr.hpp>
@@ -38,6 +39,8 @@
 #include <Moof/Aabb.hh>
 #include <Moof/Drawable.hh>
 #include <Moof/Entity.hh>
+#include <Moof/Frustum.hh>
+#include <Moof/Log.hh>
 #include <Moof/Math.hh>
 #include <Moof/Sphere.hh>
 	
@@ -45,309 +48,336 @@
 namespace Mf {
 
 
-struct OctreeNode : public Entity
+struct OctreeInsertable
 {
-	std::list<EntityPtr> objects;
+	virtual ~OctreeInsertable() {}
 
-	OctreeNode()
-	{
-		aabb_.min = Vector3(-1.0, -1.0, -1.0);
-		aabb_.max = Vector3(1.0, 1.0, 1.0);
-		sphere_.init(Vector3(0.0, 0.0, 0.0), 1.41421);
-	}
+	virtual bool isInsideAabb(const Aabb& aabb) const = 0;
+	virtual int getOctant(const Aabb& aabb) const = 0;
+};
 
-	OctreeNode(const Aabb& aabb)
-	{
-		aabb_ = aabb;
-		sphere_.point = aabb.getCenter();
-		sphere_.radius = (aabb.min - sphere_.point).length();
-	}
 
-	void draw(Scalar alpha) const
+template <class T>
+class Octree : public Entity
+{
+	typedef boost::shared_ptr<T> InsertableP;
+
+	struct Node : public Entity
 	{
-		std::list<EntityPtr>::const_iterator it;
+		std::list<InsertableP> objects;
 
-		for (it = objects.begin(); it != objects.end(); ++it)
+		Aabb aabb;
+		Sphere sphere;
+
+		Node(const Aabb& box) :
+			aabb(box)
 		{
-			(*it)->draw(alpha);
+			sphere.point = aabb.getCenter();
+			sphere.radius = (aabb.min - sphere.point).length();
 		}
-		if (!objects.empty())
-			aabb_.draw(); // temporary
-	}
 
-	void drawIfVisible(Scalar alpha, const Camera& cam) const
-	{
-		std::list<EntityPtr>::const_iterator it;
+		void draw(Scalar alpha) const
+		{
+			aabb.draw(alpha);
+		}
 
-		for (it = objects.begin(); it != objects.end(); ++it)
+		void drawIfVisible(Scalar alpha, const Frustum& frustum) const
 		{
-			(*it)->drawIfVisible(alpha, cam);
+			if (isVisible(frustum))
+			{
+				aabb.draw(alpha);
+			}
 		}
-		if (!objects.empty())
-			aabb_.draw();
-	}
 
+		void printSize()
+		{
+			logDebug("size of node %d", objects.size());
+		}
 
-	bool isVisible(const Camera& cam) const
-	{
-		if (sphere_.isVisible(cam))
+		void getAll(std::list<InsertableP>& insertables) const
 		{
-			return aabb_.isVisible(cam);
+			insertables.insert(insertables.end(), objects.begin(),
+					objects.end());
 		}
 
-		return false;
-	}
-};
+		void getIfVisible(std::list<InsertableP>& insertables,
+				const Frustum& frustum) const
+		{
+			typename std::list<InsertableP>::const_iterator it;
 
+			for (it = objects.begin(); it != objects.end(); ++it)
+			{
+				if ((*it)->isVisible(frustum)) insertables.push_back(*it);
+			}
+		}
 
-class Octree;
-typedef boost::shared_ptr<Octree> OctreePtr;
 
-class Octree
-{
+		bool isVisible(const Frustum& frustum) const
+		{
+			if (sphere.isVisible(frustum))
+			{
+				return aabb.isVisible(frustum);
+			}
+
+			return false;
+		}
+	};
 
-	stlplus::ntree<OctreeNode> root_;
 
 public:
 
-	explicit Octree(const OctreeNode& rootNode)
-	{
-		root_.insert(rootNode);
-	}
+	typedef boost::shared_ptr<Octree> Ptr;
+	typedef typename stlplus::ntree<Node>::iterator NodeP;
+	
+private:
 
-	void insert(EntityPtr entity)
-	{
-		insert(root_.root(), entity);
-	}
 
-	void insert(stlplus::ntree<OctreeNode>::iterator node, EntityPtr entity)
+	NodeP insert(InsertableP entity, NodeP node)
 	{
-		Plane::Halfspace halfspace;
-		int octantNum = -1;
+		ASSERT(node.valid() && "invalid node passed");
+		ASSERT(entity && "null entity passed");
 
-		if (!node.valid())
+		if (entity->isInsideAabb(node->aabb))
 		{
-			std::cerr << "cannot insert into invalid node" << std::endl;
-			return;
+			return insert_recurse(entity, node);
 		}
-
-		Plane xy = node->getAabb().getPlaneXY();
-		halfspace = xy.intersectsSphere(entity->getSphere());
-
-		if (halfspace == Plane::POSITIVE)
+		else
 		{
-			Plane xz = node->getAabb().getPlaneXZ();
-			halfspace = xz.intersectsSphere(entity->getSphere());
-
-			if (halfspace == Plane::POSITIVE)
-			{
-				Plane yz = node->getAabb().getPlaneYZ();
-				halfspace = yz.intersectsSphere(entity->getSphere());
-
-				if (halfspace == Plane::POSITIVE)
-				{
-					octantNum = 2;
-				}
-				else if (halfspace == Plane::NEGATIVE)
-				{
-					octantNum = 3;
-				}
-			}
-			else if (halfspace == Plane::NEGATIVE)
-			{
-				Plane yz = node->getAabb().getPlaneYZ();
-				halfspace = yz.intersectsSphere(entity->getSphere());
-
-				if (halfspace == Plane::POSITIVE)
-				{
-					octantNum = 1;
-				}
-				else if (halfspace == Plane::NEGATIVE)
-				{
-					octantNum = 0;
-				}
-			}
+			node->objects.push_back(entity);
+			return node;
 		}
-		else if (halfspace == Plane::NEGATIVE)
-		{
-			Plane xz = node->getAabb().getPlaneXZ();
-			halfspace = xz.intersectsSphere(entity->getSphere());
-
-			if (halfspace == Plane::POSITIVE)
-			{
-				Plane yz = node->getAabb().getPlaneYZ();
-				halfspace = yz.intersectsSphere(entity->getSphere());
-
-				if (halfspace == Plane::POSITIVE)
-				{
-					octantNum = 6;
-				}
-				else if (halfspace == Plane::NEGATIVE)
-				{
-					octantNum = 7;
-				}
-			}
-			else if (halfspace == Plane::NEGATIVE)
-			{
-				Plane yz = node->getAabb().getPlaneYZ();
-				halfspace = yz.intersectsSphere(entity->getSphere());
+	}
 
-				if (halfspace == Plane::POSITIVE)
-				{
-					octantNum = 5;
-				}
-				else if (halfspace == Plane::NEGATIVE)
-				{
-					octantNum = 4;
-				}
-			}
-		}
+	NodeP insert_recurse(InsertableP entity, NodeP node)
+	{
+		ASSERT(node.valid() && "invalid node passed");
+		ASSERT(entity && "null entity passed");
 
+		int octantNum = entity->getOctant(node->aabb);
 		if (octantNum == -1)
 		{
-			node->objects.push_front(entity);
-			//return node;
+			node->objects.push_back(entity);
+			return node;
 		}
 		else
 		{
-			if (root_.children(node) == 0)
+			if ((int)tree_.children(node) <= octantNum)
 			{
-				addChildren(node);
+				addChild(node, octantNum);
 			}
 
-			stlplus::ntree<OctreeNode>::iterator child = root_.child(node, octantNum);
+			NodeP child = tree_.child(node, octantNum);
+			ASSERT(child.valid() && "expected valid child node");
 
-			if (child.valid())
-			{
-				return insert(child, entity);
-			}
-			else
-			{
-				std::cerr << "expected but found no child at index " << octantNum << std::endl;
-				//return stlplus::ntree<OctreeNode>::iterator();
-			}
-			//return WeakPtr();
+			return insert_recurse(entity, child);
 		}
 	}
 
-	void addChildren(stlplus::ntree<OctreeNode>::iterator node)
+	void addChild(NodeP node, int index)
 	{
-		Aabb octant;
+		ASSERT(node.valid() && "invalid node passed");
 
-		if (!node.valid())
-		{
-			std::cerr << "cannot add children to invalid node" << std::endl;
-			return;
-		}
+		Aabb octant;
 
-		for (int i = 0; i < 8; ++i)
+		for (int i = tree_.children(node); i <= index; ++i)
 		{
-			node->getAabb().getOctant(octant, i);
-			//OctreeNode octantNode(octant);
-
-			root_.append(node, octant);
+			node->aabb.getOctant(octant, i);
+			tree_.append(node, octant);
 		}
 	}
 
-	void draw(stlplus::ntree<OctreeNode>::iterator node, Scalar alpha)
+
+	void getNearbyObjects(std::list<InsertableP>& insertables,
+			const OctreeInsertable& entity, NodeP node) const
 	{
-		if (!node.valid())
-		{
-			std::cerr << "cannot draw null child node :-(" << std::endl;
-			return;
-		}
+		ASSERT(node.valid() && "invalid node passed");
 
-		node->draw(alpha);
+		node->printSize();
 
-		for (unsigned i = 0; i < root_.children(node); ++i)
+		int octantNum = entity.getOctant(node->aabb);
+		if (octantNum != -1)
 		{
-			stlplus::ntree<OctreeNode>::iterator child = root_.child(node, i);
+			node->getAll(insertables);
 
-			if (child.valid())
+			if (octantNum < (int)tree_.children(node))
 			{
-				draw(child, alpha);
-			}
-			else
-			{
-				std::cerr << "node is not a leaf, but has an invalid child" << std::endl;
-			}
+				NodeP child = tree_.child(node, octantNum);
+				ASSERT(child.valid() && "expected valid child node");
 
+				getNearbyObjects(insertables, entity, child);
+			}
+		}
+		else
+		{
+			logDebug("getting all the rest...");
+			getAll(insertables, node);
 		}
 	}
 
-	void drawIfVisible(stlplus::ntree<OctreeNode>::iterator node,
-			Scalar alpha, const Camera& cam)
+
+	void getAll(std::list<InsertableP>& insertables, NodeP node) const
 	{
-		//node.drawIfVisible(alpha, cam);
-		
-		if (!node.valid())
+		ASSERT(node.valid() && "invalid node passed");
+
+		node->getAll(insertables);
+
+		for (unsigned i = 0; i < tree_.children(node); ++i)
 		{
-			std::cerr << "invalid child while drawing :-(" << std::endl;
-			return;
+			NodeP child = tree_.child(node, i);
+			ASSERT(child.valid() && "expected valid child node");
+
+			getAll(insertables, child);
 		}
+	}
+
+	void getIfVisible(std::list<InsertableP>& insertables,
+			const Frustum& frustum, NodeP node) const
+	{
+		ASSERT(node.valid() && "invalid node passed");
 
-		Frustum::Collision collision =
-			cam.getFrustum().containsSphere(node->getSphere());
+		// try to cull by sphere
+		Frustum::Collision collision = frustum.contains(node->sphere);
 		if (collision == Frustum::OUTSIDE) return;
 
-		collision = cam.getFrustum().containsAabb(node->getAabb());
+		// try to cull by aabb
+		collision = frustum.contains(node->aabb);
 		if (collision == Frustum::OUTSIDE) return;
 
 
 		if (collision == Frustum::INSIDE)
 		{
-			node->draw(alpha);
+			node->getAll(insertables);
 		}
 		else // collision == Frustum::INTERSECT
 		{
-			node->drawIfVisible(alpha, cam);
+			node->getIfVisible(insertables, frustum);
 		}
 
-		if (root_.children(node) > 0)
+		if (tree_.children(node) > 0)
 		{
 			if (collision == Frustum::INSIDE)
 			{
-				for (unsigned i = 0; i < root_.children(node); ++i)
+				for (unsigned i = 0; i < tree_.children(node); ++i)
 				{
-					stlplus::ntree<OctreeNode>::iterator child = root_.child(node, i);
-
-					if (child.valid())
-					{
-						draw(child, alpha);
-					}
-					else
-					{
-						std::cerr << "node is not a leaf, but has an invalid child" << std::endl;
-					}
+					NodeP child = tree_.child(node, i);
+					ASSERT(child.valid() && "expected valid child node");
 
+					getAll(insertables, child);
 				}
 			}
 			else // collision == Frustum::INTERSECT
 			{
-				for (unsigned i = 0; i < root_.children(node); ++i)
+				for (unsigned i = 0; i < tree_.children(node); ++i)
 				{
-					stlplus::ntree<OctreeNode>::iterator child = root_.child(node, i);
-
-					if (child.valid())
-					{
-						drawIfVisible(child, alpha, cam);
-					}
-					else
-					{
-						std::cerr << "node is not a leaf, but has an invalid child" << std::endl;
-					}
+					NodeP child = tree_.child(node, i);
+					ASSERT(child.valid() && "expected valid child node");
+
+					getIfVisible(insertables, frustum, child);
 				}
 			}
 		}
 	}
 
-	void drawIfVisible(Scalar alpha, const Camera& cam)
+
+	mutable stlplus::ntree<Node> tree_;
+
+
+public:
+
+	void print(NodeP node)
 	{
-		drawIfVisible(root_.root(), alpha, cam);
+		logInfo("-----");
+		logInfo("depth to node: %d", tree_.depth(node));
+		logInfo("size of node: %d", tree_.size(node));
 	}
 
-};
+	static Ptr alloc(const Node& rootNode)
+	{
+		return Ptr(new Octree(rootNode));
+	}
+
+	explicit Octree(const Node& rootNode)
+	{
+		tree_.insert(rootNode);
+	}
+
+
+	NodeP insert(InsertableP entity)
+	{
+		return insert(entity, tree_.root());
+	}
+
+	void remove(InsertableP entity, NodeP node)
+	{
+		ASSERT(entity && "null entity passed");
+		ASSERT(node.valid() && "invalid node passed");
 
+		typename std::list<InsertableP>::iterator it;
+		it = std::find(node->objects.begin(), node->objects.end(), entity);
+
+		if (it != node->objects.end())
+		{
+			node->objects.erase(it);
+		}
+	}
+
+
+	NodeP reinsert(InsertableP entity, NodeP node)
+	{
+		remove(entity, node);
+		return insert(entity);
+	}
+
+
+	void draw(Scalar alpha) const
+	{
+		std::list<InsertableP> objects;
+		getAll(objects);
+
+		typename std::list<InsertableP>::const_iterator it;
+		for (it = objects.begin(); it != objects.end(); ++it)
+		{
+			(*it)->draw(alpha);
+		}
+	}
+
+	void drawIfVisible(Scalar alpha, const Frustum& frustum) const
+	{
+		std::list<InsertableP> objects;
+		//getIfVisible(objects, frustum);
+		getNearbyObjects(objects, *savedObj);
+
+		typename std::list<InsertableP>::const_iterator it;
+		for (it = objects.begin(); it != objects.end(); ++it)
+		{
+			(*it)->draw(alpha);
+		}
+	}
+
+
+	void getAll(std::list<InsertableP>& insertables) const
+	{
+		getAll(insertables, tree_.root());
+	}
+
+	void getIfVisible(std::list<InsertableP>& insertables,
+			const Frustum& frustum) const
+	{
+		getIfVisible(insertables, frustum, tree_.root());
+	}
+
+	mutable const OctreeInsertable* savedObj;
+
+
+	void getNearbyObjects(std::list<InsertableP>& insertables,
+			const OctreeInsertable& entity) const
+	{
+		logDebug("--- GETTING NEARBY");
+		getNearbyObjects(insertables, entity, tree_.root());
+		logDebug("---");
+		savedObj = &entity;
+	}
+};
 
 
 } // namespace Mf