OGS
ExtractBoundary.cpp File Reference
#include <tclap/CmdLine.h>
#include <algorithm>
#include <memory>
#include <string>
#include <vector>
#include "BaseLib/FileTools.h"
#include "BaseLib/Logging.h"
#include "BaseLib/MPI.h"
#include "BaseLib/StringTools.h"
#include "BaseLib/TCLAPArguments.h"
#include "InfoLib/GitInfo.h"
#include "MeshLib/Elements/Element.h"
#include "MeshLib/IO/VtkIO/VtuInterface.h"
#include "MeshLib/IO/readMeshFromFile.h"
#include "MeshLib/IO/writeMeshToFile.h"
#include "MeshLib/Mesh.h"
#include "MeshLib/Node.h"
#include "MeshToolsLib/MeshSurfaceExtraction.h"
Include dependency graph for ExtractBoundary.cpp:

Go to the source code of this file.

Functions

int main (int argc, char *argv[])

Function Documentation

◆ main()

int main ( int argc,
char * argv[] )

Definition at line 25 of file ExtractBoundary.cpp.

26{
27 TCLAP::CmdLine cmd(
28 "Tool extracts the boundary of the given mesh. The documentation is "
29 "available at "
30 "https://docs.opengeosys.org/docs/tools/meshing-submeshes/"
31 "extract-boundary.\n\n"
32 "OpenGeoSys-6 software, version " +
33 GitInfoLib::GitInfo::ogs_version +
34 ".\n"
35 "Copyright (c) 2012-2026, OpenGeoSys Community "
36 "(http://www.opengeosys.org)",
37 ' ', GitInfoLib::GitInfo::ogs_version);
38 TCLAP::ValueArg<std::string> mesh_in(
39 "i", "mesh-input-file",
40 "Input (.vtu). The name of the file containing the input mesh", true,
41 "", "INPUT_FILE");
42 cmd.add(mesh_in);
43 TCLAP::ValueArg<std::string> mesh_out(
44 "o", "mesh-output-file",
45 "Output (.vtu). The name of the file the surface mesh should be "
46 "written to",
47 false, "", "OUTPUT_FILE");
48 cmd.add(mesh_out);
49
50 TCLAP::SwitchArg use_ascii_arg("", "ascii-output",
51 "If the switch is set use ascii instead of "
52 "binary format for data in the vtu output.",
53 false);
54 cmd.add(use_ascii_arg);
55
56 auto log_level_arg = BaseLib::makeLogLevelArg();
57 cmd.add(log_level_arg);
58 cmd.parse(argc, argv);
59
60 BaseLib::MPI::Setup mpi_setup(argc, argv);
61 BaseLib::initOGSLogger(log_level_arg.getValue());
62
63 std::unique_ptr<MeshLib::Mesh const> mesh(MeshLib::IO::readMeshFromFile(
64 mesh_in.getValue(), true /* compute_element_neighbors */));
65
66 if (!mesh)
67 {
68 return EXIT_FAILURE;
69 }
70
71 INFO("Mesh read: {:d} nodes, {:d} elements.", mesh->getNumberOfNodes(),
72 mesh->getNumberOfElements());
73
74 // extract surface
75 std::unique_ptr<MeshLib::Mesh> surface_mesh(
76 MeshToolsLib::BoundaryExtraction::getBoundaryElementsAsMesh(
77 *mesh,
78 MeshLib::getBulkIDString(MeshLib::MeshItemType::Node),
79 MeshLib::getBulkIDString(MeshLib::MeshItemType::Cell),
80 MeshLib::getBulkIDString(MeshLib::MeshItemType::Face)));
81
82 INFO("Created surface mesh: {:d} nodes, {:d} elements.",
83 surface_mesh->getNumberOfNodes(), surface_mesh->getNumberOfElements());
84
85 std::string out_fname(mesh_out.getValue());
86 if (out_fname.empty())
87 {
88 out_fname = BaseLib::dropFileExtension(mesh_in.getValue()) + "_sfc.vtu";
89 }
90
91 auto const data_mode =
92 use_ascii_arg.getValue() ? vtkXMLWriter::Ascii : vtkXMLWriter::Binary;
93 auto const use_compression = true;
94 if (MeshLib::IO::writeMeshToFile(*surface_mesh, out_fname, {},
95 use_compression, data_mode) != 0)
96 {
97 return EXIT_FAILURE;
98 }
99 return EXIT_SUCCESS;
100}
INFO
void INFO(fmt::format_string< Args... > fmt, Args &&... args)
Definition Logging.h:28
BaseLib::makeLogLevelArg
TCLAP::ValueArg< std::string > makeLogLevelArg()
Definition TCLAPArguments.cpp:19
BaseLib::initOGSLogger
void initOGSLogger(std::string const &log_level)
Definition Logging.cpp:56
BaseLib::dropFileExtension
std::string dropFileExtension(std::string const &filename)
Definition FileTools.cpp:165
GitInfoLib::GitInfo::ogs_version
GITINFOLIB_EXPORT const std::string ogs_version
MeshLib::IO::readMeshFromFile
MeshLib::Mesh * readMeshFromFile(const std::string &file_name, bool const compute_element_neighbors)
Definition readMeshFromFile.cpp:58
MeshLib::IO::writeMeshToFile
int writeMeshToFile(const MeshLib::Mesh &mesh, std::filesystem::path const &file_path, std::set< std::string > output_variable_names, bool const use_compression, int const data_mode)
Definition writeMeshToFile.cpp:18
MeshLib::getBulkIDString
constexpr std::string_view getBulkIDString(MeshItemType mesh_item_type)
Definition MeshLib/Properties.h:183
MeshToolsLib::BoundaryExtraction::getBoundaryElementsAsMesh
std::unique_ptr< MeshLib::Mesh > getBoundaryElementsAsMesh(MeshLib::Mesh const &bulk_mesh, std::string_view subsfc_node_id_prop_name, std::string_view subsfc_element_id_prop_name, std::string_view face_id_prop_name)
Definition MeshSurfaceExtraction.cpp:467

References MeshLib::Cell, BaseLib::dropFileExtension(), MeshLib::Face, MeshToolsLib::BoundaryExtraction::getBoundaryElementsAsMesh(), MeshLib::getBulkIDString(), INFO(), BaseLib::initOGSLogger(), BaseLib::makeLogLevelArg(), MeshLib::Node, GitInfoLib::GitInfo::ogs_version, MeshLib::IO::readMeshFromFile(), and MeshLib::IO::writeMeshToFile().