using DevExpress.XtraEditors;
using NPOI.Util;
 
namespace Yw.WinFrmUI.Phart
{
    /// <summary>
    /// (待修改)
    /// </summary>
    public partial class PumpImportByOnePointCtrl : XtraUserControl
    {

        //构造函数
        public PumpImportByOnePointCtrl()
        {
            InitializeComponent();
            this.comboBoxType.SelectedIndex = 0;
            this.txtQ.Focus();
            this.txtQ.SelectAll();

        }

        double _ratedn = -1;
        PumpQueryPointViewModel _work_point = null;
        List<Yw.Geometry.Point2d> _qh_pt_list = null;
        List<Yw.Geometry.Point2d> _qe_pt_list = null;
        List<Yw.Geometry.Point2d> _qp_pt_list = null;
        List<Yw.Geometry.Point2d> _npsh_pt_list = null;


        //计算转速比
        private void CalcNs()
        {
            if (string.IsNullOrEmpty(this.txtQ.Text) ||
                string.IsNullOrEmpty(this.txtH.Text) ||
                string.IsNullOrEmpty(this.txtN.Text))
                return;

            var YLJS = Convert.ToInt32(this.spinEditYLJS.Text);
            var Q = Convert.ToDouble(this.txtQ.Text);
            var H = Convert.ToDouble(this.txtH.Text);
            var N = Convert.ToDouble(this.txtN.Text);

            double ns = 0;
            if (this.checkEditSXB.Checked)
            {
                ns = PumpCalcHelper.CalculateNs(Q / 2, H / YLJS, N);
            }
            else
            {
                ns = PumpCalcHelper.CalculateNs(Q, H / YLJS, N);
            }

            this.txtNS.EditValue = Math.Round(ns, 2);
        }

