geometry/test/algorithms/comparable_distance.cpp at develop · boostorg/geometry · GitHub

History

175 lines (131 loc) · 5.32 KB

Raw

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

164

165

166

167

168

169

170

171

172

173

174

175

// Boost.Geometry (aka GGL, Generic Geometry Library)

// Unit Test

// This file was modified by Oracle on 2014-2021.

// Contributed and/or modified by Menelaos Karavelas, on behalf of Oracle

// Contributed and/or modified by Adam Wulkiewicz, on behalf of Oracle

// Parts of Boost.Geometry are redesigned from Geodan's Geographic Library

// 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)

#include <sstream>

#include <geometry_test_common.hpp>

#include <boost/geometry/algorithms/comparable_distance.hpp>

#include <boost/geometry/algorithms/make.hpp>

#include <boost/geometry/geometries/geometries.hpp>

#include <boost/geometry/geometries/point_xy.hpp>

#include <boost/geometry/io/wkt/read.hpp>

#include <boost/geometry/strategies/strategies.hpp>

template <typename P>

void test_distance_result()

{

typedef typename bg::default_distance_result<P, P>::type distance_type;

P p1 = bg::make(0, 0);

P p2 = bg::make(3, 0);

P p3 = bg::make(0, 4);

distance_type dr12 = bg::comparable_distance(p1, p2);

distance_type dr13 = bg::comparable_distance(p1, p3);

distance_type dr23 = bg::comparable_distance(p2, p3);

BOOST_CHECK_CLOSE(dr12, 9.000, 0.001);

BOOST_CHECK_CLOSE(dr13, 16.000, 0.001);

BOOST_CHECK_CLOSE(dr23, 25.000, 0.001);

}

template <typename P>

void test_distance_point()

{

P p1;

bg::set<0>(p1, 1);

bg::set<1>(p1, 1);

P p2;

bg::set<0>(p2, 2);

bg::set<1>(p2, 2);

typename bg::coordinate_type::type d = bg::comparable_distance(p1, p2);

BOOST_CHECK_CLOSE(d, 2.0, 0.001);

}

template <typename P>

void test_distance_segment()

{

typedef typename bg::coordinate_type::type coordinate_type;

P s1 = bg::make(2, 2);

P s2 = bg::make(3, 3);

// Check points left, right, projected-left, projected-right, on segment

P p1 = bg::make(0, 0);

P p2 = bg::make(4, 4);

P p3 = bg::make(2.4, 2.6);

P p4 = bg::make(2.6, 2.4);

P p5 = bg::make(2.5, 2.5);

bg::model::referring_segment const seg(s1, s2);

coordinate_type d1 = bg::comparable_distance(p1, seg); BOOST_CHECK_CLOSE(d1, 8.0, 0.001);

coordinate_type d2 = bg::comparable_distance(p2, seg); BOOST_CHECK_CLOSE(d2, 2.0, 0.001);

coordinate_type d3 = bg::comparable_distance(p3, seg); BOOST_CHECK_CLOSE(d3, 0.02, 0.001);

coordinate_type d4 = bg::comparable_distance(p4, seg); BOOST_CHECK_CLOSE(d4, 0.02, 0.001);

coordinate_type d5 = bg::comparable_distance(p5, seg); BOOST_CHECK_CLOSE(d5, 0.0, 0.001);

// Reverse case

coordinate_type dr1 = bg::comparable_distance(seg, p1); BOOST_CHECK_CLOSE(dr1, d1, 0.001);

coordinate_type dr2 = bg::comparable_distance(seg, p2); BOOST_CHECK_CLOSE(dr2, d2, 0.001);

}

template <typename P>

void test_distance_linestring()

{

bg::model::linestring points;

points.push_back(bg::make(1, 1));

points.push_back(bg::make(3, 3));

P p = bg::make(2, 1);

typename bg::coordinate_type::type d = bg::comparable_distance(p, points);

BOOST_CHECK_CLOSE(d, 0.5, 0.001);

p = bg::make(5, 5);

d = bg::comparable_distance(p, points);

BOOST_CHECK_CLOSE(d, 8.0, 0.001);

bg::model::linestring line;

line.push_back(bg::make(1,1));

line.push_back(bg::make(2,2));

line.push_back(bg::make(3,3));

p = bg::make(5, 5);

d = bg::comparable_distance(p, line);

BOOST_CHECK_CLOSE(d, 8.0, 0.001);

// Reverse case

d = bg::comparable_distance(line, p);

BOOST_CHECK_CLOSE(d, 8.0, 0.001);

}

template <typename P>

void test_all()

{

test_distance_result();

test_distance_point();

test_distance_segment();

test_distance_linestring();

}

template <typename T>

void test_double_result_from_integer()

{

typedef bg::model::point<T, 2, bg::cs::cartesian> point_type;

point_type point;

// Check linestring

bg::model::linestring<point_type> linestring;

bg::read_wkt("POINT(2 2)", point);

bg::read_wkt("LINESTRING(4 1,1 4)", linestring);

double normal_distance = bg::distance(point, linestring);

double comparable_distance = bg::comparable_distance(point, linestring);

BOOST_CHECK_CLOSE(normal_distance, std::sqrt(0.5), 0.001);

BOOST_CHECK_CLOSE(comparable_distance, 0.5, 0.001);

// Check polygon

bg::model::polygon<point_type> polygon;

bg::read_wkt("POLYGON((0 0,1 9,8 1,0 0),(1 1,4 1,1 4,1 1))", polygon);

normal_distance = bg::distance(point, polygon);

comparable_distance = bg::comparable_distance(point, polygon);

BOOST_CHECK_CLOSE(normal_distance, std::sqrt(0.5), 0.001);

BOOST_CHECK_CLOSE(comparable_distance, 0.5, 0.001);

}

int test_main(int, char* [])

{

test_double_result_from_integer<int>();

test_double_result_from_integer<long long>();

test_all<bg::model::d2::point_xy<float> >();

test_all<bg::model::d2::point_xy<double> >();

return 0;

}

geometry/test/algorithms/comparable_distance.cpp at develop · boostorg/geometry · GitHub

Uh oh!

Expand file tree

Breadcrumbs

comparable_distance.cpp

Latest commit

History

Breadcrumbs

comparable_distance.cpp

File metadata and controls