Updated robustness test, added linear/areal (tested for 1.49)

[SVN r77129]
This commit is contained in:
Barend Gehrels 2012-02-27 19:23:05 +00:00
parent 577f86d246
commit e14580efdb
8 changed files with 860 additions and 1 deletions

View File

@ -0,0 +1,32 @@
// Boost.Geometry (aka GGL, Generic Geometry Library)
// Robustness Test
//
// Copyright (c) 2009-2012 Barend Gehrels, Amsterdam, the Netherlands.
// Use, modification and distribution is subject to the Boost Software License,
// Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
// http://www.boost.org/LICENSE_1_0.txt)
#ifndef BOOST_GEOMETRY_COMMON_SETTINGS_HPP
#define BOOST_GEOMETRY_COMMON_SETTINGS_HPP
struct common_settings
{
bool svg;
bool wkt;
bool also_difference;
double tolerance;
int field_size;
bool triangular;
common_settings()
: svg(false)
, wkt(false)
, also_difference(false)
, tolerance(1.0e-6)
, field_size(10)
, triangular(false)
{}
};
#endif // BOOST_GEOMETRY_COMMON_SETTINGS_HPP

View File

@ -0,0 +1,46 @@
// Boost.Geometry (aka GGL, Generic Geometry Library)
// Robustness Test
//
// Copyright (c) 2009-2012 Barend Gehrels, Amsterdam, the Netherlands.
// Use, modification and distribution is subject to the Boost Software License,
// Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
// http://www.boost.org/LICENSE_1_0.txt)
#ifndef BOOST_GEOMETRY_TEST_ROBUSTNESS_MAKE_SQUARE_POLYGON_HPP
#define BOOST_GEOMETRY_TEST_ROBUSTNESS_MAKE_SQUARE_POLYGON_HPP
#include <boost/geometry.hpp>
template <typename Polygon, typename Generator, typename Settings>
inline void make_square_polygon(Polygon& polygon, Generator& generator, Settings const& settings)
{
using namespace boost::geometry;
typedef typename point_type<Polygon>::type point_type;
typedef typename coordinate_type<Polygon>::type coordinate_type;
coordinate_type x, y;
x = generator();
y = generator();
typename ring_type<Polygon>::type& ring = exterior_ring(polygon);
point_type p;
set<0>(p, x); set<1>(p, y); append(ring, p);
set<0>(p, x); set<1>(p, y + 1); append(ring, p);
set<0>(p, x + 1); set<1>(p, y + 1); append(ring, p);
set<0>(p, x + 1); set<1>(p, y); append(ring, p);
set<0>(p, x); set<1>(p, y); append(ring, p);
if (settings.triangular)
{
// Remove a point, generator says which
int c = generator() % 4;
if (c >= 1 && c <= 3)
{
ring.erase(ring.begin() + c);
}
}
}
#endif // BOOST_GEOMETRY_TEST_ROBUSTNESS_MAKE_SQUARE_POLYGON_HPP

View File

@ -0,0 +1,18 @@
# Boost.Geometry (aka GGL, Generic Geometry Library)
# Robustness Test - convex_hull
#
# Copyright (c) 2012 Barend Gehrels, Amsterdam, the Netherlands.
# Use, modification and distribution is subject to the Boost Software License,
# Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
# http://www.boost.org/LICENSE_1_0.txt)
project random_multi_points
: requirements
<include>.
<library>../../../../program_options/build//boost_program_options
<link>static
;
exe random_multi_points : random_multi_points.cpp ;

View File

@ -1,5 +1,5 @@
// Boost.Geometry (aka GGL, Generic Geometry Library)
// Robustness Test
// Robustness Test - convex_hull
// Copyright (c) 2012 Barend Gehrels, Amsterdam, the Netherlands.

View File

@ -0,0 +1,19 @@
# Boost.Geometry (aka GGL, Generic Geometry Library)
# Robustness Test - overlay - linear/areal
#
# Copyright (c) 2012 Barend Gehrels, Amsterdam, the Netherlands.
# Use, modification and distribution is subject to the Boost Software License,
# Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
# http://www.boost.org/LICENSE_1_0.txt)
project recursive_polygons_linear_areal
: requirements
<include>.
<include>../..
<library>../../../../../program_options/build//boost_program_options
<link>static
;
exe recursive_polygons_linear_areal : recursive_polygons_linear_areal.cpp ;

