namespace IStation.Curve
{
    public static class CurveExpress_Extensions
    {

        #region CurveFunction

        /// <summary>
        /// 获取曲线方程
        /// </summary>
        public static string GetCurveFunction(this CurveExpress express, string xName = "Q")
        {

            StringBuilder sb = new StringBuilder();
            sb.Append(Math.Abs(express.Index0).ToString("0.00"));

            if (express.Index1 > 0)
                sb.Append(" + ");
            else
                sb.Append(" - ");
            sb.Append(Math.Abs(express.Index1).ToString("0.0000"));
            sb.Append(" * ");
            sb.Append(xName);

            if (express.Index2 > 0)
                sb.Append(" + ");
            else
                sb.Append(" - ");
            sb.Append(Math.Abs(express.Index2).ToString("0.000000000"));
            sb.Append(" * ");
            sb.Append(xName);
            sb.Append("^2 ");


            if (express.Index3 > 0)
                sb.Append("  +  ");
            else
                sb.Append("  -  ");
            sb.Append(Math.Abs(express.Index3).ToString("0.000000000000"));
            sb.Append(" * ");
            sb.Append(xName);
            sb.Append("^3 ");

            return sb.ToString();
        }

        /// <summary>
        /// 获取
        /// </summary>
        /// <param name="functionIndex"></param>
        /// <param name="functionCoeff"></param>
        /// <param name="signValue"></param>
        /// <param name="IndexValue"></param>
        private static void GetFuncScientificValue(int functionIndex, double functionCoeff, ref double signValue, ref int IndexValue)
        {
            string strList = "";
            if (functionIndex == 0)
            {
                strList = functionCoeff.ToString("0.000e+000");
            }
            else if (functionIndex == 1)
            {
                strList = functionCoeff.ToString("0.000000e+000");
            }
            else if (functionIndex == 2)
            {
                strList = functionCoeff.ToString("#.#########e+000");
            }
            else if (functionIndex == 3)
            {
                strList = functionCoeff.ToString("#.############e+000");
            }
            string[] st1 = null;
            st1 = strList.Split('e');
            signValue = Convert.ToDouble(st1[0]);
            IndexValue = Convert.ToInt32(st1[1]);
        }

        #endregion

        #region  List<CurvePoint>  Extend

        /// <summary>
        /// 获取拟合点
        /// </summary>
        public static List<CurvePoint> GetFitPointsByExtend(this CurveExpress express, double ratioExtend, int pointNumber)
        {
            if (express == null)
                return default;
            var helper = new FitHelper(express);
            return helper.GetFitPointsByXRange(express.Min, express.Max * ratioExtend, pointNumber);
        }

        /// <summary>
        /// 获取拟合点
        /// </summary>
        public static List<CurvePoint> GetFitPointsByExtend(this CurveExpress express, double ratioExtendmin, double ratioExtendmax, int pointNumber)
        {
            if (express == null)
                return default;
            var helper = new FitHelper(express);
            return helper.GetFitPointsByXRange(express.Min * ratioExtendmin, express.Max * ratioExtendmax, pointNumber);
        }

        #endregion

        #region 获得X值

        /// <summary>
        /// 获得Y值的坐标位置,多个点里取最后一个点(是线形插值,不是曲线上的点)
        /// </summary>
        /// <param name="curve"></param>
        /// <param name="y"></param>
        /// <param name="interPointNum"></param>
        /// <returns></returns>
        public static double? GetInterPointLastX(this CurveExpress curve, double y, int interPointNum = 20)
        {
            if (curve == null)
                return null;
            var points = FitHelper.GetFitPoints(curve, interPointNum);
            if (points == null)
                return null;

            //在范围内
            double maxY = (from yyy in points select yyy.Y).Max();
            double minY = (from yyy in points select yyy.Y).Min();

            if (double.IsNaN(maxY) || double.IsNaN(minY))
                return null;
            if (y < minY || y > maxY)
                return null;

            var list = points.GetPointsByInterPointX(y);
            if (list == null || list.Count < 1)
                return null;

            return list.Last().X;
        }


