d95794ec8a
As many open source projects have started doing it, we're removing the current year from the copyright notice, so that we don't need to bump it every year. It seems like only the first year of publication is technically relevant for copyright notices, and even that seems to be something that many companies stopped listing altogether (in a version controlled codebase, the commits are a much better source of date of publication than a hardcoded copyright statement). We also now list Godot Engine contributors first as we're collectively the current maintainers of the project, and we clarify that the "exclusive" copyright of the co-founders covers the timespan before opensourcing (their further contributions are included as part of Godot Engine contributors). Also fixed "cf." Frenchism - it's meant as "refer to / see".
866 lines
37 KiB
C++
866 lines
37 KiB
C++
/**************************************************************************/
|
|
/* test_geometry_2d.h */
|
|
/**************************************************************************/
|
|
/* This file is part of: */
|
|
/* GODOT ENGINE */
|
|
/* https://godotengine.org */
|
|
/**************************************************************************/
|
|
/* Copyright (c) 2014-present Godot Engine contributors (see AUTHORS.md). */
|
|
/* Copyright (c) 2007-2014 Juan Linietsky, Ariel Manzur. */
|
|
/* */
|
|
/* Permission is hereby granted, free of charge, to any person obtaining */
|
|
/* a copy of this software and associated documentation files (the */
|
|
/* "Software"), to deal in the Software without restriction, including */
|
|
/* without limitation the rights to use, copy, modify, merge, publish, */
|
|
/* distribute, sublicense, and/or sell copies of the Software, and to */
|
|
/* permit persons to whom the Software is furnished to do so, subject to */
|
|
/* the following conditions: */
|
|
/* */
|
|
/* The above copyright notice and this permission notice shall be */
|
|
/* included in all copies or substantial portions of the Software. */
|
|
/* */
|
|
/* THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, */
|
|
/* EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF */
|
|
/* MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. */
|
|
/* IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY */
|
|
/* CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT, */
|
|
/* TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE */
|
|
/* SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE. */
|
|
/**************************************************************************/
|
|
|
|
#ifndef TEST_GEOMETRY_2D_H
|
|
#define TEST_GEOMETRY_2D_H
|
|
|
|
#include "core/math/geometry_2d.h"
|
|
|
|
#include "thirdparty/doctest/doctest.h"
|
|
|
|
namespace TestGeometry2D {
|
|
|
|
TEST_CASE("[Geometry2D] Point in circle") {
|
|
CHECK(Geometry2D::is_point_in_circle(Vector2(0, 0), Vector2(0, 0), 1.0));
|
|
|
|
CHECK(Geometry2D::is_point_in_circle(Vector2(0, 0), Vector2(11.99, 0), 12));
|
|
CHECK(Geometry2D::is_point_in_circle(Vector2(-11.99, 0), Vector2(0, 0), 12));
|
|
|
|
CHECK_FALSE(Geometry2D::is_point_in_circle(Vector2(0, 0), Vector2(12.01, 0), 12));
|
|
CHECK_FALSE(Geometry2D::is_point_in_circle(Vector2(-12.01, 0), Vector2(0, 0), 12));
|
|
|
|
CHECK(Geometry2D::is_point_in_circle(Vector2(7, -42), Vector2(4, -40), 3.7));
|
|
CHECK_FALSE(Geometry2D::is_point_in_circle(Vector2(7, -42), Vector2(4, -40), 3.5));
|
|
|
|
// This tests points on the edge of the circle. They are treated as being inside the circle.
|
|
CHECK(Geometry2D::is_point_in_circle(Vector2(1.0, 0.0), Vector2(0, 0), 1.0));
|
|
CHECK(Geometry2D::is_point_in_circle(Vector2(0.0, -1.0), Vector2(0, 0), 1.0));
|
|
}
|
|
|
|
TEST_CASE("[Geometry2D] Point in triangle") {
|
|
CHECK(Geometry2D::is_point_in_triangle(Vector2(0, 0), Vector2(-1, 1), Vector2(0, -1), Vector2(1, 1)));
|
|
CHECK_FALSE(Geometry2D::is_point_in_triangle(Vector2(-1.01, 1.0), Vector2(-1, 1), Vector2(0, -1), Vector2(1, 1)));
|
|
|
|
CHECK(Geometry2D::is_point_in_triangle(Vector2(3, 2.5), Vector2(1, 4), Vector2(3, 2), Vector2(5, 4)));
|
|
CHECK(Geometry2D::is_point_in_triangle(Vector2(-3, -2.5), Vector2(-1, -4), Vector2(-3, -2), Vector2(-5, -4)));
|
|
CHECK_FALSE(Geometry2D::is_point_in_triangle(Vector2(0, 0), Vector2(1, 4), Vector2(3, 2), Vector2(5, 4)));
|
|
|
|
// This tests points on the edge of the triangle. They are treated as being outside the triangle.
|
|
// In `is_point_in_circle` and `is_point_in_polygon` they are treated as being inside, so in order the make
|
|
// the behavior consistent this may change in the future (see issue #44717 and PR #44274).
|
|
CHECK_FALSE(Geometry2D::is_point_in_triangle(Vector2(1, 1), Vector2(-1, 1), Vector2(0, -1), Vector2(1, 1)));
|
|
CHECK_FALSE(Geometry2D::is_point_in_triangle(Vector2(0, 1), Vector2(-1, 1), Vector2(0, -1), Vector2(1, 1)));
|
|
}
|
|
|
|
TEST_CASE("[Geometry2D] Point in polygon") {
|
|
Vector<Vector2> p;
|
|
CHECK_FALSE(Geometry2D::is_point_in_polygon(Vector2(0, 0), p));
|
|
|
|
p.push_back(Vector2(-88, 120));
|
|
p.push_back(Vector2(-74, -38));
|
|
p.push_back(Vector2(135, -145));
|
|
p.push_back(Vector2(425, 70));
|
|
p.push_back(Vector2(68, 112));
|
|
p.push_back(Vector2(-120, 370));
|
|
p.push_back(Vector2(-323, -145));
|
|
CHECK_FALSE(Geometry2D::is_point_in_polygon(Vector2(-350, 0), p));
|
|
CHECK_FALSE(Geometry2D::is_point_in_polygon(Vector2(-110, 60), p));
|
|
CHECK_FALSE(Geometry2D::is_point_in_polygon(Vector2(412, 96), p));
|
|
CHECK_FALSE(Geometry2D::is_point_in_polygon(Vector2(83, 130), p));
|
|
CHECK_FALSE(Geometry2D::is_point_in_polygon(Vector2(-320, -153), p));
|
|
|
|
CHECK(Geometry2D::is_point_in_polygon(Vector2(0, 0), p));
|
|
CHECK(Geometry2D::is_point_in_polygon(Vector2(-230, 0), p));
|
|
CHECK(Geometry2D::is_point_in_polygon(Vector2(130, -110), p));
|
|
CHECK(Geometry2D::is_point_in_polygon(Vector2(370, 55), p));
|
|
CHECK(Geometry2D::is_point_in_polygon(Vector2(-160, 190), p));
|
|
|
|
// This tests points on the edge of the polygon. They are treated as being inside the polygon.
|
|
int c = p.size();
|
|
for (int i = 0; i < c; i++) {
|
|
const Vector2 &p1 = p[i];
|
|
CHECK(Geometry2D::is_point_in_polygon(p1, p));
|
|
|
|
const Vector2 &p2 = p[(i + 1) % c];
|
|
Vector2 midpoint((p1 + p2) * 0.5);
|
|
CHECK(Geometry2D::is_point_in_polygon(midpoint, p));
|
|
}
|
|
}
|
|
|
|
TEST_CASE("[Geometry2D] Polygon clockwise") {
|
|
Vector<Vector2> p;
|
|
CHECK_FALSE(Geometry2D::is_polygon_clockwise(p));
|
|
|
|
p.push_back(Vector2(5, -5));
|
|
p.push_back(Vector2(-1, -5));
|
|
p.push_back(Vector2(-5, -1));
|
|
p.push_back(Vector2(-1, 3));
|
|
p.push_back(Vector2(1, 5));
|
|
CHECK(Geometry2D::is_polygon_clockwise(p));
|
|
|
|
p.reverse();
|
|
CHECK_FALSE(Geometry2D::is_polygon_clockwise(p));
|
|
}
|
|
|
|
TEST_CASE("[Geometry2D] Line intersection") {
|
|
Vector2 r;
|
|
CHECK(Geometry2D::line_intersects_line(Vector2(2, 0), Vector2(0, 1), Vector2(0, 2), Vector2(1, 0), r));
|
|
CHECK(r.is_equal_approx(Vector2(2, 2)));
|
|
|
|
CHECK(Geometry2D::line_intersects_line(Vector2(-1, 1), Vector2(1, -1), Vector2(4, 1), Vector2(-1, -1), r));
|
|
CHECK(r.is_equal_approx(Vector2(1.5, -1.5)));
|
|
|
|
CHECK(Geometry2D::line_intersects_line(Vector2(-1, 0), Vector2(-1, -1), Vector2(1, 0), Vector2(1, -1), r));
|
|
CHECK(r.is_equal_approx(Vector2(0, 1)));
|
|
|
|
CHECK_FALSE_MESSAGE(
|
|
Geometry2D::line_intersects_line(Vector2(-1, 1), Vector2(1, -1), Vector2(0, 1), Vector2(1, -1), r),
|
|
"Parallel lines should not intersect.");
|
|
}
|
|
|
|
TEST_CASE("[Geometry2D] Segment intersection") {
|
|
Vector2 r;
|
|
|
|
CHECK(Geometry2D::segment_intersects_segment(Vector2(-1, 1), Vector2(1, -1), Vector2(1, 1), Vector2(-1, -1), &r));
|
|
CHECK(r.is_equal_approx(Vector2(0, 0)));
|
|
|
|
CHECK_FALSE(Geometry2D::segment_intersects_segment(Vector2(-1, 1), Vector2(1, -1), Vector2(1, 1), Vector2(0.1, 0.1), &r));
|
|
CHECK_FALSE(Geometry2D::segment_intersects_segment(Vector2(-1, 1), Vector2(1, -1), Vector2(0.1, 0.1), Vector2(1, 1), &r));
|
|
|
|
CHECK_FALSE_MESSAGE(
|
|
Geometry2D::segment_intersects_segment(Vector2(-1, 1), Vector2(1, -1), Vector2(0, 1), Vector2(2, -1), &r),
|
|
"Parallel segments should not intersect.");
|
|
|
|
CHECK_FALSE_MESSAGE(
|
|
Geometry2D::segment_intersects_segment(Vector2(1, 2), Vector2(3, 2), Vector2(0, 2), Vector2(-2, 2), &r),
|
|
"Non-overlapping collinear segments should not intersect.");
|
|
|
|
CHECK_MESSAGE(
|
|
Geometry2D::segment_intersects_segment(Vector2(0, 0), Vector2(0, 1), Vector2(0, 0), Vector2(1, 0), &r),
|
|
"Touching segments should intersect.");
|
|
CHECK(r.is_equal_approx(Vector2(0, 0)));
|
|
|
|
CHECK_MESSAGE(
|
|
Geometry2D::segment_intersects_segment(Vector2(0, 1), Vector2(0, 0), Vector2(0, 0), Vector2(1, 0), &r),
|
|
"Touching segments should intersect.");
|
|
CHECK(r.is_equal_approx(Vector2(0, 0)));
|
|
}
|
|
|
|
TEST_CASE("[Geometry2D] Segment intersection with circle") {
|
|
real_t minus_one = -1.0;
|
|
real_t zero = 0.0;
|
|
real_t one_quarter = 0.25;
|
|
real_t three_quarters = 0.75;
|
|
real_t one = 1.0;
|
|
|
|
CHECK_MESSAGE(
|
|
Geometry2D::segment_intersects_circle(Vector2(0, 0), Vector2(4, 0), Vector2(0, 0), 1.0) == doctest::Approx(one_quarter),
|
|
"Segment from inside to outside of circle should intersect it.");
|
|
CHECK_MESSAGE(
|
|
Geometry2D::segment_intersects_circle(Vector2(4, 0), Vector2(0, 0), Vector2(0, 0), 1.0) == doctest::Approx(three_quarters),
|
|
"Segment from outside to inside of circle should intersect it.");
|
|
|
|
CHECK_MESSAGE(
|
|
Geometry2D::segment_intersects_circle(Vector2(-2, 0), Vector2(2, 0), Vector2(0, 0), 1.0) == doctest::Approx(one_quarter),
|
|
"Segment running through circle should intersect it.");
|
|
CHECK_MESSAGE(
|
|
Geometry2D::segment_intersects_circle(Vector2(2, 0), Vector2(-2, 0), Vector2(0, 0), 1.0) == doctest::Approx(one_quarter),
|
|
"Segment running through circle should intersect it.");
|
|
|
|
CHECK_MESSAGE(
|
|
Geometry2D::segment_intersects_circle(Vector2(0, 0), Vector2(1, 0), Vector2(0, 0), 1.0) == doctest::Approx(one),
|
|
"Segment starting inside the circle and ending on the circle should intersect it");
|
|
CHECK_MESSAGE(
|
|
Geometry2D::segment_intersects_circle(Vector2(1, 0), Vector2(0, 0), Vector2(0, 0), 1.0) == doctest::Approx(zero),
|
|
"Segment starting on the circle and going inwards should intersect it");
|
|
CHECK_MESSAGE(
|
|
Geometry2D::segment_intersects_circle(Vector2(1, 0), Vector2(2, 0), Vector2(0, 0), 1.0) == doctest::Approx(zero),
|
|
"Segment starting on the circle and going outwards should intersect it");
|
|
CHECK_MESSAGE(
|
|
Geometry2D::segment_intersects_circle(Vector2(2, 0), Vector2(1, 0), Vector2(0, 0), 1.0) == doctest::Approx(one),
|
|
"Segment starting outside the circle and ending on the circle intersect it");
|
|
|
|
CHECK_MESSAGE(
|
|
Geometry2D::segment_intersects_circle(Vector2(-1, 0), Vector2(1, 0), Vector2(0, 0), 2.0) == doctest::Approx(minus_one),
|
|
"Segment completely within the circle should not intersect it");
|
|
CHECK_MESSAGE(
|
|
Geometry2D::segment_intersects_circle(Vector2(1, 0), Vector2(-1, 0), Vector2(0, 0), 2.0) == doctest::Approx(minus_one),
|
|
"Segment completely within the circle should not intersect it");
|
|
CHECK_MESSAGE(
|
|
Geometry2D::segment_intersects_circle(Vector2(2, 0), Vector2(3, 0), Vector2(0, 0), 1.0) == doctest::Approx(minus_one),
|
|
"Segment completely outside the circle should not intersect it");
|
|
CHECK_MESSAGE(
|
|
Geometry2D::segment_intersects_circle(Vector2(3, 0), Vector2(2, 0), Vector2(0, 0), 1.0) == doctest::Approx(minus_one),
|
|
"Segment completely outside the circle should not intersect it");
|
|
}
|
|
|
|
TEST_CASE("[Geometry2D] Segment intersection with polygon") {
|
|
Vector<Point2> a;
|
|
|
|
a.push_back(Point2(-2, 2));
|
|
a.push_back(Point2(3, 4));
|
|
a.push_back(Point2(1, 1));
|
|
a.push_back(Point2(2, -2));
|
|
a.push_back(Point2(-1, -1));
|
|
|
|
CHECK_MESSAGE(
|
|
Geometry2D::is_segment_intersecting_polygon(Vector2(0, 2), Vector2(2, 2), a),
|
|
"Segment from inside to outside of polygon should intersect it.");
|
|
CHECK_MESSAGE(
|
|
Geometry2D::is_segment_intersecting_polygon(Vector2(2, 2), Vector2(0, 2), a),
|
|
"Segment from outside to inside of polygon should intersect it.");
|
|
|
|
CHECK_MESSAGE(
|
|
Geometry2D::is_segment_intersecting_polygon(Vector2(2, 4), Vector2(3, 3), a),
|
|
"Segment running through polygon should intersect it.");
|
|
CHECK_MESSAGE(
|
|
Geometry2D::is_segment_intersecting_polygon(Vector2(3, 3), Vector2(2, 4), a),
|
|
"Segment running through polygon should intersect it.");
|
|
|
|
CHECK_MESSAGE(
|
|
Geometry2D::is_segment_intersecting_polygon(Vector2(0, 0), Vector2(1, 1), a),
|
|
"Segment starting inside the polygon and ending on the polygon should intersect it");
|
|
CHECK_MESSAGE(
|
|
Geometry2D::is_segment_intersecting_polygon(Vector2(1, 1), Vector2(0, 0), a),
|
|
"Segment starting on the polygon and going inwards should intersect it");
|
|
CHECK_MESSAGE(
|
|
Geometry2D::is_segment_intersecting_polygon(Vector2(-2, 2), Vector2(-2, -1), a),
|
|
"Segment starting on the polygon and going outwards should intersect it");
|
|
CHECK_MESSAGE(
|
|
Geometry2D::is_segment_intersecting_polygon(Vector2(-2, 1), Vector2(-2, 2), a),
|
|
"Segment starting outside the polygon and ending on the polygon intersect it");
|
|
|
|
CHECK_FALSE_MESSAGE(
|
|
Geometry2D::is_segment_intersecting_polygon(Vector2(-1, 2), Vector2(1, -1), a),
|
|
"Segment completely within the polygon should not intersect it");
|
|
CHECK_FALSE_MESSAGE(
|
|
Geometry2D::is_segment_intersecting_polygon(Vector2(1, -1), Vector2(-1, 2), a),
|
|
"Segment completely within the polygon should not intersect it");
|
|
CHECK_FALSE_MESSAGE(
|
|
Geometry2D::is_segment_intersecting_polygon(Vector2(2, 2), Vector2(2, -1), a),
|
|
"Segment completely outside the polygon should not intersect it");
|
|
CHECK_FALSE_MESSAGE(
|
|
Geometry2D::is_segment_intersecting_polygon(Vector2(2, -1), Vector2(2, 2), a),
|
|
"Segment completely outside the polygon should not intersect it");
|
|
}
|
|
|
|
TEST_CASE("[Geometry2D] Closest point to segment") {
|
|
Vector2 s[] = { Vector2(-4, -4), Vector2(4, 4) };
|
|
CHECK(Geometry2D::get_closest_point_to_segment(Vector2(4.1, 4.1), s).is_equal_approx(Vector2(4, 4)));
|
|
CHECK(Geometry2D::get_closest_point_to_segment(Vector2(-4.1, -4.1), s).is_equal_approx(Vector2(-4, -4)));
|
|
CHECK(Geometry2D::get_closest_point_to_segment(Vector2(-1, 1), s).is_equal_approx(Vector2(0, 0)));
|
|
|
|
Vector2 t[] = { Vector2(1, -2), Vector2(1, -2) };
|
|
CHECK_MESSAGE(
|
|
Geometry2D::get_closest_point_to_segment(Vector2(-3, 4), t).is_equal_approx(Vector2(1, -2)),
|
|
"Line segment is only a single point. This point should be the closest.");
|
|
}
|
|
|
|
TEST_CASE("[Geometry2D] Closest point to uncapped segment") {
|
|
Vector2 s[] = { Vector2(-4, -4), Vector2(4, 4) };
|
|
CHECK(Geometry2D::get_closest_point_to_segment_uncapped(Vector2(-1, 1), s).is_equal_approx(Vector2(0, 0)));
|
|
CHECK(Geometry2D::get_closest_point_to_segment_uncapped(Vector2(-4, -6), s).is_equal_approx(Vector2(-5, -5)));
|
|
CHECK(Geometry2D::get_closest_point_to_segment_uncapped(Vector2(4, 6), s).is_equal_approx(Vector2(5, 5)));
|
|
}
|
|
|
|
TEST_CASE("[Geometry2D] Closest points between segments") {
|
|
Vector2 c1, c2;
|
|
Geometry2D::get_closest_points_between_segments(Vector2(2, 2), Vector2(3, 3), Vector2(4, 4), Vector2(4, 5), c1, c2);
|
|
CHECK(c1.is_equal_approx(Vector2(3, 3)));
|
|
CHECK(c2.is_equal_approx(Vector2(4, 4)));
|
|
|
|
Geometry2D::get_closest_points_between_segments(Vector2(0, 1), Vector2(-2, -1), Vector2(0, 0), Vector2(2, -2), c1, c2);
|
|
CHECK(c1.is_equal_approx(Vector2(-0.5, 0.5)));
|
|
CHECK(c2.is_equal_approx(Vector2(0, 0)));
|
|
|
|
Geometry2D::get_closest_points_between_segments(Vector2(-1, 1), Vector2(1, -1), Vector2(1, 1), Vector2(-1, -1), c1, c2);
|
|
CHECK(c1.is_equal_approx(Vector2(0, 0)));
|
|
CHECK(c2.is_equal_approx(Vector2(0, 0)));
|
|
|
|
Geometry2D::get_closest_points_between_segments(Vector2(-3, 4), Vector2(-3, 4), Vector2(-4, 3), Vector2(-2, 3), c1, c2);
|
|
CHECK_MESSAGE(
|
|
c1.is_equal_approx(Vector2(-3, 4)),
|
|
"1st line segment is only a point, this point should be the closest point to the 2nd line segment.");
|
|
CHECK_MESSAGE(
|
|
c2.is_equal_approx(Vector2(-3, 3)),
|
|
"1st line segment is only a point, this should not matter when determining the closest point on the 2nd line segment.");
|
|
|
|
Geometry2D::get_closest_points_between_segments(Vector2(-4, 3), Vector2(-2, 3), Vector2(-3, 4), Vector2(-3, 4), c1, c2);
|
|
CHECK_MESSAGE(
|
|
c1.is_equal_approx(Vector2(-3, 3)),
|
|
"2nd line segment is only a point, this should not matter when determining the closest point on the 1st line segment.");
|
|
CHECK_MESSAGE(
|
|
c2.is_equal_approx(Vector2(-3, 4)),
|
|
"2nd line segment is only a point, this point should be the closest point to the 1st line segment.");
|
|
|
|
Geometry2D::get_closest_points_between_segments(Vector2(5, -4), Vector2(5, -4), Vector2(-2, 1), Vector2(-2, 1), c1, c2);
|
|
CHECK_MESSAGE(
|
|
c1.is_equal_approx(Vector2(5, -4)),
|
|
"Both line segments are only a point. On the 1st line segment, that point should be the closest point to the 2nd line segment.");
|
|
CHECK_MESSAGE(
|
|
c2.is_equal_approx(Vector2(-2, 1)),
|
|
"Both line segments are only a point. On the 2nd line segment, that point should be the closest point to the 1st line segment.");
|
|
}
|
|
|
|
TEST_CASE("[Geometry2D] Make atlas") {
|
|
Vector<Point2i> result;
|
|
Size2i size;
|
|
|
|
Vector<Size2i> r;
|
|
r.push_back(Size2i(2, 2));
|
|
Geometry2D::make_atlas(r, result, size);
|
|
CHECK(size == Size2i(2, 2));
|
|
CHECK(result.size() == r.size());
|
|
|
|
r.clear();
|
|
result.clear();
|
|
r.push_back(Size2i(1, 2));
|
|
r.push_back(Size2i(3, 4));
|
|
r.push_back(Size2i(5, 6));
|
|
r.push_back(Size2i(7, 8));
|
|
Geometry2D::make_atlas(r, result, size);
|
|
CHECK(result.size() == r.size());
|
|
}
|
|
|
|
TEST_CASE("[Geometry2D] Polygon intersection") {
|
|
Vector<Point2> a;
|
|
Vector<Point2> b;
|
|
Vector<Vector<Point2>> r;
|
|
|
|
a.push_back(Point2(30, 60));
|
|
a.push_back(Point2(70, 5));
|
|
a.push_back(Point2(200, 40));
|
|
a.push_back(Point2(80, 200));
|
|
|
|
SUBCASE("[Geometry2D] Both polygons are empty") {
|
|
r = Geometry2D::intersect_polygons(Vector<Point2>(), Vector<Point2>());
|
|
CHECK_MESSAGE(r.is_empty(), "Both polygons are empty. The intersection should also be empty.");
|
|
}
|
|
|
|
SUBCASE("[Geometry2D] One polygon is empty") {
|
|
r = Geometry2D::intersect_polygons(a, b);
|
|
REQUIRE_MESSAGE(r.is_empty(), "One polygon is empty. The intersection should also be empty.");
|
|
}
|
|
|
|
SUBCASE("[Geometry2D] Basic intersection") {
|
|
b.push_back(Point2(200, 300));
|
|
b.push_back(Point2(90, 200));
|
|
b.push_back(Point2(50, 100));
|
|
b.push_back(Point2(200, 90));
|
|
r = Geometry2D::intersect_polygons(a, b);
|
|
REQUIRE_MESSAGE(r.size() == 1, "The polygons should intersect each other with 1 resulting intersection polygon.");
|
|
REQUIRE_MESSAGE(r[0].size() == 3, "The resulting intersection polygon should have 3 vertices.");
|
|
CHECK(r[0][0].is_equal_approx(Point2(86.52174, 191.30436)));
|
|
CHECK(r[0][1].is_equal_approx(Point2(50, 100)));
|
|
CHECK(r[0][2].is_equal_approx(Point2(160.52632, 92.63157)));
|
|
}
|
|
|
|
SUBCASE("[Geometry2D] Intersection with one polygon being completely inside the other polygon") {
|
|
b.push_back(Point2(80, 100));
|
|
b.push_back(Point2(50, 50));
|
|
b.push_back(Point2(150, 50));
|
|
r = Geometry2D::intersect_polygons(a, b);
|
|
REQUIRE_MESSAGE(r.size() == 1, "The polygons should intersect each other with 1 resulting intersection polygon.");
|
|
REQUIRE_MESSAGE(r[0].size() == 3, "The resulting intersection polygon should have 3 vertices.");
|
|
CHECK(r[0][0].is_equal_approx(b[0]));
|
|
CHECK(r[0][1].is_equal_approx(b[1]));
|
|
CHECK(r[0][2].is_equal_approx(b[2]));
|
|
}
|
|
|
|
SUBCASE("[Geometry2D] No intersection with 2 non-empty polygons") {
|
|
b.push_back(Point2(150, 150));
|
|
b.push_back(Point2(250, 100));
|
|
b.push_back(Point2(300, 200));
|
|
r = Geometry2D::intersect_polygons(a, b);
|
|
REQUIRE_MESSAGE(r.is_empty(), "The polygons should not intersect each other.");
|
|
}
|
|
|
|
SUBCASE("[Geometry2D] Intersection with 2 resulting polygons") {
|
|
a.clear();
|
|
a.push_back(Point2(70, 5));
|
|
a.push_back(Point2(140, 7));
|
|
a.push_back(Point2(100, 52));
|
|
a.push_back(Point2(170, 50));
|
|
a.push_back(Point2(60, 125));
|
|
b.push_back(Point2(70, 105));
|
|
b.push_back(Point2(115, 55));
|
|
b.push_back(Point2(90, 15));
|
|
b.push_back(Point2(160, 50));
|
|
r = Geometry2D::intersect_polygons(a, b);
|
|
REQUIRE_MESSAGE(r.size() == 2, "The polygons should intersect each other with 2 resulting intersection polygons.");
|
|
REQUIRE_MESSAGE(r[0].size() == 4, "The resulting intersection polygon should have 4 vertices.");
|
|
CHECK(r[0][0].is_equal_approx(Point2(70, 105)));
|
|
CHECK(r[0][1].is_equal_approx(Point2(115, 55)));
|
|
CHECK(r[0][2].is_equal_approx(Point2(112.894737, 51.63158)));
|
|
CHECK(r[0][3].is_equal_approx(Point2(159.509537, 50.299728)));
|
|
|
|
REQUIRE_MESSAGE(r[1].size() == 3, "The intersection polygon should have 3 vertices.");
|
|
CHECK(r[1][0].is_equal_approx(Point2(119.692307, 29.846149)));
|
|
CHECK(r[1][1].is_equal_approx(Point2(107.706421, 43.33028)));
|
|
CHECK(r[1][2].is_equal_approx(Point2(90, 15)));
|
|
}
|
|
}
|
|
|
|
TEST_CASE("[Geometry2D] Merge polygons") {
|
|
Vector<Point2> a;
|
|
Vector<Point2> b;
|
|
Vector<Vector<Point2>> r;
|
|
|
|
a.push_back(Point2(225, 180));
|
|
a.push_back(Point2(160, 230));
|
|
a.push_back(Point2(20, 212));
|
|
a.push_back(Point2(50, 115));
|
|
|
|
SUBCASE("[Geometry2D] Both polygons are empty") {
|
|
r = Geometry2D::merge_polygons(Vector<Point2>(), Vector<Point2>());
|
|
REQUIRE_MESSAGE(r.is_empty(), "Both polygons are empty. The union should also be empty.");
|
|
}
|
|
|
|
SUBCASE("[Geometry2D] One polygon is empty") {
|
|
r = Geometry2D::merge_polygons(a, b);
|
|
REQUIRE_MESSAGE(r.size() == 1, "One polygon is non-empty. There should be 1 resulting merged polygon.");
|
|
REQUIRE_MESSAGE(r[0].size() == 4, "The resulting merged polygon should have 4 vertices.");
|
|
CHECK(r[0][0].is_equal_approx(a[0]));
|
|
CHECK(r[0][1].is_equal_approx(a[1]));
|
|
CHECK(r[0][2].is_equal_approx(a[2]));
|
|
CHECK(r[0][3].is_equal_approx(a[3]));
|
|
}
|
|
|
|
SUBCASE("[Geometry2D] Basic merge with 2 polygons") {
|
|
b.push_back(Point2(180, 190));
|
|
b.push_back(Point2(60, 140));
|
|
b.push_back(Point2(160, 80));
|
|
r = Geometry2D::merge_polygons(a, b);
|
|
REQUIRE_MESSAGE(r.size() == 1, "The merged polygons should result in 1 polygon.");
|
|
REQUIRE_MESSAGE(r[0].size() == 7, "The resulting merged polygon should have 7 vertices.");
|
|
CHECK(r[0][0].is_equal_approx(Point2(174.791077, 161.350967)));
|
|
CHECK(r[0][1].is_equal_approx(Point2(225, 180)));
|
|
CHECK(r[0][2].is_equal_approx(Point2(160, 230)));
|
|
CHECK(r[0][3].is_equal_approx(Point2(20, 212)));
|
|
CHECK(r[0][4].is_equal_approx(Point2(50, 115)));
|
|
CHECK(r[0][5].is_equal_approx(Point2(81.911758, 126.852943)));
|
|
CHECK(r[0][6].is_equal_approx(Point2(160, 80)));
|
|
}
|
|
|
|
SUBCASE("[Geometry2D] Merge with 2 resulting merged polygons (outline and hole)") {
|
|
b.push_back(Point2(180, 190));
|
|
b.push_back(Point2(140, 125));
|
|
b.push_back(Point2(60, 140));
|
|
b.push_back(Point2(160, 80));
|
|
r = Geometry2D::merge_polygons(a, b);
|
|
REQUIRE_MESSAGE(r.size() == 2, "The merged polygons should result in 2 polygons.");
|
|
|
|
REQUIRE_MESSAGE(!Geometry2D::is_polygon_clockwise(r[0]), "The merged polygon (outline) should be counter-clockwise.");
|
|
REQUIRE_MESSAGE(r[0].size() == 7, "The resulting merged polygon (outline) should have 7 vertices.");
|
|
CHECK(r[0][0].is_equal_approx(Point2(174.791077, 161.350967)));
|
|
CHECK(r[0][1].is_equal_approx(Point2(225, 180)));
|
|
CHECK(r[0][2].is_equal_approx(Point2(160, 230)));
|
|
CHECK(r[0][3].is_equal_approx(Point2(20, 212)));
|
|
CHECK(r[0][4].is_equal_approx(Point2(50, 115)));
|
|
CHECK(r[0][5].is_equal_approx(Point2(81.911758, 126.852943)));
|
|
CHECK(r[0][6].is_equal_approx(Point2(160, 80)));
|
|
|
|
REQUIRE_MESSAGE(Geometry2D::is_polygon_clockwise(r[1]), "The resulting merged polygon (hole) should be clockwise.");
|
|
REQUIRE_MESSAGE(r[1].size() == 3, "The resulting merged polygon (hole) should have 3 vertices.");
|
|
CHECK(r[1][0].is_equal_approx(Point2(98.083069, 132.859421)));
|
|
CHECK(r[1][1].is_equal_approx(Point2(158.689453, 155.370377)));
|
|
CHECK(r[1][2].is_equal_approx(Point2(140, 125)));
|
|
}
|
|
}
|
|
|
|
TEST_CASE("[Geometry2D] Clip polygons") {
|
|
Vector<Point2> a;
|
|
Vector<Point2> b;
|
|
Vector<Vector<Point2>> r;
|
|
|
|
a.push_back(Point2(225, 180));
|
|
a.push_back(Point2(160, 230));
|
|
a.push_back(Point2(20, 212));
|
|
a.push_back(Point2(50, 115));
|
|
|
|
SUBCASE("[Geometry2D] Both polygons are empty") {
|
|
r = Geometry2D::clip_polygons(Vector<Point2>(), Vector<Point2>());
|
|
CHECK_MESSAGE(r.is_empty(), "Both polygons are empty. The clip should also be empty.");
|
|
}
|
|
|
|
SUBCASE("[Geometry2D] Basic clip with one result polygon") {
|
|
b.push_back(Point2(250, 170));
|
|
b.push_back(Point2(175, 270));
|
|
b.push_back(Point2(120, 260));
|
|
b.push_back(Point2(25, 80));
|
|
r = Geometry2D::clip_polygons(a, b);
|
|
REQUIRE_MESSAGE(r.size() == 1, "The clipped polygons should result in 1 polygon.");
|
|
REQUIRE_MESSAGE(r[0].size() == 3, "The resulting clipped polygon should have 3 vertices.");
|
|
CHECK(r[0][0].is_equal_approx(Point2(100.102173, 222.298843)));
|
|
CHECK(r[0][1].is_equal_approx(Point2(20, 212)));
|
|
CHECK(r[0][2].is_equal_approx(Point2(47.588089, 122.798492)));
|
|
}
|
|
|
|
SUBCASE("[Geometry2D] Polygon b completely overlaps polygon a") {
|
|
b.push_back(Point2(250, 170));
|
|
b.push_back(Point2(175, 270));
|
|
b.push_back(Point2(10, 210));
|
|
b.push_back(Point2(55, 80));
|
|
r = Geometry2D::clip_polygons(a, b);
|
|
CHECK_MESSAGE(r.is_empty(), "Polygon 'b' completely overlaps polygon 'a'. This should result in no clipped polygons.");
|
|
}
|
|
|
|
SUBCASE("[Geometry2D] Polygon a completely overlaps polygon b") {
|
|
b.push_back(Point2(150, 200));
|
|
b.push_back(Point2(65, 190));
|
|
b.push_back(Point2(80, 140));
|
|
r = Geometry2D::clip_polygons(a, b);
|
|
REQUIRE_MESSAGE(r.size() == 2, "Polygon 'a' completely overlaps polygon 'b'. This should result in 2 clipped polygons.");
|
|
REQUIRE_MESSAGE(r[0].size() == 4, "The resulting clipped polygon should have 4 vertices.");
|
|
REQUIRE_MESSAGE(!Geometry2D::is_polygon_clockwise(r[0]), "The resulting clipped polygon (outline) should be counter-clockwise.");
|
|
CHECK(r[0][0].is_equal_approx(a[0]));
|
|
CHECK(r[0][1].is_equal_approx(a[1]));
|
|
CHECK(r[0][2].is_equal_approx(a[2]));
|
|
CHECK(r[0][3].is_equal_approx(a[3]));
|
|
REQUIRE_MESSAGE(r[1].size() == 3, "The resulting clipped polygon should have 3 vertices.");
|
|
REQUIRE_MESSAGE(Geometry2D::is_polygon_clockwise(r[1]), "The resulting clipped polygon (hole) should be clockwise.");
|
|
CHECK(r[1][0].is_equal_approx(b[1]));
|
|
CHECK(r[1][1].is_equal_approx(b[0]));
|
|
CHECK(r[1][2].is_equal_approx(b[2]));
|
|
}
|
|
}
|
|
|
|
TEST_CASE("[Geometry2D] Exclude polygons") {
|
|
Vector<Point2> a;
|
|
Vector<Point2> b;
|
|
Vector<Vector<Point2>> r;
|
|
|
|
a.push_back(Point2(225, 180));
|
|
a.push_back(Point2(160, 230));
|
|
a.push_back(Point2(20, 212));
|
|
a.push_back(Point2(50, 115));
|
|
|
|
SUBCASE("[Geometry2D] Both polygons are empty") {
|
|
r = Geometry2D::exclude_polygons(Vector<Point2>(), Vector<Point2>());
|
|
CHECK_MESSAGE(r.is_empty(), "Both polygons are empty. The excluded polygon should also be empty.");
|
|
}
|
|
|
|
SUBCASE("[Geometry2D] One polygon is empty") {
|
|
r = Geometry2D::exclude_polygons(a, b);
|
|
REQUIRE_MESSAGE(r.size() == 1, "One polygon is non-empty. There should be 1 resulting excluded polygon.");
|
|
REQUIRE_MESSAGE(r[0].size() == 4, "The resulting excluded polygon should have 4 vertices.");
|
|
CHECK(r[0][0].is_equal_approx(a[0]));
|
|
CHECK(r[0][1].is_equal_approx(a[1]));
|
|
CHECK(r[0][2].is_equal_approx(a[2]));
|
|
CHECK(r[0][3].is_equal_approx(a[3]));
|
|
}
|
|
|
|
SUBCASE("[Geometry2D] Exclude with 2 resulting polygons (outline and hole)") {
|
|
b.push_back(Point2(140, 160));
|
|
b.push_back(Point2(150, 220));
|
|
b.push_back(Point2(40, 200));
|
|
b.push_back(Point2(60, 140));
|
|
r = Geometry2D::exclude_polygons(a, b);
|
|
REQUIRE_MESSAGE(r.size() == 2, "There should be 2 resulting excluded polygons (outline and hole).");
|
|
REQUIRE_MESSAGE(r[0].size() == 4, "The resulting excluded polygon should have 4 vertices.");
|
|
REQUIRE_MESSAGE(!Geometry2D::is_polygon_clockwise(r[0]), "The resulting excluded polygon (outline) should be counter-clockwise.");
|
|
CHECK(r[0][0].is_equal_approx(a[0]));
|
|
CHECK(r[0][1].is_equal_approx(a[1]));
|
|
CHECK(r[0][2].is_equal_approx(a[2]));
|
|
CHECK(r[0][3].is_equal_approx(a[3]));
|
|
REQUIRE_MESSAGE(r[1].size() == 4, "The resulting excluded polygon should have 4 vertices.");
|
|
REQUIRE_MESSAGE(Geometry2D::is_polygon_clockwise(r[1]), "The resulting excluded polygon (hole) should be clockwise.");
|
|
CHECK(r[1][0].is_equal_approx(Point2(40, 200)));
|
|
CHECK(r[1][1].is_equal_approx(Point2(150, 220)));
|
|
CHECK(r[1][2].is_equal_approx(Point2(140, 160)));
|
|
CHECK(r[1][3].is_equal_approx(Point2(60, 140)));
|
|
}
|
|
}
|
|
|
|
TEST_CASE("[Geometry2D] Intersect polyline with polygon") {
|
|
Vector<Vector2> l;
|
|
Vector<Vector2> p;
|
|
Vector<Vector<Point2>> r;
|
|
|
|
l.push_back(Vector2(100, 90));
|
|
l.push_back(Vector2(120, 250));
|
|
|
|
p.push_back(Vector2(225, 180));
|
|
p.push_back(Vector2(160, 230));
|
|
p.push_back(Vector2(20, 212));
|
|
p.push_back(Vector2(50, 115));
|
|
|
|
SUBCASE("[Geometry2D] Both line and polygon are empty") {
|
|
r = Geometry2D::intersect_polyline_with_polygon(Vector<Vector2>(), Vector<Vector2>());
|
|
CHECK_MESSAGE(r.is_empty(), "Both line and polygon are empty. The intersection line should also be empty.");
|
|
}
|
|
|
|
SUBCASE("[Geometry2D] Line is non-empty and polygon is empty") {
|
|
r = Geometry2D::intersect_polyline_with_polygon(l, Vector<Vector2>());
|
|
CHECK_MESSAGE(r.is_empty(), "The polygon is empty while the line is non-empty. The intersection line should be empty.");
|
|
}
|
|
|
|
SUBCASE("[Geometry2D] Basic intersection with 1 resulting intersection line") {
|
|
r = Geometry2D::intersect_polyline_with_polygon(l, p);
|
|
REQUIRE_MESSAGE(r.size() == 1, "There should be 1 resulting intersection line.");
|
|
REQUIRE_MESSAGE(r[0].size() == 2, "The resulting intersection line should have 2 vertices.");
|
|
CHECK(r[0][0].is_equal_approx(Vector2(105.711609, 135.692886)));
|
|
CHECK(r[0][1].is_equal_approx(Vector2(116.805809, 224.446457)));
|
|
}
|
|
|
|
SUBCASE("[Geometry2D] Complex intersection with 2 resulting intersection lines") {
|
|
l.clear();
|
|
l.push_back(Vector2(100, 90));
|
|
l.push_back(Vector2(190, 255));
|
|
l.push_back(Vector2(135, 260));
|
|
l.push_back(Vector2(57, 200));
|
|
l.push_back(Vector2(50, 170));
|
|
l.push_back(Vector2(15, 155));
|
|
r = Geometry2D::intersect_polyline_with_polygon(l, p);
|
|
REQUIRE_MESSAGE(r.size() == 2, "There should be 2 resulting intersection lines.");
|
|
REQUIRE_MESSAGE(r[0].size() == 2, "The resulting intersection line should have 2 vertices.");
|
|
CHECK(r[0][0].is_equal_approx(Vector2(129.804565, 144.641693)));
|
|
CHECK(r[0][1].is_equal_approx(Vector2(171.527084, 221.132996)));
|
|
REQUIRE_MESSAGE(r[1].size() == 4, "The resulting intersection line should have 4 vertices.");
|
|
CHECK(r[1][0].is_equal_approx(Vector2(83.15609, 220.120087)));
|
|
CHECK(r[1][1].is_equal_approx(Vector2(57, 200)));
|
|
CHECK(r[1][2].is_equal_approx(Vector2(50, 170)));
|
|
CHECK(r[1][3].is_equal_approx(Vector2(34.980492, 163.563065)));
|
|
}
|
|
}
|
|
|
|
TEST_CASE("[Geometry2D] Clip polyline with polygon") {
|
|
Vector<Vector2> l;
|
|
Vector<Vector2> p;
|
|
Vector<Vector<Point2>> r;
|
|
|
|
l.push_back(Vector2(70, 140));
|
|
l.push_back(Vector2(160, 320));
|
|
|
|
p.push_back(Vector2(225, 180));
|
|
p.push_back(Vector2(160, 230));
|
|
p.push_back(Vector2(20, 212));
|
|
p.push_back(Vector2(50, 115));
|
|
|
|
SUBCASE("[Geometry2D] Both line and polygon are empty") {
|
|
r = Geometry2D::clip_polyline_with_polygon(Vector<Vector2>(), Vector<Vector2>());
|
|
CHECK_MESSAGE(r.is_empty(), "Both line and polygon are empty. The clipped line should also be empty.");
|
|
}
|
|
|
|
SUBCASE("[Geometry2D] Polygon is empty and line is non-empty") {
|
|
r = Geometry2D::clip_polyline_with_polygon(l, Vector<Vector2>());
|
|
REQUIRE_MESSAGE(r.size() == 1, "There should be 1 resulting clipped line.");
|
|
REQUIRE_MESSAGE(r[0].size() == 2, "The resulting clipped line should have 2 vertices.");
|
|
CHECK(r[0][0].is_equal_approx(l[0]));
|
|
CHECK(r[0][1].is_equal_approx(l[1]));
|
|
}
|
|
|
|
SUBCASE("[Geometry2D] Basic clip with 1 resulting clipped line") {
|
|
r = Geometry2D::clip_polyline_with_polygon(l, p);
|
|
REQUIRE_MESSAGE(r.size() == 1, "There should be 1 resulting clipped line.");
|
|
REQUIRE_MESSAGE(r[0].size() == 2, "The resulting clipped line should have 2 vertices.");
|
|
CHECK(r[0][0].is_equal_approx(Vector2(111.908401, 223.816803)));
|
|
CHECK(r[0][1].is_equal_approx(Vector2(160, 320)));
|
|
}
|
|
|
|
SUBCASE("[Geometry2D] Complex clip with 2 resulting clipped lines") {
|
|
l.clear();
|
|
l.push_back(Vector2(55, 70));
|
|
l.push_back(Vector2(50, 190));
|
|
l.push_back(Vector2(120, 165));
|
|
l.push_back(Vector2(122, 250));
|
|
l.push_back(Vector2(160, 320));
|
|
r = Geometry2D::clip_polyline_with_polygon(l, p);
|
|
REQUIRE_MESSAGE(r.size() == 2, "There should be 2 resulting clipped lines.");
|
|
REQUIRE_MESSAGE(r[0].size() == 3, "The resulting clipped line should have 3 vertices.");
|
|
CHECK(r[0][0].is_equal_approx(Vector2(160, 320)));
|
|
CHECK(r[0][1].is_equal_approx(Vector2(122, 250)));
|
|
CHECK(r[0][2].is_equal_approx(Vector2(121.412682, 225.038757)));
|
|
REQUIRE_MESSAGE(r[1].size() == 2, "The resulting clipped line should have 2 vertices.");
|
|
CHECK(r[1][0].is_equal_approx(Vector2(53.07737, 116.143021)));
|
|
CHECK(r[1][1].is_equal_approx(Vector2(55, 70)));
|
|
}
|
|
}
|
|
|
|
TEST_CASE("[Geometry2D] Convex hull") {
|
|
Vector<Point2> a;
|
|
Vector<Point2> r;
|
|
|
|
a.push_back(Point2(-4, -8));
|
|
a.push_back(Point2(-10, -4));
|
|
a.push_back(Point2(8, 2));
|
|
a.push_back(Point2(-6, 10));
|
|
a.push_back(Point2(-12, 4));
|
|
a.push_back(Point2(10, -8));
|
|
a.push_back(Point2(4, 8));
|
|
|
|
SUBCASE("[Geometry2D] No points") {
|
|
r = Geometry2D::convex_hull(Vector<Vector2>());
|
|
|
|
CHECK_MESSAGE(r.is_empty(), "The convex hull should be empty if there are no input points.");
|
|
}
|
|
|
|
SUBCASE("[Geometry2D] Single point") {
|
|
Vector<Point2> b;
|
|
b.push_back(Point2(4, -3));
|
|
|
|
r = Geometry2D::convex_hull(b);
|
|
REQUIRE_MESSAGE(r.size() == 1, "Convex hull should contain 1 point.");
|
|
CHECK(r[0].is_equal_approx(b[0]));
|
|
}
|
|
|
|
SUBCASE("[Geometry2D] All points form the convex hull") {
|
|
r = Geometry2D::convex_hull(a);
|
|
REQUIRE_MESSAGE(r.size() == 8, "Convex hull should contain 8 points.");
|
|
CHECK(r[0].is_equal_approx(Point2(-12, 4)));
|
|
CHECK(r[1].is_equal_approx(Point2(-10, -4)));
|
|
CHECK(r[2].is_equal_approx(Point2(-4, -8)));
|
|
CHECK(r[3].is_equal_approx(Point2(10, -8)));
|
|
CHECK(r[4].is_equal_approx(Point2(8, 2)));
|
|
CHECK(r[5].is_equal_approx(Point2(4, 8)));
|
|
CHECK(r[6].is_equal_approx(Point2(-6, 10)));
|
|
CHECK(r[7].is_equal_approx(Point2(-12, 4)));
|
|
}
|
|
|
|
SUBCASE("[Geometry2D] Add extra points inside original convex hull") {
|
|
a.push_back(Point2(-4, -8));
|
|
a.push_back(Point2(0, 0));
|
|
a.push_back(Point2(0, 8));
|
|
a.push_back(Point2(-10, -3));
|
|
a.push_back(Point2(9, -4));
|
|
a.push_back(Point2(6, 4));
|
|
|
|
r = Geometry2D::convex_hull(a);
|
|
REQUIRE_MESSAGE(r.size() == 8, "Convex hull should contain 8 points.");
|
|
CHECK(r[0].is_equal_approx(Point2(-12, 4)));
|
|
CHECK(r[1].is_equal_approx(Point2(-10, -4)));
|
|
CHECK(r[2].is_equal_approx(Point2(-4, -8)));
|
|
CHECK(r[3].is_equal_approx(Point2(10, -8)));
|
|
CHECK(r[4].is_equal_approx(Point2(8, 2)));
|
|
CHECK(r[5].is_equal_approx(Point2(4, 8)));
|
|
CHECK(r[6].is_equal_approx(Point2(-6, 10)));
|
|
CHECK(r[7].is_equal_approx(Point2(-12, 4)));
|
|
}
|
|
|
|
SUBCASE("[Geometry2D] Add extra points on border of original convex hull") {
|
|
a.push_back(Point2(9, -3));
|
|
a.push_back(Point2(-2, -8));
|
|
|
|
r = Geometry2D::convex_hull(a);
|
|
REQUIRE_MESSAGE(r.size() == 8, "Convex hull should contain 8 points.");
|
|
CHECK(r[0].is_equal_approx(Point2(-12, 4)));
|
|
CHECK(r[1].is_equal_approx(Point2(-10, -4)));
|
|
CHECK(r[2].is_equal_approx(Point2(-4, -8)));
|
|
CHECK(r[3].is_equal_approx(Point2(10, -8)));
|
|
CHECK(r[4].is_equal_approx(Point2(8, 2)));
|
|
CHECK(r[5].is_equal_approx(Point2(4, 8)));
|
|
CHECK(r[6].is_equal_approx(Point2(-6, 10)));
|
|
CHECK(r[7].is_equal_approx(Point2(-12, 4)));
|
|
}
|
|
|
|
SUBCASE("[Geometry2D] Add extra points outside border of original convex hull") {
|
|
a.push_back(Point2(-11, -1));
|
|
a.push_back(Point2(7, 6));
|
|
|
|
r = Geometry2D::convex_hull(a);
|
|
REQUIRE_MESSAGE(r.size() == 10, "Convex hull should contain 10 points.");
|
|
CHECK(r[0].is_equal_approx(Point2(-12, 4)));
|
|
CHECK(r[1].is_equal_approx(Point2(-11, -1)));
|
|
CHECK(r[2].is_equal_approx(Point2(-10, -4)));
|
|
CHECK(r[3].is_equal_approx(Point2(-4, -8)));
|
|
CHECK(r[4].is_equal_approx(Point2(10, -8)));
|
|
CHECK(r[5].is_equal_approx(Point2(8, 2)));
|
|
CHECK(r[6].is_equal_approx(Point2(7, 6)));
|
|
CHECK(r[7].is_equal_approx(Point2(4, 8)));
|
|
CHECK(r[8].is_equal_approx(Point2(-6, 10)));
|
|
CHECK(r[9].is_equal_approx(Point2(-12, 4)));
|
|
}
|
|
}
|
|
|
|
TEST_CASE("[Geometry2D] Bresenham line") {
|
|
Vector<Vector2i> r;
|
|
|
|
SUBCASE("[Geometry2D] Single point") {
|
|
r = Geometry2D::bresenham_line(Point2i(0, 0), Point2i(0, 0));
|
|
|
|
REQUIRE_MESSAGE(r.size() == 1, "The Bresenham line should contain exactly one point.");
|
|
CHECK(r[0] == Vector2i(0, 0));
|
|
}
|
|
|
|
SUBCASE("[Geometry2D] Line parallel to x-axis") {
|
|
r = Geometry2D::bresenham_line(Point2i(1, 2), Point2i(5, 2));
|
|
|
|
REQUIRE_MESSAGE(r.size() == 5, "The Bresenham line should contain exactly five points.");
|
|
CHECK(r[0] == Vector2i(1, 2));
|
|
CHECK(r[1] == Vector2i(2, 2));
|
|
CHECK(r[2] == Vector2i(3, 2));
|
|
CHECK(r[3] == Vector2i(4, 2));
|
|
CHECK(r[4] == Vector2i(5, 2));
|
|
}
|
|
|
|
SUBCASE("[Geometry2D] 45 degree line from the origin") {
|
|
r = Geometry2D::bresenham_line(Point2i(0, 0), Point2i(4, 4));
|
|
|
|
REQUIRE_MESSAGE(r.size() == 5, "The Bresenham line should contain exactly five points.");
|
|
CHECK(r[0] == Vector2i(0, 0));
|
|
CHECK(r[1] == Vector2i(1, 1));
|
|
CHECK(r[2] == Vector2i(2, 2));
|
|
CHECK(r[3] == Vector2i(3, 3));
|
|
CHECK(r[4] == Vector2i(4, 4));
|
|
}
|
|
|
|
SUBCASE("[Geometry2D] Sloped line going up one unit") {
|
|
r = Geometry2D::bresenham_line(Point2i(0, 0), Point2i(4, 1));
|
|
|
|
REQUIRE_MESSAGE(r.size() == 5, "The Bresenham line should contain exactly five points.");
|
|
CHECK(r[0] == Vector2i(0, 0));
|
|
CHECK(r[1] == Vector2i(1, 0));
|
|
CHECK(r[2] == Vector2i(2, 0));
|
|
CHECK(r[3] == Vector2i(3, 1));
|
|
CHECK(r[4] == Vector2i(4, 1));
|
|
}
|
|
|
|
SUBCASE("[Geometry2D] Sloped line going up two units") {
|
|
r = Geometry2D::bresenham_line(Point2i(0, 0), Point2i(4, 2));
|
|
|
|
REQUIRE_MESSAGE(r.size() == 5, "The Bresenham line should contain exactly five points.");
|
|
CHECK(r[0] == Vector2i(0, 0));
|
|
CHECK(r[1] == Vector2i(1, 0));
|
|
CHECK(r[2] == Vector2i(2, 1));
|
|
CHECK(r[3] == Vector2i(3, 1));
|
|
CHECK(r[4] == Vector2i(4, 2));
|
|
}
|
|
|
|
SUBCASE("[Geometry2D] Long sloped line") {
|
|
r = Geometry2D::bresenham_line(Point2i(0, 0), Point2i(11, 5));
|
|
|
|
REQUIRE_MESSAGE(r.size() == 12, "The Bresenham line should contain exactly twelve points.");
|
|
CHECK(r[0] == Vector2i(0, 0));
|
|
CHECK(r[1] == Vector2i(1, 0));
|
|
CHECK(r[2] == Vector2i(2, 1));
|
|
CHECK(r[3] == Vector2i(3, 1));
|
|
CHECK(r[4] == Vector2i(4, 2));
|
|
CHECK(r[5] == Vector2i(5, 2));
|
|
CHECK(r[6] == Vector2i(6, 3));
|
|
CHECK(r[7] == Vector2i(7, 3));
|
|
CHECK(r[8] == Vector2i(8, 4));
|
|
CHECK(r[9] == Vector2i(9, 4));
|
|
CHECK(r[10] == Vector2i(10, 5));
|
|
CHECK(r[11] == Vector2i(11, 5));
|
|
}
|
|
}
|
|
} // namespace TestGeometry2D
|
|
|
|
#endif // TEST_GEOMETRY_2D_H
|