#nullable enable
using System;
using System.Collections.Generic;
using System.Runtime.CompilerServices;
#if USINGZ
namespace Clipper2ZLib
#else
namespace Clipper2Lib
#endif
{
public struct Point64
{
public long X;
public long Y;
#if USINGZ
public long Z;
#endif
public Point64(Point64 pt)
{
X = pt.X;
Y = pt.Y;
#if USINGZ
Z = pt.Z;
#endif
}
public Point64(Point64 pt, double scale)
{
X = (long) Math.Round(pt.X * scale, MidpointRounding.AwayFromZero);
Y = (long) Math.Round(pt.Y * scale, MidpointRounding.AwayFromZero);
#if USINGZ
Z = (long) Math.Round(pt.Z * scale, MidpointRounding.AwayFromZero);
#endif
}
public Point64(long x, long y
#if USINGZ
, long z = 0
#endif
) {
X = x;
Y = y;
#if USINGZ
Z = z;
#endif
}
public Point64(double x, double y
#if USINGZ
, double z = 0.0
#endif
) {
X = (long) Math.Round(x, MidpointRounding.AwayFromZero);
Y = (long) Math.Round(y, MidpointRounding.AwayFromZero);
#if USINGZ
Z = (long) Math.Round(z, MidpointRounding.AwayFromZero);
#endif
}
public Point64(PointD pt)
{
X = (long) Math.Round(pt.x, MidpointRounding.AwayFromZero);
Y = (long) Math.Round(pt.y, MidpointRounding.AwayFromZero);
#if USINGZ
Z = pt.z;
#endif
}
public Point64(PointD pt, double scale)
{
X = (long) Math.Round(pt.x * scale, MidpointRounding.AwayFromZero);
Y = (long) Math.Round(pt.y * scale, MidpointRounding.AwayFromZero);
#if USINGZ
Z = pt.z;
#endif
}
public static bool operator ==(Point64 lhs, Point64 rhs)
{
return lhs.X == rhs.X && lhs.Y == rhs.Y;
}
public static bool operator !=(Point64 lhs, Point64 rhs)
{
return lhs.X != rhs.X || lhs.Y != rhs.Y;
}
public static Point64 operator +(Point64 lhs, Point64 rhs)
{
return new Point64(lhs.X + rhs.X, lhs.Y + rhs.Y
#if USINGZ
, lhs.Z + rhs.Z
#endif
);
}
public static Point64 operator -(Point64 lhs, Point64 rhs)
{
return new Point64(lhs.X - rhs.X, lhs.Y - rhs.Y
#if USINGZ
, lhs.Z - rhs.Z
#endif
);
}
public readonly override string ToString()
{
#if USINGZ
return $"{X},{Y},{Z} ";
#else
return $"{X},{Y} ";
#endif
}
public readonly override bool Equals(object? obj)
{
if (obj != null && obj is Point64 p)
return this == p;
return false;
}
public readonly override int GetHashCode()
{
return HashCode.Combine(X, Y); }
}
public struct PointD
{
public double x;
public double y;
#if USINGZ
public long z;
#endif
public PointD(PointD pt)
{
x = pt.x;
y = pt.y;
#if USINGZ
z = pt.z;
#endif
}
public PointD(Point64 pt)
{
x = pt.X;
y = pt.Y;
#if USINGZ
z = pt.Z;
#endif
}
public PointD(Point64 pt, double scale)
{
x = pt.X * scale;
y = pt.Y * scale;
#if USINGZ
z = pt.Z;
#endif
}
public PointD(PointD pt, double scale)
{
x = pt.x * scale;
y = pt.y * scale;
#if USINGZ
z = pt.z;
#endif
}
public PointD(long x, long y
#if USINGZ
, long z = 0
#endif
) {
this.x = x;
this.y = y;
#if USINGZ
this.z = z;
#endif
}
public PointD(double x, double y
#if USINGZ
, long z = 0
#endif
) {
this.x = x;
this.y = y;
#if USINGZ
this.z = z;
#endif
}
public readonly string ToString(int precision = 2)
{
#if USINGZ
return string.Format($"{{0:F{precision}}},{{1:F{precision}}},{{2:D}}", x,y,z);
#else
return string.Format($"{{0:F{precision}}},{{1:F{precision}}}", x,y);
#endif
}
public static bool operator ==(PointD lhs, PointD rhs)
{
return InternalClipper.IsAlmostZero(lhs.x - rhs.x) &&
InternalClipper.IsAlmostZero(lhs.y - rhs.y);
}
public static bool operator !=(PointD lhs, PointD rhs)
{
return !InternalClipper.IsAlmostZero(lhs.x - rhs.x) ||
!InternalClipper.IsAlmostZero(lhs.y - rhs.y);
}
public readonly override bool Equals(object? obj)
{
if (obj != null && obj is PointD p)
return this == p;
return false;
}
public void Negate() { x = -x; y = -y; }
public readonly override int GetHashCode()
{
return HashCode.Combine(x, y); }
}
public struct Rect64
{
public long left;
public long top;
public long right;
public long bottom;
public Rect64(long l, long t, long r, long b)
{
left = l;
top = t;
right = r;
bottom = b;
}
public Rect64(bool isValid)
{
if (isValid)
{
left = 0; top = 0; right = 0; bottom = 0;
}
else
{
left = long.MaxValue; top = long.MaxValue;
right = long.MinValue; bottom = long.MinValue;
}
}
public Rect64(Rect64 rec)
{
left = rec.left;
top = rec.top;
right = rec.right;
bottom = rec.bottom;
}
public long Width
{ readonly get => right - left;
set => right = left + value;
}
public long Height
{ readonly get => bottom - top;
set => bottom = top + value;
}
public readonly bool IsEmpty()
{
return bottom <= top || right <= left;
}
public readonly bool IsValid()
{
return left < long.MaxValue;
}
public readonly Point64 MidPoint()
{
return new Point64((left + right) /2, (top + bottom)/2);
}
public readonly bool Contains(Point64 pt)
{
return pt.X > left && pt.X < right &&
pt.Y > top && pt.Y < bottom;
}
public readonly bool Contains(Rect64 rec)
{
return rec.left >= left && rec.right <= right &&
rec.top >= top && rec.bottom <= bottom;
}
public readonly bool Intersects(Rect64 rec)
{
return (Math.Max(left, rec.left) <= Math.Min(right, rec.right)) &&
(Math.Max(top, rec.top) <= Math.Min(bottom, rec.bottom));
}
public readonly Path64 AsPath()
{
Path64 result = new Path64(4)
{
new Point64(left, top),
new Point64(right, top),
new Point64(right, bottom),
new Point64(left, bottom)
};
return result;
}
}
public struct RectD
{
public double left;
public double top;
public double right;
public double bottom;
public RectD(double l, double t, double r, double b)
{
left = l;
top = t;
right = r;
bottom = b;
}
public RectD(RectD rec)
{
left = rec.left;
top = rec.top;
right = rec.right;
bottom = rec.bottom;
}
public RectD(bool isValid)
{
if (isValid)
{
left = 0; top = 0; right = 0; bottom = 0;
}
else
{
left = double.MaxValue; top = double.MaxValue;
right = -double.MaxValue; bottom = -double.MaxValue;
}
}
public double Width
{ readonly get => right - left;
set => right = left + value;
}
public double Height
{ readonly get => bottom - top;
set => bottom = top + value;
}
public readonly bool IsEmpty()
{
return bottom <= top || right <= left;
}
public readonly PointD MidPoint()
{
return new PointD((left + right) / 2, (top + bottom) / 2);
}
public readonly bool Contains(PointD pt)
{
return pt.x > left && pt.x < right &&
pt.y > top && pt.y < bottom;
}
public readonly bool Contains(RectD rec)
{
return rec.left >= left && rec.right <= right &&
rec.top >= top && rec.bottom <= bottom;
}
public readonly bool Intersects(RectD rec)
{
return (Math.Max(left, rec.left) < Math.Min(right, rec.right)) &&
(Math.Max(top, rec.top) < Math.Min(bottom, rec.bottom));
}
public readonly PathD AsPath()
{
PathD result = new PathD(4)
{
new PointD(left, top),
new PointD(right, top),
new PointD(right, bottom),
new PointD(left, bottom)
};
return result;
}
}
public class Path64 : List<Point64>
{
public Path64() : base() { }
public Path64(int capacity = 0) : base(capacity) { }
public Path64(IEnumerable<Point64> path) : base(path) { }
public override string ToString()
{
return string.Join(", ", this);
}
}
public class Paths64 : List<Path64>
{
public Paths64() : base() { }
public Paths64(int capacity = 0) : base(capacity) { }
public Paths64(IEnumerable<Path64> paths) : base(paths) { }
public override string ToString()
{
return string.Join(Environment.NewLine, this);
}
}
public class PathD : List<PointD>
{
public PathD() : base() { }
public PathD(int capacity = 0) : base(capacity) { }
public PathD(IEnumerable<PointD> path) : base(path) { }
public string ToString(int precision = 2)
{
return string.Join(", ", ConvertAll(x => x.ToString(precision)));
}
}
public class PathsD : List<PathD>
{
public PathsD() : base() { }
public PathsD(int capacity = 0) : base(capacity) { }
public PathsD(IEnumerable<PathD> paths) : base(paths) { }
public string ToString(int precision = 2)
{
return string.Join(Environment.NewLine, ConvertAll(x => x.ToString(precision)));
}
}
public enum ClipType
{
NoClip,
Intersection,
Union,
Difference,
Xor
}
public enum PathType
{
Subject,
Clip
}
public enum FillRule
{
EvenOdd,
NonZero,
Positive,
Negative
}
public static class InternalClipper
{
internal const long MaxInt64 = 9223372036854775807;
internal const long MaxCoord = MaxInt64 / 4;
internal const double max_coord = MaxCoord;
internal const double min_coord = -MaxCoord;
internal const long Invalid64 = MaxInt64;
internal const double floatingPointTolerance = 1E-12;
internal const double defaultMinimumEdgeLength = 0.1;
private static readonly string
precision_range_error = "Error: Precision is out of range.";
public static double CrossProduct(Point64 pt1, Point64 pt2, Point64 pt3)
{
return ((double) (pt2.X - pt1.X) * (pt3.Y - pt2.Y) -
(double) (pt2.Y - pt1.Y) * (pt3.X - pt2.X));
}
public static int CrossProductSign(Point64 pt1, Point64 pt2, Point64 pt3)
{
long a = pt2.X - pt1.X;
long b = pt3.Y - pt2.Y;
long c = pt2.Y - pt1.Y;
long d = pt3.X - pt2.X;
UInt128Struct ab = MultiplyUInt64((ulong) Math.Abs(a), (ulong) Math.Abs(b));
UInt128Struct cd = MultiplyUInt64((ulong) Math.Abs(c), (ulong) Math.Abs(d));
int signAB = TriSign(a) * TriSign(b);
int signCD = TriSign(c) * TriSign(d);
if (signAB == signCD)
{
int result;
if (ab.hi64 == cd.hi64)
{
if (ab.lo64 == cd.lo64) return 0;
result = (ab.lo64 > cd.lo64) ? 1 : -1;
}
else result = (ab.hi64 > cd.hi64) ? 1 : -1;
return (signAB > 0) ? result : -result;
}
return (signAB > signCD) ? 1 : -1;
}
#if USINGZ
public static Path64 SetZ(Path64 path, long Z)
{
Path64 result = new Path64(path.Count);
foreach (Point64 pt in path) result.Add(new Point64(pt.X, pt.Y, Z));
return result;
}
#endif
[MethodImpl(MethodImplOptions.AggressiveInlining)]
internal static void CheckPrecision(int precision)
{
if (precision < -8 || precision > 8)
throw new Exception(precision_range_error);
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
internal static bool IsAlmostZero(double value)
{
return (Math.Abs(value) <= floatingPointTolerance);
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
internal static int TriSign(long x) {
return (x < 0) ? -1 : (x > 0) ? 1 : 0;
}
public struct UInt128Struct
{
public ulong lo64;
public ulong hi64;
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public static UInt128Struct MultiplyUInt64(ulong a, ulong b) {
ulong x1 = (a & 0xFFFFFFFF) * (b & 0xFFFFFFFF);
ulong x2 = (a >> 32) * (b & 0xFFFFFFFF) + (x1 >> 32);
ulong x3 = (a & 0xFFFFFFFF) * (b >> 32) + (x2 & 0xFFFFFFFF);
UInt128Struct result;
result.lo64 = (x3 & 0xFFFFFFFF) << 32 | (x1 & 0xFFFFFFFF);
result.hi64 = (a >> 32) * (b >> 32) + (x2 >> 32) + (x3 >> 32);
return result;
}
internal static bool ProductsAreEqual(long a, long b, long c, long d)
{
ulong absA = (ulong) Math.Abs(a);
ulong absB = (ulong) Math.Abs(b);
ulong absC = (ulong) Math.Abs(c);
ulong absD = (ulong) Math.Abs(d);
UInt128Struct mul_ab = MultiplyUInt64(absA, absB);
UInt128Struct mul_cd = MultiplyUInt64(absC, absD);
int sign_ab = TriSign(a) * TriSign(b);
int sign_cd = TriSign(c) * TriSign(d);
return mul_ab.lo64 == mul_cd.lo64 && mul_ab.hi64 == mul_cd.hi64 && sign_ab == sign_cd;
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
internal static bool IsCollinear(Point64 pt1, Point64 sharedPt, Point64 pt2)
{
long a = sharedPt.X - pt1.X;
long b = pt2.Y - sharedPt.Y;
long c = sharedPt.Y - pt1.Y;
long d = pt2.X - sharedPt.X;
return ProductsAreEqual(a, b, c, d);
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
internal static double DotProduct(Point64 pt1, Point64 pt2, Point64 pt3)
{
return ((double) (pt2.X - pt1.X) * (pt3.X - pt2.X) +
(double) (pt2.Y - pt1.Y) * (pt3.Y - pt2.Y));
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
internal static double CrossProduct(PointD vec1, PointD vec2)
{
return (vec1.y * vec2.x - vec2.y * vec1.x);
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
internal static double DotProduct(PointD vec1, PointD vec2)
{
return (vec1.x * vec2.x + vec1.y * vec2.y);
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
internal static long CheckCastInt64(double val)
{
if ((val >= max_coord) || (val <= min_coord)) return Invalid64;
return (long)Math.Round(val, MidpointRounding.AwayFromZero);
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public static bool GetLineIntersectPt(Point64 ln1a,
Point64 ln1b, Point64 ln2a, Point64 ln2b, out Point64 ip)
{
double dy1 = (ln1b.Y - ln1a.Y);
double dx1 = (ln1b.X - ln1a.X);
double dy2 = (ln2b.Y - ln2a.Y);
double dx2 = (ln2b.X - ln2a.X);
double det = dy1 * dx2 - dy2 * dx1;
if (det == 0.0)
{
ip = new Point64();
return false;
}
double t = ((ln1a.X - ln2a.X) * dy2 - (ln1a.Y - ln2a.Y) * dx2) / det;
if (t <= 0.0) ip = ln1a;
else if (t >= 1.0) ip = ln1b;
else
{
ip.X = (long) (ln1a.X + t * dx1);
ip.Y = (long) (ln1a.Y + t * dy1);
#if USINGZ
ip.Z = 0;
#endif
}
return true;
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public static bool GetLineIntersectPt(PointD ln1a,
PointD ln1b, PointD ln2a, PointD ln2b, out PointD ip)
{
double dy1 = (ln1b.y - ln1a.y);
double dx1 = (ln1b.x - ln1a.x);
double dy2 = (ln2b.y - ln2a.y);
double dx2 = (ln2b.x - ln2a.x);
double det = dy1 * dx2 - dy2 * dx1;
if (det == 0.0)
{
ip = new PointD();
return false;
}
double t = ((ln1a.x - ln2a.x) * dy2 - (ln1a.y - ln2a.y) * dx2) / det;
if (t <= 0.0) ip = ln1a;
else if (t >= 1.0) ip = ln1b;
else
{
ip.x = (ln1a.x + t * dx1);
ip.y = (ln1a.y + t * dy1);
#if USINGZ
ip.z = 0;
#endif
}
return true;
}
internal static bool SegsIntersect(Point64 seg1a,
Point64 seg1b, Point64 seg2a, Point64 seg2b, bool inclusive = false)
{
double dy1 = (seg1b.Y - seg1a.Y);
double dx1 = (seg1b.X - seg1a.X);
double dy2 = (seg2b.Y - seg2a.Y);
double dx2 = (seg2b.X - seg2a.X);
double cp = dy1 * dx2 - dy2 * dx1;
if (cp == 0) return false;
if (inclusive)
{
double t = ((seg1a.X - seg2a.X) * dy2 - (seg1a.Y - seg2a.Y) * dx2);
if (t == 0) return true;
if (t > 0)
{
if (cp < 0 || t > cp) return false;
}
else if (cp > 0 || t < cp) return false;
t = ((seg1a.X - seg2a.X) * dy1 - (seg1a.Y - seg2a.Y) * dx1);
if (t == 0) return true;
if (t > 0) return (cp > 0 && t <= cp);
else return (cp < 0 && t >= cp); }
else
{
double t = ((seg1a.X - seg2a.X) * dy2 - (seg1a.Y - seg2a.Y) * dx2);
if (t == 0) return false;
if (t > 0)
{
if (cp < 0 || t >= cp) return false;
}
else if (cp > 0 || t <= cp) return false;
t = ((seg1a.X - seg2a.X) * dy1 - (seg1a.Y - seg2a.Y) * dx1);
if (t == 0) return false;
if (t > 0) return (cp > 0 && t < cp);
else return (cp < 0 && t > cp); }
}
public static Rect64 GetBounds(Path64 path)
{
if (path.Count == 0) return new Rect64();
Rect64 result = Clipper.InvalidRect64;
foreach (Point64 pt in path)
{
if (pt.X < result.left) result.left = pt.X;
if (pt.X > result.right) result.right = pt.X;
if (pt.Y < result.top) result.top = pt.Y;
if (pt.Y > result.bottom) result.bottom = pt.Y;
}
return result;
}
public static Point64 GetClosestPtOnSegment(Point64 offPt,
Point64 seg1, Point64 seg2)
{
if (seg1.X == seg2.X && seg1.Y == seg2.Y) return seg1;
double dx = (seg2.X - seg1.X);
double dy = (seg2.Y - seg1.Y);
double q = ((offPt.X - seg1.X) * dx +
(offPt.Y - seg1.Y) * dy) / ((dx*dx) + (dy*dy));
if (q < 0) q = 0; else if (q > 1) q = 1;
return new Point64(
seg1.X + Math.Round(q * dx, MidpointRounding.ToEven),
seg1.Y + Math.Round(q * dy, MidpointRounding.ToEven)
);
}
public static PointInPolygonResult PointInPolygon(Point64 pt, Path64 polygon)
{
int len = polygon.Count, start = 0;
if (len < 3) return PointInPolygonResult.IsOutside;
while (start < len && polygon[start].Y == pt.Y) start++;
if (start == len) return PointInPolygonResult.IsOutside;
bool isAbove = polygon[start].Y < pt.Y, startingAbove = isAbove;
int val = 0, i = start + 1, end = len;
while (true)
{
if (i == end)
{
if (end == 0 || start == 0) break;
end = start;
i = 0;
}
if (isAbove)
{
while (i < end && polygon[i].Y < pt.Y) i++;
}
else
{
while (i < end && polygon[i].Y > pt.Y) i++;
}
if (i == end) continue;
Point64 curr = polygon[i], prev;
if (i > 0) prev = polygon[i - 1];
else prev = polygon[len - 1];
if (curr.Y == pt.Y)
{
if (curr.X == pt.X || (curr.Y == prev.Y &&
((pt.X < prev.X) != (pt.X < curr.X))))
return PointInPolygonResult.IsOn;
i++;
if (i == start) break;
continue;
}
if (pt.X < curr.X && pt.X < prev.X)
{
}
else if (pt.X > prev.X && pt.X > curr.X)
{
val = 1 - val; }
else
{
int cps2 = CrossProductSign(prev, curr, pt);
if (cps2 == 0) return PointInPolygonResult.IsOn;
if ((cps2 < 0) == isAbove) val = 1 - val;
}
isAbove = !isAbove;
i++;
}
if (isAbove == startingAbove) return val == 0 ? PointInPolygonResult.IsOutside : PointInPolygonResult.IsInside;
if (i == len) i = 0;
int cps = (i == 0) ?
CrossProductSign(polygon[len - 1], polygon[0], pt) :
CrossProductSign(polygon[i - 1], polygon[i], pt);
if (cps == 0) return PointInPolygonResult.IsOn;
if ((cps < 0) == isAbove) val = 1 - val;
return val == 0 ? PointInPolygonResult.IsOutside : PointInPolygonResult.IsInside;
}
public static bool Path2ContainsPath1(Path64 path1, Path64 path2)
{
PointInPolygonResult pip = PointInPolygonResult.IsOn;
foreach (Point64 pt in path1)
{
switch (PointInPolygon(pt, path2))
{
case PointInPolygonResult.IsOutside:
if (pip == PointInPolygonResult.IsOutside) return false;
pip = PointInPolygonResult.IsOutside;
break;
case PointInPolygonResult.IsInside:
if (pip == PointInPolygonResult.IsInside) return true;
pip = PointInPolygonResult.IsInside;
break;
default: break;
}
}
Point64 mp = GetBounds(path1).MidPoint();
return InternalClipper.PointInPolygon(mp, path2) != PointInPolygonResult.IsOutside;
}
}
}