using Newtonsoft.Json;
using System;
using System.Collections.Generic;
using System.Drawing;
using System.IO;
using System.Linq;
using System.Numerics;
using System.Text;
using System.Text.RegularExpressions;
using System.Threading.Tasks;
using System.Windows.Forms;
using Yw.WinFrmUI.Q3d;
using Yw.WinFrmUI.Q3d;
using static System.Math;
using static Yw.WinFrmUI.Q3d.ObjectEnum;

namespace Yw.WinFrmUI.Q3d
{
    public partial class NetworkViewModel
    {


        public void BuildRelation()
        {
            //读取坐标
            int k1 = 0;//表示管线索引
            int k2 = 0;//表示节点索引
            Nodes.Sort((a, b) => string.Compare(a.ID, b.ID));

            //建立点线关系链表StartNode，先将管线以Node1（节点1的ID）排序，再将Nodes按ID排序，建立两个游标k1、k2，正向一次循环，建立链表关系
            //时间复杂度 O(n)
            Links.Sort((a, b) => string.Compare(a.Node1, b.Node1));
            k2 = 0;
            k1 = 0;
            while (k1 < Links.Count)
            {
                var p = Links[k1];
                var J = Nodes[k2];
                var k0 = k2;
                while (J.ID != p.Node1 && k2 < Nodes.Count)
                {
                    k2++;
                    if (k2 < Nodes.Count) J = Nodes[k2];
                }
                if (k2 == Nodes.Count)
                {
                    k2 = k0;
                    k1++;
                    p.Visible = false;
                    continue;
                    //throw new Exception($"未找到Link[{p.ID}]的起始节点[{p.Node1}]");
                }
                p.StartNode = J;
                p.Node1 = J.ID;
                if (J.MaxDiameter < p.Diameter) J.MaxDiameter = p.Diameter;
                J.Links.Add(p);
                k1++;
            }

            //建立点线关系链表StartNode，先将管线以Node2（节点1的ID）排序，再将Nodes按ID排序，建立两个游标k1、k2，正向一次循环，建立链表关系
            //时间复杂度 O(n)
            Links.Sort((a, b) => string.Compare(a.Node2, b.Node2));
            k2 = 0;
            k1 = 0;
            while (k1 < Links.Count)
            {
                var p = Links[k1];
                var J = Nodes[k2];
                var k0 = k2;
                while (J.ID != p.Node2 && k2 < Nodes.Count)
                {
                    k2++;
                    if (k2 < Nodes.Count) J = Nodes[k2];
                }
                if (k2 == Nodes.Count)
                {
                    k2 = k0;
                    k1++;
                    p.Visible = false;

                    continue;
                    //throw new Exception($"未找到Link[{p.ID}]的终止节点[{p.Node2}]");
                }
                p.EndNode = J;
                p.Node2 = J.ID;
                if (J.MaxDiameter < p.Diameter) J.MaxDiameter = p.Diameter;
                J.Links.Add(p);
                k1++;
            }


            if (dict_dataset == null)
            {
                dict_dataset = new Dictionary<string, Dataset>();
            }
            pumps.ForEach(p =>
            {
                foreach (var kp in p.Datasets)
                {
                    if (!dict_dataset.ContainsKey(kp.Value.Name))
                    {
                        dict_dataset.Add(kp.Value.Name, kp.Value);
                    }
                }
            });
            pumps.ForEach((pump) =>
            {
                if (dict_dataset.ContainsKey(pump.HeadCurve))
                {
                    Dataset ds = dict_dataset[pump.HeadCurve];

                    if (!pump.Datasets.ContainsKey("流量扬程曲线"))
                    {
                        pump.Datasets.Add("流量扬程曲线", ds);
                    }
                    else
                    {
                        pump.Datasets["流量扬程曲线"] = ds;
                    }
                    ds.pump = pump;
                }

                if (!pump.Datasets.ContainsKey("流量扬程曲线"))
                {
                    pump.Datasets.Add("流量扬程曲线", new Dataset(pump.ID + "_Head", pump));
                }
                if (pump.Datasets["流量扬程曲线"].Data.Count <= 0)
                {
                    pump.Datasets["流量扬程曲线"].Data.Add(new PointF((float)pump.额定流量, (float)pump.额定扬程));
                }
                pump.HeadCurve = pump.Datasets["流量扬程曲线"].Name;
            });

            //dict_dataset.Clear();
            //valveNodes.ForEach(v => Nodes.Remove(v)); ;
            Hash_ID = new HashSet<string>();
            Nodes.ForEach(o => Hash_ID.Add(o.ID));
            Links.ForEach(o => Hash_ID.Add(o.ID));
        }

