X-Git-Url: http://git.tdb.fi/?a=blobdiff_plain;f=source%2Fgeometry%2Fcompositeshape.h;h=5d0cd8c59532ed3b782c3cd60908c8ce96778b6c;hb=44bd1d1ab256d397be4e2169c4ca5efdd0569d31;hp=50b04e80c506e93f7dc995187140e871dec66e27;hpb=827b227ae53884109b6885ccd1424dda8c0b3725;p=libs%2Fmath.git

diff --git a/source/geometry/compositeshape.h b/source/geometry/compositeshape.h
index 50b04e8..5d0cd8c 100644
--- a/source/geometry/compositeshape.h
+++ b/source/geometry/compositeshape.h
@@ -1,6 +1,8 @@
 #ifndef MSP_GEOMETRY_COMPOSITESHAPE_H_
 #define MSP_GEOMETRY_COMPOSITESHAPE_H_
 
+#include <algorithm>
+#include <stdexcept>
 #include <vector>
 #include "shape.h"
 
@@ -18,19 +20,25 @@ protected:
 	typedef std::vector<Shape<T, D> *> ShapeArray;
 
 	ShapeArray shapes;
+	unsigned max_isect;
 
+	CompositeShape() { }
 	CompositeShape(const Shape<T, D> &, const Shape<T, D> &);
-	CompositeShape(const ShapeArray &);
+	template<typename Iter>
+	void init_from_iter_range(const Iter &, const Iter &);
+private:
+	void init();
+protected:
 	CompositeShape(const CompositeShape &);
 	CompositeShape &operator=(const CompositeShape &);
 public:
 	virtual ~CompositeShape();
 
-	virtual HyperBox<T, D> get_axis_aligned_bounding_box() const;
+	virtual BoundingBox<T, D> get_axis_aligned_bounding_box(unsigned = 0) const;
 	virtual bool contains(const LinAl::Vector<T, D> &) const;
-	virtual bool check_intersection(const Ray<T, D> &) const;
-	virtual unsigned get_max_ray_intersections() const;
+	virtual unsigned get_max_ray_intersections() const { return max_isect; }
 	virtual unsigned get_intersections(const Ray<T, D> &, SurfacePoint<T, D> *, unsigned) const;
+	virtual Coverage get_coverage(const BoundingBox<T, D> &) const;
 };
 
 template<typename T, unsigned D, typename O>
@@ -39,17 +47,49 @@ inline CompositeShape<T, D, O>::CompositeShape(const Shape<T, D> &s1, const Shap
 	shapes.reserve(2);
 	shapes.push_back(s1.clone());
 	shapes.push_back(s2.clone());
+	init();
 }
 
 template<typename T, unsigned D, typename O>
-inline CompositeShape<T, D, O>::CompositeShape(const ShapeArray &s)
+template<typename Iter>
+inline void CompositeShape<T, D, O>::init_from_iter_range(const Iter &begin, const Iter &end)
 {
-	if(s.empty())
-		throw std::invalid_argument("CompositeShape::CompositeShape");
+	if(begin==end)
+		throw std::invalid_argument("CompositeShape::init_from_iter_range");
 
-	shapes.reserve(s.size());
-	for(typename ShapeArray::const_iterator i=s.begin(); i!=s.end(); ++i)
+	for(Iter i=begin; i!=end; ++i)
 		shapes.push_back((*i)->clone());
+	init();
+}
+
+template<typename T, unsigned D, typename O>
+inline void CompositeShape<T, D, O>::init()
+{
+	max_isect = 0;
+	for(typename ShapeArray::const_iterator i=shapes.begin(); i!=shapes.end(); ++i)
+		max_isect += (*i)->get_max_ray_intersections();
+}
+
+template<typename T, unsigned D, typename O>
+inline CompositeShape<T, D, O>::CompositeShape(const CompositeShape<T, D, O> &other):
+	shapes(other.shapes),
+	max_isect(other.max_isect)
+{
+	for(typename ShapeArray::iterator i=shapes.begin(); i!=shapes.end(); ++i)
+		*i = (*i)->clone();
+}
+
+template<typename T, unsigned D, typename O>
+inline CompositeShape<T, D, O> &CompositeShape<T, D, O>::operator=(const CompositeShape<T, D, O> &other)
+{
+	for(typename ShapeArray::iterator i=shapes.begin(); i!=shapes.end(); ++i)
+		delete *i;
+
+	shapes = other.shapes;
+	for(typename ShapeArray::iterator i=shapes.begin(); i!=shapes.end(); ++i)
+		*i = (*i)->clone();
+
+	max_isect = other.max_isect;
 }
 
 template<typename T, unsigned D, typename O>
@@ -60,9 +100,12 @@ inline CompositeShape<T, D, O>::~CompositeShape()
 }
 
 template<typename T, unsigned D, typename O>
