From 51a778282a8160619db81261f932cc019016d8fb Mon Sep 17 00:00:00 2001
From: Norihiro Watanabe <norihiro.watanabe@ufz.de>
Date: Wed, 4 Jun 2014 19:12:08 +0200
Subject: [PATCH] add AppendLinesAlongPolyLine in MeshGeoToolsLib
---
MeshGeoToolsLib/AppendLinesAlongPolyline.cpp | 69 ++++++++++++++++++++
MeshGeoToolsLib/AppendLinesAlongPolyline.h | 25 +++++++
2 files changed, 94 insertions(+)
create mode 100644 MeshGeoToolsLib/AppendLinesAlongPolyline.cpp
create mode 100644 MeshGeoToolsLib/AppendLinesAlongPolyline.h
diff --git a/MeshGeoToolsLib/AppendLinesAlongPolyline.cpp b/MeshGeoToolsLib/AppendLinesAlongPolyline.cpp
new file mode 100644
index 00000000000..074983358b3
--- /dev/null
+++ b/MeshGeoToolsLib/AppendLinesAlongPolyline.cpp
@@ -0,0 +1,69 @@
+
+/**
+ * @copyright
+ * Copyright (c) 2013, 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 "AppendLinesAlongPolyline.h"
+
+// ThirdParty/logog
+#include "logog/include/logog.hpp"
+
+// MeshLib
+#include "Mesh.h"
+#include "Node.h"
+#include "Elements/Line.h"
+#include "Elements/Element.h"
+#include "MeshEnums.h"
+#include "MeshEditing/DuplicateMeshComponents.h"
+
+#include "MeshGeoToolsLib/MeshNodesAlongPolyline.h"
+
+namespace MeshGeoToolsLib
+{
+
+MeshLib::Mesh* appendLinesAlongPolylines(const MeshLib::Mesh &mesh, const GeoLib::PolylineVec &ply_vec)
+{
+
+ // copy existing nodes and elements
+ std::vector<MeshLib::Node*> vec_new_nodes = MeshLib::copyNodeVector(mesh.getNodes());
+ std::vector<MeshLib::Element*> vec_new_eles = MeshLib::copyElementVector(mesh.getElements(), vec_new_nodes);
+ unsigned max_matID = 0;
+ for (auto e : vec_new_eles)
+ max_matID = std::max(max_matID, e->getValue());
+
+ const size_t n_ply (ply_vec.size());
+ // for each polyline
+ for (size_t k(0); k < n_ply; k++)
+ {
+ const GeoLib::Polyline* ply = (*ply_vec.getVector())[k];
+
+ // search nodes on the polyline
+ MeshGeoToolsLib::MeshNodesAlongPolyline mshNodesAlongPoly(mesh, *ply, mesh.getMinEdgeLength()*0.5);
+ auto &vec_nodes_on_ply = mshNodesAlongPoly.getNodeIDs();
+ if (vec_nodes_on_ply.empty()) {
+ std::string ply_name;
+ ply_vec.getNameOfElementByID(k, ply_name);
+ INFO("No nodes found on polyline %s", ply_name.c_str());
+ continue;
+ }
+
+ // add line elements
+ for (size_t i=0; i<vec_nodes_on_ply.size()-1; i++) {
+ std::array<MeshLib::Node*, 2> element_nodes;
+ element_nodes[0] = vec_new_nodes[vec_nodes_on_ply[i]];
+ element_nodes[1] = vec_new_nodes[vec_nodes_on_ply[i+1]];
+ vec_new_eles.push_back(new MeshLib::Line(element_nodes, max_matID+k+1));
+ }
+ }
+
+ // generate a mesh
+ const std::string new_mesh_name = mesh.getName() + "_with_lines";
+ return new MeshLib::Mesh(new_mesh_name, vec_new_nodes, vec_new_eles);
+}
+
+} // MeshGeoToolsLib
+
diff --git a/MeshGeoToolsLib/AppendLinesAlongPolyline.h b/MeshGeoToolsLib/AppendLinesAlongPolyline.h
new file mode 100644
index 00000000000..3f1bae0931a
--- /dev/null
+++ b/MeshGeoToolsLib/AppendLinesAlongPolyline.h
@@ -0,0 +1,25 @@
+/**
+ * @copyright
+ * Copyright (c) 2013, OpenGeoSys Community (http://www.opengeosys.org)
+ * Distributed under a Modified BSD License.
+ * See accompanying file LICENSE.txt or
+ * http://www.opengeosys.org/LICENSE.txt
+ */
+
+#ifndef APPENDLINESALONGPOLYLINES_H_
+#define APPENDLINESALONGPOLYLINES_H_
+
+// GeoLib
+#include "PolylineVec.h"
+
+namespace MeshLib
+{
+class Mesh;
+}
+
+namespace MeshGeoToolsLib
+{
+MeshLib::Mesh* appendLinesAlongPolylines(const MeshLib::Mesh &mesh, const GeoLib::PolylineVec &ply_vec);
+}
+
+#endif
--
GitLab