// Copyright (C) Stichting Deltares 2024. All rights reserved.
// 
// This file is part of the application DAM - UI.
// 
// DAM - UI is free software: you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published by
// the Free Software Foundation, either version 3 of the License, or
// (at your option) any later version.
// 
// This program is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
// GNU General Public License for more details.
// 
// You should have received a copy of the GNU General Public License
// along with this program. If not, see <http://www.gnu.org/licenses/>.
// 
// All names, logos, and references to "Deltares" are registered trademarks of
// Stichting Deltares and remain full property of Stichting Deltares at all times.
// All rights reserved.

using System.Collections.Generic;
using System.IO;
using Deltares.Dam.Data;
using Deltares.Dam.TestHelper;
using Deltares.Standard;
using Deltares.Standard.Logging;
using NUnit.Framework;

namespace Deltares.Dam.IntegrationTests
{
    [TestFixture]
    public class FailureMechanismTests
    {
        private const double tolerance3Decimals = 0.00051;
        private const double tolerance2Decimals = 0.0051;
        private const double tolerance1Decimals = 0.051;
        
        [Test]
        public void TestStabilityOutsideWith2DStixFiles()
        {
            string projectFilename = Path.Combine(Directory.GetCurrentDirectory(), @"TestData\StabilityOutside2DFromStixTest\BishopOutside.damx");
            List<CsvExportData> allCalculationResults = ComputeStabilityOutsideProject(projectFilename, 1);
            
            Assert.That(allCalculationResults, Has.Count.EqualTo(1));
            Assert.Multiple(() =>
            {
                Assert.That(allCalculationResults[0].LocationName, Is.EqualTo("DWP_1"));
                Assert.That(allCalculationResults[0].CalculationResult, Is.EqualTo(CalculationResult.Succeeded));
                Assert.That(allCalculationResults[0].X, Is.EqualTo(99718.000).Within(tolerance3Decimals));
                Assert.That(allCalculationResults[0].Y, Is.EqualTo(437106.000).Within(tolerance3Decimals));
                Assert.That(allCalculationResults[0].AnalysisType, Is.EqualTo(AnalysisType.NoAdaption));
                Assert.That(allCalculationResults[0].IsUplift, Is.EqualTo(false));
                Assert.That(allCalculationResults[0].ProfileName, Is.EqualTo("DWP_1.stix"));
                Assert.That(allCalculationResults[0].ProfileProbability, Is.EqualTo(100).Within(tolerance1Decimals));
                Assert.That(allCalculationResults[0].StabilityModel, Is.EqualTo(StabilityModelType.Bishop));
                Assert.That(allCalculationResults[0].SafetyFactor, Is.EqualTo(1.565).Within(tolerance3Decimals));
                Assert.That(allCalculationResults[0].ShoulderHeight, Is.Null);
                Assert.That(allCalculationResults[0].ToeAtPolderX, Is.Null);
                Assert.That(allCalculationResults[0].ToeAtPolderZ, Is.Null);
                Assert.That(allCalculationResults[0].RequiredSafetyFactor, Is.EqualTo(1.170).Within(tolerance3Decimals));
                Assert.That(allCalculationResults[0].RiverLevel, Is.EqualTo(4.4).Within(tolerance1Decimals));
                Assert.That(allCalculationResults[0].RiverLevelLow, Is.EqualTo(0.63).Within(tolerance2Decimals));
                Assert.That(allCalculationResults[0].DikeTableHeight, Is.EqualTo(4.900).Within(tolerance3Decimals));
                Assert.That(allCalculationResults[0].SlopeDampingPiezometricHeightPolderSide, Is.EqualTo(0.000).Within(tolerance3Decimals));
                Assert.That(allCalculationResults[0].DikeLength, Is.EqualTo(36.150).Within(tolerance3Decimals));

                Assert.That(allCalculationResults[0].Pl3MinUplift, Is.EqualTo(null));
                Assert.That(allCalculationResults[0].Pl3HeadAdjusted, Is.EqualTo(null));
                Assert.That(allCalculationResults[0].Pl3LocalLocationXMinUplift, Is.EqualTo(null));
                Assert.That(allCalculationResults[0].Pl3LocationXMinUplift, Is.EqualTo(null));
                Assert.That(allCalculationResults[0].Pl3LocationYMinUplift, Is.EqualTo(null));
                Assert.That(allCalculationResults[0].Pl4MinUplift, Is.EqualTo(null));
                Assert.That(allCalculationResults[0].Pl4HeadAdjusted, Is.EqualTo(null));
                Assert.That(allCalculationResults[0].Pl4LocalLocationXMinUplift, Is.EqualTo(null));
                Assert.That(allCalculationResults[0].Pl4LocationXMinUplift, Is.EqualTo(null));
                Assert.That(allCalculationResults[0].Pl4LocationYMinUplift, Is.EqualTo(null));

                Assert.That(allCalculationResults[0].LocalPipingEntryPointX, Is.EqualTo(16.190).Within(tolerance3Decimals));
                Assert.That(allCalculationResults[0].PipingEntryPointX, Is.EqualTo(24.610).Within(tolerance3Decimals));
                Assert.That(allCalculationResults[0].PipingEntryPointY, Is.EqualTo(0.000).Within(tolerance3Decimals));
                Assert.That(allCalculationResults[0].LocalPipingExitPointX, Is.EqualTo(null));
                Assert.That(allCalculationResults[0].PipingExitPointX, Is.EqualTo(null));
                Assert.That(allCalculationResults[0].PipingExitPointY, Is.EqualTo(null));
                Assert.That(allCalculationResults[0].SeepageLength, Is.EqualTo(null));
                Assert.That(allCalculationResults[0].ResultMessage, Is.EqualTo(""));
                Assert.That(allCalculationResults[0].HeaveFactor, Is.EqualTo(null));
                Assert.That(allCalculationResults[0].Notes, Is.EqualTo(""));
                Assert.That(allCalculationResults[0].NumberOfIterations, Is.EqualTo(0));
                Assert.That(allCalculationResults[0].ResultEvaluation, Is.EqualTo(ResultEvaluation.NotEvaluated));
                Assert.That(allCalculationResults[0].UpliftFactor, Is.EqualTo(null));
            });
        }

