std::string name;
gmsh::model::getPhysicalName(entity_dim, physical_tag, name);
if (!name.empty())
- try
- {
- std::map<std::string, int> id_names;
- Patterns::Tools::to_value(name, id_names);
- bool throw_anyway = false;
- bool found_boundary_id = false;
- // If the above did not throw, we keep going, and retrieve
- // all the information that we know how to translate.
- for (const auto &it : id_names)
- {
- const auto &name = it.first;
- const auto &id = it.second;
- if (entity_dim == dim && name == "MaterialID")
- {
- boundary_id = static_cast<types::boundary_id>(id);
- found_boundary_id = true;
- }
- else if (entity_dim < dim && name == "BoundaryID")
- {
- boundary_id = static_cast<types::boundary_id>(id);
- found_boundary_id = true;
- }
- else if (name == "ManifoldID")
- manifold_id = static_cast<types::manifold_id>(id);
- else
- // We did not recognize one of the keys. We'll fall
- // back to setting the boundary id to the physical tag
- // after reading all strings.
- throw_anyway = true;
- }
- // If we didn't find a BoundaryID:XX or MaterialID:XX, and
- // something was found but not recognized, then we set the
- // material id or boundary id in the catch block below,
- // using directly the physical tag
- if (throw_anyway && !found_boundary_id)
- throw;
- }
- catch (...)
- {
- // When the above didn't work, we revert to the old
- // behaviour: the physical tag itself is interpreted either
- // as a material_id or a boundary_id, and no manifold id is
- // known
- boundary_id = physical_tag;
- }
+ {
+ // Patterns::Tools::to_value throws an exception, if it can not
+ // convert name to a map from string to int.
+ try
+ {
+ std::map<std::string, int> id_names;
+ Patterns::Tools::to_value(name, id_names);
+ bool found_unrecognized_tag = false;
+ bool found_boundary_id = false;
+ // If the above did not throw, we keep going, and retrieve
+ // all the information that we know how to translate.
+ for (const auto &it : id_names)
+ {
+ const auto &name = it.first;
+ const auto &id = it.second;
+ if (entity_dim == dim && name == "MaterialID")
+ {
+ boundary_id = static_cast<types::boundary_id>(id);
+ found_boundary_id = true;
+ }
+ else if (entity_dim < dim && name == "BoundaryID")
+ {
+ boundary_id = static_cast<types::boundary_id>(id);
+ found_boundary_id = true;
+ }
+ else if (name == "ManifoldID")
+ manifold_id = static_cast<types::manifold_id>(id);
+ else
+ // We did not recognize one of the keys. We'll fall
+ // back to setting the boundary id to the physical tag
+ // after reading all strings.
+ found_unrecognized_tag = true;
+ }
+ // If we didn't find a BoundaryID:XX or MaterialID:XX, and
+ // something was found but not recognized, then we set the
+ // material id or boundary using the physical tag directly.
+ if (found_unrecognized_tag && !found_boundary_id)
+ boundary_id = physical_tag;
+ }
+ catch (...)
+ {
+ // When the above didn't work, we revert to the old
+ // behaviour: the physical tag itself is interpreted either
+ // as a material_id or a boundary_id, and no manifold id is
+ // known
+ boundary_id = physical_tag;
+ }
+ }
}
// Get the mesh elements for the entity (dim, tag):
--- /dev/null
+// ------------------------------------------------------------------------
+//
+// SPDX-License-Identifier: LGPL-2.1-or-later
+// Copyright (C) 2021 - 2023 by the deal.II authors
+//
+// This file is part of the deal.II library.
+//
+// Part of the source code is dual licensed under Apache-2.0 WITH
+// LLVM-exception OR LGPL-2.1-or-later. Detailed license information
+// governing the source code and code contributions can be found in
+// LICENSE.md and CONTRIBUTING.md at the top level directory of deal.II.
+//
+// ------------------------------------------------------------------------
+
+// Check that gmsh api correctly reads and writes a mesh with manifold
+// information, in all coordinate dimensions
+
+#include <deal.II/grid/grid_generator.h>
+#include <deal.II/grid/grid_in.h>
+#include <deal.II/grid/grid_out.h>
+#include <deal.II/grid/tria.h>
+
+#include <boost/algorithm/string.hpp>
+
+#include <fstream>
+
+#include "../tests.h"
+
+/*
+ * Test that we can use the gmsh-API version of GridIn::read_msh, when the
+ * description of the "Physical name" does not contain MaterialID, BoundaryID,
+ * or ManifoldID.
+ *
+ * Create a hypercube triangulation, write it to file so that we get a valid
+ * msh-file. Read the filecontent back in and replace MaterialId with a string
+ * that will not be recognized. Write it back to the file and try to read the
+ * triangulation from the msh-file.
+ */
+template <int dim>
+void
+test()
+{
+ // Create a mesh, write it to file.
+ Triangulation<dim> original_mesh;
+ const double left = 0.;
+ const double right = 1.;
+ const bool colorize = true;
+ GridGenerator::hyper_cube(original_mesh, left, right, colorize);
+ // Set materialID to something nonzero so that we
+ // get a physical name "MaterialID: 1" in the msh-file.
+ auto cell_on_original = original_mesh.begin_active();
+ cell_on_original->set_material_id(17);
+
+ const std::string mesh_filename = "output.msh";
+ GridOut grid_out;
+ grid_out.write_msh(original_mesh, mesh_filename);
+
+ // Read the file back and change the physical group description to be
+ // unrecognizable.
+ std::ifstream input(mesh_filename, std::ifstream::in);
+ std::stringstream buffert;
+ buffert << input.rdbuf();
+ std::string filecontent = buffert.str();
+ input.close();
+
+ boost::replace_all(filecontent, "MaterialID", "Unrecognizable");
+ boost::replace_all(filecontent, "BoundaryID", "Unrecognizable");
+ boost::replace_all(filecontent, "ManifoldID", "Unrecognizable");
+
+ // Write the mesh to file and read the triangulation back in.
+ std::ofstream output(mesh_filename);
+ output << filecontent;
+ output.close();
+
+ Triangulation<dim> triangulation_from_file;
+ GridIn<dim> grid_in(triangulation_from_file);
+ grid_in.read_msh(mesh_filename);
+
+ // The function should set material_id and boundary_id based on the physical
+ // tag. Write these to log to check them.
+ const auto cell = triangulation_from_file.begin_active();
+ deallog << "MaterialID = " << cell->material_id() << std::endl;
+ for (unsigned int f : cell->face_indices())
+ deallog << "BoundaryID = " << cell->face(f)->boundary_id() << std::endl;
+}
+
+int
+main(int argc, char **argv)
+{
+ // gmsh might be build with mpi support enabled.
+ Utilities::MPI::MPI_InitFinalize mpi_initialization(argc, argv);
+ initlog();
+
+ test<2>();
+}
projects / dealii.git / commitdiff