// Copyright (C) Stichting Deltares 2018. All rights reserved.
//
// This file is part of Ringtoets.
//
// Ringtoets 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;
using Core.Common.Base.Data;
using NUnit.Framework;
using Ringtoets.Common.Data.Probabilistics;
using Ringtoets.Common.Data.TestUtil;

namespace Ringtoets.Common.Data.Test.Probabilistics
{
    [TestFixture]
    public class VariationCoefficientLogNormalDistributionDesignVariableTest
    {
        [Test]
        public void ParameteredConstructor_ValidLogNormalDistribution_ExpectedValues()
        {
            // Setup
            var logNormalDistribution = new VariationCoefficientLogNormalDistribution(2);

            // Call
            var designVariable = new VariationCoefficientLogNormalDistributionDesignVariable(logNormalDistribution);

            // Assert
            Assert.IsInstanceOf<VariationCoefficientDesignVariable<VariationCoefficientLogNormalDistribution>>(designVariable);
            Assert.AreSame(logNormalDistribution, designVariable.Distribution);
            Assert.AreEqual(0.5, designVariable.Percentile);
        }

        [Test]
        [TestCase(75, 70, 0.95, 89.50524481)]
        [TestCase(75, 70, 0.5, 74.53727185)]
        [TestCase(75, 70, 0.05, 62.07239483)]
        [TestCase(75, 123.45, 0.95, 94.53050044)]
        [TestCase(75, 1.2345, 0.95, 76.83967484)]
        [TestCase(123.45, 70, 0.95, 137.6805705)]
        [TestCase(1.2345, 70, 0.95, 4.541272847)]
        public void GetDesignValue_ValidLogNormalDistribution_ReturnExpectedValue(
            double expectedValue, double variance, double percentile,
            double expectedResult)
        {
            // Setup
            const int numberOfDecimalPlaces = 4;
            var logNormalDistribution = new VariationCoefficientLogNormalDistribution(numberOfDecimalPlaces)
            {
                Mean = (RoundedDouble) expectedValue,
                CoefficientOfVariation = (RoundedDouble) (Math.Sqrt(variance) / expectedValue)
            };

            var designVariable = new VariationCoefficientLogNormalDistributionDesignVariable(logNormalDistribution)
            {
                Percentile = percentile
            };

            // Call
            RoundedDouble result = designVariable.GetDesignValue();

            // Assert
            Assert.AreEqual(numberOfDecimalPlaces, result.NumberOfDecimalPlaces);
            Assert.AreEqual(expectedResult, result, result.GetAccuracy());
        }

        [Test]
        [TestCase(75, 20, 10, 0.95, 260.1746)]
        [TestCase(75, 70, 10, 0.5, 17.71409)]
        [TestCase(75, 70, 10, 0.95, 240.2244)]
        [TestCase(75, 70, 10, 0.05, 10.2585)]
        [TestCase(75, 70, -30, 0.95, 341.901)]
        [TestCase(75, 123.45, 10, 0.95, 225.9596)]
        [TestCase(75, 1.2345, 10, 0.95, 200.0482)]
        [TestCase(123.45, 70, 10, 0.95, 411.83019)]
        public void GetDesignValue_ValidLogNormalDistributionWithNonZeroShift_ReturnExpectedValue(
            double expectedValue, double variance, double shift, double percentile,
            double expectedResult)
        {
            // Setup
            const int numberOfDecimalPlaces = 4;
            var logNormalDistribution = new VariationCoefficientLogNormalDistribution(numberOfDecimalPlaces)
            {
                Mean = (RoundedDouble) expectedValue,
                CoefficientOfVariation = (RoundedDouble) Math.Sqrt(variance),
                Shift = (RoundedDouble) shift
            };

            var designVariable = new VariationCoefficientLogNormalDistributionDesignVariable(logNormalDistribution)
            {
                Percentile = percentile
            };

            // Call
            RoundedDouble result = designVariable.GetDesignValue();

            // Assert
            Assert.AreEqual(numberOfDecimalPlaces, result.NumberOfDecimalPlaces);
            Assert.AreEqual(expectedResult, result, result.GetAccuracy());
        }
    }
}