using Accord;
using Accord.Genetic; 
using Optimization1D;

using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;
namespace IStation
{
    public class EpanetHelper
    {  
        EPAFunctionGUI epaFunction;

        const float MAX_ENERGE = 99999999999;


        /// <summary>
        /// err问题 
        /// </summary>
        private string ReportQuestion(int err)
        {
            ErrorDefine errorDefine = new ErrorDefine();
            #region [秦白云飞修订]2021-4-30
            //新增目的:查看报错问题
            //新增

#if DEBUG
            throw (new Exception(errorDefine.getErrorTxt(err)));
#endif
            #endregion [秦白云飞修订]2021-4-30

            return errorDefine.getErrorTxt(err);
        }

        #region const_mode
        public const int EN_ELEVATION = 0;  //  { Mode parameters }
        public const int EN_BASEDEMAND = 1;
        public const int EN_PATTERN = 2;
        public const int EN_EMITTER = 3;
        public const int EN_INITQUAL = 4;
        public const int EN_SOURCEQUAL = 5;
        public const int EN_SOURCEPAT = 6;
        public const int EN_SOURCETYPE = 7;
        public const int EN_TANKLEVEL = 8;
        public const int EN_DEMAND = 9;
        public const int EN_HEAD = 10;
        public const int EN_PRESSURE = 11;
        public const int EN_QUALITY = 12;
        public const int EN_SOURCEMASS = 13;
        public const int EN_INITVOLUME = 14;
        public const int EN_MIXMODEL = 15;
        public const int EN_MIXZONEVOL = 16;
        #endregion

        #region type_model
        enum Type_EPAObject
        {
            Node = 0,
            Link = 1,
            Pump = 2,
            Tank = 3
        }
        #endregion

        #region other
        Dictionary<ModelInput.eType, Type_EPAObject> Typelist = new Dictionary<ModelInput.eType, Type_EPAObject>() { { ModelInput.eType.level, Type_EPAObject.Tank }, { ModelInput.eType.flow, Type_EPAObject.Node }, { ModelInput.eType.press_out, Type_EPAObject.Tank }, { ModelInput.eType.pump_run, Type_EPAObject.Pump }, { ModelInput.eType.press_in, Type_EPAObject.Tank }, { ModelInput.eType.valve_initstatus, Type_EPAObject.Link } };

        List<ModelOutput.eType> const_list_nodeType = new List<ModelOutput.eType> { ModelOutput.eType.EN_JUNCTION, ModelOutput.eType.EN_RESERVOIR, ModelOutput.eType.EN_TANK };
        List<ModelOutput.eType> const_list_linkType = new List<ModelOutput.eType> { ModelOutput.eType.EN_CVPIPE, ModelOutput.eType.EN_PIPE, ModelOutput.eType.EN_PUMP, ModelOutput.eType.EN_PRV, ModelOutput.eType.EN_PSV, ModelOutput.eType.EN_FCV, ModelOutput.eType.EN_TCV, ModelOutput.eType.EN_GPV };

        enum Method
        {
            Count = 0,
            Sum = 1,
            Avg = 2
        }

