1 // Boost.Geometry (aka GGL, Generic Geometry Library)
2 // Unit Test
3
4 // Copyright (c) 2012-2014 Barend Gehrels, Amsterdam, the Netherlands.
5
6 // Use, modification and distribution is subject to the Boost Software License,
7 // Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
8 // http://www.boost.org/LICENSE_1_0.txt)
9
10 #include <test_buffer.hpp>
11
12 #include <boost/geometry/algorithms/buffer.hpp>
13 #include <boost/geometry/core/coordinate_type.hpp>
14
15 #include <boost/geometry/strategies/strategies.hpp>
16
17 #include <boost/geometry/geometries/point.hpp>
18 #include <boost/geometry/geometries/box.hpp>
19 #include <test_common/test_point.hpp>
20
21
22 static std::string const simplex = "LINESTRING(0 0,4 5)";
23 static std::string const simplex_vertical = "LINESTRING(0 0,0 1)";
24 static std::string const simplex_horizontal = "LINESTRING(0 0,1 0)";
25
26 static std::string const straight = "LINESTRING(0 0,4 5,8 10)";
27 static std::string const one_bend = "LINESTRING(0 0,4 5,7 4)";
28 static std::string const two_bends = "LINESTRING(0 0,4 5,7 4,10 6)";
29
30 static std::string const bend_near_start1 = "LINESTRING(0 0,1 0,5 2)";
31 static std::string const bend_near_start2 = "LINESTRING(0 0,1 0,2 1.5,5 3)";
32
33 static std::string const overlapping = "LINESTRING(0 0,4 5,7 4,10 6, 10 2,2 2)";
34 static std::string const curve = "LINESTRING(2 7,3 5,5 4,7 5,8 7)";
35 static std::string const tripod = "LINESTRING(5 0,5 5,1 8,5 5,9 8)"; // with spike
36
37 static std::string const degenerate0 = "LINESTRING()";
38 static std::string const degenerate1 = "LINESTRING(5 5)";
39 static std::string const degenerate2 = "LINESTRING(5 5,5 5)";
40 static std::string const degenerate3 = "LINESTRING(5 5,5 5,5 5)";
41 static std::string const degenerate4 = "LINESTRING(5 5,5 5,4 4,5 5,5 5)";
42
43 static std::string const for_collinear = "LINESTRING(2 0,0 0,0 4,6 4,6 0,4 0)";
44 static std::string const for_collinear2 = "LINESTRING(2 1,2 0,0 0,0 4,6 4,6 0,4 0,4 1)";
45
46 static std::string const chained2 = "LINESTRING(0 0,1 1,2 2)";
47 static std::string const chained3 = "LINESTRING(0 0,1 1,2 2,3 3)";
48 static std::string const chained4 = "LINESTRING(0 0,1 1,2 2,3 3,4 4)";
49
50 static std::string const field_sprayer1 = "LINESTRING(76396.40464822574 410095.6795147947,76397.85016212701 410095.211865792,76401.30666443033 410095.0466387949,76405.05892643372 410096.1007777959,76409.45103273794 410098.257640797,76412.96309264141 410101.6522238015)";
51 static std::string const aimes120 = "LINESTRING(-2.505218 52.189211,-2.505069 52.189019,-2.504941 52.188854)";
52 static std::string const aimes167 = "LINESTRING(-2.378569 52.312133,-2.37857 52.312127,-2.378544 52.31209)";
53 static std::string const aimes175 = "LINESTRING(-2.3116 52.354326,-2.311555 52.35417,-2.311489 52.354145,-2.311335 52.354178)";
54 static std::string const aimes171 = "LINESTRING(-2.393161 52.265087,-2.393002 52.264965,-2.392901 52.264891)";
55 static std::string const aimes181 = "LINESTRING(-2.320686 52.43505,-2.320678 52.435016,-2.320697 52.434978,-2.3207 52.434977,-2.320741 52.434964,-2.320807 52.434964,-2.320847 52.434986,-2.320903 52.435022)";
56
57 static std::string const crossing = "LINESTRING(0 0,10 10,10 0,0 10)";
58
59 // Simplified cases from multi_linestring tesets:
60 static std::string const mikado1 = "LINESTRING(11.406143344709896325639419956133 0.75426621160409546007485914742574,12 1,11.403846153846153299582510953769 0.75)";
61
62 static std::string const mysql_report_2015_03_02a = "LINESTRING(0 0,0 5,5 5,5 0,0 0)"; // closed
63 static std::string const mysql_report_2015_03_02b = "LINESTRING(0 1,0 5,5 5,5 0,1 0)"; // not closed, 1 difference
64 static std::string const mysql_report_2015_03_02c = "LINESTRING(0 2,0 5,5 5,5 0,2 0)"; // not closed, 2 difference
65
66 static std::string const mysql_report_2015_04_01 = "LINESTRING(103 5,107 2,111 4,116 -1,115 0,112 4)";
67
68 static std::string const mysql_report_2015_04_10a = "LINESTRING(1.922421e+307 1.520384e+308, 15 42, 89 -93,-89 -22)";
69 static std::string const mysql_report_2015_04_10b = "LINESTRING(15 42, 89 -93,-89 -22, 1.922421e+307 1.520384e+308)";
70 static std::string const mysql_report_2015_04_10c = "LINESTRING(15 42, 1.922421e+307 1.520384e+308, 89 -93,-89 -22)";
71 static std::string const mysql_report_2015_04_10d = "LINESTRING(1.922421e+307 1.520384e+308, 1.923421e+307 1.521384e+308, 15 42, 89 -93,-89 -22)";
72 static std::string const mysql_report_2015_04_10e = "LINESTRING(15 42, 89 -93,-89 -22, 1.922421e+307 1.520384e+308, 1.923421e+307 1.521384e+308)";
73 static std::string const mysql_report_2015_04_10f = "LINESTRING(15 42, 1.922421e+307 1.520384e+308, 1.923421e+307 1.521384e+308, 89 -93,-89 -22)";
74 static std::string const mysql_report_2015_04_10g = "LINESTRING(15 42, 89 -93,-89 -22)";
75
76 static std::string const mysql_report_2015_06_11 = "LINESTRING("
77 "-155.9300341531310000 4.1672727531600900, -14.0079144546799000 "
78 "-12.2485554508160000, 176.9503531590800000 -3.0930641354495000, "
79 "32.6863251871831000 -17.9691125862157000, -17.7739746299451000 "
80 "41.3177973084700000, -36.0310834162082000 59.9486214620753000, "
81 "153.1574937017440000 46.3007892930418000, 172.7795126069240000 "
82 "19.5367061763707000, -85.6306040220105000 35.0128339347489000, "
83 "-61.1404987988716000 0.3278080608359490, -127.5034592987520000 "
84 "18.6202802642343000, 114.5567005754250000 -83.7227732658958000, "
85 "-66.1134822881378000 -75.2141906159065000, -93.7363999307791000 "
86 "49.3124773443269000, -8.7182702071584100 56.2071174970861000, "
87 "7.7959787229988800 60.8845281744769000, 13.0420633931840000 "
88 "58.8150539662759000, -89.9754374613871000 26.4546501154335000, "
89 "-44.5746548960799000 -88.8122262334508000, -178.4807616092640000 "
90 "10.7770331393820000, 161.8238702890570000 -42.3894892597522000, "
91 "136.2382890452810000 28.6261570633511000, 49.6788041059295000 "
92 "61.7724885566963000, 52.7876201424690000 -61.9246644395984000, "
93 "-162.7456296900030000 11.7183989853218000, 115.6208648232840000 "
94 "51.0941612539320000, -48.7772321835054000 50.4339743128205000)";
95
96 template <bool Clockwise, typename P>
test_all()97 void test_all()
98 {
99 typedef bg::model::linestring<P> linestring;
100 typedef bg::model::polygon<P, Clockwise> polygon;
101
102 bg::strategy::buffer::join_miter join_miter;
103 bg::strategy::buffer::join_round join_round(100);
104 bg::strategy::buffer::join_round_by_divide join_round_by_divide(4);
105 bg::strategy::buffer::end_flat end_flat;
106 bg::strategy::buffer::end_round end_round(100);
107
108 // For testing MySQL issues, which uses 32 by default
109 bg::strategy::buffer::end_round end_round32(32);
110 bg::strategy::buffer::join_round join_round32(32);
111
112 // Simplex (join-type is not relevant)
113 test_one<linestring, polygon>("simplex", simplex, join_miter, end_flat, 19.209, 1.5, 1.5);
114 test_one<linestring, polygon>("simplex", simplex, join_miter, end_round, 26.2733, 1.5, 1.5);
115
116 // Should be about PI + 2
117 test_one<linestring, polygon>("simplex_vertical", simplex_vertical, join_round, end_round, 5.14, 1, 1);
118 test_one<linestring, polygon>("simplex_horizontal", simplex_horizontal, join_round, end_round, 5.14, 1, 1);
119
120 // Should be a bit less than PI + 2
121 test_one<linestring, polygon>("simplex_vertical32", simplex_vertical, join_round32, end_round32, 5.12145, 1, 1);
122 test_one<linestring, polygon>("simplex_horizontal32", simplex_horizontal, join_round32, end_round32, 5.12145, 1, 1);
123
124 test_one<linestring, polygon>("simplex_asym_neg", simplex, join_miter, end_flat, 3.202, +1.5, -1.0);
125 test_one<linestring, polygon>("simplex_asym_pos", simplex, join_miter, end_flat, 3.202, -1.0, +1.5);
126 // Do not work yet:
127 // test_one<linestring, polygon>("simplex_asym_neg", simplex, join_miter, end_round, 3.202, +1.5, -1.0);
128 // test_one<linestring, polygon>("simplex_asym_pos", simplex, join_miter, end_round, 3.202, -1.0, +1.5);
129
130 // Generates (initially) a reversed polygon, with a negative area, which is reversed afterwards in assign_parents
131 test_one<linestring, polygon>("simplex_asym_neg_rev", simplex, join_miter, end_flat, 3.202, +1.0, -1.5);
132 test_one<linestring, polygon>("simplex_asym_pos_rev", simplex, join_miter, end_flat, 3.202, -1.5, +1.0);
133
134 test_one<linestring, polygon>("straight", straight, join_round, end_flat, 38.4187, 1.5, 1.5);
135 test_one<linestring, polygon>("straight", straight, join_miter, end_flat, 38.4187, 1.5, 1.5);
136
137 // One bend/two bends (tests join-type)
138 test_one<linestring, polygon>("one_bend", one_bend, join_round, end_flat, 28.488, 1.5, 1.5);
139 test_one<linestring, polygon>("one_bend", one_bend, join_miter, end_flat, 28.696, 1.5, 1.5);
140 test_one<linestring, polygon>("one_bend", one_bend, join_round_by_divide, end_flat, 28.488, 1.5, 1.5);
141
142 test_one<linestring, polygon>("one_bend", one_bend, join_round, end_round, 35.5603, 1.5, 1.5);
143 test_one<linestring, polygon>("one_bend", one_bend, join_miter, end_round, 35.7601, 1.5, 1.5);
144
145 test_one<linestring, polygon>("two_bends", two_bends, join_round, end_flat, 39.235, 1.5, 1.5);
146 test_one<linestring, polygon>("two_bends", two_bends, join_round_by_divide, end_flat, 39.235, 1.5, 1.5);
147 test_one<linestring, polygon>("two_bends", two_bends, join_miter, end_flat, 39.513, 1.5, 1.5);
148 test_one<linestring, polygon>("two_bends_left", two_bends, join_round, end_flat, 20.028, 1.5, 0.0);
149 test_one<linestring, polygon>("two_bends_left", two_bends, join_miter, end_flat, 20.225, 1.5, 0.0);
150 test_one<linestring, polygon>("two_bends_right", two_bends, join_round, end_flat, 19.211, 0.0, 1.5);
151 test_one<linestring, polygon>("two_bends_right", two_bends, join_miter, end_flat, 19.288, 0.0, 1.5);
152
153 test_one<linestring, polygon>("bend_near_start1", bend_near_start1, join_round, end_flat, 109.2625, 9.0, 9.0);
154 test_one<linestring, polygon>("bend_near_start2", bend_near_start2, join_round, end_flat, 142.8709, 9.0, 9.0);
155
156 // Next (and all similar cases) which a offsetted-one-sided buffer has to be fixed. TODO
157 //test_one<linestring, polygon>("two_bends_neg", two_bends, join_miter, end_flat, 99, +1.5, -1.0);
158 //test_one<linestring, polygon>("two_bends_pos", two_bends, join_miter, end_flat, 99, -1.5, +1.0);
159 //test_one<linestring, polygon>("two_bends_neg", two_bends, join_round, end_flat,99, +1.5, -1.0);
160 //test_one<linestring, polygon>("two_bends_pos", two_bends, join_round, end_flat, 99, -1.5, +1.0);
161
162 test_one<linestring, polygon>("overlapping150", overlapping, join_round, end_flat, 65.6786, 1.5, 1.5);
163 test_one<linestring, polygon>("overlapping150", overlapping, join_miter, end_flat, 68.140, 1.5, 1.5);
164
165 // Different cases with intersection points on flat and (left/right from line itself)
166 test_one<linestring, polygon>("overlapping_asym_150_010", overlapping, join_round, end_flat, 48.308, 1.5, 0.25);
167 test_one<linestring, polygon>("overlapping_asym_150_010", overlapping, join_miter, end_flat, 50.770, 1.5, 0.25);
168 test_one<linestring, polygon>("overlapping_asym_150_075", overlapping, join_round, end_flat, 58.506, 1.5, 0.75);
169 test_one<linestring, polygon>("overlapping_asym_150_075", overlapping, join_miter, end_flat, 60.985, 1.5, 0.75);
170 test_one<linestring, polygon>("overlapping_asym_150_100", overlapping, join_round, end_flat, 62.514, 1.5, 1.0);
171 test_one<linestring, polygon>("overlapping_asym_150_100", overlapping, join_miter, end_flat, 64.984, 1.5, 1.0);
172
173 // Having flat end
174 test_one<linestring, polygon>("for_collinear", for_collinear, join_round, end_flat, 68.561, 2.0, 2.0);
175 test_one<linestring, polygon>("for_collinear", for_collinear, join_miter, end_flat, 72, 2.0, 2.0);
176 #if defined(BOOST_GEOMETRY_BUFFER_INCLUDE_FAILING_TESTS)
177 test_one<linestring, polygon>("for_collinear2", for_collinear2, join_round, end_flat, 74.387, 2.0, 2.0);
178 test_one<linestring, polygon>("for_collinear2", for_collinear2, join_miter, end_flat, 78.0, 2.0, 2.0);
179 #endif
180
181 test_one<linestring, polygon>("curve", curve, join_round, end_flat, 58.1944, 5.0, 3.0);
182 test_one<linestring, polygon>("curve", curve, join_miter, end_flat, 58.7371, 5.0, 3.0);
183
184 test_one<linestring, polygon>("tripod", tripod, join_miter, end_flat, 74.25, 3.0);
185 test_one<linestring, polygon>("tripod", tripod, join_miter, end_round, 116.6336, 3.0);
186
187 test_one<linestring, polygon>("chained2", chained2, join_round, end_flat, 11.3137, 2.5, 1.5);
188 test_one<linestring, polygon>("chained3", chained3, join_round, end_flat, 16.9706, 2.5, 1.5);
189 test_one<linestring, polygon>("chained4", chained4, join_round, end_flat, 22.6274, 2.5, 1.5);
190
191 test_one<linestring, polygon>("field_sprayer1", field_sprayer1, join_round, end_flat, 324.3550, 16.5, 6.5);
192 test_one<linestring, polygon>("field_sprayer1", field_sprayer1, join_round, end_round, 718.761877, 16.5, 6.5);
193 test_one<linestring, polygon>("field_sprayer1", field_sprayer1, join_miter, end_round, 718.939628, 16.5, 6.5);
194
195 test_one<linestring, polygon>("degenerate0", degenerate0, join_round, end_round, 0.0, 3.0);
196 test_one<linestring, polygon>("degenerate1", degenerate1, join_round, end_round, 28.25, 3.0);
197 test_one<linestring, polygon>("degenerate2", degenerate2, join_round, end_round, 28.2503, 3.0);
198 test_one<linestring, polygon>("degenerate3", degenerate3, join_round, end_round, 28.2503, 3.0);
199 test_one<linestring, polygon>("degenerate4", degenerate4, join_round, end_round, 36.7410, 3.0);
200 test_one<linestring, polygon>("degenerate4", degenerate4, join_round, end_flat, 8.4853, 3.0);
201
202 {
203 // These tests do test behaviour in end_round strategy:
204 // -> it should generate closed pieces, also for an odd number of points.
205 // It also tests behaviour in join_round strategy:
206 // -> it should generate e.g. 4 points for a full circle,
207 // so a quarter circle does not get points in between
208 using bg::strategy::buffer::join_round;
209 using bg::strategy::buffer::end_round;
210
211 double const d10 = 1.0;
212
213 test_one<linestring, polygon>("mysql_report_2015_03_02a_3", mysql_report_2015_03_02a, join_round(3), end_round(3), 38.000, d10);
214 test_one<linestring, polygon>("mysql_report_2015_03_02a_4", mysql_report_2015_03_02a, join_round(4), end_round(4), 38.000, d10);
215 test_one<linestring, polygon>("mysql_report_2015_03_02a_5", mysql_report_2015_03_02a, join_round(5), end_round(5), 38.790, d10);
216 test_one<linestring, polygon>("mysql_report_2015_03_02a_6", mysql_report_2015_03_02a, join_round(6), end_round(6), 38.817, d10);
217 test_one<linestring, polygon>("mysql_report_2015_03_02a_7", mysql_report_2015_03_02a, join_round(7), end_round(7), 38.851, d10);
218 test_one<linestring, polygon>("mysql_report_2015_03_02b_3", mysql_report_2015_03_02b, join_round(3), end_round(3), 36.500, d10);
219 test_one<linestring, polygon>("mysql_report_2015_03_02b_4", mysql_report_2015_03_02b, join_round(4), end_round(4), 36.500, d10);
220 test_one<linestring, polygon>("mysql_report_2015_03_02b_5", mysql_report_2015_03_02b, join_round(5), end_round(5), 37.346, d10);
221 test_one<linestring, polygon>("mysql_report_2015_03_02b_6", mysql_report_2015_03_02b, join_round(6), end_round(6), 37.402, d10);
222 test_one<linestring, polygon>("mysql_report_2015_03_02b_7", mysql_report_2015_03_02b, join_round(7), end_round(7), 37.506, d10);
223 test_one<linestring, polygon>("mysql_report_2015_03_02c_3", mysql_report_2015_03_02c, join_round(2), end_round(3), 32.500, d10);
224 test_one<linestring, polygon>("mysql_report_2015_03_02c_4", mysql_report_2015_03_02c, join_round(4), end_round(4), 32.500, d10);
225 test_one<linestring, polygon>("mysql_report_2015_03_02c_5", mysql_report_2015_03_02c, join_round(5), end_round(5), 33.611, d10);
226 test_one<linestring, polygon>("mysql_report_2015_03_02c_6", mysql_report_2015_03_02c, join_round(6), end_round(6), 33.719, d10);
227 test_one<linestring, polygon>("mysql_report_2015_03_02c_7", mysql_report_2015_03_02c, join_round(7), end_round(7), 33.901, d10);
228
229 // Testing the asymmetric end caps with odd number of points
230 double const d15 = 1.5;
231 test_one<linestring, polygon>("mysql_report_2015_03_02c_asym1", mysql_report_2015_03_02c, join_round(7), end_round(7), 39.714, d10, d15);
232 test_one<linestring, polygon>("mysql_report_2015_03_02c_asym2", mysql_report_2015_03_02c, join_round(7), end_round(7), 46.116, d15, d10);
233
234 #if defined(BOOST_GEOMETRY_BUFFER_USE_SIDE_OF_INTERSECTION)
235 double const d100 = 10;
236 test_one<linestring, polygon>("mysql_report_2015_04_01", mysql_report_2015_04_01, join_round(32), end_round(32), 632.234, d100);
237 #endif
238 }
239
240
241 {
242 double tolerance = 1.0e-10;
243
244 test_one<linestring, polygon>("aimes120", aimes120, join_miter, end_flat, 1.62669948622351512e-08, 0.000018, 0.000018, false, tolerance);
245 test_one<linestring, polygon>("aimes120", aimes120, join_round, end_round, 1.72842078427493107e-08, 0.000018, 0.000018, true, tolerance);
246
247 test_one<linestring, polygon>("aimes167", aimes167, join_miter, end_flat, 1.88900628472765675e-09, 0.000018, 0.000018, true, tolerance);
248 test_one<linestring, polygon>("aimes167", aimes167, join_round, end_round, 2.85734813587623648e-09, 0.000018, 0.000018, true, tolerance);
249
250 test_one<linestring, polygon>("aimes175", aimes175, join_miter, end_flat, 2.81111809385947709e-08, 0.000036, 0.000036, true, tolerance);
251 test_one<linestring, polygon>("aimes175", aimes175, join_round, end_round, 3.21215765097804251e-08, 0.000036, 0.000036, true, tolerance);
252
253 test_one<linestring, polygon>("aimes171", aimes171, join_miter, end_flat, 1.1721873249825876e-08, 0.000018, 0.000018, true, tolerance);
254 test_one<linestring, polygon>("aimes171", aimes171, join_round, end_round, 1.2739093335767393e-08, 0.000018, 0.000018, true, tolerance);
255 test_one<linestring, polygon>("aimes171", aimes171, join_round_by_divide, end_round, 1.2739093335767393e-08, 0.000018, 0.000018, true, tolerance);
256
257 test_one<linestring, polygon>("aimes181", aimes181, join_miter, end_flat, 2.1729405830228643e-08, 0.000036, 0.000036, true, tolerance);
258 test_one<linestring, polygon>("aimes181", aimes181, join_round, end_round, 2.57415564419716247e-08, 0.000036, 0.000036, true, tolerance);
259 test_one<linestring, polygon>("aimes181", aimes181, join_round_by_divide, end_round, 2.57415564419716247e-08, 0.000036, 0.000036, true, tolerance);
260 }
261
262 // On one compiler 1702.56530051454502 2140.78725663358819
263 // so we increase tolerance
264 test_one<linestring, polygon>("crossing", crossing, join_round32, end_flat, 1702.1, 20.0, 20.0, true, 0.5);
265 test_one<linestring, polygon>("crossing", crossing, join_round32, end_round32, 2140.4, 20.0, 20.0, true, 0.5);
266
267 test_one<linestring, polygon>("mikado1", mikado1, join_round32, end_round32, 5441135039.0979, 41751.0);
268
269 test_one<linestring, polygon>("mysql_report_2015_06_11",
270 mysql_report_2015_06_11, join_round32, end_round32,
271 27862.733459829971,
272 5.9518403867035365);
273 }
274
275 template <bool Clockwise, typename P>
test_invalid()276 void test_invalid()
277 {
278 typedef bg::model::linestring<P> linestring;
279 typedef bg::model::polygon<P, Clockwise> polygon;
280
281 bg::strategy::buffer::end_round end_round32(32);
282 bg::strategy::buffer::join_round join_round32(32);
283
284 // Linestring contains extreme differences causing numeric errors. Empty geometry is returned
285 test_one<linestring, polygon>("mysql_report_2015_04_10a", mysql_report_2015_04_10a, join_round32, end_round32, 0.0, 100.0);
286 test_one<linestring, polygon>("mysql_report_2015_04_10b", mysql_report_2015_04_10b, join_round32, end_round32, 0.0, 100.0);
287 test_one<linestring, polygon>("mysql_report_2015_04_10c", mysql_report_2015_04_10c, join_round32, end_round32, 0.0, 100.0);
288 test_one<linestring, polygon>("mysql_report_2015_04_10d", mysql_report_2015_04_10d, join_round32, end_round32, 0.0, 100.0);
289 test_one<linestring, polygon>("mysql_report_2015_04_10e", mysql_report_2015_04_10e, join_round32, end_round32, 0.0, 100.0);
290 test_one<linestring, polygon>("mysql_report_2015_04_10f", mysql_report_2015_04_10f, join_round32, end_round32, 0.0, 100.0);
291
292 // The equivalent, valid, case
293 test_one<linestring, polygon>("mysql_report_2015_04_10g", mysql_report_2015_04_10g, join_round32, end_round32, 86527.871, 100.0);
294 }
295
296 #ifdef HAVE_TTMATH
297 #include <ttmath_stub.hpp>
298 #endif
299
300
test_main(int,char * [])301 int test_main(int, char* [])
302 {
303 test_all<true, bg::model::point<double, 2, bg::cs::cartesian> >();
304 test_all<false, bg::model::point<double, 2, bg::cs::cartesian> >();
305 //test_all<bg::model::point<tt, 2, bg::cs::cartesian> >();
306
307 test_invalid<true, bg::model::point<double, 2, bg::cs::cartesian> >();
308 // test_invalid<true, bg::model::point<long double, 2, bg::cs::cartesian> >();
309 return 0;
310 }
311