        //产生曲线
        private bool CreateCurveCore()
        {
            double? k0 = 0;
            if (double.TryParse(this.txtK0.Text, out double value))
            {
                k0 = Convert.ToDouble(txtK0.Text);
                if (k0 < 1.001)
                {
                    XtraMessageBox.Show("关死点系数不能小于1");
                    return false;
                }
            }

            if (string.IsNullOrEmpty(this.txtQ.Text) || string.IsNullOrEmpty(this.txtH.Text))
            {
                MessageBox.Show("流量扬程不能为空");
                return false;
            }
            if (string.IsNullOrEmpty(this.txtN.Text))
            {
                MessageBox.Show("转速不能为空");
                return false;
            }

            _ratedn = double.Parse(this.txtN.Text);
            if (_ratedn < 10)
            {
                MessageBox.Show("转速有这么小吗?");
                return false;
            }

            _work_point = new PumpQueryPointViewModel();
            _work_point.Q = double.Parse(this.txtQ.Text);
            _work_point.H = double.Parse(this.txtH.Text);
            if (_work_point.H < 1 || _work_point.Q < 1)
            {
                MessageBox.Show("流量扬程有这么小吗?");
                return false;
            }

            if (string.IsNullOrWhiteSpace(this.txtE.Text))
            {
                double ns = PumpCalcHelper.CalculateNs(_work_point.Q, _work_point.H, _ratedn);
                double eff = PumpCalcHelper.CalculateEByPumpType(_work_point.Q, _work_point.H, _ratedn, comboBoxType.SelectedIndex + 1);

                _work_point.E = eff;
                this.txtE.EditValue = Math.Round(eff, 1);
            }
            else
            {
                _work_point.E = double.Parse(this.txtE.Text);
            }
            if (_work_point.E > 99)
            {
                MessageBox.Show("效率有这么大吗?");
                return false;
            }
            _work_point.P = PumpCalcHelper.CalculateP(_work_point.Q, _work_point.H, _work_point.E.Value);
            if (string.IsNullOrWhiteSpace(this.txtP.Text))
            {
                this.txtP.EditValue = Math.Round(_work_point.P.Value, 1);
            }
            //汽蚀
            double NPSHr = 0;
            if (!double.TryParse(this.txtNPSH.Text, out NPSHr))
            {
                NPSHr = PumpCalcHelper.CalculateNPSHrByPumpType(_work_point.Q, _work_point.H, _ratedn, comboBoxType.SelectedIndex + 1);//双吸离心泵
                NPSHr = Math.Round(NPSHr, 1);
                this.txtNPSH.EditValue = NPSHr;
            }
            //_work_point.NPSH = NPSHr;

            //得到点位置

            bool iOK = PumpDimensionlessCurvesHelper.CalcPoints(_work_point, _ratedn, k0,
                ref _qh_pt_list, ref _qe_pt_list, ref _qp_pt_list, ref _npsh_pt_list, checkEditSXB.Checked, Convert.ToInt32(spinEditYLJS.Value));
            if (!iOK)
            {
                XtraMessageBox.Show("Error:265");
                return false;
            }

            //消除驼峰
            if (k0 == null || k0 < 0.1)
            {
                double spaceH = _work_point.H / 200;//200为预估值

                for (int i = _qh_pt_list.Count - 2; i >= 0; i--)
                {
                    if (_qh_pt_list[i].Y < _qh_pt_list[i + 1].Y)
                    {
                        _qh_pt_list[i].Y = _qh_pt_list[i + 1].Y + spaceH;
                    }
                }
            }

            if (!string.IsNullOrEmpty(txtMaxRatio.Text))
            {
                double qre = 0;
                if (double.TryParse(txtMaxRatio.Text, out qre))
                {
                    if (qre > 1.01 && qre < 1.8)
                    {
                        var ex_q = _work_point.Q * qre;
                        var max_q = _qh_pt_list.Last().X;
                        var ratio = ex_q / max_q;
                        if (qre > 1.35)
                        {
                            var pointInfoQH2 = GetFitPointListByExtend(_qh_pt_list, ratio);
                            var max_powert_E = GetLineInsert(_qe_pt_list[_qe_pt_list.Count - 2].X, _qe_pt_list[_qe_pt_list.Count - 1].X, _qe_pt_list[_qe_pt_list.Count - 2].Y, _qe_pt_list[_qe_pt_list.Count - 1].Y, qre * _work_point.Q);
                            _qe_pt_list.Add(new Yw.Geometry.Point2d(qre * _work_point.Q, max_powert_E));
                            _qh_pt_list = pointInfoQH2;
                            _qe_pt_list = _qe_pt_list.GetFitPointList(20);
                            _qp_pt_list = CalculateP(_qh_pt_list, _qe_pt_list);
                        }
                        else
                        {
                            _qh_pt_list = GetFitPointListByExtend(_qh_pt_list, ratio);
                            _qe_pt_list = GetFitPointListByExtend(_qe_pt_list, ratio);
                            _qp_pt_list = GetFitPointListByExtend(_qp_pt_list, ratio);
                        }
                    }
                }
            }

            var qh_pt_list = _qh_pt_list.ToList();
            var qe_pt_list = _qe_pt_list.ToList();
            var qp_pt_list = _qp_pt_list.ToList();

            var qhCurve = new Yw.Geometry.CubicSpline2d(qh_pt_list);
            var qeCurve = new Yw.Geometry.CubicSpline2d(qe_pt_list);
            var qPCurve = new Yw.Geometry.CubicSpline2d(qp_pt_list);

            this.xtrPerform2dChart1.SetBindingData(qhCurve, qeCurve, qPCurve);
            return true;
        }


        public List<Yw.Geometry.Point2d> GetFitPointListByExtend(List<Yw.Geometry.Point2d> points, double ratioExtend, int pointNumber = 20)
        {
            if (points == null)
                return default;

            var express = new Yw.Geometry.CubicSpline2d(points);
            return express.GetPointListByXRatioRange(1, ratioExtend);
        }