        string GetTypeValue(int typeCode, int ObjType, Method method, out float Value, int valueCode = -1)
        {
            Value = -1;
            List<float> listValue = new List<float>();
            float result = 0;
            int err = 0;
            if (ObjType == 0)
            {
                int NLink1 = 0;
                err = EpanetBase.ENgetcount(EpanetBase.Const_class.Const_Component.EN_NODECOUNT, ref NLink1);
                if (err > 9) return ReportQuestion(err);
                for (int i = 1; i < NLink1; i++)
                {
                    int type = -1;
                    err = EpanetBase.ENgetnodetype(i, ref type);
                    if (err > 9) return ReportQuestion(err);
                    if (type == typeCode)
                    {
                        if (method == Method.Count)
                            result++;
                        else
                        {
                            float value = -1;
                            err = EpanetBase.ENgetnodevalue(i, valueCode, ref value);
                            if (err > 9) return ReportQuestion(err);
                            listValue.Add(value);
                        }
                    }
                }
            }
            else
            {
                int NLink1 = 0; //int pump_count = 0;
                err = EpanetBase.ENgetcount(EpanetBase.Const_class.Const_Component.EN_LINKCOUNT, ref NLink1);
                if (err > 9) return ReportQuestion(err);
                for (int i = 1; i < NLink1; i++)
                {
                    int type = -1;
                    err = EpanetBase.ENgetlinktype(i, ref type);
                    if (err > 9) return ReportQuestion(err);

                    if (type == typeCode)
                    {
                        if (method == Method.Count)
                            result++;
                        else
                        {
                            float value = -1;
                            err = EpanetBase.ENgetlinkvalue(i, valueCode, ref value);
                            if (err > 9) return ReportQuestion(err);
                            listValue.Add(value);
                        }
                    }
                }
            }
            if (method == Method.Sum)
            {
                result = 0;
                for (int i = 0; i < listValue.Count; i++)
                    result += listValue[i];
            }
            else if (method == Method.Avg)
            {
                result = 0;
                for (int i = 0; i < listValue.Count; i++)
                    result += listValue[i];
                result /= listValue.Count;
            }
            Value = result;

            return null;
        }
        string SetPumpValue(List<bool> isfreq, float value, List<string> Combine)
        {
            int err = 0;
            if (value != -1)
                for (int i = 0; i < Combine.Count; i++)
                {

                    if (!isfreq[i]) //工频泵
                    {
                        if (value > 0)
                        {
                            err = SetItemValue(Combine[i], 1, Type_EPAObject.Pump);
                            if (err > 9) return ReportQuestion(err);
                            return null;


                        }
                        else
                        {
                            err = SetItemValue(Combine[i], 0, Type_EPAObject.Pump);
                            if (err > 9) return ReportQuestion(err);
                            return null;
                        }
                    }
                    err = SetItemValue(Combine[i], value / 50.0, Type_EPAObject.Pump);
                    if (err > 9) return ReportQuestion(err);
                }

            return null;
        }
        string SetPumpValue(List<bool> isfreq, double[] value, List<string> Combine)
        {
            int err = 0;
            for (int i = 0; i < Combine.Count; i++)
            {

                if (!isfreq[i]) //工频泵
                {
                    if (value[i] > 0)
                    {
                        err = SetItemValue(Combine[i], 1, Type_EPAObject.Pump);
                        if (err > 9) return ReportQuestion(err);
                        return null;


                    }
                    else
                    {
                        err = SetItemValue(Combine[i], 0, Type_EPAObject.Pump);
                        if (err > 9) return ReportQuestion(err);
                        return null;
                    }
                }
                err = SetItemValue(Combine[i], value[i] / 50.0, Type_EPAObject.Pump);
                if (err > 9) return ReportQuestion(err);
            }

            return null;
        }
        int SetTypeValue(List<ModelInput> InputData, ModelInput.eType type, float value = -1, int Num = 0)
        {
            int err = 0;
            var InputData_type = (InputData as List<ModelInput>).FindAll(m => (m.Type == type));
            if (value == -1)
                for (int i = 0; i < InputData_type.Count; i++)
                {
                    err = SetItemValue(InputData_type[i].ObjectID, InputData_type[i].Value, Typelist[type]);
                    if (err > 9) return err;
                }
            if (value != -1)
                for (int i = 1; i <= Num; i++)
                {
                    err = SetItemValue(type + "_" + i.ToString(), value, Typelist[type]);
                    if (err > 9) return err;
                }
            return err;
        }
        int SetItemValue(string ID, dynamic value, Type_EPAObject type_Obj)
        {
            int Index = -1;
            int err = 0;
            switch (type_Obj)
            {
                case (Type_EPAObject.Tank):

                    err = EpanetBase.ENgetnodeindex(ID, ref Index);
                    err = EpanetBase.ENsetnodevalue(Index, EN_TANKLEVEL, value);  //源头压力基数
                    break;
                case (Type_EPAObject.Node):
                    err = EpanetBase.ENgetnodeindex(ID, ref Index);
                    err = EpanetBase.ENsetnodevalue(Index, EN_BASEDEMAND, value);  //源头压力基数
                    break;
                case (Type_EPAObject.Pump):
                    float speed = (float.Parse(value.ToString()));
                    if (speed < 0.1) speed = 0;
                    if (speed > 0)
                        err = EpanetBase.ENsetlinkvalue(Index, 4, 1);
                    else
                        err = EpanetBase.ENsetlinkvalue(Index, 4, 0);
                    //float speed = float.Parse(value.ToString());
                    err = EpanetBase.ENgetlinkindex(ID, ref Index);
                    //err = EpanetBase.ENsetlinkvalue(Index, 5, speed);
                    err = EpanetBase.ENsetinistatus(ID, speed.ToString());
                    //err = EpanetBase.ENsaveinpfile("C:\\EPA\\修改后" + DateTime.Now.ToString("yyyyMMddHHmmss") + ".inp");
                    if (err > 9) return err;


                    break;
                case (Type_EPAObject.Link):
                    err = EpanetBase.ENgetlinkindex(ID, ref Index);
                    err = EpanetBase.ENsetlinkvalue(Index, EpanetBase.Const_class.Const_Link.EN_INITSTATUS, value);
                    break;
                default:
                    err = -1;
                    break;
            }

            return err;




        }

