27int main(
int argc,
char* argv[])
30 "Query mesh information.\n\n"
31 "OpenGeoSys-6 software, version " +
34 "Copyright (c) 2012-2025, OpenGeoSys Community "
35 "(http://www.opengeosys.org)",
37 TCLAP::UnlabeledValueArg<std::string> mesh_arg(
38 "mesh-file",
"input mesh file",
true,
"",
"string");
40 TCLAP::MultiArg<std::size_t> eleId_arg(
"e",
"element-id",
"element ID",
43 TCLAP::MultiArg<std::size_t> nodeId_arg(
"n",
"node-id",
"node ID",
false,
46 TCLAP::SwitchArg showNodeWithMaxEle_arg(
47 "",
"show-node-with-max-elements",
48 "show a node having the max number of connected elements",
false);
49 cmd.add(showNodeWithMaxEle_arg);
51 cmd.parse(argc, argv);
55 const std::string filename(mesh_arg.getValue());
66 std::vector<std::size_t> selected_node_ids;
67 if (showNodeWithMaxEle_arg.getValue())
69 auto itr = std::max_element(
70 mesh->getNodes().begin(), mesh->getNodes().end(),
73 return mesh->getElementsConnectedToNode(*i).size() <
74 mesh->getElementsConnectedToNode(*j).size();
76 if (itr != mesh->getNodes().end())
79 selected_node_ids.push_back(node->
getID());
82 selected_node_ids.insert(selected_node_ids.end(),
83 nodeId_arg.getValue().begin(),
84 nodeId_arg.getValue().end());
86 auto const materialIds = materialIDs(*mesh);
87 for (
auto ele_id : eleId_arg.getValue())
89 std::stringstream out;
90 out << std::scientific
91 << std::setprecision(std::numeric_limits<double>::max_digits10);
92 out <<
"--------------------------------------------------------"
94 auto* ele = mesh->getElement(ele_id);
95 out <<
"# Element " << ele->getID() << std::endl;
96 out <<
"Type : " << CellType2String(ele->getCellType()) << std::endl;
99 out <<
"Mat ID : " << (*materialIds)[ele_id] << std::endl;
101 out <<
"Nodes: " << std::endl;
102 for (
unsigned i = 0; i < ele->getNumberOfNodes(); i++)
104 out << ele->getNode(i)->getID() <<
" " << *ele->getNode(i)
107 out <<
"Content: " << ele->getContent() << std::endl;
108 out <<
"Neighbors: ";
109 for (
unsigned i = 0; i < ele->getNumberOfNeighbors(); i++)
111 if (ele->getNeighbor(i))
113 out << ele->getNeighbor(i)->getID() <<
" ";
121 INFO(
"{:s}", out.str());
125 for (
auto node_id : selected_node_ids)
127 std::stringstream out;
128 out << std::scientific
129 << std::setprecision(std::numeric_limits<double>::max_digits10);
130 out <<
"--------------------------------------------------------"
133 out <<
"# Node " << node->
getID() << std::endl;
134 out <<
"Coordinates: " << *node << std::endl;
135 out <<
"Connected elements ("
136 << mesh->getElementsConnectedToNode(*node).size() <<
"): ";
137 for (
auto ele : mesh->getElementsConnectedToNode(*node))
139 out << ele->getID() <<
" ";
142 out <<
"Connected nodes (" << connections[node->
getID()].size()
144 for (
auto nd : connections[node->
getID()])
146 out << nd->getID() <<
" ";
149 INFO(
"{:s}", out.str());