30 std::string
const& id_area_fname, std::string
const& csv_fname,
31 std::vector<std::pair<std::size_t, double>>
const& ids_and_areas,
32 std::vector<MeshLib::Node*>
const& mesh_nodes)
34 std::ofstream ids_and_area_out(id_area_fname);
35 if (!ids_and_area_out)
37 OGS_FATAL(
"Unable to open the file '{:s}' - aborting.", id_area_fname);
39 std::ofstream csv_out(csv_fname);
42 OGS_FATAL(
"Unable to open the file '{:s}' - aborting.", csv_fname);
45 ids_and_area_out.precision(std::numeric_limits<double>::max_digits10);
46 csv_out.precision(std::numeric_limits<double>::max_digits10);
48 csv_out <<
"ID x y z area node_id\n";
49 for (std::size_t k(0); k < ids_and_areas.size(); k++)
51 ids_and_area_out << ids_and_areas[k].first <<
" "
52 << ids_and_areas[k].second <<
"\n";
53 csv_out << k <<
" " << *(mesh_nodes[k]) << ids_and_areas[k].second
54 <<
" " << ids_and_areas[k].first <<
"\n";
56 ids_and_area_out <<
"\n";
60int main(
int argc,
char* argv[])
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 " +
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,
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(
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 "
89 false,
"",
"BASE_FILENAME_OUTPUT");
90 cmd.add(out_base_fname);
92 cmd.parse(argc, argv);
96 std::unique_ptr<MeshLib::Mesh> surface_mesh(
98 INFO(
"Mesh read: {:d} nodes, {:d} elements.",
99 surface_mesh->getNumberOfNodes(), surface_mesh->getNumberOfElements());
105 if (!surface_mesh->getProperties().existsPropertyVector<std::size_t>(
106 id_prop_name.getValue()))
109 surface_mesh->getProperties().createNewPropertyVector<std::size_t>(
111 surface_mesh->getNumberOfNodes(), 1);
114 ERR(
"Fatal error: could not create property.");
117 std::iota(orig_node_ids->
begin(), orig_node_ids->
end(), 0);
122 surface_mesh->getProperties().getPropertyVector<std::size_t>(
123 id_prop_name.getValue());
126 std::vector<double> areas(
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&>);
137 std::string path(out_base_fname.getValue());
142 std::string
const id_and_area_fname(path +
".txt");
143 std::string
const csv_fname(path +
".csv");
145 writeToFile(id_and_area_fname, csv_fname, ids_and_areas,
146 surface_mesh->getNodes());
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)