        /// <summary>
        /// The safety factors calculated in version 19.2.1 are different from the actual calculated safety factors.
        /// This can be explained by the Sigma-Tau table model which is differently implemented:
        /// in version 19.2.1 Sigma is the effective vertical stress whereas in the current version it is the effective normal stress.
        /// </summary>
        [Test]
        public void TestStabilityOutsideWith1DProfiles()
        {
            string projectFilename = Path.Combine(Directory.GetCurrentDirectory(), @"TestData\StabilityOutside1DProfileTest\Actualisatie.damx");
            List<CsvExportData> allCalculationResults = ComputeStabilityOutsideProject(projectFilename, 59);
            
            Assert.That(allCalculationResults, Has.Count.EqualTo(36));
            Assert.Multiple(() =>
            {
                Assert.That(allCalculationResults[1].LocationName, Is.EqualTo("RK210-560"));
                Assert.That(allCalculationResults[1].CalculationResult, Is.EqualTo(CalculationResult.Succeeded));
                Assert.That(allCalculationResults[1].X, Is.EqualTo(88218.047).Within(tolerance3Decimals));
                Assert.That(allCalculationResults[1].Y, Is.EqualTo(448871.233).Within(tolerance3Decimals));
                Assert.That(allCalculationResults[1].AnalysisType, Is.EqualTo(AnalysisType.NoAdaption));
                Assert.That(allCalculationResults[1].IsUplift, Is.EqualTo(false));
                Assert.That(allCalculationResults[1].ProfileName, Is.EqualTo("Segment_1_1D2"));
                Assert.That(allCalculationResults[1].ProfileProbability, Is.EqualTo(2.43));
                Assert.That(allCalculationResults[1].StabilityModel, Is.EqualTo(StabilityModelType.Bishop));
                // In version 19.2.1, the safety factor was 1.079
                Assert.That(allCalculationResults[1].SafetyFactor, Is.EqualTo(1.129).Within(tolerance3Decimals));
                Assert.That(allCalculationResults[1].ShoulderHeight, Is.EqualTo(null));
                Assert.That(allCalculationResults[1].ToeAtPolderX, Is.EqualTo(null));
                Assert.That(allCalculationResults[1].ToeAtPolderZ, Is.EqualTo(null));
                Assert.That(allCalculationResults[1].RequiredSafetyFactor, Is.EqualTo(0.900));
                Assert.That(allCalculationResults[1].RiverLevel, Is.EqualTo(-2.090));
                Assert.That(allCalculationResults[1].RiverLevelLow, Is.EqualTo(-2.520));
                Assert.That(allCalculationResults[1].DikeTableHeight, Is.EqualTo(-2.000).Within(tolerance3Decimals));
                Assert.That(allCalculationResults[1].SlopeDampingPiezometricHeightPolderSide, Is.EqualTo(0.000));
                Assert.That(allCalculationResults[1].DikeLength, Is.EqualTo(14.12).Within(tolerance3Decimals * 10));

                Assert.That(allCalculationResults[1].Pl3MinUplift, Is.EqualTo(null));
                Assert.That(allCalculationResults[1].Pl3HeadAdjusted, Is.EqualTo(null));
                Assert.That(allCalculationResults[1].Pl3LocalLocationXMinUplift, Is.EqualTo(null));
                Assert.That(allCalculationResults[1].Pl3LocationXMinUplift, Is.EqualTo(null));
                Assert.That(allCalculationResults[1].Pl3LocationYMinUplift, Is.EqualTo(null));
                Assert.That(allCalculationResults[1].Pl4MinUplift, Is.EqualTo(null));
                Assert.That(allCalculationResults[1].Pl4HeadAdjusted, Is.EqualTo(null));
                Assert.That(allCalculationResults[1].Pl4LocalLocationXMinUplift, Is.EqualTo(null));
                Assert.That(allCalculationResults[1].Pl4LocationXMinUplift, Is.EqualTo(null));
                Assert.That(allCalculationResults[1].Pl4LocationYMinUplift, Is.EqualTo(null));

                Assert.That(allCalculationResults[1].LocalPipingEntryPointX, Is.EqualTo(9.798).Within(tolerance3Decimals));
                Assert.That(allCalculationResults[1].PipingEntryPointX, Is.EqualTo(88216.844).Within(tolerance3Decimals));
                Assert.That(allCalculationResults[1].PipingEntryPointY, Is.EqualTo(448870.287).Within(tolerance3Decimals));
                Assert.That(allCalculationResults[1].LocalPipingExitPointX, Is.EqualTo(null));
                Assert.That(allCalculationResults[1].PipingExitPointX, Is.EqualTo(null));
                Assert.That(allCalculationResults[1].PipingExitPointY, Is.EqualTo(null));
                Assert.That(allCalculationResults[1].SeepageLength, Is.EqualTo(null));
                Assert.That(allCalculationResults[1].ResultMessage, Is.EqualTo(""));
                Assert.That(allCalculationResults[1].HeaveFactor, Is.EqualTo(null));
                Assert.That(allCalculationResults[1].Notes, Is.EqualTo(""));
                Assert.That(allCalculationResults[1].NumberOfIterations, Is.EqualTo(0));
                Assert.That(allCalculationResults[1].ResultEvaluation, Is.EqualTo(ResultEvaluation.NotEvaluated));
                Assert.That(allCalculationResults[1].UpliftFactor, Is.EqualTo(null));

                // Check safety factors of few other soil profiles
                // In version 19.2.1, safety factors were: [4]=1.079, [8]=0.967, [10]=1.046, [11]=1.225 [12]=1.046
                Assert.That(allCalculationResults[4].ProfileName, Is.EqualTo("Segment_1_1D5"));
                Assert.That(allCalculationResults[4].SafetyFactor, Is.EqualTo(1.129).Within(tolerance3Decimals));
                Assert.That(allCalculationResults[8].ProfileName, Is.EqualTo("Segment_1_1D9"));
                Assert.That(allCalculationResults[8].SafetyFactor, Is.EqualTo(1.049).Within(tolerance3Decimals));
                Assert.That(allCalculationResults[10].ProfileName, Is.EqualTo("Segment_1_1D11"));
                Assert.That(allCalculationResults[10].SafetyFactor, Is.EqualTo(1.111).Within(tolerance3Decimals));
                Assert.That(allCalculationResults[11].ProfileName, Is.EqualTo("Segment_1_1D12"));
                Assert.That(allCalculationResults[11].SafetyFactor, Is.EqualTo(1.253).Within(tolerance3Decimals));
                Assert.That(allCalculationResults[12].ProfileName, Is.EqualTo("Segment_1_1D13"));
                Assert.That(allCalculationResults[12].SafetyFactor, Is.EqualTo(1.111).Within(tolerance3Decimals));
            });
        }

        /// <summary>
        /// Compute stability output project
        /// </summary>
        /// <param name="projectFilename"></param>
        /// <param name="expectedLocations"></param>
        /// <returns></returns>
        private static List<CsvExportData> ComputeStabilityOutsideProject(string projectFilename, int expectedLocations)
        {
            return ComputeHelper.ComputeStabilityProject(projectFilename, StabilityModelType.Bishop, false, expectedLocations, 1, out List<LogMessage> _);
        }
    }
}