introduced interface for handling element quality

MeshLib/MeshQuality/ElementQualityInterface.h

+/**
+ * \file   ElementQualityInterface.h
+ * \author Karsten Rink
+ * \date   2015-03-24
+ * \brief  Definition of the ElementQualityInterface class.
+ *
+ * \copyright
+ * Copyright (c) 2012-2015, OpenGeoSys Community (http://www.opengeosys.org)
+ *            Distributed under a Modified BSD License.
+ *              See accompanying file LICENSE.txt or
+ *              http://www.opengeosys.org/project/license
+ *
+ */
+
+#ifndef ELEMENTQUALITYINTERFACE_H_
+#define ELEMENTQUALITYINTERFACE_H_
+
+#include <vector>
+
+#include "BaseLib/Histogram.h"
+
+#include "MeshLib/Mesh.h"
+#include "MeshLib/MeshQuality/ElementQualityMetric.h"
+#include "MeshLib/MeshQuality/EdgeRatioMetric.h"
+#include "MeshLib/MeshQuality/AreaMetric.h"
+#include "MeshLib/MeshQuality/VolumeMetric.h"
+#include "MeshLib/MeshQuality/AngleSkewMetric.h"
+#include "MeshLib/MeshQuality/RadiusEdgeRatioMetric.h"
+
+namespace MeshLib
+{
+
+/**
+ * Interface class for handling mesh element quality metrics
+ */
+class ElementQualityInterface
+{
+public:
+	ElementQualityInterface(MeshLib::Mesh const& mesh, MeshQualityType t)
+	: _type(t), _mesh(mesh), _quality_tester(nullptr)
+	{
+		init(_mesh, _type);
+	}
+
+	~ElementQualityInterface() 
+	{ 
+		delete _quality_tester; 
+	}
+
+	std::vector<double> const& getQualityVector() const
+	{ 
+		return _quality_tester->getElementQuality(); 
+	}
+
+	int writeHistogram(std::string const& file_name) const
+	{
+		if (_quality_tester == nullptr)
+			return 1;
+
+		double const n_bins (1 + 3.3 * log (static_cast<float>(_mesh.getNElements())));
+		BaseLib::Histogram<double> const histogram (_quality_tester->getHistogram(static_cast<size_t>(n_bins)));
+		histogram.write(file_name, _mesh.getName(), MeshQualityType2String(_type));
+		return 0;
+	}
+	
+private:
+	void init(MeshLib::Mesh const& mesh, MeshQualityType t)
+	{
+		if (t == MeshQualityType::EDGERATIO)
+			_quality_tester = new MeshLib::EdgeRatioMetric(mesh);
+		else if (t == MeshQualityType::AREA)
+			_quality_tester = new MeshLib::AreaMetric(mesh);
+		else if (t == MeshQualityType::VOLUME)
+			_quality_tester = new MeshLib::VolumeMetric(mesh);
+		else if (t == MeshQualityType::EQUIANGLESKEW)
+			_quality_tester = new MeshLib::AngleSkewMetric(mesh);
+		else if (t == MeshQualityType::RADIUSEDGERATIO)
+			_quality_tester = new MeshLib::RadiusEdgeRatioMetric(mesh);
+		else
+		{
+			ERR("VtkVisPipeline::checkMeshQuality(): Unknown MeshQualityType.");
+			return;
+		}
+		_quality_tester->calculateQuality();
+	}
+
+	MeshQualityType const _type;
+	MeshLib::Mesh const& _mesh;
+	MeshLib::ElementQualityMetric* _quality_tester;
+};
+
+}
+
+#endif /* ELEMENTQUALITYINTERFACE_H_ */