        /// <summary>
        /// 获得Y值的坐标位置,可能有多个点(是线形插值,不是曲线上的点):会计算延长部分:而且要求points点按照从小到大排序
        /// </summary>
        /// <param name="curve"></param>
        /// <param name="y"></param>
        /// <param name="interPointNum"></param>
        /// <returns></returns>
        public static List<CurvePoint> GetInterPointExX(this CurveExpress curve, double y, int interPointNum)
        {
            if (curve == null)
                return null;
            var points = FitHelper.GetFitPoints(curve, interPointNum);
            if (points == null)
                return null;

            //在范围内
            double maxY = (from yyy in points select yyy.Y).Max();
            double minY = (from yyy in points select yyy.Y).Min();

            if (double.IsNaN(maxY) || double.IsNaN(minY))
                return null;
            if (y > minY && y < maxY)
            {
                return points.GetPointsByInterPointX(y);
            }

            List<CurvePoint> ptList = new List<CurvePoint>();

            //判断是否在前面的延长部分
            double x_start = CurveLineHelper.GetXbyY(points[0].X, points[0].Y, points[1].X, points[1].Y, y);
            if (x_start <= points[0].X)
            {
                if (x_start > 0 && x_start > points[0].X * 0.9)//放一点余量
                    ptList.Add(new CurvePoint(x_start, y));
            }

            //判断是否在后面的延长部分:U型曲线 是前后都有
            double x_end = CurveLineHelper.GetXbyY(points[points.Count - 2].X, points[points.Count - 2].Y, points[points.Count - 1].X, points[points.Count - 1].Y, y);
            if (x_end >= points[points.Count - 1].X && x_end < points[points.Count - 1].X * 1.1)//放一点余量
            {
                ptList.Add(new CurvePoint(x_end, y));
            }

            if (ptList.Count == 0)
                return null;
            return ptList;
        }



        /// <summary>
        /// 获得Y值的坐标位置,可能有多个点(是线形插值,不是曲线上的点):会计算延长部分:而且要求points点按照从小到大排序
        /// </summary>
        /// <param name="curve"></param>
        /// <param name="y"></param>
        /// <param name="interPointNum"></param>
        /// <param name="extendRatio"></param>
        /// <returns></returns>
        public static List<CurvePoint> GetInterPointExX(this CurveExpress curve, double y, int interPointNum, double extendRatio)
        {
            if (curve == null)
                return null;
            var points = GetFitPointsByExtend(curve, extendRatio, interPointNum);
            if (points == null)
                return null;

            //在范围内
            double maxY = (from yyy in points select yyy.Y).Max();
            double minY = (from yyy in points select yyy.Y).Min();
            if (y > minY && y < maxY)
            {
                return points.GetPointsByInterPointX(y);
            }

            List<CurvePoint> ptList = new List<CurvePoint>();
            //判断是否在前面的延长部分
            double x_start = CurveLineHelper.GetXbyY(points[0].X, points[0].Y, points[1].X, points[1].Y, y);
            if (x_start <= points[0].X)
            {
                if (x_start > 0 && x_start > points[0].X * 0.9)//放一点余量
                    ptList.Add(new CurvePoint(x_start, y));
            }

            //判断是否在后面的延长部分:U型曲线 是前后都有
            double x_end = CurveLineHelper.GetXbyY(points[points.Count - 2].X, points[points.Count - 2].Y, points[points.Count - 1].X, points[points.Count - 1].Y, y);
            if (x_end >= points[points.Count - 1].X && x_end < points[points.Count - 1].X * extendRatio)//放一点余量
            {
                ptList.Add(new CurvePoint(x_end, y));
            }

            if (ptList.Count == 0)
                return null;
            return ptList;
        }

        #endregion 获得X值

        #region 获得Y值

        /// <summary>
        /// 获取Y值最大值的点
        /// </summary> 
        public static CurvePoint GetMaxPoint(this CurveExpress express, int pointNumber = 200)
        {
            if (express == null)
                return null;
            if (express.FitType <= 0 && express.DefinePoints != null)
            {
                return express.DefinePoints.GetPointByMaxY();
            }
            else if (express.FitType == eFitType.ThroughPoint)
            {
                return express.DefinePoints.GetPointByMaxY();
            }
            else
            {
                var points = express.GetFitPoints(pointNumber);
                return points.GetPointByMaxY();
            }
        }

        /// <summary>
        /// 获取Y值最大值的点
        /// </summary> 
        public static CurvePoint GetMaxPointByRange(this CurveExpress express, double rangeMin, double rangeMax, int pointNumber = 200)
        {
            if (express == null)
                return null;
            var points = express.GetFitPointsByXRange(rangeMin, rangeMax, pointNumber);
            return points.GetPointByMaxY();
        }


