private static VNode ProcessDataEvent(VDataEvent e, VNode Root, VoronoiGraph VG, double ys, out VDataNode[] CircleCheckList)
        {
            if (Root == null)
            {
                Root = new VDataNode(e.DataPoint);
                CircleCheckList = new VDataNode[] { (VDataNode)Root };
                return Root;
            }
            //Знаходимо вузол для заміщення
            VNode C = VNode.FindDataNode(Root, ys, e.DataPoint[0]);
            //Створюємо піддерево з одним ребром, але двома VEdgeNodes
            VoronoiEdge VE = new VoronoiEdge();
            VE.LeftData = ((VDataNode)C).DataPoint;
            VE.RightData = e.DataPoint;
            VE.VVertexA = Fortune.VVUnkown;
            VE.VVertexB = Fortune.VVUnkown;
            VG.Edges.Add(VE);

            VNode SubRoot;
            if (Math.Abs(VE.LeftData[1] - VE.RightData[1]) < 1e-10)
            {
                if (VE.LeftData[0] < VE.RightData[0])
                {
                    SubRoot = new VEdgeNode(VE, false);
                    SubRoot.Left = new VDataNode(VE.LeftData);
                    SubRoot.Right = new VDataNode(VE.RightData);
                }
                else
                {
                    SubRoot = new VEdgeNode(VE, true);
                    SubRoot.Left = new VDataNode(VE.RightData);
                    SubRoot.Right = new VDataNode(VE.LeftData);
                }
                CircleCheckList = new VDataNode[] { (VDataNode)SubRoot.Left, (VDataNode)SubRoot.Right };
            }
            else
            {
                SubRoot = new VEdgeNode(VE, false);
                SubRoot.Left = new VDataNode(VE.LeftData);
                SubRoot.Right = new VEdgeNode(VE, true);
                SubRoot.Right.Left = new VDataNode(VE.RightData);
                SubRoot.Right.Right = new VDataNode(VE.LeftData);
                CircleCheckList = new VDataNode[] { (VDataNode)SubRoot.Left, (VDataNode)SubRoot.Right.Left, (VDataNode)SubRoot.Right.Right };
            }

            //"Застосовуємо" піддерево
            if (C.Parent == null)
                return SubRoot;
            C.Parent.Replace(C, SubRoot);
            return Root;
        }
 public static Edge Pack_into_Edge(VoronoiEdge VE)
 {
     int x1 = 0;
     int y1 = 0;
     int x2 = 0;
     int y2 = 0;
     if (VE.IsInfinite)
     {
         x1 = (int)(1000 * VE.DirectionVector[0] + VE.FixedPoint[0]);
         y1 = (int)(1000 * VE.DirectionVector[1] + VE.FixedPoint[1]);
         x2 = (int)(-1000 * VE.DirectionVector[0] + VE.FixedPoint[0]);
         y2 = (int)(-1000 * VE.DirectionVector[1] + VE.FixedPoint[1]);
     }
     else if (VE.IsPartlyInfinite)
     {
         x1 = (int)VE.FixedPoint[0];
         y1 = (int)VE.FixedPoint[1];
         x2 = (int)(1000 * VE.DirectionVector[0] + VE.FixedPoint[0]);
         y2 = (int)(1000 * VE.DirectionVector[1] + VE.FixedPoint[1]);
     }
     else
     {
         x1 = (int)VE.VVertexA[0];
         y1 = (int)VE.VVertexA[1];
         x2 = (int)VE.VVertexB[0];
         y2 = (int)VE.VVertexB[1];
     }
     Point p1 = new Point(x1, y1);
     Point p2 = new Point(x2, y2);
     Edge e = new Edge(p1, p2);
     return e;
 }
        private static VNode ProcessCircleEvent(VCircleEvent e, VNode Root, VoronoiGraph VG, double ys, out VDataNode[] CircleCheckList)
        {
            VDataNode a, b, c;
            VEdgeNode e1, e2;
            b = e.NodeN;
            a = VNode.LeftDataNode(b);
            c = VNode.RightDataNode(b);
            if (a == null || b.Parent == null || c == null || !a.DataPoint.Equals(e.NodeL.DataPoint) || !c.DataPoint.Equals(e.NodeR.DataPoint))
            {
                CircleCheckList = new VDataNode[] { };
                return Root; // повертаємось, бо графік змінився
            }
            e1 = (VEdgeNode)b.Parent;
            CircleCheckList = new VDataNode[] { a, c };
            //Створюємо нову вершину
            Vector VNew = new Vector(e.Center[0], e.Center[1]);
            VG.Vertizes.Add(VNew);
            //2. виясняємо, чи а або с знаходяться у віддаленій частині дерева (інший - брат b), і призначаємо нову вершину
            if (e1.Left == b) // c - брат
            {
                e2 = VNode.EdgeToRightDataNode(a);
                // замінюємо e1 правим нащадком
                e1.Parent.Replace(e1, e1.Right);
            }
            else // a - брат
            {
                e2 = VNode.EdgeToRightDataNode(b);

                // замінюємо e1 лівим нащадком
                e1.Parent.Replace(e1, e1.Left);
            }
            e1.Edge.AddVertex(VNew);
            e2.Edge.AddVertex(VNew);

            //Замінюємо e2 новим ребром
            VoronoiEdge VE = new VoronoiEdge();
            VE.LeftData = a.DataPoint;
            VE.RightData = c.DataPoint;
            VE.AddVertex(VNew);
            VG.Edges.Add(VE);

            VEdgeNode VEN = new VEdgeNode(VE, false);
            VEN.Left = e2.Left;
            VEN.Right = e2.Right;
            if (e2.Parent == null)
                return VEN;
            e2.Parent.Replace(e2, VEN);
            return Root;
        }
 private void DrawVoronoiEdge(VoronoiEdge VE)
 {
     int x1 = 0;
     int y1 = 0;
     int x2 = 0;
     int y2 = 0;
     if (VE.IsInfinite)
     {
         x1 = (int)(1000 * VE.DirectionVector[0] + VE.FixedPoint[0]);
         y1 = (int)(1000 * VE.DirectionVector[1] + VE.FixedPoint[1]);
         x2 = (int)(-1000 * VE.DirectionVector[0] + VE.FixedPoint[0]);
         y2 = (int)(-1000 * VE.DirectionVector[1] + VE.FixedPoint[1]);
     }
     else if (VE.IsPartlyInfinite)
     {
         x1 = (int)VE.FixedPoint[0];
         y1 = (int)VE.FixedPoint[1];
         x2 = (int)(1000 * VE.DirectionVector[0] + VE.FixedPoint[0]);
         y2 = (int)(1000 * VE.DirectionVector[1] + VE.FixedPoint[1]);
     }
     else
     {
         x1 = (int)VE.VVertexA[0];
         y1 = (int)VE.VVertexA[1];
         x2 = (int)VE.VVertexB[0];
         y2 = (int)VE.VVertexB[1];
     }
     GL.LineWidth(1);
     GL.Begin(BeginMode.Lines);
     GL.Color3(Color.Black);
     GL.Vertex2(x1, y1);
     GL.Vertex2(x2, y2);
     GL.End();
 }