        /// <summary>
        /// 
        /// </summary>
        /// <returns></returns>
        public int GetOutputString(out List<ModelOutputNode> ListOutputNode, out List<ModelOutputLink> ListOutputLink)
        {
            int err = 0;
            //string result = null;
            ListOutputNode = new List<ModelOutputNode>();
            ListOutputLink = new List<ModelOutputLink>();
            int Njunction = 0;
            err = EpanetBase.ENgetcount(EpanetBase.Const_class.Const_Component.EN_NODECOUNT, ref Njunction);
            if (err > 9) return err;

            int NLink = 0;
            err = EpanetBase.ENgetcount(EpanetBase.Const_class.Const_Component.EN_LINKCOUNT, ref NLink);
            if (err > 9) return err;
            List<float> output = new List<float>() { Njunction, NLink };

            for (int i = 1; i <= Njunction; i++)
            {
                ModelOutputNode mOutput = new ModelOutputNode();
                StringBuilder ID = new StringBuilder();

                float Value = 0; int type = -1;
                err = EpanetBase.ENgetnodeid(i, ID);
                if (err > 9) return err;

                mOutput.ObjectID = ID.ToString();
                err = EpanetBase.ENgetnodetype(i, ref type);
                if (err > 9) return err;
                mOutput.Type = (ModelOutputNode.eType)type;
                var list_ValueIndex = mOutput.Get_list_ValueIndex();
                for (int j = 0; j < list_ValueIndex.Length; j++)
                {

                    err = EpanetBase.ENgetnodevalue(i, list_ValueIndex[j], ref Value);  //源头压力基数
                    if (err > 9) return err;
                    //if (Value <= 0.0001) Value = 0;
                    mOutput.SetValue(j, Value);


                }
                ListOutputNode.Add(mOutput);
            }
            for (int i = 1; i <= NLink; i++)
            {
                ModelOutputLink mOutput = new ModelOutputLink();
                StringBuilder idBuilder = new StringBuilder();
                //mOutput.Value = new List<float>();
                float Value = 0; int type = -1;
                err = EpanetBase.ENgetlinkid(i, idBuilder);  //源头压力基数
                if (err > 9) return err;
                mOutput.ObjectID = idBuilder.ToString();
                err = EpanetBase.ENgetlinktype(i, ref type);
                if (err > 9) return err;
                mOutput.Type = (ModelOutputLink.eType)(type + 10);
                var list_ValueIndex = mOutput.Get_list_ValueIndex();
                for (int j = 0; j < list_ValueIndex.Length; j++)
                {
                    err = EpanetBase.ENgetlinkvalue(i, list_ValueIndex[j], ref Value);  //源头压力基数
                    if (err > 9) return err;
                    //if (Value <= 0.0001) Value = 0;
                    mOutput.SetValue(j, Value);
                }
                ListOutputLink.Add(mOutput);
            }
            return err;
        }
        private string SetPumpStatus()
        {
            int NLink = 0; int err = 0;
            err = EpanetBase.ENgetcount(EpanetBase.Const_class.Const_Component.EN_LINKCOUNT, ref NLink);
            if (err > 9) return ReportQuestion(err);


            for (int i = 1; i <= NLink; i++)
            {
                //ModelOutputLink mOutput = new ModelOutputLink();
                //StringBuilder idBuilder = new StringBuilder();
                ////mOutput.Value = new List<float>();
                //float Value = 0; 
                //err = EpanetBase.ENgetlinkid(i, idBuilder);  //源头压力基数
                //if (err > 9) return ReportQuestion(err);
                //mOutput.ObjectID = idBuilder.ToString();
                int type = -1;
                err = EpanetBase.ENgetlinktype(i, ref type);
                if (err > 9) return ReportQuestion(err);
                if (type == EpanetBase.Const_class.Const_Link_types.EN_PUMP)
                {

                    err = EpanetBase.ENsetlinkvalue(i, EpanetBase.Const_class.Const_Link.EN_INITSTATUS, 0);

                }

            }
            return null;
        }
        #endregion


