Skip to content
Snippets Groups Projects
Node.cpp 4.31 KiB
Newer Older
#include <limits>
#include <cmath>
#include <iostream>
Thomas Steinreiter's avatar
Thomas Steinreiter committed
#include <algorithm>
#include "Node.hpp"

namespace nbody {
	Node::Node(Tree* _tree):tree(_tree) {}
	Box Node::getBB() const {
	void Node::setBB(const Box& bb_) {
		bb = bb_;
	//check if node needs to be splitted during tree build
	bool Node::isSplitable() const {
		if (bodies.size() <= tree->maxLeafBodies) {
		//this is to prevent errors with collocated particles
		if (bb.volume() <= std::numeric_limits<float>::epsilon()) {
			result = false;
		}
		return result;
	}

	void Node::extendBBforBodies() {
		bb.extendForBodies(bodies);
	}

	void Node::extendBBtoCube() {
	std::vector<Body> Node::getBodies() const {
	}

	void Node::insertBefore(Node* node) {
		node->next = this;
		node->prev = this->prev;
		this->prev->next = node;
		this->prev = node;
	}

	void Node::insertAfter(Node* node) {
		this->next->insertBefore(node);
	}

	void Node::unlink() {
		this->next->prev = this->prev;
		this->prev->next = this->next;
		this->next = this;
		this->prev = this;
	}

	bool Node::isCorrect() const {
		if (afterSubtree == nullptr) {
Thomas Steinreiter's avatar
Thomas Steinreiter committed
			std::cerr << "after subtree null\n";
		if (!bb.isCorrectBox()) {
Thomas Steinreiter's avatar
Thomas Steinreiter committed
			std::cerr << "bb wrong\n";
		for (std::size_t i = 0; i < 3; i++) {
			if (bb.min[i] > bb.max[i]) {
				std::cerr << "bb " << i << " min " << bb.min[i] << " max " << bb.max[i] << '\n';
		if (std::any_of(std::begin(bodies), std::end(bodies), [&](const Body& b) {return !isContained(b, bb); })) {
Thomas Steinreiter's avatar
Thomas Steinreiter committed
			std::cerr << "bb out of bounds\n";
			return false;
		if (!leaf) {
			Node* current = next;
			std::size_t children = 0;
			while (current != nullptr && current != afterSubtree) {
				current = current->afterSubtree;
				children++;
			}
			if (current == nullptr) {
Thomas Steinreiter's avatar
Thomas Steinreiter committed
				std::cerr << "afterSubtree null\n";
			if (children != tree->numberOfChildren()) {
Thomas Steinreiter's avatar
Thomas Steinreiter committed
				std::cerr << "wrong number of children " << children << '\n';
			current = next;
			for (std::size_t i = 0; i < tree->numberOfChildren(); i++) {
				current = current->afterSubtree;
			}
			if (current != afterSubtree) {
Thomas Steinreiter's avatar
Thomas Steinreiter committed
				std::cerr << "last sibling afterSubtree inconsistent\n";
		if (!leaf && !bodies.empty()) {
Thomas Steinreiter's avatar
Thomas Steinreiter committed
			std::cerr << "non-empty inner node\n";
		if (leaf && nextSibling != nullptr && next != nextSibling) {
Thomas Steinreiter's avatar
Thomas Steinreiter committed
			std::cerr << "wrong next sibling\n";
	//update overall node information
	void Node::update() {
		double position[3] = {0.0, 0.0, 0.0};
		double mass = 0.0;

			for (const auto& body : bodies) {
				mass += body.mass;
				for (std::size_t i = 0; i < 3; i++) {
					position[i] += body.position[i] * body.mass;
			for (Node* node = next; node != nullptr; node = node->nextSibling) {
				mass += node->representative.mass;
			}
		}
		for (std::size_t i = 0; i < 3; i++) {
			representative.position[i] = position[i] / mass;
		representative.mass = mass;
	//get criterion to check if node is sufficient for force evaluation
	double Node::getL() const {
		return bb.maxSidelength();
	void Node::print(std::size_t parallelId) const {
		for (const auto& body : bodies) {
			std::cout << "  ";
			body.print(parallelId);
	//check if node is sufficient for force evaluation
	bool Node::sufficientForBody(const Body& body) const {
		double distance = 0.0;

		for (std::size_t i = 0; i < 3; i++) {
			distance += (representative.position[i] - body.position[i]) * (representative.position[i] - body.position[i]);
		return sqrt(distance) > getL(); // TODO(steinret): do squared comparison
	//check if node is sufficient for force evaluation for all bodies in box
	bool Node::sufficientForBox(const Box& box) const {
		return bb.distanceToBox(box) > getL();
	void Node::setBodies(const std::vector<Body>& bodies_) {
		bodies = bodies_;
	void Node::setBodies(std::vector<Body>&& bodies_) {
		bodies = std::move(bodies_);
	//get local bodies
Thomas Steinreiter's avatar
Thomas Steinreiter committed
	void Node::extractLocalBodiesTo(std::vector<Body>& result) {
		std::copy_if(std::begin(bodies), std::end(bodies), std::back_inserter(result), [](const Body& b) {return !b.refinement; });
		bodies.clear();
} // namespace nbody