public bool Contains(Point2LR point) { if (point.X > LongRational.Max(this.start.X, this.end.X) || point.X < LongRational.Min(this.start.X, this.end.X) || (point.Y > LongRational.Max(this.start.Y, this.end.Y) || point.Y < LongRational.Min(this.start.Y, this.end.Y))) { return(false); } if (this.start == this.end) { return(true); } Vector2LR delta = this.GetDelta(); Vector2LR u = point - this.start; Vector2LR v = new Vector2LR(-delta.Y, delta.X); if (!Vector2LR.DotProduct(u, v).IsZero) { return(false); } LongRational longRational = Vector2LR.DotProduct(u, delta); if (longRational.IsNegative) { return(false); } LongRational lengthSquared = delta.GetLengthSquared(); return(!(longRational > lengthSquared)); }
public Point2LR GetClosestPoint(Point2LR point) { Vector2LR u = point - this.start; Vector2LR v = this.end - this.start; LongRational lengthSquared = v.GetLengthSquared(); LongRational longRational = Vector2LR.DotProduct(u, v); return(longRational.IsNegative || !(longRational <= lengthSquared) ? (longRational.IsNegative ? this.start : this.end) : longRational / lengthSquared * v + this.start); }
public LongRational GetDistanceSquared(Point2LR point) { Vector2LR u = point - this.start; Vector2LR v = this.end - this.start; LongRational lengthSquared = v.GetLengthSquared(); LongRational longRational = Vector2LR.DotProduct(u, v); return(longRational.IsNegative || !(longRational <= lengthSquared) ? (!longRational.IsNegative ? (this.end - point).GetLengthSquared() : (this.start - point).GetLengthSquared()) : (longRational / lengthSquared * v + this.start - point).GetLengthSquared()); }
public LongRational GetNormalizedProjection(Point2LR point) { if (this.start == this.end) { return(LongRational.NaN); } Vector2LR delta = this.GetDelta(); return(Vector2LR.DotProduct(point - this.start, delta) / delta.GetLengthSquared()); }
public bool ContainsProjection(Point2LR point) { Vector2LR delta = this.GetDelta(); LongRational longRational = Vector2LR.DotProduct(point - this.start, delta); if (longRational.IsNegative) { return(false); } LongRational lengthSquared = delta.GetLengthSquared(); return(!(longRational.Square() > lengthSquared)); }
public static bool GetIntersectionCoefficients( Ray2LR a, Segment2LR b, out LongRational?p, out LongRational?q, out bool areParallel) { if (a.direction.IsZero) { areParallel = false; p = new LongRational?(); bool flag; if (flag = b.Contains(a.origin)) { Vector2LR delta = b.GetDelta(); q = new LongRational?(Vector2LR.DotProduct(a.origin - b.Start, delta) / delta.GetLengthSquared()); } else { q = new LongRational?(); } return(flag); } Vector2LR v = new Vector2LR(-a.direction.Y, a.direction.X); Vector2LR u = b.Start - a.Origin; LongRational longRational1 = Vector2LR.DotProduct(u, v); LongRational longRational2 = Vector2LR.DotProduct(b.End - a.Origin, v); if (longRational1.IsPositive && longRational2.IsPositive || longRational1.IsNegative && longRational2.IsNegative) { p = new LongRational?(); q = new LongRational?(); areParallel = false; return(false); } Vector2LR delta1 = b.GetDelta(); if (delta1.IsZero) { areParallel = false; p = new LongRational?(Vector2LR.DotProduct(u, a.direction) / a.direction.GetLengthSquared()); q = new LongRational?(LongRational.Zero); return(true); } LongRational longRational3 = a.Direction.X * delta1.Y - a.Direction.Y * delta1.X; if (longRational3.IsZero) { p = new LongRational?(); q = new LongRational?(); areParallel = true; return(true); } areParallel = false; LongRational longRational4 = Vector2LR.CrossProduct(u, delta1); if (longRational3.IsPositive) { if (longRational4.IsNegative) { p = new LongRational?(); q = new LongRational?(); return(false); } } else if (longRational4.IsPositive) { p = new LongRational?(); q = new LongRational?(); return(false); } LongRational longRational5 = Vector2LR.CrossProduct(u, a.Direction); if (longRational3.IsPositive) { if (longRational5.IsNegative || longRational5 > longRational3) { p = new LongRational?(); q = new LongRational?(); return(false); } } else if (longRational5.IsPositive || longRational5 < longRational3) { p = new LongRational?(); q = new LongRational?(); return(false); } p = new LongRational?(longRational4 / longRational3); q = new LongRational?(longRational5 / longRational3); return(true); }