        private bool _isInitialSuccess = false;

        /// <summary>
        /// 初始化
        /// </summary>
        /// <param name="file"></param>
        /// <returns></returns>
        public string Initial(string file)
        {

            if (!System.IO.File.Exists(file))
                return "inp 不存在";
            var err = EpanetBase.ENopen(file, "0.rpt", "0.bin");
            if (err == 0)
            {
                _isInitialSuccess = true;
                return null;
            }
            return ReportQuestion(err);
        }

        public string GetResultByFlowAndPumpStatus(List<ModelInput> InputData, out List<ModelOutputNode> modelOutputNodes, out List<ModelOutputLink> modelOutputLinks)
        {
            #region 初始判断

            modelOutputLinks = null;
            modelOutputNodes = null;
            int err = 0;
            //string result = null;
            if (!_isInitialSuccess)
                return "初始化不成功";
            if (InputData == null)
                return "输入参数为空";
            #endregion
            #region 传参计算
            err = SetTypeValue(InputData, ModelInput.eType.press_in);
            if (err > 9) return ReportQuestion(err);
            err = SetTypeValue(InputData, ModelInput.eType.flow);
            if (err > 9) return ReportQuestion(err);
            err = SetTypeValue(InputData, ModelInput.eType.pump_run);
            if (err > 9) return ReportQuestion(err);
            err = SetTypeValue(InputData, ModelInput.eType.level);
            if (err > 9) return ReportQuestion(err);
            err = SetTypeValue(InputData, ModelInput.eType.valve_initstatus);
            if (err > 9) return ReportQuestion(err);

            err = EpanetBase.ENsolveH();
            if (err > 9) return ReportQuestion(err);
            #endregion
          /*  var path = System.IO.Directory.GetCurrentDirectory() + @"\model\out.inp";
            err = EpanetBase.ENsaveinpfile(path);
            if (err > 9) return ReportQuestion(err);*/
            #region 返回结果
            err = GetOutputString(out modelOutputNodes, out modelOutputLinks);
            if (err > 9) return ReportQuestion(err);
            return null;
            #endregion 

        }

        public string GetResultByPressAndPumpStatus(List<ModelInput> InputData, out List<ModelOutputNode> modelOutputNodes, out List<ModelOutputLink> modelOutputLinks)
        {


            #region 初始判断

            modelOutputLinks = null;
            modelOutputNodes = null;
            //string result = null;
            int err = 0;
            if (!_isInitialSuccess)
                return "初始化不成功";
            if (InputData == null)
                return "输入参数为空";
            #endregion

            #region 传参计算
            err = SetTypeValue(InputData, ModelInput.eType.press_in);
            if (err > 9) return ReportQuestion(err);
            err = SetTypeValue(InputData, ModelInput.eType.press_out);
            if (err > 9) return ReportQuestion(err);
            err = SetTypeValue(InputData, ModelInput.eType.pump_run);
            if (err > 9) return ReportQuestion(err);
            err = SetTypeValue(InputData, ModelInput.eType.level);
            if (err > 9) return ReportQuestion(err);
            err = SetTypeValue(InputData, ModelInput.eType.valve_initstatus);
            if (err > 9) return ReportQuestion(err);

            err = EpanetBase.ENsolveH();
            if (err > 9) return ReportQuestion(err);
            #endregion
           /* var path = System.IO.Directory.GetCurrentDirectory() + @"\model\out.inp";
            err = EpanetBase.ENsaveinpfile(path);
            if (err > 9) return ReportQuestion(err);*/
            #region 返回结果
            err = GetOutputString(out modelOutputNodes, out modelOutputLinks);
            if (err > 9) return ReportQuestion(err);
            return null;
            #endregion
        }