View File

@ -0,0 +1,501 @@
// Boost.Geometry (aka GGL, Generic Geometry Library)
// Robustness Test
// Copyright (c) 2012 Barend Gehrels, Amsterdam, the Netherlands.
// Use, modification and distribution is subject to the Boost Software License,
// Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
// http://www.boost.org/LICENSE_1_0.txt)
#if defined(_MSC_VER)
# pragma warning( disable : 4244 )
# pragma warning( disable : 4267 )
#endif
#include <fstream>
#include <sstream>
#include <boost/foreach.hpp>
#include <boost/program_options.hpp>
#include <boost/random/linear_congruential.hpp>
#include <boost/random/uniform_int.hpp>
#include <boost/random/uniform_real.hpp>
#include <boost/random/variate_generator.hpp>
#include <boost/timer.hpp>
#include <boost/geometry.hpp>
#include <boost/geometry/geometries/geometries.hpp>
#include <boost/geometry/geometries/point_xy.hpp>
#include <boost/geometry/multi/geometries/multi_geometries.hpp>
#include <boost/geometry/extensions/io/svg/svg_mapper.hpp>
#include <boost/geometry/extensions/algorithms/midpoints.hpp>
#include <common/common_settings.hpp>
#include <common/make_square_polygon.hpp>
namespace bg = boost::geometry;
template <typename Geometry1, typename Geometry2, typename Geometry3>
void create_svg(std::string const& filename
, Geometry1 const& mp
, Geometry2 const& ls
, Geometry3 const& difference
, Geometry3 const& intersection
)
{
typedef typename boost::geometry::point_type<Geometry1>::type point_type;
std::ofstream svg(filename.c_str());
boost::geometry::svg_mapper<point_type> mapper(svg, 800, 800);
boost::geometry::model::box<point_type> box;
bg::envelope(mp, box);
bg::buffer(box, box, 1.0);
mapper.add(box);
bg::envelope(ls, box);
bg::buffer(box, box, 1.0);
mapper.add(box);
mapper.map(mp, "fill-opacity:0.5;fill:rgb(153,204,0);stroke:rgb(153,204,0);stroke-width:3");
mapper.map(ls, "stroke-opacity:0.9;stroke:rgb(0,0,0);stroke-width:1");
mapper.map(intersection,"opacity:0.6;stroke:rgb(0,128,0);stroke-width:5");
mapper.map(difference, "opacity:0.6;stroke:rgb(255,0,255);stroke-width:5"); //;stroke-dasharray:1,7;stroke-linecap:round
}
template <typename Linestring, typename Generator, typename Settings>
inline void make_random_linestring(Linestring& line, Generator& generator, Settings const& settings)
{
using namespace boost::geometry;
typedef typename point_type<Linestring>::type point_type;
typedef typename coordinate_type<Linestring>::type coordinate_type;
coordinate_type x, y;
x = generator();
y = generator();
int count = 3 + generator() % 6;
int d = 0; // direction
for (int i = 0; i <= count && x <= settings.field_size; i++, x++, d = 1 - d)
{
append(line, make<point_type>(x, y + d));
append(line, make<point_type>(x, y + 1 - d));
}
if (d == 0 && generator() % 4 < 3 && y >= 2)
{
d = 1 - d;
x--;
y -= 2;
count = 3 + generator() % 6;
for (int i = 0; i <= count && x >= 0; i++, x--, d = 1 - d)
{
append(line, make<point_type>(x, y + d));
append(line, make<point_type>(x, y + 1 - d));
}
}
//if (settings.triangular)
//{
// // Remove a point, generator says which
// int c = generator() % 4;
// if (c >= 1 && c <= 3)
// {
// ring.erase(ring.begin() + c);
// }
//}
}
template<typename Geometry>
class inside_check
{
Geometry const& m_geo;
bool& m_result;
public :
inside_check(Geometry const& geo, bool& result)
: m_geo(geo)
, m_result(result)
{}
inline inside_check<Geometry> operator=(inside_check<Geometry> const& input)
{
return inside_check<Geometry>(input.m_geo, input.m_result);
}
template <typename Point>
inline void operator()(Point const& p)
{
if (! bg::covered_by(p, m_geo))
{
if (m_result)
{
std::cout << "Does not fulfill inside check" << std::endl;
}
m_result = false;
}
}
};
template<typename Geometry>
class outside_check
{
Geometry const& m_geo;
bool& m_result;
public :
outside_check(Geometry const& geo, bool& result)
: m_geo(geo)
, m_result(result)
{}
inline outside_check<Geometry> operator=(outside_check<Geometry> const& input)
{
return outside_check<Geometry>(input.m_geo, input.m_result);
}
template <typename Point>
inline void operator()(Point const& p)
{
if (bg::within(p, m_geo))
{
if (m_result)
{
std::cout << "Does not fulfill outside check" << std::endl;
}
m_result = false;
}
}
};
template<typename Segment>
class border2_check
{
Segment const& m_segment;
bool& m_result;
public :
border2_check(Segment const& seg, bool& result)
: m_segment(seg)
, m_result(result)
{}
inline border2_check<Segment> operator=(border2_check<Segment> const& input)
{
return border2_check<Segment>(input.m_segment, input.m_result);
}
template <typename Segment2>
inline void operator()(Segment2 const& segment)
{
// Create copies (TODO: find out why referring_segment does not compile)
typedef typename bg::point_type<Segment2>::type pt;
typedef bg::model::segment<pt> segment_type;
typedef bg::strategy::intersection::relate_cartesian_segments
<
bg::policies::relate::segments_intersection_points
<
segment_type,
segment_type,
bg::segment_intersection_points<pt>
>
> policy;
segment_type seg1, seg2;
bg::convert(m_segment, seg1);
bg::convert(segment, seg2);
bg::segment_intersection_points<pt> is = policy::apply(seg1, seg2);
if (is.count == 2)
{
if (m_result)
{
std::cout << "Does not fulfill border2 check" << std::endl;
}
m_result = true;
}
}
};
template<typename Geometry>
class border_check
{
Geometry const& m_geo;
bool& m_result;
public :
border_check(Geometry const& geo, bool& result)
: m_geo(geo)
, m_result(result)
{}
inline border_check<Geometry> operator=(border_check<Geometry> const& input)
{
return border_check<Geometry>(input.m_geo, input.m_result);
}
template <typename Segment>
inline void operator()(Segment const& s)
{
bool on_border = false;
border2_check<Segment> checker(s, on_border);
bg::for_each_segment(m_geo, checker);
if (on_border)
{
if (m_result)
{
std::cout << "Does not fulfill border check" << std::endl;
}
m_result = false;
}
}
};
template <typename MultiPolygon, typename Linestring, typename Settings>
bool verify(std::string const& caseid, MultiPolygon const& mp, Linestring const& ls, Settings const& settings)
{
bg::model::multi_linestring<Linestring> difference, intersection;
bg::difference(ls, mp, difference);
bg::intersection(ls, mp, intersection);
//typedef typename bg::length_result_type<point_type>::type length_type;
bool result = true;
// 1) Check by length
typedef double length_type;
length_type len_input = bg::length(ls);
length_type len_difference = bg::length(difference);
length_type len_intersection = bg::length(intersection);
if (! bg::math::equals(len_input, len_difference + len_intersection))
{
std::cout << "Input: " << len_input
<< " difference: " << len_difference
<< " intersection: " << len_intersection
<< std::endl;
std::cout << "Does not fulfill length check" << std::endl;
result = false;
}
// 2) Check by within and covered by
inside_check<MultiPolygon> ic(mp, result);
bg::for_each_point(intersection, ic);
outside_check<MultiPolygon> oc(mp, result);
bg::for_each_point(difference, oc);
border_check<MultiPolygon> bc(mp, result);
bg::for_each_segment(difference, bc);
// 3) check also the mid-points from the difference to remove false positives
BOOST_FOREACH(Linestring const& d, difference)
{
Linestring difference_midpoints;
bg::midpoints(d, false, std::back_inserter(difference_midpoints));
outside_check<MultiPolygon> ocm(mp, result);
bg::for_each_point(difference_midpoints, ocm);
}
bool svg = settings.svg;
bool wkt = settings.wkt;
if (! result)
{
std::cout << "ERROR " << caseid << std::endl;
std::cout << bg::wkt(mp) << std::endl;
std::cout << bg::wkt(ls) << std::endl;
svg = true;
wkt = true;
}
if (svg)
{
std::ostringstream filename;
filename << caseid << "_"
<< typeid(typename bg::coordinate_type<Linestring>::type).name()
<< ".svg";
create_svg(filename.str(), mp, ls, difference, intersection);
}
if (wkt)
{
std::ostringstream filename;
filename << caseid << "_"
<< typeid(typename bg::coordinate_type<Linestring>::type).name()
<< ".wkt";
std::ofstream stream(filename.str().c_str());
stream << bg::wkt(mp) << std::endl;
stream << bg::wkt(ls) << std::endl;
}
return result;
}
template <typename MultiPolygon, typename Generator>
bool test_linear_areal(MultiPolygon& result, int& index,
Generator& generator,
int level, common_settings const& settings)
{
MultiPolygon p, q;
// Generate two boxes
if (level == 0)
{
p.resize(1);
q.resize(1);
make_square_polygon(p.front(), generator, settings);
make_square_polygon(q.front(), generator, settings);
bg::correct(p);
bg::correct(q);
}
else
{
bg::correct(p);
bg::correct(q);
if (! test_linear_areal(p, index, generator, level - 1, settings)
|| ! test_linear_areal(q, index, generator, level - 1, settings))
{
return false;
}
}
typedef typename boost::range_value<MultiPolygon>::type polygon;
MultiPolygon mp;
bg::detail::union_::union_insert
<
polygon
>(p, q, std::back_inserter(mp));
bg::unique(mp);
bg::simplify(mp, result, 0.01);
bg::correct(mp);
// Generate a linestring
typedef typename bg::point_type<MultiPolygon>::type point_type;
typedef bg::model::linestring<point_type> linestring_type;
linestring_type ls;
make_random_linestring(ls, generator, settings);
std::ostringstream out;
out << "recursive_la_" << index++ << "_" << level;
return verify(out.str(), mp, ls, settings);
}
template <typename T, bool Clockwise, bool Closed>
void test_all(int seed, int count, int level, common_settings const& settings)
{
boost::timer t;
typedef boost::minstd_rand base_generator_type;
base_generator_type generator(seed);
boost::uniform_int<> random_coordinate(0, settings.field_size - 1);
boost::variate_generator<base_generator_type&, boost::uniform_int<> >
coordinate_generator(generator, random_coordinate);
typedef bg::model::polygon
<
bg::model::d2::point_xy<T>, Clockwise, Closed
> polygon;
typedef bg::model::multi_polygon<polygon> mp;
int index = 0;
for(int i = 0; i < count; i++)
{
mp p;
test_linear_areal<mp>(p, index, coordinate_generator, level, settings);
}
std::cout
<< "geometries: " << index
<< " type: " << typeid(T).name()
<< " time: " << t.elapsed() << std::endl;
}
int main(int argc, char** argv)
{
try
{
namespace po = boost::program_options;
po::options_description description("=== recursive_polygons_linear_areal ===\nAllowed options");
int count = 1;
int seed = static_cast<unsigned int>(std::time(0));
int level = 3;
bool ccw = false;
bool open = false;
common_settings settings;
std::string form = "box";
description.add_options()
("help", "Help message")
("seed", po::value<int>(&seed), "Initialization seed for random generator")
("count", po::value<int>(&count)->default_value(1), "Number of tests")
("diff", po::value<bool>(&settings.also_difference)->default_value(false), "Include testing on difference")
("level", po::value<int>(&level)->default_value(3), "Level to reach (higher->slower)")
("form", po::value<std::string>(&form)->default_value("box"), "Form of the polygons (box, triangle)")
("ccw", po::value<bool>(&ccw)->default_value(false), "Counter clockwise polygons")
("open", po::value<bool>(&open)->default_value(false), "Open polygons")
("size", po::value<int>(&settings.field_size)->default_value(10), "Size of the field")
("wkt", po::value<bool>(&settings.wkt)->default_value(false), "Create a WKT of the inputs, for all tests")
("svg", po::value<bool>(&settings.svg)->default_value(false), "Create a SVG for all tests")
;
po::variables_map varmap;
po::store(po::parse_command_line(argc, argv, description), varmap);
po::notify(varmap);
if (varmap.count("help")
|| (form != "box" && form != "triangle"))
{
std::cout << description << std::endl;
return 1;
}
settings.triangular = form != "box";
if (ccw && open)
{
test_all<double, false, false>(seed, count, level, settings);
}
else if (ccw)
{
test_all<double, false, true>(seed, count, level, settings);
}
else if (open)
{
test_all<double, true, false>(seed, count, level, settings);
}
else
{
test_all<double, true, true>(seed, count, level, settings);
}
#if defined(HAVE_TTMATH)
// test_all<ttmath_big, true, true>(seed, count, max, svg, level);
#endif
}
catch(std::exception const& e)
{
std::cout << "Exception " << e.what() << std::endl;
}
catch(...)
{
std::cout << "Other exception" << std::endl;
}
return 0;
}