-inline HyperBox<T, D> CompositeShape<T, D, O>::get_axis_aligned_bounding_box() const
+inline BoundingBox<T, D> CompositeShape<T, D, O>::get_axis_aligned_bounding_box(unsigned detail) const
 {
-	HyperBox<T, D> aabb;
+	if(detail)
+		return this->bisect_axis_aligned_bounding_box(detail);
+
+	BoundingBox<T, D> aabb;
 	for(typename ShapeArray::const_iterator i=shapes.begin(); i!=shapes.end(); ++i)
 	{
 		if(i==shapes.begin())
@@ -76,63 +119,98 @@ inline HyperBox<T, D> CompositeShape<T, D, O>::get_axis_aligned_bounding_box() c
 template<typename T, unsigned D, typename O>
 inline bool CompositeShape<T, D, O>::contains(const LinAl::Vector<T, D> &point) const
 {
-	bool inside = Ops::init_inside();
+	bool inside = false;
 	for(typename ShapeArray::const_iterator i=shapes.begin(); i!=shapes.end(); ++i)
 	{
-		inside = Ops::combine_inside(inside, (*i)->contains(point));
-		if(Ops::is_inside_decided(inside))
-			break;
+		inside = (*i)->contains(point);
+		if(Ops::shortcircuit(inside))
+			return inside;
 	}
 	return inside;
 }
 
-template<typename T, unsigned D, typename O>
-inline bool CompositeShape<T, D, O>::check_intersection(const Ray<T, D> &ray) const
-{
-	return get_intersections(ray, 0, 1);
-}
-
-template<typename T, unsigned D, typename O>
-inline unsigned CompositeShape<T, D, O>::get_max_ray_intersections() const
-{
-	unsigned max_isect = 0;
-	for(typename ShapeArray::const_iterator i=shapes.begin(); i!=shapes.end(); ++i)
-		max_isect += (*i)->get_max_ray_intersections();
-	return max_isect;
-}
-
 template<typename T, unsigned D, typename O>
 inline unsigned CompositeShape<T, D, O>::get_intersections(const Ray<T, D> &ray, SurfacePoint<T, D> *points, unsigned size) const
 {
+	SurfacePoint<T, D> *buffer = points;
+	unsigned buf_size = size;
+	if(!points && !Ops::shortcircuit(true))
+	{
+		buffer = new SurfacePoint<T, D>[max_isect];
+		buf_size = max_isect;
+	}
+
+	int start_nesting = 0;
 	unsigned n = 0;
-	for(typename ShapeArray::const_iterator i=shapes.begin(); (n<size && i!=shapes.end()); ++i)
+	for(typename ShapeArray::const_iterator i=shapes.begin(); (n<buf_size && i!=shapes.end()); ++i)
 	{
 		unsigned base = n;
-		unsigned count = (*i)->get_intersections(ray, points+base, size-base);
-		for(unsigned j=0; (n<size && j<count); ++j)
+		unsigned count = (*i)->get_intersections(ray, buffer+base, buf_size-base);
+		bool start_inside = (*i)->contains(ray.get_start());
+
+		if(!count && !start_inside)
 		{
-			SurfacePoint<T, D> &pt = points[base+j];
+			if(Ops::shortcircuit(false))
+				return 0;
+			continue;
+		}
 
-			bool surface = Ops::init_surface();
-			for(typename ShapeArray::const_iterator k=shapes.begin(); k!=shapes.end(); ++k)
-				if(k!=i)
-					surface = Ops::combine_surface(surface, (*k)->contains(pt.position));
+		start_nesting += start_inside;
 
-			if(surface)
-			{
-				if(points && base+j!=n)
-					points[n] = pt;
+		if(!base)
+		{
+			n = count;
+			continue;
+		}
 
-				++n;
+		sort_points(buffer, base+count);
+
+		int nesting = start_nesting;
+		unsigned k = 0;
+		for(unsigned j=0; j<base+count; ++j)
+		{
+			if(Ops::shortcircuit(nesting+buffer[j].entry<2))
+			{
+				if(j!=k)
+					buffer[k] = buffer[j];
+				++k;
 			}
+
+			nesting += buffer[j].entry*2-1;
 		}
+
+		if(!k && Ops::shortcircuit(false))
+			return 0;
+
+		n = k;
+
+		if(i!=shapes.begin())
+			start_nesting = (start_nesting>!Ops::shortcircuit(true));
 	}
 
-	sort_points(points, n);
+	if(buffer!=points)
+		delete[] buffer;
 
 	return n;
 }
 
+template<typename T, unsigned D, typename O>
+inline Coverage CompositeShape<T, D, O>::get_coverage(const BoundingBox<T, D> &bbox) const
+{
+	Coverage coverage = NO_COVERAGE;
+	for(typename ShapeArray::const_iterator i=shapes.begin(); i!=shapes.end(); ++i)
+	{
+		Coverage c = (*i)->get_coverage(bbox);
+		if(i==shapes.begin() || Ops::shortcircuit(c>coverage))
+			coverage = c;
+
+		if(coverage!=PARTIAL_COVERAGE && Ops::shortcircuit(coverage==FULL_COVERAGE))
+			break;
+	}
+
+	return coverage;
+}
+
 } // namespace Geometry
 } // namespace Msp