        /// <summary>
        /// 获取水泵转速比组合
        /// </summary>
        /// <param name="InputData"></param>
        /// <param name="List_Combine"></param>
        /// <param name="List_modelOutputNodes"></param>
        /// <param name="List_modelOutputLinks"></param>
        /// <param name="opt_combines"></param>
        /// <param name="opt_energes"></param>
        /// <param name="populationSize"></param>
        /// <param name="iterations"></param>
        /// <param name="range"></param>
        /// <param name="accuracy"></param>
        /// <returns></returns>
        public string GetResultByFlowAndPress(List<ModelInput> InputData, List<List<string>> List_Combine,
            out List<List<ModelOutputNode>> List_modelOutputNodes, out List<List<ModelOutputLink>> List_modelOutputLinks,
            out List<List<double>> opt_combines, out List<double> opt_energes, double rangeMin = 0.8, double rangeMax = 1,
            int populationSize = 40, int iterations = 20, double accuracy = 0.1)
        {

            #region 初始判断
            List_modelOutputNodes = new List<List<ModelOutputNode>>();
            List_modelOutputLinks = new List<List<ModelOutputLink>>();
            opt_combines = new List<List<double>>();
            opt_energes = new List<double>();
            string result = null;
            var range = new DoubleRange(rangeMin, rangeMax);

            if (!_isInitialSuccess)
                return "初始化不成功";
            if (InputData == null)
                return "输入参数为空";
            #endregion

            //float min_energe = MAX_ENERGE;
            for (int i = 0; i < List_Combine.Count; i++)
            {
                List<ModelOutputNode> modelOutputNodes;
                List<ModelOutputLink> modelOutputLinks;
                List<double> opt_combine = new List<double>();
                double energe;
                calcParms parms;
                var dr = new DoubleRange(range.Min * 50, range.Max * 50);
                result = InputInit(InputData, List_Combine[i], out parms, dr.Min, dr.Max);
                if (result != null)
                    return result;
                string result0 = _GetResultByFlowAndPress(InputData, List_Combine[i], out modelOutputNodes, out modelOutputLinks, out opt_combine, out energe, dr.Min, dr.Max, populationSize, iterations, accuracy);
                List_modelOutputNodes.Add(modelOutputNodes);
                List_modelOutputLinks.Add(modelOutputLinks);
                for (var j = 0; j < opt_combine.Count; j++) opt_combine[j] /= 50.0;
                opt_combines.Add(opt_combine);
                opt_energes.Add(energe);
            }
            return result;
        }

        #region other
        public void GetResult(double[] list_bestresult)
        {
            List<ModelOutputNode> nodes = null;
            List<ModelOutputLink> links = null;
            List<double> pumpflows = null;
            //double out_energe = 999999999f;
            double[] li = new double[list_bestresult.Length];
            for (int i = 0; i < li.Length; i++) li[i] = list_bestresult[i];
            this.epaFunction.GetResult(li, out nodes, out links, out pumpflows, true);
        }
        private string _GetResultByFlowAndPress(List<ModelInput> InputData, List<string> combine, out List<ModelOutputNode> modelOutputNodes, out List<ModelOutputLink> modelOutputLinks, out List<double> List_Hz, out double energe, double rangeMin = 0.8, double rangeMax = 1, int populationSize = 40, int iterations = 20, double accuracy = 0.1)
        {


            #region 初始判断
            modelOutputNodes = null;
            modelOutputLinks = null;
            //string result = null;
            //int err = 0;
            List_Hz = new List<double>();
            energe = 9999999999;
            #endregion
            #region 初始设置
            SetPumpStatus();
            #endregion
            var range = new DoubleRange(rangeMin, rangeMax);
            SearchSolution(combine, out List_Hz, out modelOutputNodes, out modelOutputLinks, out energe, rangeMin, rangeMax, populationSize, iterations, accuracy);
            return "";
        }

