add reviseMesh tool

Applications/Utils/MeshEdit/CMakeLists.txt

@@ -103,3 +103,18 @@ ADD_CATALYST_DEPENDENCY(checkMesh)
 
 SET_TARGET_PROPERTIES(checkMesh
 	PROPERTIES FOLDER Utilities)
+
+ADD_EXECUTABLE( reviseMesh reviseMesh.cpp)
+TARGET_LINK_LIBRARIES ( reviseMesh
+	BaseLib
+	FileIO
+	MathLib
+	MeshLib
+	InSituLib
+	${CATALYST_LIBRARIES}
+)
+ADD_CATALYST_DEPENDENCY(reviseMesh)
+
+SET_TARGET_PROPERTIES(reviseMesh
+	PROPERTIES FOLDER Utilities)
+	
\ No newline at end of file

Applications/Utils/MeshEdit/reviseMesh.cpp

+/**
+ * @copyright
+ * Copyright (c) 2012-2014, OpenGeoSys Community (http://www.opengeosys.org)
+ *            Distributed under a Modified BSD License.
+ *              See accompanying file LICENSE.txt or
+ *              http://www.opengeosys.org/LICENSE.txt
+ *
+ */
+
+#include <array>
+#include <string>
+
+#include "logog/include/logog.hpp"
+#include "tclap/CmdLine.h"
+
+#include "BaseLib/BuildInfo.h"
+#include "BaseLib/StringTools.h"
+#include "BaseLib/LogogSimpleFormatter.h"
+#include "BaseLib/FileTools.h"
+
+#include "MeshLib/Node.h"
+#include "MeshLib/Elements/Element.h"
+#include "MeshLib/Mesh.h"
+#include "MeshLib/MeshEditing/MeshRevision.h"
+
+#include "FileIO/readMeshFromFile.h"
+#include "FileIO/writeMeshToFile.h"
+
+int main(int argc, char *argv[])
+{
+	LOGOG_INITIALIZE();
+	logog::Cout* logog_cout (new logog::Cout);
+	BaseLib::LogogSimpleFormatter *custom_format (new BaseLib::LogogSimpleFormatter);
+	logog_cout->SetFormatter(*custom_format);
+
+	try {
+		TCLAP::CmdLine cmd("Mesh revision tool", ' ', BaseLib::BuildInfo::git_version_sha1);
+		TCLAP::ValueArg<std::string> input_arg("i", "input-mesh-file","input mesh file",true,"","string");
+		cmd.add( input_arg );
+		TCLAP::ValueArg<std::string> output_arg("o", "output-mesh-file","output mesh file",true,"","string");
+		cmd.add( output_arg );
+		TCLAP::SwitchArg simplify_arg("s","simplify","simplify the mesh (removing duplicated nodes)");
+		cmd.add( simplify_arg );
+		TCLAP::ValueArg<double> eps_arg("e", "eps","Minimum distance for nodes not to be collapsed",
+										false, std::numeric_limits<double>::epsilon(),"float");
+		cmd.add( eps_arg );
+		TCLAP::ValueArg<unsigned> minDim_arg("d", "min-ele-dim","Minimum dimension of elements to be inserted into new mesh",
+										false, 1, "unsigned");
+		cmd.add( minDim_arg );
+
+		cmd.parse( argc, argv );
+
+		// read a mesh file
+		MeshLib::Mesh* org_mesh (FileIO::readMeshFromFile(input_arg.getValue()));
+		if (!org_mesh)
+			return 0;
+		INFO("Mesh read: %d nodes, %d elements.", org_mesh->getNNodes(), org_mesh->getNElements());
+
+		// revise the mesh
+		MeshLib::Mesh* new_mesh = nullptr;
+		if (simplify_arg.getValue()) {
+			INFO("Simplifying the mesh...");
+			MeshLib::MeshRevision rev(const_cast<MeshLib::Mesh&>(*org_mesh));
+			unsigned int minDim = (minDim_arg.isSet() ? minDim_arg.getValue() : org_mesh->getDimension());
+			new_mesh = rev.simplifyMesh("revised_mesh", eps_arg.getValue(), minDim);
+		}
+
+		// write into a file
+		if (new_mesh) {
+			INFO("Revised mesh: %d nodes, %d elements.", new_mesh->getNNodes(), new_mesh->getNElements());
+			FileIO::writeMeshToFile(*new_mesh, output_arg.getValue());
+		}
+
+		delete org_mesh;
+		delete new_mesh;
+	} catch (TCLAP::ArgException &e) {
+		ERR("error: %s for arg %s", e.error().c_str(), e.argId().c_str());
+	}
+
+	delete custom_format;
+	delete logog_cout;
+	LOGOG_SHUTDOWN();
+
+	return 0;
+}