Merge pull request #588 from endJunction/DebuggingOfElements

[MeL] Simple Element::operator<<() for debugging.

MeshLib/Elements/Element.cpp

@@ -186,4 +186,22 @@ bool Element::isBoundaryElement() const
         [](MeshLib::Element const*const e){ return e == nullptr; });
 }
 
+#ifndef NDEBUG
+std::ostream& operator<<(std::ostream& os, Element const& e)
+{
+	os << "Element #" << e._id << " value " << e._value << " @ " << &e
+		<< " with " << e.getNNeighbors() << " neighbours\n";
+
+	unsigned const nnodes = e.getNNodes();
+	MeshLib::Node* const* const nodes = e.getNodes();
+	os << "MeshElemType: "
+		<< static_cast<std::underlying_type<MeshElemType>::type>(e.getGeomType())
+		<< " with " << nnodes << " nodes: { ";
+	for (unsigned n = 0; n < nnodes; ++n)
+		os << nodes[n]->getID() << " @ " << nodes[n] << "  ";
+	os << "}\n";
+	return os;
+}
+#endif  // NDEBUG
+
 }

MeshLib/Elements/Element.h

@@ -207,6 +207,9 @@ public:
 	 */
 	virtual bool testElementNodeOrder() const = 0;
 
+#ifndef NDEBUG
+	friend std::ostream& operator<<(std::ostream& os, Element const& e);
+#endif  // NDEBUG
 
 protected:
 	/// Constructor for a generic mesh element without an array of mesh nodes.