using System.Collections.Generic;
using System.Linq;
using System.Xml.Linq;
using Deltares.DamEngine.Calculators.KernelWrappers.DamMacroStability;
using Deltares.DamEngine.Calculators.KernelWrappers.DamMacroStability.Assemblers;
using Deltares.DamEngine.Data.Design;
using Deltares.DamEngine.Data.General;
using Deltares.DamEngine.Data.Geometry;
using Deltares.DamEngine.Data.Geotechnics;
using Deltares.DamEngine.Data.Standard.Logging;
using NUnit.Framework;

namespace Deltares.DamEngine.Calculators.Tests.KernelWrappers.DamMacroStability
{
    [TestFixture]
    public class MStabXmlDocTests
    {
        [Test, Category("Work_In_Progress")]
        public void TestCreateMStabXmlDoc()
        {
            var soilGeometry2DName = "testgeom2D";
            var riverLevel = 1.1;
            var mstabDesignEmbankment = new MStabDesignEmbankment();
            var trafficLoad = 2.2;
            var requiredSafetyFactor = 1.2;
            List<LogMessage> errorMessages;

            var stabilityProjectFilename = "testproject";
            var soilDbName = "soilmaterials.mdb";
            var scenario = new DesignScenario();
            scenario.Location = new Location();
            scenario.Location.SoildatabaseName = soilDbName;
            // ToDo zant_mxl I think a Location needs Segments for location.GetMostProbableProfile while creating Pl lines
            //#Bka : Esther. to get this working you must add a lot more data otherwise this test will never work
            scenario.Location.Segment = new Segment();
            scenario.Location.Segment.SoilProfileProbabilities.Add(new SoilGeometryProbability()); //#Bka SoilprofProb should be filled!
            var soilProfile = CreatePipingSellmeijerProfileWithOneSandlayer();
            var surfaceLine = CreateSurfaceLineTutorial1();

            XDocument mstabXml = MStabXmlDoc.CreateMStabXmlDoc(stabilityProjectFilename, scenario, soilProfile, soilGeometry2DName,
                riverLevel, mstabDesignEmbankment, surfaceLine, trafficLoad, requiredSafetyFactor, out errorMessages);
            mstabXml.Save(stabilityProjectFilename + ".xml");

            XElement inputElement = (from element in mstabXml.Root.Descendants()
                                     where element.Name.LocalName == DamMStabAssembler.XmlElementNameInput
                                     select element).Single();

            XAttribute attribute = inputElement.Attribute(DamMStabAssembler.XmlAttributeMStabFileName);
            Assert.AreEqual(stabilityProjectFilename, attribute.Value);

            attribute = inputElement.Attribute(DamMStabAssembler.XmlAttributeSoilDBName);
            Assert.AreEqual(soilDbName, attribute.Value);

            XElement surfaceLineElement = (from element in inputElement.Descendants()
                                           where element.Name.LocalName == DamMStabAssembler.XmlElementSurfaceLine
                                           select element).Single();
            Assert.AreEqual(10, surfaceLineElement.Elements().Count());

            XElement charPoints = (from element in inputElement.Descendants()
                                   where element.Name.LocalName == DamMStabAssembler.XmlElementCharacteristicPoints
                                   select element).Single();
            Assert.AreEqual(10, charPoints.Elements().Count());
        }
        private static SoilProfile1D CreatePipingSellmeijerProfileWithOneSandlayer()
        {
            SoilProfile1D soilProfile1D = new SoilProfile1D();
            SoilLayer1D soilLayer1D1 = new SoilLayer1D();
            soilLayer1D1.Name = "L0";
            soilLayer1D1.TopLevel = 10.0;
            soilLayer1D1.Soil = new Soil("Topmaterial", 1.0, 1.0);
            soilLayer1D1.Soil.PermeabKx = 0.0003;
            soilLayer1D1.Soil.DiameterD70 = 0.0003;
            soilLayer1D1.Soil.WhitesConstant = 0.5;
            soilLayer1D1.Soil.BeddingAngle = 57.0;
            soilLayer1D1.IsAquifer = false;
            soilProfile1D.Layers.Add(soilLayer1D1);

            SoilLayer1D soilLayer1D2 = new SoilLayer1D();
            soilLayer1D2.Name = "L1";
            soilLayer1D2.TopLevel = -2.0;
            soilLayer1D2.Soil = new Soil("Sand", 22.0, 20.0);
            soilLayer1D2.Soil.PermeabKx = 0.0001;
            soilLayer1D2.Soil.DiameterD70 = 0.0002;
            soilLayer1D2.Soil.WhitesConstant = 0.25;
            soilLayer1D2.Soil.BeddingAngle = 37.0;
            soilLayer1D2.IsAquifer = true;
            soilProfile1D.Layers.Add(soilLayer1D2);
            soilProfile1D.BottomLevel = -10.0;
            return soilProfile1D;
        }
        private static SurfaceLine2 CreateSurfaceLineTutorial1(bool includingTraffic = false)
        {
            SurfaceLine2 surfaceLine2 = new SurfaceLine2();
            surfaceLine2.Name = "Tutorial1";
            surfaceLine2.Geometry = new GeometryPointString();
            surfaceLine2.CharacteristicPoints.GeometryMustContainPoint = true;
            SurfaceLine2 line = surfaceLine2;
            line.EnsurePointOfType(0.0, 0.0, CharacteristicPointType.SurfaceLevelOutside);
            line.EnsurePointOfType(10.0, 0.0, CharacteristicPointType.DikeToeAtRiver);
            line.EnsurePointOfType(34.5, 5.0, CharacteristicPointType.DikeTopAtRiver);
            if (includingTraffic)
            {
                line.EnsurePointOfType(35.0, 5.0, CharacteristicPointType.TrafficLoadOutside);
                line.EnsurePointOfType(38.5, 5.0, CharacteristicPointType.TrafficLoadInside);
            }
            line.EnsurePointOfType(40.5, 5.0, CharacteristicPointType.DikeTopAtPolder);
            line.EnsurePointOfType(50.5, 0.0, CharacteristicPointType.DikeToeAtPolder);
            line.EnsurePointOfType(58.5, 0.0, CharacteristicPointType.DitchDikeSide);
            line.EnsurePointOfType(59.5, -2.0, CharacteristicPointType.BottomDitchDikeSide);
            line.EnsurePointOfType(61.5, -2.0, CharacteristicPointType.BottomDitchPolderSide);
            line.EnsurePointOfType(61.5, 0.0, CharacteristicPointType.DitchPolderSide);
            line.EnsurePointOfType(75.0, 0.0, CharacteristicPointType.SurfaceLevelInside);
            line.Geometry.SyncCalcPoints();
            return line;
        }
    }
}