// Copyright (C) Stichting Deltares 2016. 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 . // // 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 Core.Common.TestUtil; using NUnit.Framework; using Ringtoets.Common.Data.Probabilistics; using Ringtoets.Common.Data.TestUtil; namespace Ringtoets.Common.Data.Test.Probabilistics { [TestFixture] public class LogNormalDistributionTest { [Test] public void DefaultConstructor_ExpectedValues() { // Call var distribution = new LogNormalDistribution(); // Assert Assert.IsInstanceOf(distribution); int numberOfDecimalPlaces = RoundedDouble.MaximumNumberOfDecimalPlaces; double expectedAccuracy = Math.Pow(10.0, -numberOfDecimalPlaces); Assert.AreEqual(Math.Exp(-0.5), distribution.Mean, expectedAccuracy); Assert.AreEqual(numberOfDecimalPlaces, distribution.Mean.NumberOfDecimalPlaces); Assert.AreEqual(Math.Sqrt((Math.Exp(1) - 1) * Math.Exp(1)), distribution.StandardDeviation, expectedAccuracy); Assert.AreEqual(numberOfDecimalPlaces, distribution.StandardDeviation.NumberOfDecimalPlaces); Assert.AreEqual(0.0, distribution.Shift, expectedAccuracy); Assert.AreEqual(numberOfDecimalPlaces, distribution.Shift.NumberOfDecimalPlaces); } [Test] [TestCase(1)] [TestCase(9)] [TestCase(15)] public void Constructor_WithParameter_ExpectedValues(int numberOfDecimalPlaces) { // Call var distribution = new LogNormalDistribution(numberOfDecimalPlaces); // Assert Assert.IsInstanceOf(distribution); double expectedAccuracy = Math.Pow(10.0, -numberOfDecimalPlaces); Assert.AreEqual(Math.Exp(-0.5), distribution.Mean, expectedAccuracy); Assert.AreEqual(numberOfDecimalPlaces, distribution.Mean.NumberOfDecimalPlaces); Assert.AreEqual(Math.Sqrt((Math.Exp(1) - 1)*Math.Exp(1)), distribution.StandardDeviation, expectedAccuracy); Assert.AreEqual(numberOfDecimalPlaces, distribution.StandardDeviation.NumberOfDecimalPlaces); Assert.AreEqual(0.0, distribution.Shift, expectedAccuracy); Assert.AreEqual(numberOfDecimalPlaces, distribution.Shift.NumberOfDecimalPlaces); } [Test] public void Constructor_InvalidNumberOfDecimalPlaces_ThrowArgumentOutOfRangeException() { // Call TestDelegate call = () => new LogNormalDistribution(0); // Assert const string expectedMessage = "Value must be in range [1, 15]."; TestHelper.AssertThrowsArgumentExceptionAndTestMessage(call, expectedMessage); } [Test] public void Mean_SetNewValue_GetValueRoundedToGivenNumberOfDecimalPlaces() { // Setup var value = 1.23456789; var numberOfDecimalPlaces = 2; var distribution = new LogNormalDistribution(numberOfDecimalPlaces); // Call distribution.Mean = (RoundedDouble) value; // Assert Assert.AreEqual(numberOfDecimalPlaces, distribution.Mean.NumberOfDecimalPlaces); Assert.AreEqual(new RoundedDouble(numberOfDecimalPlaces, value), distribution.Mean); } [Test] [TestCase(0)] [TestCase(0 + 1e-3, Description = "Invalid mean due to rounding to 0.0")] [TestCase(-123.45)] public void Mean_SettingToLessThanOrEqualTo0_ThrowArgumentOutOfRangeException(double mean) { // Setup var distribution = new LogNormalDistribution(2); // Call TestDelegate call = () => distribution.Mean = (RoundedDouble) mean; // Assert const string expectedMessage = "Gemiddelde moet groter zijn dan 0."; TestHelper.AssertThrowsArgumentExceptionAndTestMessage(call, expectedMessage); } [Test] [TestCase(1.23456789)] [TestCase(0 - 1e-3, Description = "Valid standard deviation due to rounding to 0.0")] public void StandardDeviation_SetNewValue_GetValueRoundedToGivenNumberOfDecimalPlaces(double standardDeviation) { // Setup var numberOfDecimalPlaces = 2; var distribution = new LogNormalDistribution(numberOfDecimalPlaces); // Call distribution.StandardDeviation = (RoundedDouble) standardDeviation; // Assert Assert.AreEqual(numberOfDecimalPlaces, distribution.StandardDeviation.NumberOfDecimalPlaces); Assert.AreEqual(new RoundedDouble(numberOfDecimalPlaces, standardDeviation), distribution.StandardDeviation); } [Test] [TestCase(-4)] [TestCase(0 - 1e-2)] public void StandardDeviation_SettingToLessThan0_ThrowArgumentOutOfRangeException(double standardDeviation) { // Setup var distribution = new LogNormalDistribution(2); // Call TestDelegate call = () => distribution.StandardDeviation = (RoundedDouble) standardDeviation; // Assert const string expectedMessage = "Standaardafwijking (\u03C3) moet groter zijn dan of gelijk zijn aan 0."; TestHelper.AssertThrowsArgumentExceptionAndTestMessage(call, expectedMessage); } [Test] [TestCase(1, 5.6)] [TestCase(3, 5.647)] [TestCase(4, 5.6473)] [TestCase(15, 5.647300000000000)] public void Shift_SetNewValue_GetValueRoundedToGivenNumberOfDecimalPlaces(int numberOfDecimalPlaces, double expectedStandardDeviation) { // Setup var distribution = new LogNormalDistribution(numberOfDecimalPlaces) { Mean = new RoundedDouble(2, 10.0) }; // Call distribution.Shift = new RoundedDouble(4, 5.6473); // Assert Assert.AreEqual(numberOfDecimalPlaces, distribution.Shift.NumberOfDecimalPlaces); Assert.AreEqual(expectedStandardDeviation, distribution.Shift.Value); } [Test] public void Shift_SetIllegalValue_ThrowArgumentOutOfRangeException() { // Setup var distribution = new LogNormalDistribution(2) { Mean = new RoundedDouble(2, 10.0), StandardDeviation = new RoundedDouble(2, 1.0) }; // Call TestDelegate call = () => distribution.Shift = new RoundedDouble(2, 100.0); // Assert const string expectedMessage = "De verschuiving mag niet groter zijn dan de verwachtingswaarde."; TestHelper.AssertThrowsArgumentExceptionAndTestMessage(call, expectedMessage); } [Test] public void Clone_Always_ReturnNewInstanceWithCopiedValues() { // Setup var random = new Random(21); var distribution = new LogNormalDistribution(random.Next(1, 16)) { Mean = random.NextRoundedDouble(), StandardDeviation = random.NextRoundedDouble() }; // Call object clone = distribution.Clone(); // Assert Assert.IsInstanceOf(clone); var clonedDistribution = (LogNormalDistribution)clone; DistributionAssert.AreEqual(distribution, clonedDistribution); } } }