OGS
MoveGeometry.cpp File Reference

Detailed Description

A small tool to translate geometries.

Author
Karsten Rink
Date
2015-08-04

Definition in file MoveGeometry.cpp.

#include <tclap/CmdLine.h>
#include <mpi.h>
#include "GeoLib/GEOObjects.h"
#include "GeoLib/IO/XmlIO/Boost/BoostXmlGmlInterface.h"
#include "InfoLib/GitInfo.h"
Include dependency graph for MoveGeometry.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 25 of file MoveGeometry.cpp.

26{
27 TCLAP::CmdLine cmd(
28 "Moves the points of a geometry by a given displacement vector\n\n"
29 "OpenGeoSys-6 software, version " +
31 ".\n"
32 "Copyright (c) 2012-2024, OpenGeoSys Community "
33 "(http://www.opengeosys.org)",
35 TCLAP::ValueArg<double> z_arg("z", "z", "displacement in z direction",
36 false, 0.0, "z-displacement");
37 cmd.add(z_arg);
38 TCLAP::ValueArg<double> y_arg("y", "y", "displacement in y direction",
39 false, 0.0, "y-displacement");
40 cmd.add(y_arg);
41 TCLAP::ValueArg<double> x_arg("x", "x", "displacement in x direction",
42 false, 0.0, "x-displacement");
43 cmd.add(x_arg);
44 TCLAP::ValueArg<std::string> geo_output_arg(
45 "o", "output", "output geometry file (*.gml)", true, "", "output file");
46 cmd.add(geo_output_arg);
47 TCLAP::ValueArg<std::string> geo_input_arg(
48 "i", "input", "input geometry file (*.gml)", true, "", "input file");
49 cmd.add(geo_input_arg);
50 cmd.parse(argc, argv);
51
52#ifdef USE_PETSC
53 MPI_Init(&argc, &argv);
54#endif
55
56 GeoLib::GEOObjects geo_objects;
57 GeoLib::IO::BoostXmlGmlInterface xml(geo_objects);
58 try
59 {
60 if (!xml.readFile(geo_input_arg.getValue()))
61 {
62#ifdef USE_PETSC
63 MPI_Finalize();
64#endif
65 return EXIT_FAILURE;
66 }
67 }
68 catch (std::runtime_error const& err)
69 {
70 ERR("Failed to read file `{:s}'.", geo_input_arg.getValue());
71 ERR("{:s}", err.what());
72#ifdef USE_PETSC
73 MPI_Finalize();
74#endif
75 return EXIT_FAILURE;
76 }
77
78 Eigen::Vector3d displacement = Eigen::Vector3d::Zero();
79 if (x_arg.isSet())
80 {
81 displacement[0] = x_arg.getValue();
82 }
83 if (y_arg.isSet())
84 {
85 displacement[1] = y_arg.getValue();
86 }
87 if (z_arg.isSet())
88 {
89 displacement[2] = z_arg.getValue();
90 }
91
92 auto const geo_name = geo_objects.getGeometryNames()[0];
93
94 std::vector<GeoLib::Point*> const* point_vec =
95 geo_objects.getPointVec(geo_name);
96 std::size_t const n_points = point_vec->size();
97 for (std::size_t i = 0; i < n_points; ++i)
98 {
99 for (std::size_t c = 0; c < 3; ++c)
100 {
101 (*(*point_vec)[i])[c] += displacement[c];
102 }
103 }
104
105 xml.export_name = geo_name;
106 BaseLib::IO::writeStringToFile(xml.writeToString(),
107 geo_output_arg.getValue());
108
109#ifdef USE_PETSC
110 MPI_Finalize();
111#endif
112 return EXIT_SUCCESS;
113}
void ERR(fmt::format_string< Args... > fmt, Args &&... args)
Definition Logging.h:45
Container class for geometric objects.
Definition GEOObjects.h:57
std::vector< std::string > getGeometryNames() const
Returns the names of all geometry vectors.
const std::vector< Point * > * getPointVec(const std::string &name) const
int writeStringToFile(std::string_view content, std::filesystem::path const &file_path)
Definition Writer.cpp:45
GITINFOLIB_EXPORT const std::string ogs_version

References ERR(), BaseLib::IO::XMLInterface::export_name, GeoLib::GEOObjects::getGeometryNames(), GeoLib::GEOObjects::getPointVec(), GitInfoLib::GitInfo::ogs_version, GeoLib::IO::BoostXmlGmlInterface::readFile(), BaseLib::IO::writeStringToFile(), and BaseLib::IO::Writer::writeToString().