        public bool BuildFromInp(string filePath)
        {
            if (filePath == null || !File.Exists(filePath)) return false;
            List<Coor> points = new List<Coor>();
            dict_dataset = new Dictionary<string, Dataset>();
            StreamReader sr = new StreamReader(filePath);
            //try
            {
                Nodes = new NodeViewModelList();
                Links = new LinkViewModelList();

                string line;
                string section = "";
                while ((line = sr.ReadLine()) != null)
                {
                    if (line.Trim().StartsWith("["))
                    {
                        section = line.TrimStart('[').TrimEnd(']');
                    }
                    else
                    {
                        string s = line.Trim('\t').Trim(' ');
                        if (s.Length == 0 || s[0] == ';') continue;
                        Parts parts = new Parts(line);
                        section = section.ToUpper();
                        IBaseViewModel o = null;
                        switch (section)
                        {
                            case "JUNCTIONS":
                                {
                                    NodeViewModel j = null;
                                    if (parts[-2] == "Meter")
                                    {
                                        var m = new MeterViewModel();
                                        //m.DlTemplateID = parts.ToString(-4, "");
                                        j = m;
                                    }
                                    else if (parts[-2] == "Nozzle")
                                    {
                                        var n = new NozzleViewModel();
                                        j = n;
                                        n.FlowCoefficient = parts.ToFloat(-4, 0);
                                    }
                                    else
                                        j = new JunctionViewModel();

                                    o = j;
                                    j.ID = parts.ToString(0, null);
                                    j.Z = parts.ToFloat(1, 0);
                                    j.Demand = parts.ToFloat(2, 0);
                                    j.PatternID = parts.ToString(3, null);
                                    if (j.PatternID == "NONE") j.PatternID = null;
                                    j.Level = parts.ToInt(-1, 0);

                                    Nodes.Add(j);
                                }
                                break;

                            case "RESERVOIRS":
                                {
                                    ReservoirViewModel r = new ReservoirViewModel();
                                    o = r;
                                    r.ID = parts.ToString(0, null);
                                    r.Head = parts.ToFloat(1, 0);
                                    r.PatternID = parts.ToString(2, null);
                                    r.Level = parts.ToInt(-1, 0);
                                    r.Z = parts.ToFloat(-2, 0);
                                    Nodes.Add(r);
                                }
                                break;

                            case "TANKS":
                                {
                                    TankViewModel tank = new TankViewModel();
                                    o = tank;
                                    tank.ID = parts.ToString(0, null);
                                    tank.Z = parts.ToFloat(1, 0);
                                    tank.InitLevel = parts.ToFloat(2, 0);
                                    tank.MinLevel = parts.ToFloat(3, 0);
                                    tank.MaxLevel = parts.ToFloat(4, 0);
                                    tank.Diameter = parts.ToFloat(5, 0);
                                    tank.MinVol = parts.ToFloat(6, 0);
                                    tank.VolCurve = parts.ToString(7, null);
                                    tank.Overflow = "";
                                    tank.Level = parts.ToInt(-1, 0);
                                    Nodes.Add(tank);
                                }
                                break;


                            case "VALVES":
                                {
                                    ValveViewModel valve = new ValveViewModel();
                                    o = valve;
                                    valve.ID = parts.ToString(0, null);

                                    // 取出Node1和Node2中的字母部分，例如“S201326593”被取出为“201326593”
                                    valve.Node1 = parts.ToString(1, null); // Regex.Replace(parts[1], "[^0-9]", "");
                                    valve.Node2 = parts.ToString(2, null); // Regex.Replace(parts[2], "[^0-9]", "");
                                    valve.Diameter = parts.ToFloat(3, 0);
                                    valve.Type = parts.ToString(4, null);
                                    valve.Setting = parts.ToString(5, null);
                                    valve.MinorLoss = parts.ToFloat(6, 0);
                                    valve.Level = parts.ToInt(-1, 0);
                                    Links.Add(valve);
                                }
                                break;

                            case "PUMPS":
                                {
                                    PumpViewModel pump = new PumpViewModel();
                                    o = pump;
                                    pump.ID = parts.ToString(0, null);

                                    // 取出Node1和Node2中的字母部分，例如“S201326593”被取出为“201326593”
                                    pump.Node1 = parts.ToString(1, null); // Regex.Replace(parts[1], "[^0-9]", "");
                                    pump.Node2 = parts.ToString(2, null); // Regex.Replace(parts[2], "[^0-9]", "");
                                    int index = 3;
                                    string label = null;
                                    while ((label = parts.ToString(index, null)) != null)
                                    {
                                        label = label.ToUpper();
                                        switch (label)
                                        {
                                            case "HEAD":
                                                pump.HeadCurve = parts.ToString(index + 1, "PumpDefault");
                                                break;

                                            case "SPEED":
                                                pump.当前转速 = parts.ToFloat(index + 1, 0);
                                                break;
                                        }
                                        index += 2;
                                    }
                                    //pump.Diameter = parts.ToFloat(3, 0);

                                    pump.Level = parts.ToInt(-1, 0);
                                    Links.Add(pump);
                                }
                                break;

                            case "CURVES":
                                {
                                    string ID = parts.ToString(0, "");
                                    if (!dict_dataset.ContainsKey(ID))
                                    {
                                        var ds = new Dataset(ID, null);
                                        dict_dataset.Add(ID, ds);
                                    }
                                    dict_dataset[ID]._data.Add(new PointF(parts.ToFloat(1, 0), parts.ToFloat(2, 0)));
                                }
                                break;

                            case "COORDINATES":
                                {
                                    string id = parts[0];
                                    float x;
                                    float y;
                                    if (float.TryParse(parts[1], out x) && float.TryParse(parts[2], out y))
                                    {
                                        points.Add(new Coor(id, new PointF(x, y)));
                                    }
                                }
                                break;

                            case "STATUS":
                                {
                                    LinkViewModel link = Links.Find(l => l.ID == parts.ToString(0, null));
                                    if (link != null)
                                        link.Status = StringToStatus(parts.ToString(1, "OPEN"));
                                }
                                break;
                        }
                        if (o != null)
                        {
                            // o.Tags = new TagList(parts.ToString(-3, null));
                        }
                    }
                }
                sr.Close();

                if (!dict_dataset.ContainsKey("GPVDefault"))
                {
                    var data = new Dataset("GPVDefault", null);
                    data.Data = new List<PointF>()
                    {
                        new PointF(0,0),
                        new PointF(100,0),
                    };
                    dict_dataset.Add("GPVDefault", data);
                }
                if (!dict_dataset.ContainsKey("PumpDefault"))
                {
                    var data = new Dataset("PumpDefault", null);
                    data.Data = new List<PointF>()
                    {
                        new PointF(0f            ,   45.38f),
                        new PointF(83.33333333f     ,45.25f),
                        new PointF(111.1111111f     ,45.12f),
                        new PointF(138.8888889f     ,44.96f),
                        new PointF(166.6666667f     ,44.76f),
                        new PointF(194.4444444f     ,44.52f),
                        new PointF(222.2222222f     ,44.24f),
                        new PointF(250f             ,43.92f),
                        new PointF(277.7777778f     ,43.56f),
                        new PointF(305.5555556f     ,43.17f),
                        new PointF(333.3333333f     ,42.73f),
                        new PointF(361.1111111f     ,42.25f),
                        new PointF(388.8888889f     ,41.74f),
                        new PointF(416.6666667f     ,41.18f),
                        new PointF(444.4444444f     ,40.58f),
                        new PointF(472.2222222f     ,39.95f),
                        new PointF(500f             ,39.28f),
                        new PointF(527.7777778f     ,38.56f),
                        new PointF(555.5555556f     ,37.81f),
                        new PointF(583.3333333f     ,37.02f),
                        new PointF(611.1111111f     ,36.19f),
                        new PointF(638.8888889f     ,35.32f),
                        new PointF(666.6666667f     ,34.41f),
                        new PointF(694.4444444f     ,33.46f),
                        new PointF(722.2222222f     ,32.47f),
                        new PointF(750f             ,31.44f),
                        new PointF(777.7777778f     ,30.37f),
                        new PointF(805.5555556f     ,29.27f),
                    };
                    dict_dataset.Add("PumpDefault", data);
                }

                //读取坐标
                int k1 = 0;
                int k2 = 0;
                Nodes.Sort((a, b) => string.Compare(a.ID, b.ID));
                points.Sort((a, b) => string.Compare(a.ID, b.ID));
                k1 = 0;
                k2 = 0;
                while (k1 < Nodes.Count)
                {
                    var J = Nodes[k1];
                    var coor = points[k2];

                    while (J.ID != coor.ID && k2 < points.Count)
                    {
                        k2++;
                        if (k2 < points.Count) coor = points[k2];
                    }
                    if (k2 == points.Count)
                    {
                        throw new Exception($"未找到Node[{J.ID}]的坐标");
                    }
                    J.X = coor.Position.X;
                    J.Y = coor.Position.Y;

                    k1++;
                }

                BuildRelation();

                return true;
            }
        }