View File

@ -0,0 +1,20 @@

Microsoft Visual Studio Solution File, Format Version 9.00
# Visual C++ Express 2005
Project("{8BC9CEB8-8B4A-11D0-8D11-00A0C91BC942}") = "recursive_polygons_linear_areal", "recursive_polygons_linear_areal.vcproj", "{1E269699-9450-4DD6-ACC3-C6A9AD2FC6E0}"
EndProject
Global
GlobalSection(SolutionConfigurationPlatforms) = preSolution
Debug|Win32 = Debug|Win32
Release|Win32 = Release|Win32
EndGlobalSection
GlobalSection(ProjectConfigurationPlatforms) = postSolution
{1E269699-9450-4DD6-ACC3-C6A9AD2FC6E0}.Debug|Win32.ActiveCfg = Debug|Win32
{1E269699-9450-4DD6-ACC3-C6A9AD2FC6E0}.Debug|Win32.Build.0 = Debug|Win32
{1E269699-9450-4DD6-ACC3-C6A9AD2FC6E0}.Release|Win32.ActiveCfg = Release|Win32
{1E269699-9450-4DD6-ACC3-C6A9AD2FC6E0}.Release|Win32.Build.0 = Release|Win32
EndGlobalSection
GlobalSection(SolutionProperties) = preSolution
HideSolutionNode = FALSE
EndGlobalSection
EndGlobal

