using ErrorAnalysis.Repository; using System; using System.Collections.Generic; using System.Linq; using System.Text; using System.Threading.Tasks; using static System.Runtime.InteropServices.JavaScript.JSType; namespace ErrorAnalysis.Service { public class COMergeCalcService { //Ncoef=1e8*2.54*2.54*2*6*11*2*77*0.2/0.6*speed; //Fcoef=1e8*2.54*2.54*4*6*2*11*77*0.2/0.6*speed; const double _gg = (9.0 - 0.0) / 255.0;//0.03529 const double _offset = 0.0; const int _cWinStartIndex = (int)(3.3 / _gg + _offset); //C窗口index const int _cWinEndIndex = (int)(4.78 / _gg + _offset); // width= 4.78-3.3 = 1.48 const int _oWinStartIndex = (int)(4.88 / _gg + _offset); // org const int _oWinEndIndex = (int)(6.36 / _gg + _offset); // width = 6.36-4.88=1.48 private static double GetFarInterplolateResult(COFarResultRepository repos, string modelId, double porosity, double sw, double speed, double depth) { } private static double GetNearInterplolateResult(COFarResultRepository repos, string modelId, double porosity, double sw, double speed, double depth) { } public static double GetFarMergeCResult(string connectionString, string modelId, double porosity, double sw, double speed, double depth) { if (porosity > 40) throw new InvalidDataException("Porosity value out of range!"); var repository = new COFarResultRepository(connectionString); double result = 0; if (porosity % 5 == 0) { var cResult = repository.GetCOFarResult(modelId, (int)porosity, 0); if (cResult == null) throw new InvalidDataException("COFarResult not found"); var cArr = cResult.InelasticSpec?.Split(',').Select(v => Convert.ToDouble(v)).ToArray(); if (cArr == null) throw new InvalidDataException("COFarResult InelasticSpec is null"); var length = _cWinEndIndex - _cWinStartIndex + 1; var originC = cArr.Skip(_cWinStartIndex).Take(length).Sum(); var coef = originC * 1e8 * 2.54 * 2.54 * 4 * 6 * 2 * 11 * 77 * 0.2 / 0.6 * speed * depth * 0.07; if (sw > 0) result = Utility.Interpolate(sw, 100, 0, 0, coef); } else { double ceilingResult = 0; var ceilingPorosity = Math.Ceiling(porosity / 5) * 5; var ceilingCResult = repository.GetCOFarResult(modelId, (int)ceilingPorosity, 0); if (ceilingCResult == null) throw new InvalidDataException("COFarResult not found"); var ceilingCArr = ceilingCResult.InelasticSpec?.Split(',').Select(v => Convert.ToDouble(v)).ToArray(); if (ceilingCArr == null) throw new InvalidDataException("COFarResult InelasticSpec is null"); var ceilingLength = _cWinEndIndex - _cWinStartIndex + 1; var ceilingOriginC = ceilingCArr.Skip(_cWinStartIndex).Take(ceilingLength).Sum(); var ceilingCoef = ceilingOriginC * 1e8 * 2.54 * 2.54 * 4 * 6 * 2 * 11 * 77 * 0.2 / 0.6 * speed; if (sw > 0) ceilingResult = Utility.Interpolate(sw, 100, 0, 0, ceilingCoef); double floorResult = 0; var floorPorosity = Math.Floor(porosity / 5) * 5; var floorCResult = repository.GetCOFarResult(modelId, (int)floorPorosity, 0); if (floorCResult == null) throw new InvalidDataException("COFarResult not found"); var floorCArr = floorCResult.InelasticSpec?.Split(',').Select(v => Convert.ToDouble(v)).ToArray(); if (floorCArr == null) throw new InvalidDataException("COFarResult InelasticSpec is null"); var floorLength = _cWinEndIndex - _cWinStartIndex + 1; var floorOriginC = floorCArr.Skip(_cWinStartIndex).Take(floorLength).Sum(); var floorCoef = floorOriginC * 1e8 * 2.54 * 2.54 * 4 * 6 * 2 * 11 * 77 * 0.2 / 0.6 * speed; if (sw > 0) floorResult = Utility.Interpolate(sw, 100, 0, 0, floorCoef); result = Utility.Interpolate(porosity, floorPorosity, floorResult, ceilingPorosity, ceilingResult); } return result; } public static double GetFarMergeOResult(string connectionString, string modelId, double porosity, double sw, double speed) { var repository = new COFarResultRepository(connectionString); var cResult = repository.GetCOFarResult(modelId, porosity, 100); if (cResult == null) throw new InvalidDataException("COFarResult not found"); var cArr = cResult.InelasticSpec?.Split(',').Select(v => Convert.ToDouble(v)).ToArray(); if (cArr == null) throw new InvalidDataException("COFarResult InelasticSpec is null"); var length = _oWinEndIndex - _oWinStartIndex + 1; var result = cArr.Skip(_oWinStartIndex).Take(length).Sum(); return result; } public static double GetNearMergeCResult(string connectionString, string modelId, double porosity, double sw, double speed) { var repository = new CONearResultRepository(connectionString); var cResult = repository.GetCONearResult(modelId, porosity, 0); if (cResult == null) throw new InvalidDataException("CONearResult not found"); var cArr = cResult.InelasticSpec?.Split(',').Select(v => Convert.ToDouble(v)).ToArray(); if (cArr == null) throw new InvalidDataException("CONearResult InelasticSpec is null"); var length = _cWinEndIndex - _cWinStartIndex + 1; var result = cArr.Skip(_cWinStartIndex).Take(length).Sum(); return result; } public static double GetNearMergeOResult(string connectionString, string modelId, double porosity, double sw, double speed) { var repository = new CONearResultRepository(connectionString); var cResult = repository.GetCONearResult(modelId, porosity, 100); if (cResult == null) throw new InvalidDataException("CONearResult not found"); var cArr = cResult.InelasticSpec?.Split(',').Select(v => Convert.ToDouble(v)).ToArray(); if (cArr == null) throw new InvalidDataException("CONearResult InelasticSpec is null"); var length = _oWinEndIndex - _oWinStartIndex + 1; var result = cArr.Skip(_oWinStartIndex).Take(length).Sum(); return result; } } }