892 lines
38 KiB
C++
892 lines
38 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;
|
|
// Basis Path Testing suite
|
|
SUBCASE("[Geometry2D] Both segments degenerate to a point") {
|
|
Geometry2D::get_closest_points_between_segments(Vector2(0, 0), Vector2(0, 0), Vector2(0, 0), Vector2(0, 0), c1, c2);
|
|
CHECK(c1.is_equal_approx(Vector2(0, 0)));
|
|
CHECK(c2.is_equal_approx(Vector2(0, 0)));
|
|
}
|
|
|
|
SUBCASE("[Geometry2D] Closest point on second segment trajectory is above [0,1]") {
|
|
Geometry2D::get_closest_points_between_segments(Vector2(50, -25), Vector2(50, -10), Vector2(-50, 10), Vector2(-40, 10), c1, c2);
|
|
CHECK(c1.is_equal_approx(Vector2(50, -10)));
|
|
CHECK(c2.is_equal_approx(Vector2(-40, 10)));
|
|
}
|
|
|
|
SUBCASE("[Geometry2D] Parallel segments") {
|
|
Geometry2D::get_closest_points_between_segments(Vector2(2, 1), Vector2(4, 3), Vector2(2, 3), Vector2(4, 5), c1, c2);
|
|
CHECK(c1.is_equal_approx(Vector2(3, 2)));
|
|
CHECK(c2.is_equal_approx(Vector2(2, 3)));
|
|
}
|
|
|
|
SUBCASE("[Geometry2D] Closest point on second segment trajectory is within [0,1]") {
|
|
Geometry2D::get_closest_points_between_segments(Vector2(2, 4), Vector2(2, 3), Vector2(1, 1), Vector2(4, 4), c1, c2);
|
|
CHECK(c1.is_equal_approx(Vector2(2, 3)));
|
|
CHECK(c2.is_equal_approx(Vector2(2.5, 2.5)));
|
|
}
|
|
|
|
SUBCASE("[Geometry2D] Closest point on second segment trajectory is below [0,1]") {
|
|
Geometry2D::get_closest_points_between_segments(Vector2(-20, -20), Vector2(-10, -40), Vector2(10, 25), Vector2(25, 40), c1, c2);
|
|
CHECK(c1.is_equal_approx(Vector2(-20, -20)));
|
|
CHECK(c2.is_equal_approx(Vector2(10, 25)));
|
|
}
|
|
|
|
SUBCASE("[Geometry2D] Second segment degenerates to a point") {
|
|
Geometry2D::get_closest_points_between_segments(Vector2(1, 2), Vector2(2, 1), Vector2(3, 3), Vector2(3, 3), c1, c2);
|
|
CHECK(c1.is_equal_approx(Vector2(1.5, 1.5)));
|
|
CHECK(c2.is_equal_approx(Vector2(3, 3)));
|
|
}
|
|
|
|
SUBCASE("[Geometry2D] First segment degenerates to a point") {
|
|
Geometry2D::get_closest_points_between_segments(Vector2(1, 1), Vector2(1, 1), Vector2(2, 2), Vector2(4, 4), c1, c2);
|
|
CHECK(c1.is_equal_approx(Vector2(1, 1)));
|
|
CHECK(c2.is_equal_approx(Vector2(2, 2)));
|
|
}
|
|
// End Basis Path Testing suite
|
|
|
|
SUBCASE("[Geometry2D] Segments are equal vectors") {
|
|
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)));
|
|
}
|
|
|
|
SUBCASE("[Geometry2D] Standard case") {
|
|
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)));
|
|
}
|
|
|
|
SUBCASE("[Geometry2D] Segments intersect") {
|
|
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)));
|
|
}
|
|
}
|
|
|
|
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(121.412682, 225.038757)));
|
|
CHECK(r[0][1].is_equal_approx(Vector2(122, 250)));
|
|
CHECK(r[0][2].is_equal_approx(Vector2(160, 320)));
|
|
REQUIRE_MESSAGE(r[1].size() == 2, "The resulting clipped line should have 2 vertices.");
|
|
CHECK(r[1][0].is_equal_approx(Vector2(55, 70)));
|
|
CHECK(r[1][1].is_equal_approx(Vector2(53.07737, 116.143021)));
|
|
}
|
|
}
|
|
|
|
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
|