#ifndef KDTREE_H

#define KDTREE_H

#include <iostream>

#include <fstream>

#include "container.h"

#include "vector.h"

#include "scene.h"

using namespace std;

class KdNode

{

	Float split;

	short axis; /* 0,1,2 => x,y,z; 3 => leaf */

public:

	union {

		KdNode *children;

		ShapeList *shapes;

	};

	KdNode() : axis(3) { shapes = new ShapeList(); };

	~KdNode();

	void setAxis(short aAxis) { axis = aAxis; };

	short getAxis() { return axis; };

	void setSplit(Float aSplit) { split = aSplit; };

	Float getSplit() { return split; };

	void setLeaf() { axis = 3; };

	bool isLeaf() { return axis == 3; };

	KdNode *getLeftChild() { return children; };

	KdNode *getRightChild() { return children+1; };

	void addShape(Shape* aShape) { shapes->push_back(aShape); };

	void subdivide(BBox bbox, int depth);

};

class KdTree: public Container

{

	KdNode *root;

	bool built;

	int max_depth;

public:

	KdTree() : Container(), root(NULL), built(false), max_depth(32) {};

	~KdTree() { if (root) delete root; };

	void addShape(Shape* aShape) { Container::addShape(aShape); built = false; };

	Shape *nearest_intersection(const Shape *origin_shape, const Ray &ray,

		Float &nearest_distance);

	void optimize() { build(); };

	void build();

	void save(ostream &str, KdNode *node = NULL);

	void load(istream &str, KdNode *node = NULL);

	void setMaxDepth(int md) { max_depth = md; };

};

#endif