View File

@ -0,0 +1,223 @@
<?xml version="1.0" encoding="Windows-1252"?>
<VisualStudioProject
ProjectType="Visual C++"
Version="8.00"
Name="recursive_polygons_linear_areal"
ProjectGUID="{1E269699-9450-4DD6-ACC3-C6A9AD2FC6E0}"
RootNamespace="recursive_polygons_linear_areal"
Keyword="Win32Proj"
>
<Platforms>
<Platform
Name="Win32"
/>
</Platforms>
<ToolFiles>
</ToolFiles>
<Configurations>
<Configuration
Name="Debug|Win32"
OutputDirectory="$(SolutionDir)$(ConfigurationName)"
IntermediateDirectory="$(ConfigurationName)\recursive_polygons_linear_areal"
ConfigurationType="1"
InheritedPropertySheets="..\..\..\boost.vsprops"
CharacterSet="1"
>
<Tool
Name="VCPreBuildEventTool"
/>
<Tool
Name="VCCustomBuildTool"
/>
<Tool
Name="VCXMLDataGeneratorTool"
/>
<Tool
Name="VCWebServiceProxyGeneratorTool"
/>
<Tool
Name="VCMIDLTool"
/>
<Tool
Name="VCCLCompilerTool"
Optimization="0"
AdditionalIncludeDirectories="../../../../../..;.;../.."
PreprocessorDefinitions="WIN32;_DEBUG;_WINDOWS;BOOST_ALL_NO_LIB"
BasicRuntimeChecks="3"
RuntimeLibrary="1"
UsePrecompiledHeader="0"
WarningLevel="2"
Detect64BitPortabilityProblems="true"
DebugInformationFormat="1"
/>
<Tool
Name="VCManagedResourceCompilerTool"
/>
<Tool
Name="VCResourceCompilerTool"
/>
<Tool
Name="VCPreLinkEventTool"
/>
<Tool
Name="VCLinkerTool"
AdditionalLibraryDirectories=""
GenerateDebugInformation="true"
SubSystem="1"
TargetMachine="1"
/>
<Tool
Name="VCALinkTool"
/>
<Tool
Name="VCManifestTool"
/>
<Tool
Name="VCXDCMakeTool"
/>
<Tool
Name="VCBscMakeTool"
/>
<Tool
Name="VCFxCopTool"
/>
<Tool
Name="VCAppVerifierTool"
/>
<Tool
Name="VCWebDeploymentTool"
/>
<Tool
Name="VCPostBuildEventTool"
/>
</Configuration>
<Configuration
Name="Release|Win32"
OutputDirectory="$(SolutionDir)$(ConfigurationName)"
IntermediateDirectory="$(ConfigurationName)\recursive_polygons_linear_areal"
ConfigurationType="1"
InheritedPropertySheets="..\..\..\boost.vsprops"
CharacterSet="1"
WholeProgramOptimization="0"
>
<Tool
Name="VCPreBuildEventTool"
/>
<Tool
Name="VCCustomBuildTool"
/>
<Tool
Name="VCXMLDataGeneratorTool"
/>
<Tool
Name="VCWebServiceProxyGeneratorTool"
/>
<Tool
Name="VCMIDLTool"
/>
<Tool
Name="VCCLCompilerTool"
AdditionalIncludeDirectories="../../../../../..;.;../.."
PreprocessorDefinitions="WIN32;NDEBUG;_WINDOWS;BOOST_ALL_NO_LIB"
RuntimeLibrary="0"
UsePrecompiledHeader="0"
WarningLevel="3"
Detect64BitPortabilityProblems="true"
/>
<Tool
Name="VCManagedResourceCompilerTool"
/>
<Tool
Name="VCResourceCompilerTool"
/>
<Tool
Name="VCPreLinkEventTool"
/>
<Tool
Name="VCLinkerTool"
SubSystem="1"
OptimizeReferences="2"
EnableCOMDATFolding="2"
TargetMachine="1"
/>
<Tool
Name="VCALinkTool"
/>
<Tool
Name="VCManifestTool"
/>
<Tool
Name="VCXDCMakeTool"
/>
<Tool
Name="VCBscMakeTool"
/>
<Tool
Name="VCFxCopTool"
/>
<Tool
Name="VCAppVerifierTool"
/>
<Tool
Name="VCWebDeploymentTool"
/>
<Tool
Name="VCPostBuildEventTool"
/>
</Configuration>
</Configurations>
<References>
</References>
<Files>
<Filter
Name="program options"
>
<File
RelativePath="..\..\..\..\..\program_options\src\cmdline.cpp"
>
</File>
<File
RelativePath="..\..\..\..\..\program_options\src\config_file.cpp"
>
</File>
<File
RelativePath="..\..\..\..\..\program_options\src\convert.cpp"
>
</File>
<File
RelativePath="..\..\..\..\..\program_options\src\options_description.cpp"
>
</File>
<File
RelativePath="..\..\..\..\..\program_options\src\parsers.cpp"
>
</File>
<File
RelativePath="..\..\..\..\..\program_options\src\positional_options.cpp"
>
</File>
<File
RelativePath="..\..\..\..\..\program_options\src\split.cpp"
>
</File>
<File
RelativePath="..\..\..\..\..\program_options\src\utf8_codecvt_facet.cpp"
>
</File>
<File
RelativePath="..\..\..\..\..\program_options\src\value_semantic.cpp"
>
</File>
<File
RelativePath="..\..\..\..\..\program_options\src\variables_map.cpp"
>
</File>
</Filter>
<File
RelativePath=".\recursive_polygons_linear_areal.cpp"
>
</File>
</Files>
<Globals>
</Globals>
</VisualStudioProject>