        /// <summary>
        /// 获取最大最小Y
        /// </summary>
        /// <param name=""></param> 
        public static bool GetMinMaxPointY(this CurveExpress express, out double maxY, out double minY)
        {
            if (express == null)
            {
                maxY = minY = 0;
                return false;
            }
            var fitPoints = express.GetFitPoints(50);
            if (fitPoints == null)
            {
                maxY = minY = 0;
                return false;
            }
            maxY = fitPoints.Max(x => x.Y);
            minY = fitPoints.Min(x => x.Y);
            return true;
        }

        #endregion 取Y的最大最小值

        #region SimuValue

        /// <summary>
        /// 根据相似点设置相似曲线,相似点(选择点)的x:流量y:杨程
        /// </summary> 
        public static double GetSimuValue(this CurveExpress express, CurvePoint simularPoint, double originValue, double extendRatio = 1)
        {
            if (express == null)
                return -3;
            if (simularPoint.X < 0.1 || simularPoint.Y < 0.1)
                return -2;
            if (simularPoint.X > express.Max * extendRatio * 1.5)
                return -4;

            CurvePoint sectPoint = GetSectPoint(express, simularPoint, extendRatio);
            if (sectPoint == null || sectPoint.IsZeroPoint())
                return -5;

            //计算相似点的转速/直径
            return PumpCalculateHelper.CalculateSimuByH(originValue, sectPoint.Y, simularPoint.Y);
        }

        //ratioIgnore:作用 当simularPoint超出曲线范围时,曲线扩大的倍数
        public static CurvePoint GetSectPoint(this CurveExpress CurveExpress, CurvePoint simularPoint, double ratioIgnore)
        {
            return GetSectPointGeneral(CurveExpress, simularPoint, 2, ratioIgnore);
        }

        #endregion

        #region 曲线与H=K*Q^i的交点  


        //通过点simularPoint和点(0,0)的直线,与曲线Curve的交点(没有,返回Point(0,0))
        public static CurvePoint GetSectPointLine(List<CurvePoint> CurvePoints, CurvePoint simularPoint)
        {
            CurvePoint sectPoint = new CurvePoint(0, 0);
            if (CurvePoints == null || CurvePoints.Count < 2)
                return sectPoint;

            //计算直线的K
            if (simularPoint.X < 1)
                return sectPoint;
            double a = simularPoint.Y / simularPoint.X;
            if (a < 0.0001)
                return sectPoint;

            //与2点连成直线的交点,判断交点是否在2点之间,即可是曲线的交点
            double b, c;
            double x;
            for (int i = 0; i < CurvePoints.Count - 1; i++)
            {
                CurveLineHelper.GetKandB(CurvePoints[i], CurvePoints[i + 1], out b, out c);

                /*解方程
                 * y=ax
                 * y=bx+c
                 */
                if (Math.Abs(a - b) < 0.001)
                    continue;

                x = c / (a - b);
                if (IStation.Curve.UtilsHelper.IsMiddle(CurvePoints[i].X, CurvePoints[i + 1].X, x))
                {
                    sectPoint.X = x;
                    sectPoint.Y = a * x;
                    return sectPoint;
                }
            }

            return sectPoint;
        }

        //通过点simularPoint和点(0,0)的抛物线,与曲线Curve的交点(没有,返回Point(0,0))
        //曲线公式:H=K*Q^2
        public static CurvePoint GetSectPointParabola(List<CurvePoint> CurvePoints, CurvePoint simularPoint)
        {
            return ParabolaCurveHelper.GetSectPoint(CurvePoints, simularPoint, 0);
        }

        public static CurvePoint GetSectPointParabola(CurveExpress CurveExpress, CurvePoint simularPoint)
        {
            CurvePoint sectPoint = new CurvePoint(0, 0);
            if (CurveExpress == null)
                return sectPoint;

            return GetSectPointParabola(CurveExpress.GetFitPointsByExtend(1.2, 50), simularPoint);
        }

        public static CurvePoint GetSectPointLine(this CurveExpress CurveExpress, CurvePoint simularPoint)
        {
            CurvePoint sectPoint = new CurvePoint(0, 0);
            if (CurveExpress == null)
                return sectPoint;

            //计算直线的K
            if (simularPoint.X < 1)
                return sectPoint;
            double a = simularPoint.Y / simularPoint.X;
            if (a < 0.0001)
                return sectPoint;

            //点越多越精确
            return GetSectPointLine(CurveExpress.GetFitPoints(100), simularPoint);
        }