        private StatusType StringToStatus(string status)
        {
            switch (status)
            {
                case "CLOESD":
                case "0":
                    return StatusType.CLOSED;
                    break;

                case "OPEN":
                case "1":
                    return StatusType.OPEN;
                    break;

                case "ACTIVE":
                    return StatusType.ACTIVE;
                    break;

                default:
                    return StatusType.DEFAULT;
                    break;
            }
        }


    }

    public class Vector3D
    {
        public double X { get; set; }
        public double Y { get; set; }
        public double Z { get; set; }

        public Vector3D(double x, double y, double z)
        {
            X = x;
            Y = y;
            Z = z;
        }

        public static double CalculateAngle(Vector3D v1, Vector3D v2)
        {
            // 计算两个向量的点积
            double dotProduct = v1.X * v2.X + v1.Y * v2.Y + v1.Z * v2.Z;

            // 计算两个向量的模长
            double v1Magnitude = Math.Sqrt(v1.X * v1.X + v1.Y * v1.Y + v1.Z * v1.Z);
            double v2Magnitude = Math.Sqrt(v2.X * v2.X + v2.Y * v2.Y + v2.Z * v2.Z);

            // 计算夹角（弧度）
            double angleInRadians = Math.Acos(dotProduct / (v1Magnitude * v2Magnitude));

            // 将弧度转换为度数
            //double angleInDegrees = angleInRadians * (180 / Math.PI);

            return angleInRadians;
        }

        public double Length
        {
            get
            {
                return Math.Sqrt(X * X + Y * Y + Z * Z);
            }
        }
    }
}