using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;


namespace IStation.Model
{
    /// <summary>
    /// 样条曲线（三次拟合）：使用MathNet
    /// </summary>
    public class TwinRelateCurve3M : IFitCurvePoint
    {
        /// <summary>
        /// 
        /// </summary>
        /// <param name="express"></param>
        public TwinRelateCurve3M(CurveExpress express)
        {
            _express = express;
        }
        /// <summary>
        /// 
        /// </summary>
        /// <param name="points"></param>
        /// <exception cref="Exception"></exception>
        public TwinRelateCurve3M(List<CurvePoint> points)
        {
            if (points == null || points.Count < 4)
            {
                throw new Exception("出现数值计算问题 ，无法解出三次样条曲线！");
            }

            _express = BuildCurveExpress(points);
        }

        private CurveExpress? _express = null;
 


        /// <summary>
        /// 创建曲线表达式
        /// </summary>
        public static CurveExpress BuildCurveExpress(List<CurvePoint> points)
        {
            if (points == null || points.Count() < 3)
                return null;
            Model.CurveExpress express = new CurveExpress();
            express.DefinePoints =  points ;

            var X = (from x in points select x.X).ToArray();
            var Y = (from x in points select x.Y).ToArray();

            var min = X.Min();
            var max = X.Max();
            var res = IStation.Numerics.Fit.Polynomial(X, Y, 3);

            //曲线类型
            express.FitPow = 3;
            express.FitType = eCurveFitType.CubicCurve;
            express.DefinePoints = points.Select(x => new CurvePoint(x)).ToList();

            //范围
            express.Min = X.Min(); //最小值
            express.Max = X.Max();//最大值


            //0次方
            express.Index0 = res[0];
            //1次方
            express.Index1 = res[1];
            //2次方
            express.Index2 = res[2];
            //3次方
            express.Index3 = res[3];

            return express;
        }

        /// <summary>
        /// 获取拟合点Y
        /// </summary>
        public double GetFitPointY(double x)
        {
            if (_express == null)
                return default;
            var dx = x;// - _express.Min;  //不用减掉最小值, 和alglib有点不一样
            return _express.Index3 * dx * dx * dx + _express.Index2 * dx * dx + _express.Index1 * dx + _express.Index0;
        }

        /// <summary>
        /// 获取拟合点列表
        /// </summary>
        /// <param name="pointNumber"></param>
        /// <returns></returns>
        public List<IStation.Model.CurvePoint>? GetFitPoints(int pointNumber)
        {
            if (_express == null)
                return null;
            
            if (pointNumber < 1)
                pointNumber = 12;
            if (pointNumber > 10000)
                pointNumber = 10000;
            double space = (_express.Max - _express.Min) / (pointNumber - 1);
            if (space < 0.0001)
                return null;
            List< IStation.Model.CurvePoint> points = new List<IStation.Model.CurvePoint>();
            for (int i = 0; i < pointNumber; i++)
            {
                double x = space * i + _express.Min;
                double y = GetFitPointY(x);
                points.Add(new IStation.Model.CurvePoint(x, y));
            }
            return points;
        }

        /// <summary>
        /// 通过区间获取拟合点列表
        /// </summary>
        /// <param name="x_min"></param>
        /// <param name="x_max"></param>
        /// <param name="pointNumber"></param>
        /// <returns></returns>
        public List<CurvePoint>? GetFitPointsByRange(double x_min, double x_max, int pointNumber)
        {
            if (pointNumber < 1)
                pointNumber = 12;
            if (pointNumber > 10000)
                pointNumber = 10000;
            if(x_min >= x_max)
            {//交换大小
                var x_m = x_max;
                x_max = x_min;
                x_min = x_m;
            }
            double space = (x_max - x_min) / (pointNumber - 1);
            if (space < 0.0001)
                return default;
     
            if (_express == null)
                return null;

            if (pointNumber < 1)
                pointNumber = 12;
            if (pointNumber > 10000)
                pointNumber = 10000;
    
            if (space < 0.0001)
                return null;
            List<IStation.Model.CurvePoint> points = new List<IStation.Model.CurvePoint>();
            for (int i = 0; i < pointNumber; i++)
            {
                double x = space * i + x_min;
                double y = GetFitPointY(x);
                points.Add(new IStation.Model.CurvePoint(x, y));
            }
            return points;
        }

        /// <summary>
        /// 获取拟合点Y
        /// </summary>
        public static double GetFitPointY(CurveExpress express, double x)
        {
            if (express == null)
                return default;
            if (express.IsNull)
                return default;
            return new TwinRelateCurve3M(express).GetFitPointY(x);
        }

        /// <summary>
        /// 获取拟合点Y
        /// </summary>
        public static double GetFitPointY(List<CurvePoint> points, double x)
        {
            if (points == null || points.Count < 4)
                return default;
            return new TwinRelateCurve3M(points).GetFitPointY(x);
        }

        /// <summary>
        /// 获取拟合点
        /// </summary>
        public static List<Model.CurvePoint>? GetFitPoints(CurveExpress express, int pointNumber)
        {
            if (express == null)
                return default;
            return GetFitPointsByRange(express, express.Min, express.Max, pointNumber);
        }

        /// <summary>
        /// 通过区间获取拟合点列表
        /// </summary>
        public static List<Model.CurvePoint>? GetFitPointsByRange(CurveExpress express, double x_min, double x_max, int pointNumber)
        {
            if (express == null)
                return default;
            return new TwinRelateCurve3M(express).GetFitPointsByRange(x_min, x_max, pointNumber);
        }


 
 





    }
}