1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
|
/* Copyright (c) 2017, 2025, Oracle and/or its affiliates.
This program is free software; you can redistribute it and/or modify
it under the terms of the GNU General Public License, version 2.0,
as published by the Free Software Foundation.
This program is designed to work with certain software (including
but not limited to OpenSSL) that is licensed under separate terms,
as designated in a particular file or component or in included license
documentation. The authors of MySQL hereby grant you an additional
permission to link the program and your derivative works with the
separately licensed software that they have either included with
the program or referenced in the documentation.
This program is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
GNU General Public License, version 2.0, for more details.
You should have received a copy of the GNU General Public License
along with this program; if not, write to the Free Software
Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA */
#ifndef UNITTEST_GUNIT_GIS_TEST_GEOM_H_INCLUDED
#define UNITTEST_GUNIT_GIS_TEST_GEOM_H_INCLUDED
#include <cmath> // M_PI, M_PI_2
#include "sql/gis/geometries_cs.h" // gis::{Cartesian,Geographic}*
// Header-only definitions of sample geometries useful for testing
// Note that the internal point representation for geographic coordinates is
// (eastwards (-pi, pi], northwards [pi/2, pi/2]) in radians.
namespace {
namespace test_geom {
using Cpt = gis::Cartesian_point;
using Cls = gis::Cartesian_linestring;
using Cgc = gis::Cartesian_geometrycollection;
using Gpt = gis::Geographic_point;
using Gls = gis::Geographic_linestring;
using Ggc = gis::Geographic_geometrycollection;
inline Cpt point_origin() { return {0.0, 0.0}; }
inline Cls linestring_line() {
// Slanted line
Cls ls;
ls.push_back(Cpt{-1, -1});
ls.push_back(Cpt{+1, +1});
return ls;
}
inline Cls linestring_right_angle() {
// _|-shape
Cls ls;
ls.push_back(Cpt{-1, -1});
ls.push_back(Cpt{+1, -1});
ls.push_back(Cpt{+1, +1});
return ls;
}
inline Cls linestring_triangle() {
// Closed _|-shape
Cls ls;
ls.push_back(Cpt{-1, -1});
ls.push_back(Cpt{+1, -1});
ls.push_back(Cpt{+1, +1});
ls.push_back(Cpt{-1, -1});
return ls;
}
inline Cls linestring_selfintersecting() {
// Open bow / fish-shape, like: ><|
Cls ls;
ls.push_back(Cpt{-1, +1});
ls.push_back(Cpt{+1, -1});
ls.push_back(Cpt{+1, +1});
ls.push_back(Cpt{-1, -1});
return ls;
}
inline Cls linestring_selfintersecting_closed() {
// Bow-shape, like |><|
Cls ls;
ls.push_back(Cpt{-1, +1});
ls.push_back(Cpt{+1, -1});
ls.push_back(Cpt{+1, +1});
ls.push_back(Cpt{-1, -1});
ls.push_back(Cpt{-1, +1});
return ls;
}
inline Cls linestring_selftangential() {
Cls ls;
ls.push_back(Cpt{-1.0, -0.75});
ls.push_back(Cpt{+1.0, -0.75});
ls.push_back(Cpt{+0.0, +0.00}); // Tangent point
ls.push_back(Cpt{+0.5, -0.50});
ls.push_back(Cpt{-0.5, -0.50}); // Tangent line
ls.push_back(Cpt{+0.5, +0.50}); // Touch at midpoint
ls.push_back(Cpt{-1.0, -0.75});
return ls;
}
inline Gls linestring_pole_alias() {
// Line between two coordinates both mapping to the north-pole
Gls ls;
ls.push_back(Gpt{0, M_PI_2});
ls.push_back(Gpt{1, M_PI_2});
return ls;
}
inline Cgc geometrycollection_empty() { return {}; }
inline Cgc geometrycollection_crossintersecting() {
Cgc gc;
{
Cls ls;
ls.push_back(Cpt{-1, -1});
ls.push_back(Cpt{+1, +1});
gc.push_back(ls);
}
{
Cls ls;
ls.push_back(Cpt{-1, +1});
ls.push_back(Cpt{+1, -1});
gc.push_back(ls);
}
return gc;
}
inline Cgc geometrycollection_subintersecting() {
Cgc gc;
gc.push_back(linestring_selfintersecting());
return gc;
}
inline Ggc bug26495068() {
// Regression test for Bug#26495068
// Certain geometries crossing the anti-meridian give incorrect results.
auto deg = [](double x) { return x * M_PI / 180; };
Ggc gc;
{
Gls ls;
ls.push_back(Gpt{deg(135), deg(0)});
ls.push_back(Gpt{deg(-150), deg(36)});
gc.push_back(ls);
}
{
Gls ls;
ls.push_back(Gpt{deg(-101), deg(0)});
ls.push_back(Gpt{deg(-178), deg(30)});
gc.push_back(ls);
}
return gc;
}
} // namespace test_geom
} // namespace
#endif // include guard
|
