int main (int argc, char* argv[])
{
LOGOG_INITIALIZE();
logog::Cout* logog_cout (new logog::Cout);
BaseLib::LogogSimpleFormatter *const custom_format (new BaseLib::LogogSimpleFormatter);
logog_cout->SetFormatter(*custom_format);
TCLAP::CmdLine cmd("Converting a mesh in FEFLOW file format (ASCII, version 5.4) to a vtk unstructured grid file (new OGS file format) or to the old OGS file format - see options.", ' ', "0.1");
TCLAP::ValueArg<std::string> ogs_mesh_arg(
"o",
"out",
"filename for output mesh (if extension is msh, old OGS fileformat is written)",
true,
"",
"filename as string");
cmd.add(ogs_mesh_arg);
TCLAP::ValueArg<std::string> feflow_mesh_arg(
"i",
"in",
"FEFLOW input file (*.fem)",
true,
"",
"filename as string");
cmd.add(feflow_mesh_arg);
cmd.parse(argc, argv);
// *** read mesh
INFO("Reading %s.", feflow_mesh_arg.getValue().c_str());
#ifndef WIN32
BaseLib::MemWatch mem_watch;
unsigned long mem_without_mesh (mem_watch.getVirtMemUsage());
#endif
BaseLib::RunTime run_time;
run_time.start();
FileIO::FEFLOWInterface feflowIO(nullptr);
MeshLib::Mesh const*const mesh(feflowIO.readFEFLOWFile(feflow_mesh_arg.getValue()));
if (mesh == nullptr) {
INFO("Could not read mesh from %s.", feflow_mesh_arg.getValue().c_str());
return -1;
}
#ifndef WIN32
unsigned long mem_with_mesh (mem_watch.getVirtMemUsage());
INFO("Mem for mesh: %i MB", (mem_with_mesh - mem_without_mesh)/(1024*1024));
#endif
run_time.stop();
INFO("Time for reading: %f seconds.", run_time.elapsed());
INFO("Read %d nodes and %d elements.", mesh->getNNodes(), mesh->getNElements());
// *** write mesh in new format
std::string ogs_mesh_fname(ogs_mesh_arg.getValue());
if (BaseLib::getFileExtension(ogs_mesh_fname).compare("msh") == 0) {
INFO("Writing %s.", ogs_mesh_fname.c_str());
FileIO::Legacy::MeshIO mesh_io;
mesh_io.setMesh(mesh);
mesh_io.writeToFile(ogs_mesh_fname);
} else {
if (BaseLib::getFileExtension(ogs_mesh_fname).compare("vtu") != 0) {
ogs_mesh_fname += ".vtu";
}
INFO("Writing %s.", ogs_mesh_fname.c_str());
FileIO::BoostVtuInterface mesh_io;
mesh_io.setMesh(mesh);
mesh_io.writeToFile(ogs_mesh_fname);
}
INFO("\tDone.");
delete mesh;
}
int main (int argc, char* argv[])
{
LOGOG_INITIALIZE();
logog::Cout* logog_cout (new logog::Cout);
BaseLib::LogogSimpleFormatter *const custom_format (new BaseLib::LogogSimpleFormatter);
logog_cout->SetFormatter(*custom_format);
TCLAP::CmdLine cmd("Converting meshes in gmsh file format (ASCII, version 2.2) to a vtk unstructured grid file (new OGS file format) or to the old OGS file format - see options.", ' ', "0.1");
TCLAP::ValueArg<std::string> ogs_mesh_arg(
"o",
"out",
"filename for output mesh (if extension is msh, old OGS fileformat is written)",
true,
"",
"filename as string");
cmd.add(ogs_mesh_arg);
TCLAP::ValueArg<std::string> gmsh_mesh_arg(
"i",
"in",
"gmsh input file",
true,
"",
"filename as string");
cmd.add(gmsh_mesh_arg);
TCLAP::SwitchArg exclude_lines_arg("e", "exclude-lines",
"if set, lines will not be written to the ogs mesh");
cmd.add(exclude_lines_arg);
cmd.parse(argc, argv);
// *** read mesh
INFO("Reading %s.", gmsh_mesh_arg.getValue().c_str());
#ifndef WIN32
BaseLib::MemWatch mem_watch;
unsigned long mem_without_mesh (mem_watch.getVirtMemUsage());
#endif
BaseLib::RunTime run_time;
run_time.start();
MeshLib::Mesh * mesh(FileIO::GMSHInterface::readGMSHMesh(gmsh_mesh_arg.getValue()));
if (mesh == nullptr) {
INFO("Could not read mesh from %s.", gmsh_mesh_arg.getValue().c_str());
return -1;
}
#ifndef WIN32
unsigned long mem_with_mesh (mem_watch.getVirtMemUsage());
INFO("Mem for mesh: %i MB", (mem_with_mesh - mem_without_mesh)/(1024*1024));
#endif
INFO("Time for reading: %f seconds.", run_time.elapsed());
INFO("Read %d nodes and %d elements.", mesh->getNNodes(), mesh->getNElements());
// *** remove line elements on request
if (exclude_lines_arg.getValue()) {
auto ex = MeshLib::ElementSearch(*mesh);
ex.searchByElementType(MeshLib::MeshElemType::LINE);
auto m = MeshLib::removeElements(*mesh, ex.getSearchedElementIDs(), mesh->getName()+"-withoutLines");
if (m != nullptr) {
INFO("Removed %d lines.", mesh->getNElements() - m->getNElements());
std::swap(m, mesh);
delete m;
} else {
INFO("Mesh does not contain any lines.");
}
}
// *** write mesh in new format
std::string ogs_mesh_fname(ogs_mesh_arg.getValue());
if (BaseLib::getFileExtension(ogs_mesh_fname).compare("msh") == 0) {
INFO("Writing %s.", ogs_mesh_fname.c_str());
FileIO::Legacy::MeshIO mesh_io;
mesh_io.setMesh(mesh);
mesh_io.writeToFile(ogs_mesh_fname);
} else {
if (BaseLib::getFileExtension(ogs_mesh_fname).compare("vtu") != 0) {
ogs_mesh_fname += ".vtu";
}
INFO("Writing %s.", ogs_mesh_fname.c_str());
FileIO::VtuInterface mesh_io(mesh);
mesh_io.writeToFile(ogs_mesh_fname);
}
INFO("\tDone.");
delete mesh;
}