        private void SearchSolution(List<string> combine, out List<double> list_BestHzs, out List<ModelOutputNode> modelOutputNodes, out List<ModelOutputLink> modelOutputLinks, out double energe, double rangeMin = 0.8, double rangeMax = 1, int populationSize = 40, int iterations = 20, double accuracy = 0.1)
        {
            var list_BestEnerge = new List<double>();
            list_BestHzs = new List<double>();
            var list_AvgHz = new List<double>();

            modelOutputNodes = null;
            modelOutputLinks = null;
            energe = 999999999f;


            var range = new DoubleRange(rangeMin, rangeMax);
            var num_Total = range.Length / accuracy;


            var LENGTH_PUMP = 0;
            var value_LENGTH_PUMP = 1;
            var i = 0;
            while (value_LENGTH_PUMP < num_Total && i < 1000)
            {
                value_LENGTH_PUMP *= 2;
                LENGTH_PUMP += 1;
                i++;
            }

            var chromosomeLength = combine.Count * LENGTH_PUMP;
            BinaryChromosome ancestor = new BinaryChromosome(chromosomeLength);



            EPAFunctionGUI fitnessFunction = this.epaFunction;
            var IntselectionMethod = 0;
            object obj;
            if (IntselectionMethod != 0)
            {
                if (IntselectionMethod != 1)
                {
                    ISelectionMethod selectionMethod = new RouletteWheelSelection();
                    obj = selectionMethod;
                }
                else
                {
                    obj = new RankSelection();
                }
            }
            else
            {
                obj = new EliteSelection();
            }

            var optimizationMode = 1;
            epaFunction.Mode = ((optimizationMode != 0) ? OptimizationFunctionXD.Modes.Minimization : OptimizationFunctionXD.Modes.Maximization);
            Population population = new Population(populationSize, ancestor, fitnessFunction, (ISelectionMethod)obj);

            int num = 1;
            //double[,] array = new double[showOnlyBest ? 1 : populationSize, 2];
            List<ResultCombine> list_result = new List<ResultCombine>();
            //float[,] array = new double[showOnlyBest ? 1 : populationSize,combine.Count+1];
            bool needToStop = false;
            bool showOnlyBest = true;
            while (!needToStop)
            {
                population.RunEpoch();
                list_result.Clear();
                if (showOnlyBest)
                {
                    list_result.Add(new ResultCombine());
                    list_result[0].arrayHz = epaFunction.Translate(population.BestChromosome);
                    list_result[0].totalEnerge = epaFunction.OptimizationFunction(list_result[0].arrayHz);
                    Console.WriteLine($"迭代次数{num}，总能耗{list_result[0].totalEnerge}");

                }
                else
                {
                    for (int j = 0; j < populationSize; j++)
                    {
                        list_result.Add(new ResultCombine());
                        list_result[j].arrayHz = epaFunction.Translate(population[j]);
                        list_result[j].totalEnerge = epaFunction.OptimizationFunction(list_result[j].arrayHz);
                    }
                }
                //chart.UpdateDataSeries("avg", array);
                //SetText(currentIterationBox, num.ToString());

                double[] list_bestresult = (epaFunction.Translate(population.BestChromosome));
                double bestEnerge = epaFunction.OptimizationFunction(list_bestresult);




                //SetText(currentValueBox, bestEnerge.ToString());
                List<double> pumpflows = null;

                this.epaFunction.GetResult(list_bestresult, out modelOutputNodes, out modelOutputLinks, out pumpflows, false);


                double sum = 0;

                energe = bestEnerge;

                sum /= list_bestresult.Length;
                list_BestEnerge.Add(bestEnerge);
                list_AvgHz.Add(sum);
                list_BestHzs = list_bestresult.ToList();



                num++;
                if (iterations != 0 && num > iterations)
                {

                    break;
                }

            }
        }

        public string InputInit(List<ModelInput> InputData, List<string> combine, out calcParms parms, double rangeMin = 0.8, double rangeMax = 1)
        {
            string result = null;
            #region 初始判断
            parms = null;
            //string result = null;
            int err = 0;

            var doubleRange = new DoubleRange(rangeMin, rangeMax);

            #endregion



            err = SetTypeValue(InputData, ModelInput.eType.press_in);
            if (err > 9) return ReportQuestion(err);
            err = SetTypeValue(InputData, ModelInput.eType.flow);
            if (err > 9) return ReportQuestion(err);
            err = SetTypeValue(InputData, ModelInput.eType.level);
            if (err > 9) return ReportQuestion(err);
            err = SetTypeValue(InputData, ModelInput.eType.valve_initstatus);
            if (err > 9) return ReportQuestion(err);

            List<bool> list_isfrequence = new List<bool>();
            combine.ForEach(ID => { bool isFre = InputData.Find(data => data.ObjectID == ID).Value == 1; list_isfrequence.Add(isFre); });
            var pressOutPoit = InputData.FindAll((m) => m.Type == ModelInput.eType.press_out).Select(p => (double)p.Value).ToList();
            var flowPointID = InputData.FindAll((m) => m.Type == ModelInput.eType.flow).Select(f => f.ObjectID).ToList();
            List<int> flowPointIndex = new List<int>();


            flowPointID.ForEach(FID =>
            {
                if (result != null) return;
                int index = -1;
                err = EpanetBase.ENgetnodeindex(FID, ref index);
                flowPointIndex.Add(index);
                if (err > 9) result = ReportQuestion(err);
            });
            if (result != null) return result;

            parms = new calcParms
            {
                list_isfrequence = list_isfrequence,
                flowPointIndex = flowPointIndex,
                demandPress = pressOutPoit

            };
            #region 计算单泵的长度

            var range = doubleRange;

            double accuracy = 0.1;
            var num_Total = range.Length / accuracy;


            var LENGTH_PUMP = 0;
            var value_LENGTH_PUMP = 1;
            var i = 0;
            while (value_LENGTH_PUMP < num_Total && i < 1000)
            {
                value_LENGTH_PUMP *= 2;
                LENGTH_PUMP += 1;
                i++;
            }

            #endregion
            List<DoubleRange> drs = new List<DoubleRange>();
            combine.ForEach(m => { drs.Add(range); });

            epaFunction = new EPAFunctionGUI(drs.ToArray(), LENGTH_PUMP, parms, this, combine);
            return null;
        }

