int main(int argc, char *argv[])
{
ApplicationsLib::LogogSetup logog_setup;
TCLAP::CmdLine cmd(
"Query mesh information.\n\n"
"OpenGeoSys-6 software, version " +
BaseLib::BuildInfo::git_describe +
".\n"
"Copyright (c) 2012-2019, OpenGeoSys Community "
"(http://www.opengeosys.org)",
' ', BaseLib::BuildInfo::git_describe);
TCLAP::UnlabeledValueArg<std::string> mesh_arg("mesh-file","input mesh file",true,"","string");
cmd.add( mesh_arg );
TCLAP::MultiArg<std::size_t> eleId_arg("e","element-id","element ID",false,"number");
cmd.add( eleId_arg );
TCLAP::MultiArg<std::size_t> nodeId_arg("n","node-id","node ID",false,"number");
cmd.add( nodeId_arg );
TCLAP::SwitchArg showNodeWithMaxEle_arg("", "show-node-with-max-elements", "show a node having the max number of connected elements", false);
cmd.add( showNodeWithMaxEle_arg );
cmd.parse( argc, argv );
const std::string filename(mesh_arg.getValue());
// read the mesh file
auto const mesh = std::unique_ptr<MeshLib::Mesh>(
MeshLib::IO::readMeshFromFile(filename));
if (!mesh)
{
return EXIT_FAILURE;
}
std::vector<std::size_t> selected_node_ids;
if (showNodeWithMaxEle_arg.getValue())
{
auto itr = std::max_element(mesh->getNodes().begin(), mesh->getNodes().end(),
[](MeshLib::Node* i, MeshLib::Node* j) { return i->getNumberOfElements() < j->getNumberOfElements(); });
if (itr != mesh->getNodes().end())
{
MeshLib::Node* node = *itr;
selected_node_ids.push_back(node->getID());
}
}
selected_node_ids.insert(selected_node_ids.end(), nodeId_arg.getValue().begin(), nodeId_arg.getValue().end());
auto const materialIds = materialIDs(*mesh);
for (auto ele_id : eleId_arg.getValue())
{
std::stringstream out;
out << std::scientific
<< std::setprecision(std::numeric_limits<double>::digits10);
out << "--------------------------------------------------------" << std::endl;
auto* ele = mesh->getElement(ele_id);
out << "# Element " << ele->getID() << std::endl;
out << "Type : " << CellType2String(ele->getCellType()) << std::endl;
if (materialIds)
{
out << "Mat ID : " << (*materialIds)[ele_id] << std::endl;
}
out << "Nodes: " << std::endl;
for (unsigned i = 0; i < ele->getNumberOfNodes(); i++)
{
out << ele->getNode(i)->getID() << " " << *ele->getNode(i)
<< std::endl;
}
out << "Content: " << ele->getContent() << std::endl;
out << "Neighbors: ";
for (unsigned i=0; i<ele->getNumberOfNeighbors(); i++)
{
if (ele->getNeighbor(i))
{
out << ele->getNeighbor(i)->getID() << " ";
}
else
{
out << "none ";
}
}
out << std::endl;
INFO("%s", out.str().c_str());
}
for (auto node_id : selected_node_ids)
{
std::stringstream out;
out << std::scientific
<< std::setprecision(std::numeric_limits<double>::digits10);
out << "--------------------------------------------------------" << std::endl;
MeshLib::Node const* node = mesh->getNode(node_id);
out << "# Node " << node->getID() << std::endl;
out << "Coordinates: " << *node << std::endl;
out << "Connected elements (" << node->getNumberOfElements() << "): ";
for (auto ele : node->getElements())
{
out << ele->getID() << " ";
}
out << std::endl;
out << "Connected nodes (" << node->getConnectedNodes().size() << "): ";
for (auto nd : node->getConnectedNodes())
{
out << nd->getID() << " ";
}
out << std::endl;
INFO("%s", out.str().c_str());
}
}
int main(int argc, char *argv[])
{
ApplicationsLib::LogogSetup logog_setup;
TCLAP::CmdLine cmd("Query mesh information", ' ', BaseLib::BuildInfo::git_describe);
TCLAP::UnlabeledValueArg<std::string> mesh_arg("mesh-file","input mesh file",true,"","string");
cmd.add( mesh_arg );
TCLAP::MultiArg<std::size_t> eleId_arg("e","element-id","element ID",false,"number");
cmd.add( eleId_arg );
TCLAP::MultiArg<std::size_t> nodeId_arg("n","node-id","node ID",false,"number");
cmd.add( nodeId_arg );
cmd.parse( argc, argv );
const std::string filename(mesh_arg.getValue());
// read the mesh file
auto const mesh = std::unique_ptr<MeshLib::Mesh>(
MeshLib::IO::readMeshFromFile(filename));
if (!mesh)
return EXIT_FAILURE;
auto materialIds = mesh->getProperties().getPropertyVector<int>("MaterialIDs");
for (auto ele_id : eleId_arg.getValue())
{
std::stringstream out;
out << std::scientific
<< std::setprecision(std::numeric_limits<double>::digits10);
out << "--------------------------------------------------------" << std::endl;
auto* ele = mesh->getElement(ele_id);
out << "# Element " << ele->getID() << std::endl;
out << "Type : " << CellType2String(ele->getCellType()) << std::endl;
if (materialIds)
out << "Mat ID : " << (*materialIds)[ele_id] << std::endl;
out << "Nodes: " << std::endl;
for (unsigned i=0; i<ele->getNNodes(); i++)
out << ele->getNode(i)->getID() << " " << *ele->getNode(i) << std::endl;
out << "Content: " << ele->getContent() << std::endl;
out << "Neighbors: ";
for (unsigned i=0; i<ele->getNNeighbors(); i++)
{
if (ele->getNeighbor(i))
out << ele->getNeighbor(i)->getID() << " ";
else
out << "none ";
}
out << std::endl;
INFO("%s", out.str().c_str());
}
for (auto node_id : nodeId_arg.getValue())
{
std::stringstream out;
out << std::scientific
<< std::setprecision(std::numeric_limits<double>::digits10);
out << "--------------------------------------------------------" << std::endl;
auto* node = mesh->getNode(node_id);
out << "# Node" << node->getID() << std::endl;
out << "Coordinates: " << *node << std::endl;
out << "Connected elements: " ;
for (unsigned i=0; i<node->getNElements(); i++)
out << node->getElement(i)->getID() << " ";
out << std::endl;
INFO("%s", out.str().c_str());
}
}
int main(int argc, char *argv[])
{
LOGOG_INITIALIZE();
logog::Cout* logog_cout (new logog::Cout);
BaseLib::LogogSimpleFormatter *custom_format (new BaseLib::LogogSimpleFormatter);
logog_cout->SetFormatter(*custom_format);
TCLAP::CmdLine cmd("Query mesh information", ' ', BaseLib::BuildInfo::git_describe);
TCLAP::UnlabeledValueArg<std::string> mesh_arg("mesh-file","input mesh file",true,"","string");
cmd.add( mesh_arg );
TCLAP::MultiArg<std::size_t> eleId_arg("e","element-id","element ID",false,"number");
cmd.add( eleId_arg );
TCLAP::MultiArg<std::size_t> nodeId_arg("n","node-id","node ID",false,"number");
cmd.add( nodeId_arg );
cmd.parse( argc, argv );
const std::string filename(mesh_arg.getValue());
// read the mesh file
const MeshLib::Mesh* mesh = FileIO::readMeshFromFile(filename);
if (!mesh)
return 1;
std::cout << std::scientific << std::setprecision(12);
for (auto ele_id : eleId_arg.getValue())
{
std::cout << "--------------------------------------------------------" << std::endl;
auto* ele = mesh->getElement(ele_id);
std::cout << "# Element " << ele->getID() << std::endl;
std::cout << "Type : " << CellType2String(ele->getCellType()) << std::endl;
std::cout << "Mat ID : " << ele->getValue() << std::endl;
std::cout << "Nodes: " << std::endl;
for (unsigned i=0; i<ele->getNNodes(); i++)
std::cout << ele->getNode(i)->getID() << " " << *ele->getNode(i) << std::endl;
std::cout << "Content: " << ele->getContent() << std::endl;
std::cout << "Neighbors: ";
for (unsigned i=0; i<ele->getNNeighbors(); i++)
{
if (ele->getNeighbor(i))
std::cout << ele->getNeighbor(i)->getID() << " ";
else
std::cout << "none ";
}
std::cout << std::endl;
}
for (auto node_id : nodeId_arg.getValue())
{
std::cout << "--------------------------------------------------------" << std::endl;
auto* node = mesh->getNode(node_id);
std::cout << "# Node" << node->getID() << std::endl;
std::cout << "Coordinates: " << *node << std::endl;
std::cout << "Connected elements: " ;
for (unsigned i=0; i<node->getNElements(); i++)
std::cout << node->getElement(i)->getID() << " ";
std::cout << std::endl;
}
delete mesh;
delete custom_format;
delete logog_cout;
LOGOG_SHUTDOWN();
}
