OGS
queryMesh.cpp File Reference

Detailed Description

Definition in file queryMesh.cpp.

#include <tclap/CmdLine.h>
#include <array>
#include <iomanip>
#include <memory>
#include <sstream>
#include <string>
#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/readMeshFromFile.h"
#include "MeshLib/Mesh.h"
#include "MeshLib/Node.h"
Include dependency graph for queryMesh.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 29 of file queryMesh.cpp.

30{
31 TCLAP::CmdLine cmd(
32 "Query mesh information.\n\n"
33 "OpenGeoSys-6 software, version " +
35 ".\n"
36 "Copyright (c) 2012-2025, OpenGeoSys Community "
37 "(http://www.opengeosys.org)",
39 TCLAP::UnlabeledValueArg<std::string> mesh_arg(
40 "mesh-file", "Input (.vtu) mesh file", true, "", "INPUT_FILE");
41 cmd.add(mesh_arg);
42 TCLAP::MultiArg<std::size_t> eleId_arg("e", "element-id", "element ID",
43 false, "ELEMENT_ID");
44 cmd.add(eleId_arg);
45 TCLAP::MultiArg<std::size_t> nodeId_arg("n", "node-id", "node ID", false,
46 "NODE_ID");
47 cmd.add(nodeId_arg);
48 TCLAP::SwitchArg showNodeWithMaxEle_arg(
49 "", "show-node-with-max-elements",
50 "show a node having the max number of connected elements", false);
51 cmd.add(showNodeWithMaxEle_arg);
52
53 auto log_level_arg = BaseLib::makeLogLevelArg();
54 cmd.add(log_level_arg);
55 cmd.parse(argc, argv);
56
57 BaseLib::MPI::Setup mpi_setup(argc, argv);
58 BaseLib::initOGSLogger(log_level_arg.getValue());
59
60 const std::string filename(mesh_arg.getValue());
61
62 // read the mesh file
63 auto const mesh =
64 std::unique_ptr<MeshLib::Mesh>(MeshLib::IO::readMeshFromFile(
65 filename, true /* compute_element_neighbors */));
66 if (!mesh)
67 {
68 return EXIT_FAILURE;
69 }
70
71 std::vector<std::size_t> selected_node_ids;
72 if (showNodeWithMaxEle_arg.getValue())
73 {
74 auto itr = std::max_element(
75 mesh->getNodes().begin(), mesh->getNodes().end(),
76 [&mesh](MeshLib::Node* i, MeshLib::Node* j)
77 {
78 return mesh->getElementsConnectedToNode(*i).size() <
79 mesh->getElementsConnectedToNode(*j).size();
80 });
81 if (itr != mesh->getNodes().end())
82 {
83 MeshLib::Node* node = *itr;
84 selected_node_ids.push_back(node->getID());
85 }
86 }
87 selected_node_ids.insert(selected_node_ids.end(),
88 nodeId_arg.getValue().begin(),
89 nodeId_arg.getValue().end());
90
91 auto const materialIds = materialIDs(*mesh);
92 for (auto ele_id : eleId_arg.getValue())
93 {
94 std::stringstream out;
95 out << std::scientific
96 << std::setprecision(std::numeric_limits<double>::max_digits10);
97 out << "--------------------------------------------------------"
98 << std::endl;
99 auto* ele = mesh->getElement(ele_id);
100 out << "# Element " << ele->getID() << std::endl;
101 out << "Type : " << CellType2String(ele->getCellType()) << std::endl;
102 if (materialIds)
103 {
104 out << "Mat ID : " << (*materialIds)[ele_id] << std::endl;
105 }
106 out << "Nodes: " << std::endl;
107 for (unsigned i = 0; i < ele->getNumberOfNodes(); i++)
108 {
109 out << ele->getNode(i)->getID() << " " << *ele->getNode(i)
110 << std::endl;
111 }
112 out << "Content: " << ele->getContent() << std::endl;
113 out << "Neighbors: ";
114 for (unsigned i = 0; i < ele->getNumberOfNeighbors(); i++)
115 {
116 if (ele->getNeighbor(i))
117 {
118 out << ele->getNeighbor(i)->getID() << " ";
119 }
120 else
121 {
122 out << "none ";
123 }
124 }
125 out << std::endl;
126 INFO("{:s}", out.str());
127 }
128
129 auto const& connections = MeshLib::calculateNodesConnectedByElements(*mesh);
130 for (auto node_id : selected_node_ids)
131 {
132 std::stringstream out;
133 out << std::scientific
134 << std::setprecision(std::numeric_limits<double>::max_digits10);
135 out << "--------------------------------------------------------"
136 << std::endl;
137 MeshLib::Node const* node = mesh->getNode(node_id);
138 out << "# Node " << node->getID() << std::endl;
139 out << "Coordinates: " << *node << std::endl;
140 out << "Connected elements ("
141 << mesh->getElementsConnectedToNode(*node).size() << "): ";
142 for (auto ele : mesh->getElementsConnectedToNode(*node))
143 {
144 out << ele->getID() << " ";
145 }
146 out << std::endl;
147 out << "Connected nodes (" << connections[node->getID()].size()
148 << "): ";
149 for (auto nd : connections[node->getID()])
150 {
151 out << nd->getID() << " ";
152 }
153 out << std::endl;
154 INFO("{:s}", out.str());
155 }
156 return EXIT_SUCCESS;
157}
void INFO(fmt::format_string< Args... > fmt, Args &&... args)
Definition Logging.h:36
std::string getValue(std::string const &line, std::string const &val_name, bool is_string)
std::size_t getID() const
TCLAP::ValueArg< std::string > makeLogLevelArg()
void initOGSLogger(std::string const &log_level)
Definition Logging.cpp:64
GITINFOLIB_EXPORT const std::string ogs_version
MeshLib::Mesh * readMeshFromFile(const std::string &file_name, bool const compute_element_neighbors)
std::vector< std::vector< Node * > > calculateNodesConnectedByElements(Mesh const &mesh)
Definition Mesh.cpp:309
PropertyVector< int > const * materialIDs(Mesh const &mesh)
Definition Mesh.cpp:269
std::string CellType2String(const CellType t)
Given a MeshElemType this returns the appropriate string.

References MeshLib::calculateNodesConnectedByElements(), MathLib::Point3dWithID::getID(), INFO(), BaseLib::initOGSLogger(), BaseLib::makeLogLevelArg(), GitInfoLib::GitInfo::ogs_version, and MeshLib::IO::readMeshFromFile().