Exemplo n.º 1
0
        private void UpdateCircumcircle()
        {
            var p0 = this.Vertices[0];
            var p1 = this.Vertices[1];
            var p2 = this.Vertices[2];
            var dA = p0.X * p0.X + p0.Y * p0.Y;
            var dB = p1.X * p1.X + p1.Y * p1.Y;
            var dC = p2.X * p2.X + p2.Y * p2.Y;

            var aux1 = dA * (p2.Y - p1.Y)
                       + dB * (p0.Y - p2.Y)
                       + dC * (p1.Y - p0.Y);
            var aux2 = -(dA * (p2.X - p1.X)
                         + dB * (p0.X - p2.X)
                         + dC * (p1.X - p0.X));
            var div = 2 * (p0.X * (p2.Y - p1.Y)
                           + p1.X * (p0.Y - p2.Y)
                           + p2.X * (p1.Y - p0.Y));

            if (Math.Abs(div) < Settings.Tolerance)
            {
                throw new Exception("Divisor too small");
            }

            var center = new DelaunayPoint(aux1 / div, aux2 / div);

            this.Circumcenter  = center;
            this.RadiusSquared = (center.X - p0.X) * (center.X - p0.X)
                                 + (center.Y - p0.Y) * (center.Y - p0.Y);
        }
Exemplo n.º 2
0
        /// <summary>
        ///     Initializes a new instance of the <see cref="DelaunayTriangle" /> class.
        /// </summary>
        /// <param name="point1">Point A.</param>
        /// <param name="point2">Point B.</param>
        /// <param name="point3">Point C.</param>
        public DelaunayTriangle(DelaunayPoint point1, DelaunayPoint point2, DelaunayPoint point3)
        {
            this.Vertices = new List <DelaunayPoint>();
            if (!IsCounterClockwise(point1, point2, point3))
            {
                this.Vertices.Add(point1);
                this.Vertices.Add(point3);
                this.Vertices.Add(point2);
            }
            else
            {
                this.Vertices.Add(point1);
                this.Vertices.Add(point2);
                this.Vertices.Add(point3);
            }

            this.Vertices[0].AdjacentTriangles.Add(this);
            this.Vertices[1].AdjacentTriangles.Add(this);
            this.Vertices[2].AdjacentTriangles.Add(this);
            this.UpdateCircumcircle();
        }
Exemplo n.º 3
0
 /// <summary>
 ///     Initializes a new instance of the <see cref="DelaunayEdge" /> class.
 /// </summary>
 /// <param name="startPoint">Start point.</param>
 /// <param name="endPoint">End point.</param>
 public DelaunayEdge(DelaunayPoint startPoint, DelaunayPoint endPoint) : base(startPoint, endPoint)
 {
 }