        public double GetLineInsert(double X1, double X2, double Y1, double Y2, double X)
        {
            if (Math.Abs(X2 - X1) < 1E-06)
            {
                return 0.0;
            }

            if (Math.Abs(Y2 - Y1) < 1E-06)
            {
                return 0.0;
            }

            return Y1 + (Y2 - Y1) / (X2 - X1) * (X - X1);
        }



        public void Set(double flow, double head, double speed, double? eff, double level, bool is_sxp)
        {
            this.txtQ.EditValue = flow;
            this.txtH.EditValue = head;
            this.txtN.EditValue = speed;
            this.txtE.EditValue = eff;
            this.spinEditYLJS.EditValue = level;
            this.checkEditSXB.Checked = is_sxp;
        }

        /// <summary>
        /// 获取数据
        /// </summary>
        /// <returns></returns>
        public bool Get(out string other_name, out Yw.Pump.CurveQH qh, out Yw.Pump.CurveQE qe, out Yw.Pump.CurveQP qp)
        {
            other_name = this.txtCurveCode.Text.Trim();
            qh = null;
            qe = null;
            qp = null;

            if (!CreateCurveCore())
                return false;

            if (string.IsNullOrEmpty(other_name))
            {
                XtraMessageBox.Show("请输入曲线名称!");
                return false;
            }
            qh = new Pump.CurveQH(Yw.Pump.eFeatType.Cubic, _qh_pt_list);
            qe = new Pump.CurveQE(Yw.Pump.eFeatType.Cubic, _qe_pt_list);
            qp = new Pump.CurveQP(Yw.Pump.eFeatType.Cubic, _qp_pt_list);

            return true;
        }


        #region 事件

        private void txtQ_EditValueChanged(object sender, EventArgs e)
        {
            CalcNs();
        }

        private void txtH_EditValueChanged(object sender, EventArgs e)
        {
            CalcNs();
        }

        private void txtN_EditValueChanged(object sender, EventArgs e)
        {
            CalcNs();
        }

        private void spinEditYLJS_EditValueChanged(object sender, EventArgs e)
        {
            CalcNs();
        }

        private void labReckon_Click(object sender, EventArgs e)
        {
            //水力参数
            try
            {
                if (string.IsNullOrEmpty(txtH.Text) || string.IsNullOrEmpty(txtQ.Text))
                {
                    MessageBox.Show("流量扬程不能为空");
                    return;
                }
                //转速
                if (string.IsNullOrEmpty(txtN.Text))
                {
                    MessageBox.Show("转速不能为空");
                    return;
                }
                _ratedn = double.Parse(txtN.Text);
                if (_ratedn < 10)
                {
                    MessageBox.Show("转速有这么小吗?");
                    return;
                }
                _work_point = new();
                _work_point.Q = double.Parse(txtQ.Text);
                _work_point.H = double.Parse(txtH.Text);
                if (_work_point.H < 1 || _work_point.Q < 1)
                {
                    MessageBox.Show("流量扬程有这么小吗?");
                    return;
                }

            }
            catch
            {
                MessageBox.Show("请输入数字");
                return;
            }


            double eff = PumpCalcHelper.CalculateEByPumpType(_work_point.Q, _work_point.H, _ratedn, comboBoxType.SelectedIndex + 1);
            txtE.Text = string.Format("{0:0.0}", eff);
            _work_point.E = eff;

            double NPSHr = PumpCalcHelper.CalculateNPSHrByPumpType(_work_point.Q, _work_point.H, _ratedn, comboBoxType.SelectedIndex + 1);//双吸离心泵

            NPSHr = Math.Round(NPSHr, 1);
            txtNPSH.Text = NPSHr.ToString();
        }

