* working on threading

bh_tree_mpi/datastructures/BarnesHutTree.cpp

@@ -58,7 +58,7 @@ namespace nbody {
 	}
 
 	void BarnesHutTree::splitNode(Node* current) {
-		if (current->isSplitable(this->control.maxLeafBodies)) {
+		if (current->isSplitable(&this->control)) {
 			vector<Box> subboxes = this->splitBB(current);
 			current->leaf = false;
 			Node* after = current->next;

bh_tree_mpi/datastructures/Node.cpp

@@ -34,14 +34,18 @@ namespace nbody {
 		this->bb = bb;
 	}
 
-	bool Node::isSplitable(unsigned int maxBodies) {
-		if (this->bodies.size() < maxBodies + 1) {
-			return false;
+	bool Node::isSplitable(Control* control) {
+		bool result = true;
+
+		pthread_rwlock_rdlock(&control->lock);
+		if (this->bodies.size() < control->maxLeafBodies + 1) {
+			result = false;
 		}
 		if (this->bb.volume() <= FLT_EPSILON) {
-			return false;
+			result = false;
 		}
-		return true;
+		pthread_rwlock_unlock(&control->lock);
+		return result;
 	}
 
 	void Node::extendBBforBodies() {

bh_tree_mpi/datastructures/Node.hpp

@@ -29,7 +29,7 @@ namespace nbody {
 	public:
 		Node(Tree* tree);
 		virtual ~Node();
-		virtual bool isSplitable(unsigned int maxBodies);
+		virtual bool isSplitable(Control* control);
 		virtual void extendBBforBodies();
 		virtual void extendBBtoCube();
 		virtual Box getBB();

bh_tree_mpi/parallelization/MpiSimulation.cpp

@@ -209,6 +209,9 @@ namespace nbody {
 		this->tree.computeForces();
 		this->tree.advance();
 		this->domains[this->mpiRank] = this->tree.getLocalBB();
+		if (!this->tree.isCorrect()) {
+			cout << "tree not correct" << endl;
+		}
 	}
 
 }