int main (int argc, char* argv[])
{
ApplicationsLib::LogogSetup logog_setup;
TCLAP::CmdLine cmd(
"Append line elements into a mesh.\n\n"
"OpenGeoSys-6 software, version " +
BaseLib::BuildInfo::ogs_version +
".\n"
"Copyright (c) 2012-2019, OpenGeoSys Community "
"(http://www.opengeosys.org)",
' ', BaseLib::BuildInfo::ogs_version);
TCLAP::ValueArg<std::string> mesh_in("i", "mesh-input-file",
"the name of the file containing the input mesh", true,
"", "file name of input mesh");
cmd.add(mesh_in);
TCLAP::ValueArg<std::string> mesh_out("o", "mesh-output-file",
"the name of the file the mesh will be written to", true,
"", "file name of output mesh");
cmd.add(mesh_out);
TCLAP::ValueArg<std::string> geoFileArg("g", "geo-file",
"the name of the geometry file which contains polylines", true, "", "the name of the geometry file");
cmd.add(geoFileArg);
TCLAP::ValueArg<std::string> gmsh_path_arg("g", "gmsh-path",
"the path to the gmsh binary",
false, "", "path as string");
cmd.add(gmsh_path_arg);
// parse arguments
cmd.parse(argc, argv);
// read GEO objects
GeoLib::GEOObjects geo_objs;
FileIO::readGeometryFromFile(geoFileArg.getValue(), geo_objs,
gmsh_path_arg.getValue());
std::vector<std::string> geo_names;
geo_objs.getGeometryNames (geo_names);
if (geo_names.empty ())
{
ERR("No geometries found.");
return EXIT_FAILURE;
}
const GeoLib::PolylineVec* ply_vec (geo_objs.getPolylineVecObj(geo_names[0]));
if (!ply_vec)
{
ERR("Could not find polylines in geometry '%s'.",
geo_names.front().c_str());
return EXIT_FAILURE;
}
// read a mesh
MeshLib::Mesh const*const mesh (MeshLib::IO::readMeshFromFile(mesh_in.getValue()));
if (!mesh)
{
ERR("Mesh file '%s' not found", mesh_in.getValue().c_str());
return EXIT_FAILURE;
}
INFO("Mesh read: %d nodes, %d elements.", mesh->getNumberOfNodes(), mesh->getNumberOfElements());
// add line elements
std::unique_ptr<MeshLib::Mesh> new_mesh =
MeshGeoToolsLib::appendLinesAlongPolylines(*mesh, *ply_vec);
INFO("Mesh created: %d nodes, %d elements.", new_mesh->getNumberOfNodes(), new_mesh->getNumberOfElements());
MeshLib::IO::writeMeshToFile(*new_mesh, mesh_out.getValue());
return EXIT_SUCCESS;
}
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("Append line elements into a mesh.", ' ', "0.1");
TCLAP::ValueArg<std::string> mesh_in("i", "mesh-input-file",
"the name of the file containing the input mesh", true,
"", "file name of input mesh");
cmd.add(mesh_in);
TCLAP::ValueArg<std::string> mesh_out("o", "mesh-output-file",
"the name of the file the mesh will be written to", true,
"", "file name of output mesh");
cmd.add(mesh_out);
TCLAP::ValueArg<std::string> geoFileArg("g", "geo-file",
"the name of the geometry file which contains polylines", true, "", "the name of the geometry file");
cmd.add(geoFileArg);
// parse arguments
cmd.parse(argc, argv);
// read GEO objects
GeoLib::GEOObjects geo_objs;
FileIO::BoostXmlGmlInterface xml(geo_objs);
xml.readFile(geoFileArg.getValue());
std::vector<std::string> geo_names;
geo_objs.getGeometryNames (geo_names);
if (geo_names.empty ())
{
std::cout << "no geometries found" << std::endl;
return -1;
}
const GeoLib::PolylineVec* ply_vec (geo_objs.getPolylineVecObj(geo_names[0]));
if (!ply_vec)
{
std::cout << "could not found polylines" << std::endl;
return -1;
}
// read a mesh
MeshLib::Mesh const*const mesh (FileIO::readMeshFromFile(mesh_in.getValue()));
if (!mesh)
{
ERR("Mesh file %s not found", mesh_in.getValue().c_str());
return 1;
}
INFO("Mesh read: %d nodes, %d elements.", mesh->getNNodes(), mesh->getNElements());
// add line elements
std::unique_ptr<MeshLib::Mesh> new_mesh =
MeshGeoToolsLib::appendLinesAlongPolylines(*mesh, *ply_vec);
INFO("Mesh created: %d nodes, %d elements.", new_mesh->getNNodes(), new_mesh->getNElements());
// write into a file
FileIO::Legacy::MeshIO meshIO;
meshIO.setMesh(new_mesh.get());
meshIO.writeToFile(mesh_out.getValue());
delete custom_format;
delete logog_cout;
LOGOG_SHUTDOWN();
return 1;
}
