61{
62 TCLAP::CmdLine cmd(
63 "The tool computes the area per node of the surface mesh and writes "
64 "the information as txt and csv data.\n\n"
65 "OpenGeoSys-6 software, version " +
67 ".\n"
68 "Copyright (c) 2012-2025, OpenGeoSys Community "
69 "(http://www.opengeosys.org)",
71 TCLAP::ValueArg<std::string> mesh_in(
72 "i", "mesh-input-file",
73 "Input (.vtu). The name of the file containing the input mesh", true,
74 "", "INPUT_FILE");
75 cmd.add(mesh_in);
76
77 std::vector<std::string> allowed_types_vector{
78 "bulk_node_ids", "bulk_element_ids", "bulk_edge_ids", "bulk_face_ids"};
79 TCLAP::ValuesConstraint<std::string> allowed_types(allowed_types_vector);
80 TCLAP::ValueArg<std::string> id_prop_name(
81 "", "id-prop-name",
82 "the name of the property containing the id information", false,
83 "bulk_node_ids", &allowed_types);
84 cmd.add(id_prop_name);
85 TCLAP::ValueArg<std::string> out_base_fname(
86 "p", "output-base-name",
87 "Output (.csv | .txt). The path and base file name the output will be "
88 "written to",
89 false, "", "BASE_FILENAME_OUTPUT");
90 cmd.add(out_base_fname);
91
92 cmd.parse(argc, argv);
93
95
96 std::unique_ptr<MeshLib::Mesh> surface_mesh(
98 INFO(
"Mesh read: {:d} nodes, {:d} elements.",
99 surface_mesh->getNumberOfNodes(), surface_mesh->getNumberOfElements());
100
101
103
104
105 if (!surface_mesh->getProperties().existsPropertyVector<std::size_t>(
106 id_prop_name.getValue()))
107 {
108 orig_node_ids =
109 surface_mesh->getProperties().createNewPropertyVector<std::size_t>(
111 surface_mesh->getNumberOfNodes(), 1);
112 if (!orig_node_ids)
113 {
114 ERR(
"Fatal error: could not create property.");
115 return EXIT_FAILURE;
116 }
117 std::iota(orig_node_ids->begin(), orig_node_ids->end(), 0);
118 }
119 else
120 {
121 orig_node_ids =
122 surface_mesh->getProperties().getPropertyVector<std::size_t>(
123 id_prop_name.getValue());
124 }
125
126 std::vector<double> areas(
128 *surface_mesh));
129
130
131 std::vector<std::pair<std::size_t, double>> ids_and_areas;
132 std::transform(orig_node_ids->cbegin(), orig_node_ids->cend(),
133 areas.cbegin(), std::back_inserter(ids_and_areas),
134 std::make_pair<std::size_t const&, double const&>);
135
136
137 std::string path(out_base_fname.getValue());
138 if (path.empty())
139 {
141 }
142 std::string const id_and_area_fname(path + ".txt");
143 std::string const csv_fname(path + ".csv");
144
145 writeToFile(id_and_area_fname, csv_fname, ids_and_areas,
146 surface_mesh->getNodes());
147
148 return EXIT_SUCCESS;
149}
static void writeToFile(std::string const &id_area_fname, std::string const &csv_fname, std::vector< std::pair< std::size_t, double > > const &ids_and_areas, std::vector< MeshLib::Node * > const &mesh_nodes)
void INFO(fmt::format_string< Args... > fmt, Args &&... args)
void ERR(fmt::format_string< Args... > fmt, Args &&... args)
std::string dropFileExtension(std::string const &filename)
GITINFOLIB_EXPORT const std::string ogs_version
MeshLib::Mesh * readMeshFromFile(const std::string &file_name, bool const compute_element_neighbors)