        float accur = 0.1f;
        // <summary>
        /// 二分查找
        /// </summary>
        /// <param name="arr"></param>
        /// <param name="low">开始索引 0</param>
        /// <param name="high">结束索引 </param>
        /// <param name="key">要查找的对象</param>
        /// <returns></returns>
        private float BinarySearch(List<string> combine, List<bool> list_isfrequence, int flowPointIndex, float lowHz, float highHz, float pressure)
        {

            string result = null;
            float mid = (lowHz + highHz) / 2;
            if (lowHz > highHz)
                return -1;
            else
            {
                if (highHz - lowHz <= accur)
                {

                    return mid;
                }

            }
            float press_out;
            float energe;
            result = getPressByPumps(combine, list_isfrequence, mid, flowPointIndex, out press_out, out energe);
            if (result != null) return -1;


            if (press_out > pressure)
                return BinarySearch(combine, list_isfrequence, flowPointIndex, lowHz, mid, pressure);
            else
                return BinarySearch(combine, list_isfrequence, flowPointIndex, mid, highHz, pressure);

        }

        private string getPressByPumps(List<string> combine, List<bool> isfrequence, float Hz, int flowPointIndex, out float press, out float energe)
        {
            int err = 0;

            #region 初始判断

            press = 0;
            energe = 99999999f;
            string result = null;
            #endregion



            var PressNodeID = -1;
            result = SetPumpValue(isfrequence, 0, combine);
            if (result != null) return result;

            result = SetPumpValue(isfrequence, ((float)Hz), combine);
            if (result != null) return result;



            //err = EpanetBase.ENsaveinpfile(@"D:\Desktop\C程序设计\项目源码\水力建模\eptools\out.inp");

            float PressValue = -1;
            int t = 0;
            //err = EpanetBase.ENinitH(0);
            //err = EpanetBase.ENrunH(ref t);
            err = EpanetBase.ENsolveH();
            if (err > 9) return ReportQuestion(err);
            err = EpanetBase.ENgetnodevalue(flowPointIndex, EpanetBase.Const_class.Const_Node.EN_HEAD, ref PressValue);
            if (err > 9) return ReportQuestion(err);
            press = PressValue;

            float TotalEnerge = 9999999999f;

            result = GetTypeValue(EpanetBase.Const_class.Const_Link_types.EN_PUMP, 1, Method.Sum, out TotalEnerge, EpanetBase.Const_class.Const_Link.EN_ENERGY);
            if (result != null) return result;
            energe = TotalEnerge;



            return result;
        }