        private void labCalcuE_Click(object sender, EventArgs e)
        {
            #region 得到参数
            double Q, H, P;
            try
            {
                if (string.IsNullOrEmpty(txtH.Text) || string.IsNullOrEmpty(txtQ.Text) || string.IsNullOrEmpty(txtE.Text))
                {
                    MessageBox.Show("流量,扬程和功率,不能为空");
                    return;
                }
                H = double.Parse(txtH.Text);
                Q = double.Parse(txtQ.Text);
                P = double.Parse(txtP.Text);
                if (H < 1 || Q < 1 || P < 1)
                {
                    MessageBox.Show("参数有这么小吗");
                    return;
                }
            }
            catch
            {
                MessageBox.Show("请输入数字");
                return;
            }
            #endregion

            double E = PumpCalcHelper.CalculateE(Q, H, P);
            txtE.Text = string.Format("{0:0.0}", E);

            if (E > 99)
            {
                MessageBox.Show("功率是否输入有误?");
                return;
            }
        }

        private void labCalcuP_Click(object sender, EventArgs e)
        {
            #region 得到参数
            double Q, H, E;
            try
            {
                if (string.IsNullOrEmpty(txtH.Text) || string.IsNullOrEmpty(txtQ.Text) || string.IsNullOrEmpty(txtE.Text))
                {
                    MessageBox.Show("流量,扬程和功率,不能为空");
                    return;
                }
                H = double.Parse(txtH.Text);
                Q = double.Parse(txtQ.Text);
                E = double.Parse(txtE.Text);
                if (H < 0.1 || Q < 0.1 || E < 1)
                {
                    MessageBox.Show("参数有这么小吗");
                    return;
                }
            }
            catch
            {
                MessageBox.Show("请输入数字");
                return;
            }
            #endregion

            double P = PumpCalcHelper.CalculateP(Q, H, E);
            if (P < 1)
            {
                txtP.Text = string.Format("{0:0.000}", P);
            }
            else if (P < 10)
            {
                txtP.Text = string.Format("{0:0.00}", P);
            }
            else if (P < 100)
            {
                txtP.Text = string.Format("{0:0.0}", P);
            }
            else
            {
                txtP.Text = string.Format("{0:N0}", P);
            }
        }

        //生成曲线
        private void btnCreate_Click(object sender, EventArgs e)
        {
            CreateCurveCore();
        }

        #endregion


        public static List<Yw.Geometry.Point2d> CalculateP(List<Yw.Geometry.Point2d> pointQH, List<Yw.Geometry.Point2d> pointQE)
        {
            if (pointQH == null || pointQH.Count() <= 2)
                return null;
            bool isFromZero = true;
            if (pointQH.First().X > 500 || pointQH.First().X > pointQH.Last().X * 0.2)
                isFromZero = false;
            return CalculateP_RefZero(pointQH, pointQE, Constant.WaterDensity, -1, isFromZero);
        }

        //有时 会把零流量点扬程传入 ref_zero_power
        public static List<Yw.Geometry.Point2d> CalculateP_RefZero(
            List<Yw.Geometry.Point2d> pointQH,
            List<Yw.Geometry.Point2d> pointQE,
            double midu,
            double ref_zero_power,
            bool isFromZero)
        {
            if (pointQH == null || pointQH.Count < 3)
                return null;
            if (pointQE == null || pointQE.Count < 3)
                return null;

            var expressQH = new Yw.Geometry.CubicSpline2d(pointQH);
            var expressQE = new Yw.Geometry.CubicSpline2d(pointQE);

            var expressQP = CalculateP_RefZero(expressQH, expressQE, midu, ref_zero_power, isFromZero, false);

            return expressQP.GetPointList(pointQE.Count());
        }