        //曲线H=K*Q^i 与曲线Curve的交点(没有,返回Point(0,0))
        public static CurvePoint GetSectPointGeneral(List<CurvePoint> CurvePoints, CurvePoint simularPoint, double index)
        {
            CurvePoint sectPoint = new CurvePoint(0, 0);
            if (CurvePoints == null || CurvePoints.Count < 2)
                return sectPoint;

            if (simularPoint.X < 0.1)
                return sectPoint;

            if (Math.Abs(index - 1) < 0.01)
                return GetSectPointLine(CurvePoints, simularPoint);
            if (Math.Abs(index - 2) < 0.01)
                return GetSectPointParabola(CurvePoints, simularPoint);

            //计算系数K
            double fixK = simularPoint.Y / Math.Pow(simularPoint.X, index);
            if (fixK < 0.000001)
                return sectPoint;

            //思路是从simularPoint开始逐个增加0.1,直到k值最接近fixK
            double space = (CurvePoints.Last().X - simularPoint.X) / 1200;
            double x = simularPoint.X;
            double y, k;
            do
            {
                x = x + space;
                y = 0;
                if (!CurvePoints.GetInterPointY(x, out y))
                {
                    break;
                }

                k = y / Math.Pow(x, index);
            } while (k > fixK);

            sectPoint.X = x;
            sectPoint.Y = y;
            return sectPoint;
        }

        //曲线H=K*Q^i 与曲线Curve的交点(没有,返回Point(0,0))
        public static CurvePoint GetSectPointGeneral(this CurveExpress express, CurvePoint simularPoint, double index)
        {
            CurvePoint sectPoint = new CurvePoint(0, 0);
            if (express == null)
                return sectPoint;

            if (simularPoint.X < 1)
                return sectPoint;

            //检查是否在曲线的区域范围内
            double maxQ = express.Max;
            double maxH = express.GetFitPointY(maxQ);
            double k1 = maxH / Math.Pow(maxQ, index);
            double k2 = simularPoint.Y / Math.Pow(simularPoint.X, index);
            if (k1 > k2)
            {
                return sectPoint;
            }

            if (Math.Abs(index - 1) < 0.01)
                return GetSectPointLine(express, simularPoint);
            if (Math.Abs(index - 2) < 0.01)
                return GetSectPointParabola(express, simularPoint);

            //计算系数K
            double fixK = simularPoint.Y / Math.Pow(simularPoint.X, index);
            if (fixK < 0.000001)
                return sectPoint;

            //思路是从simularPoint开始逐个增加0.1,直到k值最接近fixK
            double space = (express.Max - simularPoint.X) / 1000;
            double x = simularPoint.X;
            double y, k;
            do
            {
                x = x + space;
                y = express.GetFitPointY(x);
                k = y / Math.Pow(x, index);
            } while (k > fixK);

            sectPoint.X = x;
            sectPoint.Y = y;
            return sectPoint;
        }

        //ratioIgnore:作用 当simularPoint超出范围时,扩大的倍数
        public static CurvePoint GetSectPointGeneral(this CurveExpress curveExpress, CurvePoint simularPoint, double index, double ratioIgnore)
        {
            CurvePoint sectPoint = new CurvePoint(0, 0);
            if (curveExpress == null)
                return sectPoint;

            if (simularPoint.X < 1)
                return sectPoint;
            var CurveExpressEx = curveExpress.Clone();
            //检查是否在曲线的区域范围内
            double maxQ = CurveExpressEx.Max;
            double maxH = CurveExpressEx.GetFitPointY(maxQ);
            double k1 = maxH / Math.Pow(maxQ, index);
            double k2 = simularPoint.Y / Math.Pow(simularPoint.X, index);
            if (k1 > k2)
            {
                CurveExpressEx.Max = CurveExpressEx.Max * ratioIgnore;//放大1.2倍
            }

            if (Math.Abs(index - 1) < 0.01)
                return GetSectPointLine(CurveExpressEx, simularPoint);
            if (Math.Abs(index - 2) < 0.01)
                return GetSectPointParabola(CurveExpressEx, simularPoint);

            //计算系数K
            double fixK = simularPoint.Y / Math.Pow(simularPoint.X, index);
            if (fixK < 0.000001)
                return sectPoint;

            //思路是从simularPoint开始逐个增加0.1,直到k值最接近fixK
            double space = (CurveExpressEx.Max - simularPoint.X) / 1000;
            double x = simularPoint.X;
            double y, k;
            do
            {
                x = x + space;
                y = CurveExpressEx.GetFitPointY(x);
                k = y / Math.Pow(x, index);
            } while (k > fixK);

            sectPoint.X = x;
            sectPoint.Y = y;
            return sectPoint;
        }

        #endregion 曲线与H=K*Q^i的交点 protect类型,给子类调用,怎么覆盖GetSectPoint
    }
}