Index: experimental/Intersection/LineUtilities.cpp |
diff --git a/experimental/Intersection/LineUtilities.cpp b/experimental/Intersection/LineUtilities.cpp |
deleted file mode 100644 |
index 53319c4560823bc168976b125e2083a6bc69931d..0000000000000000000000000000000000000000 |
--- a/experimental/Intersection/LineUtilities.cpp |
+++ /dev/null |
@@ -1,133 +0,0 @@ |
-/* |
- * Copyright 2012 Google Inc. |
- * |
- * Use of this source code is governed by a BSD-style license that can be |
- * found in the LICENSE file. |
- */ |
-#include "CurveIntersection.h" |
-#include "LineUtilities.h" |
- |
-bool implicitLine(const _Line& line, double& slope, double& axisIntercept) { |
- _Point delta; |
- tangent(line, delta); |
- bool moreHorizontal = fabs(delta.x) > fabs(delta.y); |
- if (moreHorizontal) { |
- slope = delta.y / delta.x; |
- axisIntercept = line[0].y - slope * line[0].x; |
- } else { |
- slope = delta.x / delta.y; |
- axisIntercept = line[0].x - slope * line[0].y; |
- } |
- return moreHorizontal; |
-} |
- |
-int reduceOrder(const _Line& line, _Line& reduced) { |
- reduced[0] = line[0]; |
- int different = line[0] != line[1]; |
- reduced[1] = line[different]; |
- return 1 + different; |
-} |
- |
-void sub_divide(const _Line& line, double t1, double t2, _Line& dst) { |
- _Point delta; |
- tangent(line, delta); |
- dst[0].x = line[0].x - t1 * delta.x; |
- dst[0].y = line[0].y - t1 * delta.y; |
- dst[1].x = line[0].x - t2 * delta.x; |
- dst[1].y = line[0].y - t2 * delta.y; |
-} |
- |
-// may have this below somewhere else already: |
-// copying here because I thought it was clever |
- |
-// Copyright 2001, softSurfer (www.softsurfer.com) |
-// This code may be freely used and modified for any purpose |
-// providing that this copyright notice is included with it. |
-// SoftSurfer makes no warranty for this code, and cannot be held |
-// liable for any real or imagined damage resulting from its use. |
-// Users of this code must verify correctness for their application. |
- |
-// Assume that a class is already given for the object: |
-// Point with coordinates {float x, y;} |
-//=================================================================== |
- |
-// isLeft(): tests if a point is Left|On|Right of an infinite line. |
-// Input: three points P0, P1, and P2 |
-// Return: >0 for P2 left of the line through P0 and P1 |
-// =0 for P2 on the line |
-// <0 for P2 right of the line |
-// See: the January 2001 Algorithm on Area of Triangles |
-// return (float) ((P1.x - P0.x)*(P2.y - P0.y) - (P2.x - P0.x)*(P1.y - P0.y)); |
-double is_left(const _Line& line, const _Point& pt) { |
- _Vector P0 = line[1] - line[0]; |
- _Vector P2 = pt - line[0]; |
- return P0.cross(P2); |
-} |
- |
-double t_at(const _Line& line, const _Point& pt) { |
- double dx = line[1].x - line[0].x; |
- double dy = line[1].y - line[0].y; |
- if (fabs(dx) > fabs(dy)) { |
- if (approximately_zero(dx)) { |
- return 0; |
- } |
- return (pt.x - line[0].x) / dx; |
- } |
- if (approximately_zero(dy)) { |
- return 0; |
- } |
- return (pt.y - line[0].y) / dy; |
-} |
- |
-static void setMinMax(double x, int flags, double& minX, double& maxX) { |
- if (minX > x && (flags & (kFindTopMin | kFindBottomMin))) { |
- minX = x; |
- } |
- if (maxX < x && (flags & (kFindTopMax | kFindBottomMax))) { |
- maxX = x; |
- } |
-} |
- |
-void x_at(const _Point& p1, const _Point& p2, double top, double bottom, |
- int flags, double& minX, double& maxX) { |
- if (AlmostEqualUlps(p1.y, p2.y)) { |
- // It should be OK to bail early in this case. There's another edge |
- // which shares this end point which can intersect without failing to |
- // have a slope ... maybe |
- return; |
- } |
- |
- // p2.x is always greater than p1.x -- the part of points (p1, p2) are |
- // moving from the start of the cubic towards its end. |
- // if p1.y < p2.y, minX can be affected |
- // if p1.y > p2.y, maxX can be affected |
- double slope = (p2.x - p1.x) / (p2.y - p1.y); |
- int topFlags = flags & (kFindTopMin | kFindTopMax); |
- if (topFlags && ((top <= p1.y && top >= p2.y) |
- || (top >= p1.y && top <= p2.y))) { |
- double x = p1.x + (top - p1.y) * slope; |
- setMinMax(x, topFlags, minX, maxX); |
- } |
- int bottomFlags = flags & (kFindBottomMin | kFindBottomMax); |
- if (bottomFlags && ((bottom <= p1.y && bottom >= p2.y) |
- || (bottom >= p1.y && bottom <= p2.y))) { |
- double x = p1.x + (bottom - p1.y) * slope; |
- setMinMax(x, bottomFlags, minX, maxX); |
- } |
-} |
- |
-void xy_at_t(const _Line& line, double t, double& x, double& y) { |
- double one_t = 1 - t; |
- if (&x) { |
- x = one_t * line[0].x + t * line[1].x; |
- } |
- if (&y) { |
- y = one_t * line[0].y + t * line[1].y; |
- } |
-} |
- |
-_Point xy_at_t(const _Line& line, double t) { |
- double one_t = 1 - t; |
- _Point result = { one_t * line[0].x + t * line[1].x, one_t * line[0].y + t * line[1].y }; |
- return result; |
-} |