global_execution_context.h

#ifndef GLOBAL_EXECUTION_CONTEXT_H
#define GLOBAL_EXECUTION_CONTEXT_H

#include "processing.h"

#include <CGAL/AABB_tree.h>
#if CGAL_VERSION_NR >= 1060000000
#include <CGAL/AABB_traits_3.h>
#else
#include <CGAL/AABB_traits.h>
#endif
#include <CGAL/AABB_face_graph_triangle_primitive.h>

// State that is relevant accross the different radii for which the
// program is executed. Used to map back the semantics from the original 
// IFC input to the newly constructed polyhedra.
template <typename TreeKernel>
struct global_execution_context : public execution_context {
	typedef CGAL::Polyhedron_3<TreeKernel> TreeShapeType;
	typedef CGAL::AABB_face_graph_triangle_primitive<TreeShapeType, CGAL::Default, CGAL::Tag_false> Primitive;
#if CGAL_VERSION_NR >= 1060000000
	typedef CGAL::AABB_traits_3<TreeKernel, Primitive> AAbbTraits;
#else
	typedef CGAL::AABB_traits<TreeKernel, Primitive> AAbbTraits;
#endif
	typedef typename AAbbTraits::Bounding_box Bounding_box;
	typedef CGAL::AABB_tree<AAbbTraits> AAbbTree;
	typedef typename AAbbTree::Primitive_id Primitive_id;
	
	typedef std::vector<std::pair<
		item_info*,
		std::list<cgal_shape_t::Facet_handle>>> segmentation_return_type;

	typedef std::vector<std::pair<
		item_info*,
		std::list<std::vector<std::vector<size_t>>::const_iterator>>> segmentation_return_type_2;

	AAbbTree tree;
	
	std::list<item_info*> infos;
	
	// A reference is kept to the original shapes in a std::list.
	// Later an aabb tree is used map eroded triangle centroids
	// back to the original elements to preserve semantics.
	std::list<TreeShapeType> triangulated_shape_memory;

	// A pointer only to conserve memory
	std::map<typename TreeShapeType::Facet_handle, item_info*> facet_to_info;

#ifdef GEOBIM_DEBUG
	simple_obj_writer obj_;
#endif

	global_execution_context();

	void operator()(shape_callback_item* item);

	void finalize();

	segmentation_return_type segment(const cgal_shape_t& input);

	typedef CGAL::Simple_cartesian<double>::Point_3 P3;
	std::vector<P3>* point_lookup;
	segmentation_return_type_2 segment(const non_manifold_polyhedron<CGAL::Simple_cartesian<double>>& input);

	const std::list<item_info*>& all_item_infos() const { return infos; }
};

#endif