        public string getPressByPumps(List<string> combine, List<bool> isfrequence, double[] list_Hz, List<int> flowPointIndex, out List<double> press, out double energe, out List<ModelOutputNode> modelOutputNodes, out List<ModelOutputLink> modelOutputLinks, out List<double> pumpflows, bool isNeedOutput = false, string path = null)
        {
            int err = 0;

            #region 初始判断
            modelOutputNodes = null;
            modelOutputLinks = null;
            pumpflows = null;
            press = new List<double>();

            energe = 99999999f;
            string result = null;
            #endregion



            var PressNodeID = -1;
            result = SetPumpValue(isfrequence, 0, combine);
            if (result != null) return result;

            result = SetPumpValue(isfrequence, list_Hz, combine);
            if (result != null) return result;



            //err = EpanetBase.ENsaveinpfile(@"D:\Desktop\C程序设计\项目源码\水力建模\eptools\out.inp");

            List<float> PressValue = new List<float>();
            int t = 0;
            //err = EpanetBase.ENinitH(0);
            //err = EpanetBase.ENrunH(ref t);
            //path = @"D:\Desktop\C程序设计\项目源码\水力建模\eptools_v2\Debug\model\out.inp";
            //EpanetBase.ENsaveinpfile(path);
            err = EpanetBase.ENsolveH();
            if (err > 9) return ReportQuestion(err);
            if (result != null) return result;


            var press0 = new List<double>();
            flowPointIndex.ForEach(FIndex =>
            {
                if (result != null) return;
                float PressValue0 = -1;
                err = EpanetBase.ENgetnodevalue(FIndex, EpanetBase.Const_class.Const_Node.EN_HEAD, ref PressValue0);
                press0.Add(PressValue0);
                if (err > 9) result = ReportQuestion(err);
            });
            press = press0;



            float TotalEnerge = 9999999999f;

            result = GetTypeValue(EpanetBase.Const_class.Const_Link_types.EN_PUMP, 1, Method.Sum, out TotalEnerge, EpanetBase.Const_class.Const_Link.EN_ENERGY);
            if (result != null) return result;

            bool isCannotDeliver = false;

            combine.ForEach(pump_id =>
            {
                if (isCannotDeliver) return;
                int pumpindex = -1;
                if ((err = EpanetBase.ENgetlinkindex(pump_id, ref pumpindex)) > 9) { result = ReportQuestion(err); return; }
                float pumpstatus = -1;
                if ((err = EpanetBase.ENgetlinkvalue(pumpindex, EpanetBase.Const_class.Const_Link.EN_STATUS, ref pumpstatus)) > 9) { result = ReportQuestion(err); return; }
                float pumpFlow = -1;
                err = EpanetBase.ENgetlinkvalue(pumpindex, EpanetBase.Const_class.Const_Link.EN_FLOW, ref pumpFlow);
                if (pumpstatus == 0 || pumpFlow <= 50) isCannotDeliver = true;
            });
            if (isCannotDeliver) return null;
            energe = TotalEnerge;


            if (isNeedOutput)
            {
                var pf = new List<double>();
                combine.ForEach(pump_id =>
                {

                    int pumpindex = -1;
                    err = EpanetBase.ENgetlinkindex(pump_id, ref pumpindex); if (err > 9) ;
                    float pumpFlow = -1;
                    err = EpanetBase.ENgetlinkvalue(pumpindex, EpanetBase.Const_class.Const_Link.EN_FLOW, ref pumpFlow);
                    pf.Add(pumpFlow);
                });
                pumpflows = pf;
                EpanetBase.ENsaveinpfile(path);
                GetOutputString(out modelOutputNodes, out modelOutputLinks);
            }
            return result;
        }

        public string GetPumpList(out List<Pump> pumps)
        {
            #region 初始判断

            //modelOutputLinks = null;
            //modelOutputNodes = null;
            string result = null;
            int err = 0;
            pumps = new List<Pump>();
            if (!_isInitialSuccess)
                return "初始化不成功";
            //if (InputData == null)
            //    return "输入参数为空";
            #endregion



            #region 获取

            int NLink1 = 0;
            err = EpanetBase.ENgetcount(EpanetBase.Const_class.Const_Component.EN_NODECOUNT, ref NLink1);
            if (err > 9) return ReportQuestion(err);
            for (int i = 1; i <= NLink1; i++)
            {
                int type = -1;
                err = EpanetBase.ENgetnodetype(i, ref type);
                if (err > 9) return ReportQuestion(err);
                if (type == EpanetBase.Const_class.Const_Link_types.EN_PUMP)
                {
                    Pump pump = new Pump();

                    err = EpanetBase.ENgetlinkvalue(i, EpanetBase.Const_class.Const_Link.EN_STATUS, ref pump.value);
                    if (err > 9) return ReportQuestion(err);

                    StringBuilder s = new StringBuilder();
                    err = EpanetBase.ENgetlinkid(i, s);
                    pump.ID = s.ToString();
                    if (err > 9) return ReportQuestion(err);

                    pump.index = i;

                    pumps.Add(pump);

                }
            }




            return null;





            #endregion
        }
    }
    public class Pump
    {
        public string ID = null;
        public int index;
        public float value;

    }
    public class calcParms
    {
        public List<bool> list_isfrequence = null;
        public List<int> flowPointIndex = null;
        public List<double> demandPress = null;
    }
    public class ResultCombine
    {
        public double[] arrayHz ;
        public double totalEnerge;
    }
    #endregion

}