OGS
checkMesh.cpp File Reference

Detailed Description

Definition in file checkMesh.cpp.

#include <spdlog/fmt/ranges.h>
#include <tclap/CmdLine.h>
#include <array>
#include <string>
#include "BaseLib/FileTools.h"
#include "BaseLib/Logging.h"
#include "BaseLib/MPI.h"
#include "BaseLib/MemWatch.h"
#include "BaseLib/RunTime.h"
#include "BaseLib/StringTools.h"
#include "BaseLib/TCLAPArguments.h"
#include "GeoLib/AABB.h"
#include "InfoLib/GitInfo.h"
#include "MeshLib/Elements/Element.h"
#include "MeshLib/IO/readMeshFromFile.h"
#include "MeshLib/Mesh.h"
#include "MeshLib/Node.h"
#include "MeshToolsLib/MeshInformation.h"
#include "MeshToolsLib/MeshQuality/MeshValidation.h"
Include dependency graph for checkMesh.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 32 of file checkMesh.cpp.

33{
34 TCLAP::CmdLine cmd(
35 "Checks mesh properties.\n\n"
36 "OpenGeoSys-6 software, version " +
38 ".\n"
39 "Copyright (c) 2012-2025, OpenGeoSys Community "
40 "(http://www.opengeosys.org)",
42 TCLAP::UnlabeledValueArg<std::string> mesh_arg(
43 "mesh-file", "Input (.vtu | .msh). Mesh file", true, "", "INPUT_FILE");
44 cmd.add(mesh_arg);
45 TCLAP::SwitchArg valid_arg("v", "validation", "validate the mesh");
46 cmd.add(valid_arg);
47 TCLAP::SwitchArg print_properties_arg(
48 "p", "print_properties", "print properties stored in the mesh");
49 cmd.add(print_properties_arg);
50 auto log_level_arg = BaseLib::makeLogLevelArg();
51 cmd.add(log_level_arg);
52
53 cmd.parse(argc, argv);
54
55 BaseLib::MPI::Setup mpi_setup(argc, argv);
56 BaseLib::initOGSLogger(log_level_arg.getValue());
57
58 // read the mesh file
59 BaseLib::MemWatch mem_watch;
60 const unsigned long mem_without_mesh(mem_watch.getVirtMemUsage());
61 BaseLib::RunTime run_time;
62 run_time.start();
63 std::unique_ptr<MeshLib::Mesh> mesh(MeshLib::IO::readMeshFromFile(
64 mesh_arg.getValue(), true /* compute_element_neighbors */));
65 if (!mesh)
66 {
67 return EXIT_FAILURE;
68 }
69
70 const unsigned long mem_with_mesh(mem_watch.getVirtMemUsage());
71 if (mem_with_mesh > 0)
72 {
73 INFO("Memory size: {} MiB",
74 (mem_with_mesh - mem_without_mesh) / (1024 * 1024));
75 (void)mem_with_mesh;
76 }
77 INFO("Time for reading: {:g} s", run_time.elapsed());
78
79 // Geometric information
80 const GeoLib::AABB aabb(
82 INFO("Axis aligned bounding box: {}", aabb);
83
84 auto const [min, max] = minMaxEdgeLength(mesh->getElements());
85 INFO("Min/max edge lengths: [{:g}, {:g}]", min, max);
86
87 // Element information
88
90
91 if (print_properties_arg.isSet())
92 {
94 INFO("MaterialID-list: [{}]",
96 ", "));
97 }
98
99 if (valid_arg.isSet())
100 {
101 // MeshValidation outputs error messages
102 // Remark: MeshValidation can modify the original mesh
104 }
105 return EXIT_SUCCESS;
106}
void INFO(fmt::format_string< Args... > fmt, Args &&... args)
Definition Logging.h:36
unsigned long getVirtMemUsage()
Definition MemWatch.cpp:59
Count the running time.
Definition RunTime.h:29
double elapsed() const
Get the elapsed time in seconds.
Definition RunTime.h:42
void start()
Start the timer.
Definition RunTime.h:32
Class AABB is an axis aligned bounding box around a given set of geometric points of (template) type ...
Definition AABB.h:56
static GeoLib::AABB getBoundingBox(const MeshLib::Mesh &mesh)
Returns the bounding box of the mesh.
static void writeMeshValidationResults(MeshLib::Mesh &mesh)
static std::vector< int > getMaterialIDs(const MeshLib::Mesh &mesh)
writes out a list of all material IDs that occur in the mesh.
static void writeAllNumbersOfElementTypes(const MeshLib::Mesh &mesh)
writes all numbers of element types
static void writePropertyVectorInformation(const MeshLib::Mesh &mesh)
writes out property vector information
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::pair< double, double > minMaxEdgeLength(std::vector< Element * > const &elements)
Returns the minimum and maximum edge length for given elements.
Definition Mesh.cpp:190

References BaseLib::RunTime::elapsed(), MeshToolsLib::MeshInformation::getBoundingBox(), MeshToolsLib::MeshInformation::getMaterialIDs(), BaseLib::MemWatch::getVirtMemUsage(), INFO(), BaseLib::initOGSLogger(), BaseLib::makeLogLevelArg(), GitInfoLib::GitInfo::ogs_version, MeshLib::IO::readMeshFromFile(), BaseLib::RunTime::start(), MeshToolsLib::MeshInformation::writeAllNumbersOfElementTypes(), MeshToolsLib::MeshInformation::writeMeshValidationResults(), and MeshToolsLib::MeshInformation::writePropertyVectorInformation().