ceph/src/boost/libs/geometry/test/algorithms/buffer/buffer_linestring.cpp

   1 // Boost.Geometry (aka GGL, Generic Geometry Library)
   2 // Unit Test
   3
   4 // Copyright (c) 2012-2019 Barend Gehrels, Amsterdam, the Netherlands.
   5
   6 // This file was modified by Oracle on 2016-2021.
   7 // Modifications copyright (c) 2016-2021, Oracle and/or its affiliates.
   8 // Contributed and/or modified by Adam Wulkiewicz, on behalf of Oracle
   9
  10 // Use, modification and distribution is subject to the Boost Software License,
  11 // Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
  12 // http://www.boost.org/LICENSE_1_0.txt)
  13
  14 #include "test_buffer.hpp"
  15
  16 #include <boost/geometry/algorithms/buffer.hpp>
  17 #include <boost/geometry/core/coordinate_type.hpp>
  18
  19 #include <boost/geometry/strategies/strategies.hpp>
  20
  21 #include <boost/geometry/geometries/point.hpp>
  22 #include <boost/geometry/geometries/box.hpp>
  23 #include "test_common/test_point.hpp"
  24
  25
  26 static std::string const simplex = "LINESTRING(0 0,4 5)";
  27 static std::string const simplex_vertical = "LINESTRING(0 0,0 1)";
  28 static std::string const simplex_horizontal = "LINESTRING(0 0,1 0)";
  29
  30 static std::string const straight = "LINESTRING(0 0,4 5,8 10)";
  31 static std::string const one_bend = "LINESTRING(0 0,4 5,7 4)";
  32 static std::string const two_bends = "LINESTRING(0 0,4 5,7 4,10 6)";
  33
  34 static std::string const bend_near_start1 = "LINESTRING(0 0,1 0,5 2)";
  35 static std::string const bend_near_start2 = "LINESTRING(0 0,1 0,2 1.5,5 3)";
  36
  37 static std::string const overlapping = "LINESTRING(0 0,4 5,7 4,10 6, 10 2,2 2)";
  38 static std::string const curve = "LINESTRING(2 7,3 5,5 4,7 5,8 7)";
  39 static std::string const tripod = "LINESTRING(5 0,5 5,1 8,5 5,9 8)"; // with spike
  40
  41 static std::string const degenerate0 = "LINESTRING()";
  42 static std::string const degenerate1 = "LINESTRING(5 5)";
  43 static std::string const degenerate2 = "LINESTRING(5 5,5 5)";
  44 static std::string const degenerate3 = "LINESTRING(5 5,5 5,5 5)";
  45 static std::string const degenerate4 = "LINESTRING(5 5,5 5,4 4,5 5,5 5)";
  46
  47 static std::string const for_collinear = "LINESTRING(2 0,0 0,0 4,6 4,6 0,4 0)";
  48 static std::string const for_collinear2 = "LINESTRING(2 1,2 0,0 0,0 4,6 4,6 0,4 0,4 1)";
  49
  50 static std::string const chained2 = "LINESTRING(0 0,1 1,2 2)";
  51 static std::string const chained3 = "LINESTRING(0 0,1 1,2 2,3 3)";
  52 static std::string const chained4 = "LINESTRING(0 0,1 1,2 2,3 3,4 4)";
  53
  54 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)";
  55 static std::string const aimes120 = "LINESTRING(-2.505218 52.189211,-2.505069 52.189019,-2.504941 52.188854)";
  56 static std::string const aimes167 = "LINESTRING(-2.378569 52.312133,-2.37857 52.312127,-2.378544 52.31209)";
  57 static std::string const aimes175 = "LINESTRING(-2.3116 52.354326,-2.311555 52.35417,-2.311489 52.354145,-2.311335 52.354178)";
  58 static std::string const aimes171 = "LINESTRING(-2.393161 52.265087,-2.393002 52.264965,-2.392901 52.264891)";
  59 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)";
  60
  61 static std::string const crossing = "LINESTRING(0 0,10 10,10 0,0 10)";
  62
  63 // Simplified cases from multi_linestring tesets:
  64 static std::string const mikado1 = "LINESTRING(11.406143344709896325639419956133 0.75426621160409546007485914742574,12 1,11.403846153846153299582510953769 0.75)";
  65
  66 static std::string const mysql_report_2015_03_02a = "LINESTRING(0 0,0 5,5 5,5 0,0 0)"; // closed
  67 static std::string const mysql_report_2015_03_02b = "LINESTRING(0 1,0 5,5 5,5 0,1 0)"; // not closed, 1 difference
  68 static std::string const mysql_report_2015_03_02c = "LINESTRING(0 2,0 5,5 5,5 0,2 0)"; // not closed, 2 difference
  69
  70 static std::string const mysql_report_2015_04_01 = "LINESTRING(103 5,107 2,111 4,116 -1,115 0,112 4)";
  71
  72 static std::string const mysql_report_2015_04_10a = "LINESTRING(1.922421e+307 1.520384e+308, 15 42, 89 -93,-89 -22)";
  73 static std::string const mysql_report_2015_04_10b = "LINESTRING(15 42, 89 -93,-89 -22, 1.922421e+307 1.520384e+308)";
  74 static std::string const mysql_report_2015_04_10c = "LINESTRING(15 42, 1.922421e+307 1.520384e+308, 89 -93,-89 -22)";
  75 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)";
  76 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)";
  77 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)";
  78 static std::string const mysql_report_2015_04_10g = "LINESTRING(15 42, 89 -93,-89 -22)";
  79
  80 static std::string const mysql_report_2015_06_11 = "LINESTRING("
  81     "-155.9300341531310000 4.1672727531600900, -14.0079144546799000 "
  82     "-12.2485554508160000, 176.9503531590800000 -3.0930641354495000, "
  83     "32.6863251871831000 -17.9691125862157000, -17.7739746299451000 "
  84     "41.3177973084700000, -36.0310834162082000 59.9486214620753000, "
  85     "153.1574937017440000 46.3007892930418000, 172.7795126069240000 "
  86     "19.5367061763707000, -85.6306040220105000 35.0128339347489000, "
  87     "-61.1404987988716000 0.3278080608359490, -127.5034592987520000 "
  88     "18.6202802642343000, 114.5567005754250000 -83.7227732658958000, "
  89     "-66.1134822881378000 -75.2141906159065000, -93.7363999307791000 "
  90     "49.3124773443269000, -8.7182702071584100 56.2071174970861000, "
  91     "7.7959787229988800 60.8845281744769000, 13.0420633931840000 "
  92     "58.8150539662759000, -89.9754374613871000 26.4546501154335000, "
  93     "-44.5746548960799000 -88.8122262334508000, -178.4807616092640000 "
  94     "10.7770331393820000, 161.8238702890570000 -42.3894892597522000, "
  95     "136.2382890452810000 28.6261570633511000, 49.6788041059295000 "
  96     "61.7724885566963000, 52.7876201424690000 -61.9246644395984000, "
  97     "-162.7456296900030000 11.7183989853218000, 115.6208648232840000 "
  98     "51.0941612539320000, -48.7772321835054000 50.4339743128205000)";
  99
 100 static std::string const mysql_report_2015_09_08a = "LINESTRING(1 1, 2 1, 1.765258e+308 4, -1 1, 10 4)";
 101 static std::string const mysql_report_2015_09_08b = "LINESTRING(2199023255556 16777218, 32770 8194, 1.417733e+308 7.823620e+307, -8 -9, 2147483649 20)";
 102 static std::string const mysql_report_2015_09_08c = "LINESTRING(-5 -8, 2 8, 2.160023e+307 1.937208e+307, -4 -3, -5 -4, 8796093022208 281474976710653)";
 103
 104 static std::string const mysql_23023665 = "LINESTRING(0 0, 0 5, 5 5, 5 0, 0 0)";
 105
 106 static std::string const mysql_25662426 = "LINESTRING(170 4756, 168 4756, 168 4759, 168 4764, 171 4764, 171 4700)";
 107 static std::string const mysql_25662426a = "LINESTRING(170 4756, 168 4756, 168 4759, 168 4764, 171 4764, 171 4750)";
 108
 109 static std::string const issue_596 = "LINESTRING(292979.660 6688731.370, 292979.600 6688733.420, 292979.540 6688735.440, 292979.540 6688735.500)";
 110 static std::string const issue_803 = "LINESTRING(2773.6899360413681 -17.493356405074767,2767.9084041267429 -9.0288291122330797,2765.7403296587586 -1.6929054585683376,2765.3273630934282 6.6613381477509392e-16,2762.8495637014462 10.157432751410026,2759.9587977441333 24.829280059133669,2756.3453402974928 38.372523727680374,2752.0091913615238 45.708447381345117)";
 111
 112 template <bool Clockwise, typename P>
 113 void test_all()
 114 {
 115     typedef bg::model::linestring<P> linestring;
 116     typedef bg::model::polygon<P, Clockwise> polygon;
 117     typedef typename bg::coordinate_type<P>::type coor_type;
 118
 119     bg::strategy::buffer::join_miter join_miter;
 120     bg::strategy::buffer::join_round join_round(100);
 121     bg::strategy::buffer::join_round_by_divide join_round_by_divide(4);
 122     bg::strategy::buffer::end_flat end_flat;
 123     bg::strategy::buffer::end_round end_round(100);
 124
 125     // For testing MySQL issues, which uses 32 by default
 126     bg::strategy::buffer::end_round end_round32(32);
 127     bg::strategy::buffer::join_round join_round32(32);
 128
 129     ut_settings const settings;
 130     ut_settings const specific_settings
 131             = BOOST_GEOMETRY_CONDITION((std::is_same<coor_type, long double>::value))
 132               ? ut_settings(0.02) : settings;
 133
 134     // Simplex (join-type is not relevant)
 135     test_one<linestring, polygon>("simplex", simplex, join_miter, end_flat, 19.209, 1.5);
 136     test_one<linestring, polygon>("simplex", simplex, join_miter, end_round, 26.2733, 1.5);
 137
 138     // Should be about PI + 2
 139     test_one<linestring, polygon>("simplex_vertical", simplex_vertical, join_round, end_round, 5.14, 1);
 140     test_one<linestring, polygon>("simplex_horizontal", simplex_horizontal, join_round, end_round, 5.14, 1);
 141
 142     // Should be a bit less than PI + 2
 143     test_one<linestring, polygon>("simplex_vertical32", simplex_vertical, join_round32, end_round32, 5.12145, 1);
 144     test_one<linestring, polygon>("simplex_horizontal32", simplex_horizontal, join_round32, end_round32, 5.12145, 1);
 145
 146     test_one<linestring, polygon>("simplex_asym_neg", simplex, join_miter, end_flat, 3.2016, +1.5, settings, -1.0);
 147     test_one<linestring, polygon>("simplex_asym_pos", simplex, join_miter, end_flat, 3.2016, -1.0, settings, +1.5);
 148     // Do not work yet:
 149     //    test_one<linestring, polygon>("simplex_asym_neg", simplex, join_miter, end_round, 3.202, +1.5, settings, -1.0);
 150     //    test_one<linestring, polygon>("simplex_asym_pos", simplex, join_miter, end_round, 3.202, -1.0, settings, +1.5);
 151
 152     // Generates (initially) a reversed polygon, with a negative area, which is reversed afterwards in assign_parents
 153     test_one<linestring, polygon>("simplex_asym_neg_rev", simplex, join_miter, end_flat, 3.2016, +1.0, settings, -1.5);
 154     test_one<linestring, polygon>("simplex_asym_pos_rev", simplex, join_miter, end_flat, 3.2016, -1.5, settings, +1.0);
 155
 156     test_one<linestring, polygon>("straight", straight, join_round, end_flat, 38.4187, 1.5);
 157     test_one<linestring, polygon>("straight", straight, join_miter, end_flat, 38.4187, 1.5);
 158
 159     // One bend/two bends (tests join-type)
 160     test_one<linestring, polygon>("one_bend", one_bend, join_round, end_flat, 28.496, 1.5);
 161     test_one<linestring, polygon>("one_bend", one_bend, join_miter, end_flat, 28.696, 1.5);
 162     test_one<linestring, polygon>("one_bend", one_bend, join_round_by_divide, end_flat, 28.497, 1.5);
 163
 164     test_one<linestring, polygon>("one_bend", one_bend, join_round, end_round, 35.5603, 1.5);
 165     test_one<linestring, polygon>("one_bend", one_bend, join_miter, end_round, 35.7601, 1.5);
 166
 167     test_one<linestring, polygon>("two_bends", two_bends, join_round, end_round, 46.2995, 1.5);
 168     test_one<linestring, polygon>("two_bends", two_bends, join_round, end_flat, 39.235, 1.5);
 169     test_one<linestring, polygon>("two_bends", two_bends, join_round_by_divide, end_flat, 39.235, 1.5);
 170     test_one<linestring, polygon>("two_bends", two_bends, join_miter, end_flat, 39.513, 1.5);
 171     test_one<linestring, polygon>("two_bends_left", two_bends, join_round, end_flat, 20.025, 1.5, settings, 0.0);
 172     test_one<linestring, polygon>("two_bends_left", two_bends, join_miter, end_flat, 20.225, 1.5, settings, 0.0);
 173     test_one<linestring, polygon>("two_bends_right", two_bends, join_round, end_flat, 19.211, 0.0, settings, 1.5);
 174     test_one<linestring, polygon>("two_bends_right", two_bends, join_miter, end_flat, 19.288, 0.0, settings, 1.5);
 175
 176     test_one<linestring, polygon>("bend_near_start1", bend_near_start1, join_round, end_flat, 109.2625, 9.0);
 177     test_one<linestring, polygon>("bend_near_start2", bend_near_start2, join_round, end_flat, 142.8709, 9.0);
 178
 179     // Next (and all similar cases) which a offsetted-one-sided buffer has to be fixed. TODO
 180     //test_one<linestring, polygon>("two_bends_neg", two_bends, join_miter, end_flat, 99, +1.5, settings, -1.0);
 181     //test_one<linestring, polygon>("two_bends_pos", two_bends, join_miter, end_flat, 99, -1.5, settings, +1.0);
 182     //test_one<linestring,  polygon>("two_bends_neg", two_bends, join_round, end_flat,99, +1.5, settings, -1.0);
 183     //test_one<linestring, polygon>("two_bends_pos", two_bends, join_round, end_flat, 99, -1.5, settings, +1.0);
 184
 185     test_one<linestring, polygon>("overlapping150", overlapping, join_round, end_flat, 65.6786, 1.5);
 186     test_one<linestring, polygon>("overlapping150", overlapping, join_miter, end_flat, 68.140, 1.5);
 187
 188     // Different cases with intersection points on flat and (left/right from line itself)
 189     test_one<linestring, polygon>("overlapping_asym_150_010", overlapping, join_round, end_flat, 48.308, 1.5, settings, 0.25);
 190     test_one<linestring, polygon>("overlapping_asym_150_010", overlapping, join_miter, end_flat, 50.770, 1.5, settings, 0.25);
 191     test_one<linestring, polygon>("overlapping_asym_150_075", overlapping, join_round, end_flat, 58.506, 1.5, settings, 0.75);
 192     test_one<linestring, polygon>("overlapping_asym_150_075", overlapping, join_miter, end_flat, 60.985, 1.5, settings, 0.75);
 193     test_one<linestring, polygon>("overlapping_asym_150_100", overlapping, join_round, end_flat, 62.514, 1.5, settings, 1.0);
 194     test_one<linestring, polygon>("overlapping_asym_150_100", overlapping, join_miter, end_flat, 64.984, 1.5, settings, 1.0);
 195
 196     // Having flat end
 197     test_one<linestring, polygon>("for_collinear", for_collinear, join_round, end_flat, 68.561, 2.0);
 198     test_one<linestring, polygon>("for_collinear", for_collinear, join_miter, end_flat, 72, 2.0);
 199 #if defined(BOOST_GEOMETRY_TEST_FAILURES)
 200     test_one<linestring, polygon>("for_collinear2", for_collinear2, join_round, end_flat, 74.387, 2.0);
 201     test_one<linestring, polygon>("for_collinear2", for_collinear2, join_miter, end_flat, 78.0, 2.0);
 202 #endif
 203
 204     test_one<linestring, polygon>("curve", curve, join_round, end_flat, 58.1944, 5.0, settings, 3.0);
 205     test_one<linestring, polygon>("curve", curve, join_miter, end_flat, 58.7371, 5.0, settings, 3.0);
 206
 207     test_one<linestring, polygon>("tripod", tripod, join_miter, end_flat, 74.25, 3.0);
 208     test_one<linestring, polygon>("tripod", tripod, join_miter, end_round, 116.6336, 3.0);
 209
 210     test_one<linestring, polygon>("chained2", chained2, join_round, end_flat, 11.3137, 2.5, settings, 1.5);
 211     test_one<linestring, polygon>("chained3", chained3, join_round, end_flat, 16.9706, 2.5, settings, 1.5);
 212     test_one<linestring, polygon>("chained4", chained4, join_round, end_flat, 22.6274, 2.5, settings, 1.5);
 213
 214     test_one<linestring, polygon>("field_sprayer1", field_sprayer1, join_round, end_flat, 324.3550, 16.5, settings, 6.5);
 215     test_one<linestring, polygon>("field_sprayer1", field_sprayer1, join_round, end_round, {718.686, 718.762}, 16.5, settings, 6.5);
 216     test_one<linestring, polygon>("field_sprayer1", field_sprayer1, join_miter, end_round, {718.845, 718.940}, 16.5, settings, 6.5);
 217
 218     test_one<linestring, polygon>("degenerate0", degenerate0, join_round, end_round, 0.0, 3.0);
 219     test_one<linestring, polygon>("degenerate1", degenerate1, join_round, end_round, 28.25, 3.0);
 220     test_one<linestring, polygon>("degenerate2", degenerate2, join_round, end_round, 28.2503, 3.0);
 221     test_one<linestring, polygon>("degenerate3", degenerate3, join_round, end_round, 28.2503, 3.0);
 222     test_one<linestring, polygon>("degenerate4", degenerate4, join_round, end_round, 36.7410, 3.0);
 223     test_one<linestring, polygon>("degenerate4", degenerate4, join_round, end_flat, 8.4853, 3.0);
 224
 225     {
 226         // These tests do test behaviour in end_round strategy:
 227         // -> it should generate closed pieces, also for an odd number of points.
 228         // It also tests behaviour in join_round strategy:
 229         // -> it should generate e.g. 4 points for a full circle,
 230         //    so a quarter circle does not get points in between
 231         using bg::strategy::buffer::join_round;
 232         using bg::strategy::buffer::end_round;
 233
 234         double const d10 = 1.0;
 235
 236         test_one<linestring, polygon>("mysql_report_2015_03_02a_3", mysql_report_2015_03_02a, join_round(3), end_round(3), 38.000, d10);
 237         test_one<linestring, polygon>("mysql_report_2015_03_02a_4", mysql_report_2015_03_02a, join_round(4), end_round(4), 38.000, d10);
 238         test_one<linestring, polygon>("mysql_report_2015_03_02a_5", mysql_report_2015_03_02a, join_round(5), end_round(5), 38.790, d10);
 239         test_one<linestring, polygon>("mysql_report_2015_03_02a_6", mysql_report_2015_03_02a, join_round(6), end_round(6), 38.817, d10);
 240         test_one<linestring, polygon>("mysql_report_2015_03_02a_7", mysql_report_2015_03_02a, join_round(7), end_round(7), 38.851, d10);
 241         test_one<linestring, polygon>("mysql_report_2015_03_02b_3", mysql_report_2015_03_02b, join_round(3), end_round(3), 36.500, d10);
 242         test_one<linestring, polygon>("mysql_report_2015_03_02b_4", mysql_report_2015_03_02b, join_round(4), end_round(4), 36.500, d10);
 243         test_one<linestring, polygon>("mysql_report_2015_03_02b_5", mysql_report_2015_03_02b, join_round(5), end_round(5), 37.346, d10);
 244         test_one<linestring, polygon>("mysql_report_2015_03_02b_6", mysql_report_2015_03_02b, join_round(6), end_round(6), 37.402, d10);
 245         test_one<linestring, polygon>("mysql_report_2015_03_02b_7", mysql_report_2015_03_02b, join_round(7), end_round(7), 37.506, d10);
 246         test_one<linestring, polygon>("mysql_report_2015_03_02c_3", mysql_report_2015_03_02c, join_round(2), end_round(3), 32.500, d10);
 247         test_one<linestring, polygon>("mysql_report_2015_03_02c_4", mysql_report_2015_03_02c, join_round(4), end_round(4), 32.500, d10);
 248         test_one<linestring, polygon>("mysql_report_2015_03_02c_5", mysql_report_2015_03_02c, join_round(5), end_round(5), 33.611, d10);
 249         test_one<linestring, polygon>("mysql_report_2015_03_02c_6", mysql_report_2015_03_02c, join_round(6), end_round(6), 33.719, d10);
 250         test_one<linestring, polygon>("mysql_report_2015_03_02c_7", mysql_report_2015_03_02c, join_round(7), end_round(7), 33.901, d10);
 251
 252         // Testing the asymmetric end caps with odd number of points
 253         double const d15 = 1.5;
 254         test_one<linestring, polygon>("mysql_report_2015_03_02c_asym1", mysql_report_2015_03_02c, join_round(7), end_round(7), 39.714, d10, settings, d15);
 255         test_one<linestring, polygon>("mysql_report_2015_03_02c_asym2", mysql_report_2015_03_02c, join_round(7), end_round(7), 46.116, d15, settings, d10);
 256
 257         double const d100 = 10;
 258         test_one<linestring, polygon>("mysql_report_2015_04_01", mysql_report_2015_04_01, join_round(32), end_round(32), 632.234, d100);
 259     }
 260
 261     if (! BOOST_GEOMETRY_CONDITION((std::is_same<coor_type, float>::value)))
 262     {
 263         ut_settings settings;
 264         settings.tolerance = 0.1;
 265         settings.use_ln_area = true;
 266
 267         test_one<linestring, polygon>("aimes120", aimes120, join_miter, end_flat, 1.62669948622351512e-08, 0.000018, settings);
 268         test_one<linestring, polygon>("aimes120", aimes120, join_round, end_round, 1.72842078427493107e-08, 0.000018, settings);
 269
 270         test_one<linestring, polygon>("aimes167", aimes167, join_miter, end_flat, 1.88900628472765675e-09, 0.000018, settings);
 271         test_one<linestring, polygon>("aimes167", aimes167, join_round, end_round, 2.85734813587623648e-09, 0.000018, settings);
 272
 273         test_one<linestring, polygon>("aimes175", aimes175, join_miter, end_flat, 2.81111809385947709e-08, 0.000036, settings);
 274         test_one<linestring, polygon>("aimes175", aimes175, join_round, end_round, 3.21215765097804251e-08, 0.000036, settings);
 275
 276         test_one<linestring, polygon>("aimes171", aimes171, join_miter, end_flat, 1.1721873249825876e-08, 0.000018, settings);
 277         test_one<linestring, polygon>("aimes171", aimes171, join_round, end_round, 1.2739093335767393e-08, 0.000018, settings);
 278         test_one<linestring, polygon>("aimes171", aimes171, join_round_by_divide, end_round, 1.2739093335767393e-08, 0.000018, settings);
 279
 280         test_one<linestring, polygon>("aimes181", aimes181, join_miter, end_flat, 2.1729405830228643e-08, 0.000036, settings);
 281         test_one<linestring, polygon>("aimes181", aimes181, join_round, end_round, 2.57415564419716247e-08, 0.000036, settings);
 282         test_one<linestring, polygon>("aimes181", aimes181, join_round_by_divide, end_round, 2.57415564419716247e-08, 0.000036, settings);
 283     }
 284
 285     {
 286         // Expectations can also be 1702.56530051454502 2140.78725663358819
 287         // so we increase tolerance
 288         ut_settings settings(0.5);
 289         test_one<linestring, polygon>("crossing", crossing, join_round32, end_flat, 1702.1, 20.0, settings);
 290         test_one<linestring, polygon>("crossing", crossing, join_round32, end_round32, 2140.4, 20.0, settings);
 291     }
 292
 293     test_one<linestring, polygon>("mikado1", mikado1, join_round32, end_round32, 5441135039.0979, 41751.0);
 294
 295     {
 296         // This case gave a spike (invalid result) up to 1.70
 297         // Fixed by general form for line/line intersection in miter calculation
 298         using bg::strategy::buffer::join_round;
 299         using bg::strategy::buffer::end_round;
 300         test_one<linestring, polygon>("issue_596", issue_596, join_round(12), end_round(12), 0.12462779, 0.015);
 301         test_one<linestring, polygon>("issue_596", issue_596, join_miter, end_round(12), 0.12462807, 0.015);
 302     }
 303
 304     {
 305         using bg::strategy::buffer::join_round;
 306         using bg::strategy::buffer::end_round;
 307         test_one<linestring, polygon>("issue_803", issue_803, join_miter, end_round(36), 1664.3714, 10.0);
 308         test_one<linestring, polygon>("issue_803", issue_803, join_round(36), end_round(36), 1664.0528, 10.0);
 309     }
 310
 311     test_one<linestring, polygon>("mysql_report_2015_06_11",
 312             mysql_report_2015_06_11, join_round32, end_round32,
 313             27862.733459829971,
 314             5.9518403867035365);
 315
 316     if (BOOST_GEOMETRY_CONDITION((std::is_same<coor_type, double>::value)))
 317     {
 318         test_one<linestring, polygon>("mysql_report_2015_09_08a", mysql_report_2015_09_08a, join_round32, end_round32, 0.0, 1.0);
 319         test_one<linestring, polygon>("mysql_report_2015_09_08b", mysql_report_2015_09_08b, join_round32, end_round32, 0.0, 1099511627778.0);
 320         test_one<linestring, polygon>("mysql_report_2015_09_08c", mysql_report_2015_09_08c, join_round32, end_round32, 0.0, 0xbe);
 321     }
 322
 323     test_one<linestring, polygon>("mysql_23023665_1", mysql_23023665, join_round32, end_flat, 459.1051, 10);
 324     test_one<linestring, polygon>("mysql_23023665_2", mysql_23023665, join_round32, end_flat, 6877.7097, 50);
 325
 326     test_one<linestring, polygon>("mysql_25662426", mysql_25662426, join_round32, end_round32, 1, 0, 1660.6673, 10);
 327
 328     // Test behaviour with different buffer sizes, generating internally turns on different locations
 329     test_one<linestring, polygon>("mysql_25662426a_05", mysql_25662426a, join_round32, end_round32, 27.6156, 0.5);
 330 #if defined(BOOST_GEOMETRY_USE_RESCALING) || defined(BOOST_GEOMETRY_TEST_FAILURES)
 331     // Fails without rescaling
 332     test_one<linestring, polygon>("mysql_25662426a_1", mysql_25662426a, join_round32, end_round32, 54.9018, 1.0);
 333 #endif
 334     test_one<linestring, polygon>("mysql_25662426a_2", mysql_25662426a, join_round32, end_round32, 103.6072, 2.0);
 335     test_one<linestring, polygon>("mysql_25662426a_3", mysql_25662426a, join_round32, end_round32, 152.1163, 3.0);
 336     test_one<linestring, polygon>("mysql_25662426a_4", mysql_25662426a, join_round32, end_round32, 206.4831, 4.0);
 337     test_one<linestring, polygon>("mysql_25662426a_5", mysql_25662426a, join_round32, end_round32, 266.8505, 5.0);
 338     test_one<linestring, polygon>("mysql_25662426a_10", mysql_25662426a, join_round32, end_round32, 660.7355, 10.0);
 339
 340     test_one<linestring, polygon>("mysql_25662426a_05", mysql_25662426a, join_round32, end_flat, 26.8352, 0.5);
 341     test_one<linestring, polygon>("mysql_25662426a_1", mysql_25662426a, join_round32, end_flat, 53.3411, 1.0);
 342     test_one<linestring, polygon>("mysql_25662426a_2", mysql_25662426a, join_round32, end_flat, 97.3644, 2.0);
 343     test_one<linestring, polygon>("mysql_25662426a_3", mysql_25662426a, join_round32, end_flat, 138.0697, 3.0);
 344     test_one<linestring, polygon>("mysql_25662426a_4", mysql_25662426a, join_round32, end_flat, 181.5115, 4.0);
 345     test_one<linestring, polygon>("mysql_25662426a_5", mysql_25662426a, join_round32, end_flat, 227.8325, 5.0);
 346     test_one<linestring, polygon>("mysql_25662426a_10", mysql_25662426a, join_round32, end_flat, 534.1083, 10.0, specific_settings);
 347
 348     // Asymmetric buffers
 349     // Mostly left
 350     test_one<linestring, polygon>("mysql_25662426a_mostly_left_05", mysql_25662426a, join_round32, end_round32, 15.6638, 0.5, settings, 0.05);
 351     test_one<linestring, polygon>("mysql_25662426a_mostly_left_1", mysql_25662426a, join_round32, end_round32, 32.9553, 1.0, settings, 0.1);
 352     test_one<linestring, polygon>("mysql_25662426a_mostly_left_2", mysql_25662426a, join_round32, end_round32, 72.1159, 2.0, settings, 0.2);
 353     test_one<linestring, polygon>("mysql_25662426a_mostly_left_3", mysql_25662426a, join_round32, end_round32, 116.3802, 3.0, settings, 0.3);
 354     test_one<linestring, polygon>("mysql_25662426a_mostly_left_4", mysql_25662426a, join_round32, end_round32, 165.9298, 4.0, specific_settings, 0.4);
 355     test_one<linestring, polygon>("mysql_25662426a_mostly_left_5", mysql_25662426a, join_round32, end_round32, 220.8054, 5.0, specific_settings, 0.5);
 356     test_one<linestring, polygon>("mysql_25662426a_mostly_left_10", mysql_25662426a, join_round32, end_round32, 577.3742, 10.0, specific_settings, 1.0);
 357
 358     // Mostly right
 359     test_one<linestring, polygon>("mysql_25662426a_mostly_right_05", mysql_25662426a, join_round32, end_round32, 14.3419, 0.05, settings, 0.5);
 360 #if defined(BOOST_GEOMETRY_USE_RESCALING) || defined(BOOST_GEOMETRY_TEST_FAILURES)
 361     // Fails without rescaling
 362     test_one<linestring, polygon>("mysql_25662426a_mostly_right_1", mysql_25662426a, join_round32, end_round32, 27.1955, 0.1, settings, 1.0);
 363 #endif
 364     test_one<linestring, polygon>("mysql_25662426a_mostly_right_2", mysql_25662426a, join_round32, end_round32, 43.1821, 0.2, settings, 2.0);
 365     test_one<linestring, polygon>("mysql_25662426a_mostly_right_3", mysql_25662426a, join_round32, end_round32, 54.4337, 0.3, settings, 3.0);
 366     test_one<linestring, polygon>("mysql_25662426a_mostly_right_4", mysql_25662426a, join_round32, end_round32, 75.6376, 0.4, settings, 4.0);
 367     test_one<linestring, polygon>("mysql_25662426a_mostly_right_5", mysql_25662426a, join_round32, end_round32, 102.6952, 0.5, settings, 5.0);
 368     test_one<linestring, polygon>("mysql_25662426a_mostly_right_10", mysql_25662426a, join_round32, end_round32, 260.1850, 1.0, settings, 10.0);
 369
 370     // Left
 371     test_one<linestring, polygon>("mysql_25662426a_left_05", mysql_25662426a, join_round32, end_round32, 14.2803, 0.5, settings, 0.0);
 372     test_one<linestring, polygon>("mysql_25662426a_left_1", mysql_25662426a, join_round32, end_round32, 30.1214, 1.0, settings, 0.0);
 373     test_one<linestring, polygon>("mysql_25662426a_left_2", mysql_25662426a, join_round32, end_round32, 66.4858, 2.0, ut_settings(0.01, false), 0.0); // It has a self touching point
 374     test_one<linestring, polygon>("mysql_25662426a_left_3", mysql_25662426a, join_round32, end_round32, 108.3305, 3.0, settings, 0.0);
 375     test_one<linestring, polygon>("mysql_25662426a_left_4", mysql_25662426a, join_round32, end_round32, 155.5128, 4.0, specific_settings, 0.0);
 376     test_one<linestring, polygon>("mysql_25662426a_left_5", mysql_25662426a, join_round32, end_round32, 208.1289, 5.0, specific_settings, 0.0);
 377     test_one<linestring, polygon>("mysql_25662426a_left_10", mysql_25662426a, join_round32, end_round32, 554.8818, 10.0, specific_settings, 0.0);
 378
 379     // Right
 380     test_one<linestring, polygon>("mysql_25662426a_right_05", mysql_25662426a, join_round32, end_round32, 12.9451, 0.0, settings, 0.5);
 381     test_one<linestring, polygon>("mysql_25662426a_right_1", mysql_25662426a, join_round32, end_round32, 24.3902, 0.0, settings, 1.0);
 382     test_one<linestring, polygon>("mysql_25662426a_right_2", mysql_25662426a, join_round32, end_round32, 37.5607, 0.0, settings, 2.0);
 383     test_one<linestring, polygon>("mysql_25662426a_right_3", mysql_25662426a, join_round32, end_round32, 46.2741, 0.0, settings, 3.0);
 384     test_one<linestring, polygon>("mysql_25662426a_right_4", mysql_25662426a, join_round32, end_round32, 70.2429, 0.0, settings, 4.0);
 385     test_one<linestring, polygon>("mysql_25662426a_right_5", mysql_25662426a, join_round32, end_round32, 95.7545, 0.0, settings, 5.0);
 386     test_one<linestring, polygon>("mysql_25662426a_right_10", mysql_25662426a, join_round32, end_round32, 244.5337, 0.0, settings, 10.0);
 387
 388     // Flat end for one-sided buffer
 389     test_one<linestring, polygon>("mysql_25662426a_left_05", mysql_25662426a, join_round32, end_flat, 14.0853, 0.5, settings, 0.0);
 390     test_one<linestring, polygon>("mysql_25662426a_right_05", mysql_25662426a, join_round32, end_flat, 12.7500, 0.0, settings, 0.5);
 391     test_one<linestring, polygon>("mysql_25662426a_left_25", mysql_25662426a, join_round32, end_flat, 82.1318, 2.5, settings, 0.0);
 392     test_one<linestring, polygon>("mysql_25662426a_right_25", mysql_25662426a, join_round32, end_flat, 39.0000, 0.0, settings, 2.5);
 393 }
 394
 395 template <bool Clockwise, typename P>
 396 void test_invalid()
 397 {
 398     typedef typename bg::coordinate_type<P>::type coor_type;
 399     if (! BOOST_GEOMETRY_CONDITION((std::is_same<coor_type, double>::value)))
 400     {
 401         return;
 402     }
 403     typedef bg::model::linestring<P> linestring;
 404     typedef bg::model::polygon<P, Clockwise> polygon;
 405
 406     bg::strategy::buffer::end_round end_round32(32);
 407     bg::strategy::buffer::join_round join_round32(32);
 408
 409     // Linestring contains extreme differences causing numeric errors. Empty geometry is returned
 410     test_one<linestring, polygon>("mysql_report_2015_04_10a", mysql_report_2015_04_10a, join_round32, end_round32, 0.0, 100.0);
 411     test_one<linestring, polygon>("mysql_report_2015_04_10b", mysql_report_2015_04_10b, join_round32, end_round32, 0.0, 100.0);
 412     test_one<linestring, polygon>("mysql_report_2015_04_10c", mysql_report_2015_04_10c, join_round32, end_round32, 0.0, 100.0);
 413     test_one<linestring, polygon>("mysql_report_2015_04_10d", mysql_report_2015_04_10d, join_round32, end_round32, 0.0, 100.0);
 414     test_one<linestring, polygon>("mysql_report_2015_04_10e", mysql_report_2015_04_10e, join_round32, end_round32, 0.0, 100.0);
 415     test_one<linestring, polygon>("mysql_report_2015_04_10f", mysql_report_2015_04_10f, join_round32, end_round32, 0.0, 100.0);
 416
 417     // The equivalent, valid, case
 418     test_one<linestring, polygon>("mysql_report_2015_04_10g", mysql_report_2015_04_10g, join_round32, end_round32, 86527.871, 100.0);
 419 }
 420
 421
 422 int test_main(int, char* [])
 423 {
 424     BoostGeometryWriteTestConfiguration();
 425
 426     test_all<true, bg::model::point<default_test_type, 2, bg::cs::cartesian> >();
 427
 428 #if ! defined(BOOST_GEOMETRY_TEST_ONLY_ONE_ORDER)
 429     test_all<false, bg::model::point<default_test_type, 2, bg::cs::cartesian> >();
 430 #endif
 431
 432     test_invalid<true, bg::model::point<default_test_type, 2, bg::cs::cartesian> >();
 433
 434 #if ! defined(BOOST_GEOMETRY_TEST_ONLY_ONE_TYPE)
 435     test_invalid<true, bg::model::point<long double, 2, bg::cs::cartesian> >();
 436 #endif
 437
 438 #if defined(BOOST_GEOMETRY_TEST_FAILURES)
 439     BoostGeometryWriteExpectedFailures(2, 4, 11, 3);
 440 #endif
 441
 442     return 0;
 443 }