53 int const number_of_subdivisions,
54 std::string&& polyline_name,
55 std::string&& geometry_name,
59 point0, point1, number_of_subdivisions);
60 std::vector<GeoLib::Point*> points(intermediate_points.begin(),
61 intermediate_points.end());
62 geometry.
addPointVec(std::move(points), geometry_name,
66 std::vector<std::size_t> polyline_point_ids(point_vec.getVector().size());
67 std::iota(begin(polyline_point_ids), end(polyline_point_ids), 0);
71 std::move(polyline_name));
78 std::array<GeoLib::Point, 4>
const& points,
79 std::array<int, 4>
const& number_of_subdivisions_per_edge)
81 std::vector<GeoLib::Point*> quad_points;
83 auto addPointsOnLine = [&quad_points](
auto const& begin,
auto const& end,
84 auto const number_of_subdivisions)
87 begin, end, number_of_subdivisions);
88 quad_points.insert(quad_points.end(), intermediate_points.begin(),
89 --intermediate_points.end());
90 delete intermediate_points.back();
94 addPointsOnLine(points[0], points[1], number_of_subdivisions_per_edge[0]);
95 addPointsOnLine(points[1], points[2], number_of_subdivisions_per_edge[1]);
96 addPointsOnLine(points[2], points[3], number_of_subdivisions_per_edge[2]);
97 addPointsOnLine(points[3], points[0], number_of_subdivisions_per_edge[3]);
104 int const number_of_subdivisions_first_x,
105 int const number_of_subdivisions_second_x,
106 int const number_of_subdivisions_first_y,
107 int const number_of_subdivisions_second_y,
108 int const number_of_subdivisions_first_z,
109 int const number_of_subdivisions_second_z,
110 std::string&& quad_name, std::string&& geometry_name,
113 std::array<GeoLib::Point, 4> edge_points;
114 edge_points[0] = point0;
115 edge_points[2] = point1;
116 std::array<int, 4> number_of_subdivisions;
117 if (point0[0] != point1[0] && point0[1] != point1[1] &&
118 point0[2] == point1[2])
125 number_of_subdivisions = {
126 number_of_subdivisions_first_x, number_of_subdivisions_first_y,
127 number_of_subdivisions_second_x, number_of_subdivisions_second_y};
129 else if (point0[0] != point1[0] && point0[1] == point1[1] &&
130 point0[2] != point1[2])
137 number_of_subdivisions = {
138 number_of_subdivisions_first_x, number_of_subdivisions_first_z,
139 number_of_subdivisions_second_x, number_of_subdivisions_second_z};
141 else if (point0[0] == point1[0] && point0[1] != point1[1] &&
142 point0[2] != point1[2])
149 number_of_subdivisions = {
150 number_of_subdivisions_first_y, number_of_subdivisions_first_z,
151 number_of_subdivisions_second_y, number_of_subdivisions_second_z};
155 ERR(
"Input coordinates don't describe an axis aligned polyline or "
165 std::vector<std::size_t> polyline_point_ids(point_vec.getVector().size() +
167 std::iota(begin(polyline_point_ids), end(polyline_point_ids), 0);
168 polyline_point_ids.back() = polyline_point_ids.front();
172 std::move(quad_name));
175 std::move(name_map));
179int main(
int argc,
char* argv[])
182 "Generate point, axis parallel polyline or axis parallel quad "
184 "OpenGeoSys-6 software, version " +
187 "Copyright (c) 2012-2025, OpenGeoSys Community "
188 "(http://www.opengeosys.org)",
190 TCLAP::ValueArg<double> z0(
"",
"z0",
"z coordinate of the first point",
193 TCLAP::ValueArg<double> y0(
"",
"y0",
"y coordinate of the first point",
196 TCLAP::ValueArg<double> x0(
"",
"x0",
"x coordinate of the first point",
199 TCLAP::ValueArg<double> z1(
"",
"z1",
"z coordinate of the first point",
202 TCLAP::ValueArg<double> y1(
"",
"y1",
"y coordinate of the first point",
205 TCLAP::ValueArg<double> x1(
"",
"x1",
"x coordinate of the first point",
208 TCLAP::ValueArg<int> nx(
"",
"nx",
"number of subdivisions in x direction",
211 TCLAP::ValueArg<int> nx1(
"",
"nx1",
"number of subdivisions in x direction",
214 TCLAP::ValueArg<int> ny(
"",
"ny",
"number of subdivisions in y direction",
217 TCLAP::ValueArg<int> ny1(
"",
"ny1",
"number of subdivisions in y direction",
220 TCLAP::ValueArg<int> nz(
"",
"nz",
"number of subdivisions in z direction",
221 false, 0,
"number of subdivisions in z direction");
223 TCLAP::ValueArg<int> nz1(
"",
"nz1",
"number of subdivisions in z direction",
226 TCLAP::ValueArg<std::string> geometry_name(
227 "",
"geometry_name",
"name of the generated geometry",
false,
228 "conceptual model",
"GEOMETRY_NAME");
229 cmd.add(geometry_name);
230 TCLAP::ValueArg<std::string> polyline_name(
231 "",
"polyline_name",
"name of the generated polyline",
false,
232 "polyline",
"POLYLINE_NAME");
233 cmd.add(polyline_name);
234 TCLAP::ValueArg<std::string> geo_output_arg(
235 "o",
"output",
"Output (.gml) geometry file",
true,
"",
"OUTPUT_FILE");
236 cmd.add(geo_output_arg);
238 cmd.add(log_level_arg);
239 cmd.parse(argc, argv);
244 auto const p0 =
GeoLib::Point{x0.getValue(), y0.getValue(), z0.getValue()};
245 auto const p1 =
GeoLib::Point{x1.getValue(), y1.getValue(), z1.getValue()};
248 auto constexpr eps = std::numeric_limits<double>::epsilon();
249 if (p1[0] - p0[0] < eps && p1[1] - p0[1] < eps && p1[2] - p0[2] < eps)
252 std::string(geometry_name.getValue()),
255 else if ((p1[0] - p0[0] >= eps && p1[1] - p0[1] < eps &&
256 p1[2] - p0[2] < eps) ||
257 (p1[0] - p0[0] < eps && p1[1] - p0[1] >= eps &&
258 p1[2] - p0[2] < eps) ||
259 (p1[0] - p0[0] < eps && p1[1] - p0[1] < eps &&
260 p1[2] - p0[2] >= eps))
263 p0, p1, nx.getValue(), std::string(polyline_name.getValue()),
264 std::string(geometry_name.getValue()), geometry);
268 auto eval = [](
int v,
int v1)
280 p0, p1, nx.getValue(), eval(nx.getValue(), nx1.getValue()),
281 ny.getValue(), eval(ny.getValue(), ny1.getValue()),
282 nz.getValue(), eval(nz.getValue(), nz1.getValue()),
283 std::string(polyline_name.getValue()),
284 std::string(geometry_name.getValue()),
285 geometry) == EXIT_FAILURE)
294 geo_output_arg.getValue());
int generateQuadGeometry(GeoLib::Point const &point0, GeoLib::Point const &point1, int const number_of_subdivisions_first_x, int const number_of_subdivisions_second_x, int const number_of_subdivisions_first_y, int const number_of_subdivisions_second_y, int const number_of_subdivisions_first_z, int const number_of_subdivisions_second_z, std::string &&quad_name, std::string &&geometry_name, GeoLib::GEOObjects &geometry)