        public static Yw.Geometry.CubicSpline2d CalculateP_RefZero(
         Yw.Geometry.CubicSpline2d expressQH,
         Yw.Geometry.CubicSpline2d expressQE,
         double midu, double ref_zero_power, bool isFromZero, bool isSavePointInExpress)
        {
            if (expressQH == null)
                return null;
            if (expressQE == null)
                return null;
            int point_num = 16;//点不能太多
            var pointQH = expressQH.GetPointList(point_num);
            if (pointQH == null || pointQH.Count < 3)
                return null;
            var pointQE = expressQE.GetPointList(point_num);
            if (pointQE == null || pointQE.Count < 3)
                return null;

            List<Yw.Geometry.Point2d> pointQP = new List<Yw.Geometry.Point2d>();
            //间距
            double min_flow = Math.Max(pointQH.First().X, pointQE.First().X);
            double max_flow = Math.Min(pointQH.Last().X, pointQE.Last().X);

            //
            double Q, H, E, P;

            double space = (max_flow - min_flow) / (point_num - 1);

            var fitCurveQE = expressQE.Copy();
            var fitCurveQH = expressQH.Copy();

            if (isFromZero)
            {
                #region 保证从0开始:  保证Q=0时,P不是计算出来的;
                for (int i = 5; i < point_num; i++)//前面几个点不要
                {
                    Q = pointQE.First().X + space * i;
                    E = fitCurveQE.GetPointY(Q);
                    H = fitCurveQH.GetPointY(Q);
                    if (H < 0.1 || E < 0.5)
                        continue;

                    P = PumpCalcHelper.CalculateP(Q, H, E, midu);
                    if (P < 2)
                    {
                        P = Math.Round(P, 3);
                    }
                    else if (P < 30)
                    {
                        P = Math.Round(P, 2);
                    }
                    else if (P < 100)
                    {
                        P = Math.Round(P, 1);
                    }
                    else
                    {
                        P = Math.Round(P, 0);
                    }
                    pointQP.Add(new Yw.Geometry.Point2d(Q, P));
                }

                pointQP = pointQP.GetFitPointList(point_num);
                if (pointQP == null)
                    return null;



                if (ref_zero_power > 0.1)
                {
                    pointQP.Insert(0, new Yw.Geometry.Point2d(0, ref_zero_power));
                }
                else
                {
                    if (pointQP[0].Y >= pointQP[1].Y)
                    {
                        pointQP[0].Y = pointQP[1].Y * 0.95;
                    }
                    double startP = Yw.Geometry.LineHelper.GetYbyX(pointQP[0].X, pointQP[0].Y, pointQP[1].X, pointQP[1].Y, 0);
                    if (startP < 0.001)
                        startP = pointQP[0].Y;

                    if (startP < 2)
                    {
                        startP = Math.Round(startP, 3);
                    }
                    else if (startP < 30)
                    {
                        startP = Math.Round(startP, 2);
                    }
                    else if (startP < 100)
                    {
                        startP = Math.Round(startP, 1);
                    }
                    else
                    {
                        startP = Math.Round(startP, 0);
                    }

                    pointQP.Insert(0, new Yw.Geometry.Point2d(0, startP));
                }


                return new Yw.Geometry.CubicSpline2d(pointQP);
                #endregion
            }
            else
            {
                #region 保证从0开始:  保证Q=0时,P不是计算出来的;
                for (int i = 0; i < point_num; i++)//前面几个点不要
                {
                    Q = pointQE.First().X + space * i;
                    E = fitCurveQE.GetPointY(Q);
                    H = fitCurveQH.GetPointY(Q);
                    if (H < 0.1 || E < 0.5)
                        continue;

                    P = PumpCalcHelper.CalculateP(Q, H, E, midu);
                    if (P < 2)
                    {
                        P = Math.Round(P, 3);
                    }
                    else if (P < 30)
                    {
                        P = Math.Round(P, 2);
                    }
                    else if (P < 100)
                    {
                        P = Math.Round(P, 1);
                    }
                    else
                    {
                        P = Math.Round(P, 0);
                    }
                    pointQP.Add(new Yw.Geometry.Point2d(Q, P));
                }

                pointQP = pointQP.GetFitPointList(point_num);
                if (pointQP == null)
                    return null;

                return new Yw.Geometry.CubicSpline2d(pointQP);
                #endregion
            }

        }

    }
}