// 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 System.Linq;
using Core.Common.Base.Geometry;
using Core.Common.TestUtil;
using NUnit.Framework;
using Rhino.Mocks;
using Ringtoets.ClosingStructures.Data;
using Ringtoets.ClosingStructures.Data.TestUtil;
using Ringtoets.Common.Data.AssessmentSection;
using Ringtoets.Common.Data.Contribution;
using Ringtoets.Common.Data.DikeProfiles;
using Ringtoets.Common.Data.FailureMechanism;
using Ringtoets.HydraRing.Calculation.Calculator.Factory;
using Ringtoets.HydraRing.Calculation.Data;
using Ringtoets.HydraRing.Calculation.Data.Input.Structures;
using Ringtoets.HydraRing.Calculation.TestUtil;
using Ringtoets.HydraRing.Calculation.TestUtil.Calculator;
using Ringtoets.HydraRing.Data;
namespace Ringtoets.ClosingStructures.Service.Test
{
[TestFixture]
public class ClosingStructuresCalculationServiceTest
{
private static readonly string testDataPath = TestHelper.GetTestDataPath(TestDataPath.Ringtoets.Integration.Service, "HydraRingCalculation");
[Test]
public void Calculate_InvalidInFlowModelType_ThrowsNotSupportedException()
{
// Setup
var closingStructuresFailureMechanism = new ClosingStructuresFailureMechanism();
var mockRepository = new MockRepository();
var assessmentSectionStub = mockRepository.Stub();
mockRepository.ReplayAll();
closingStructuresFailureMechanism.AddSection(new FailureMechanismSection("test section", new[]
{
new Point2D(0, 0),
new Point2D(1, 1)
}));
FailureMechanismSection failureMechanismSection = closingStructuresFailureMechanism.Sections.First();
ClosingStructuresCalculation calculation = new TestClosingStructuresCalculation()
{
InputParameters =
{
InflowModelType = (ClosingStructureInflowModelType) 100
}
};
var service = new ClosingStructuresCalculationService();
// Call
using (new HydraRingCalculatorFactoryConfig())
{
var calculator = ((TestHydraRingCalculatorFactory) HydraRingCalculatorFactory.Instance).StructuresClosureCalculator;
TestDelegate call = () => service.Calculate(calculation, assessmentSectionStub, failureMechanismSection,
closingStructuresFailureMechanism.GeneralInput, closingStructuresFailureMechanism.Contribution,
testDataPath);
StructuresClosureCalculationInput[] calculationInputs = calculator.ReceivedInputs.ToArray();
Assert.AreEqual(0, calculationInputs.Length);
var exception = Assert.Throws(call);
Assert.AreEqual("ClosingStructureInflowModelType", exception.Message);
}
mockRepository.VerifyAll();
}
[Test]
[TestCase(true, false)]
[TestCase(true, true)]
[TestCase(false, false)]
public void Calculate_VariousVerticalWallCalculations_InputPropertiesCorrectlySentToCalculator(bool useForeshore, bool useBreakWater)
{
// Setup
var closingStructuresFailureMechanism = new ClosingStructuresFailureMechanism();
var mockRepository = new MockRepository();
var assessmentSectionStub = CreateAssessmentSectionStub(closingStructuresFailureMechanism, mockRepository);
mockRepository.ReplayAll();
closingStructuresFailureMechanism.AddSection(new FailureMechanismSection("test section", new[]
{
new Point2D(0, 0),
new Point2D(1, 1)
}));
ClosingStructuresCalculation calculation = new TestClosingStructuresCalculation()
{
InputParameters =
{
HydraulicBoundaryLocation = assessmentSectionStub.HydraulicBoundaryDatabase.Locations.First(hl => hl.Id == 1300001)
}
};
if (useForeshore)
{
calculation.InputParameters.ForeshoreProfile = new ForeshoreProfile(new Point2D(0, 0),
new[]
{
new Point2D(1, 1),
new Point2D(2, 2)
},
useBreakWater ? new BreakWater(BreakWaterType.Wall, 3.0) : null,
new ForeshoreProfile.ConstructionProperties());
}
FailureMechanismSection failureMechanismSection = closingStructuresFailureMechanism.Sections.First();
using (new HydraRingCalculatorFactoryConfig())
{
var calculator = ((TestHydraRingCalculatorFactory) HydraRingCalculatorFactory.Instance).StructuresClosureCalculator;
// Call
new ClosingStructuresCalculationService().Calculate(calculation,
assessmentSectionStub,
failureMechanismSection,
closingStructuresFailureMechanism.GeneralInput,
closingStructuresFailureMechanism.Contribution,
testDataPath);
// Assert
StructuresClosureCalculationInput[] calculationInputs = calculator.ReceivedInputs.ToArray();
Assert.AreEqual(1, calculationInputs.Length);
Assert.AreEqual(testDataPath, calculator.HydraulicBoundaryDatabaseDirectory);
Assert.AreEqual(assessmentSectionStub.Id, calculator.RingId);
GeneralClosingStructuresInput generalInput = closingStructuresFailureMechanism.GeneralInput;
ClosingStructuresInput input = calculation.InputParameters;
var expectedInput = new StructuresClosureVerticalWallCalculationInput(
1300001,
new HydraRingSection(1, failureMechanismSection.GetSectionLength(), input.StructureNormalOrientation),
useForeshore ? input.ForeshoreGeometry.Select(c => new HydraRingForelandPoint(c.X, c.Y)) : new HydraRingForelandPoint[0],
useBreakWater ? new HydraRingBreakWater((int) input.BreakWater.Type, input.BreakWater.Height) : null,
generalInput.GravitationalAcceleration,
input.FactorStormDurationOpenStructure,
input.FailureProbabilityOpenStructure,
input.FailureProbabilityReparation,
input.IdenticalApertures,
input.AllowedLevelIncreaseStorage.Mean, input.AllowedLevelIncreaseStorage.StandardDeviation,
generalInput.ModelFactorStorageVolume.Mean, generalInput.ModelFactorStorageVolume.StandardDeviation,
input.StorageStructureArea.Mean, input.StorageStructureArea.CoefficientOfVariation,
generalInput.ModelFactorInflowVolume,
input.FlowWidthAtBottomProtection.Mean, input.FlowWidthAtBottomProtection.StandardDeviation,
input.CriticalOvertoppingDischarge.Mean, input.CriticalOvertoppingDischarge.CoefficientOfVariation,
input.FailureProbabilityStructureWithErosion,
input.StormDuration.Mean, input.StormDuration.CoefficientOfVariation,
input.ProbabilityOpenStructureBeforeFlooding,
generalInput.ModelFactorOvertoppingFlow.Mean, generalInput.ModelFactorOvertoppingFlow.StandardDeviation,
input.StructureNormalOrientation,
input.ModelFactorSuperCriticalFlow.Mean, input.ModelFactorSuperCriticalFlow.StandardDeviation,
input.LevelCrestStructureNotClosing.Mean, input.LevelCrestStructureNotClosing.StandardDeviation,
input.WidthFlowApertures.Mean, input.WidthFlowApertures.CoefficientOfVariation,
input.DeviationWaveDirection);
StructuresClosureVerticalWallCalculationInput actualInput = (StructuresClosureVerticalWallCalculationInput) calculationInputs[0];
HydraRingDataEqualityHelper.AreEqual(expectedInput, actualInput);
Assert.IsFalse(calculator.IsCanceled);
}
mockRepository.VerifyAll();
}
[Test]
[TestCase(true, false)]
[TestCase(true, true)]
[TestCase(false, false)]
public void Calculate_VariousLowSillCalculations_InputPropertiesCorrectlySentToCalculator(bool useForeshore, bool useBreakWater)
{
// Setup
var closingStructuresFailureMechanism = new ClosingStructuresFailureMechanism();
var mockRepository = new MockRepository();
var assessmentSectionStub = CreateAssessmentSectionStub(closingStructuresFailureMechanism, mockRepository);
mockRepository.ReplayAll();
closingStructuresFailureMechanism.AddSection(new FailureMechanismSection("test section", new[]
{
new Point2D(0, 0),
new Point2D(1, 1)
}));
ClosingStructuresCalculation calculation = new TestClosingStructuresCalculation()
{
InputParameters =
{
HydraulicBoundaryLocation = assessmentSectionStub.HydraulicBoundaryDatabase.Locations.First(hl => hl.Id == 1300001),
InflowModelType = ClosingStructureInflowModelType.LowSill
}
};
if (useForeshore)
{
calculation.InputParameters.ForeshoreProfile = new ForeshoreProfile(new Point2D(0, 0),
new[]
{
new Point2D(1, 1),
new Point2D(2, 2)
},
useBreakWater ? new BreakWater(BreakWaterType.Wall, 3.0) : null,
new ForeshoreProfile.ConstructionProperties());
}
FailureMechanismSection failureMechanismSection = closingStructuresFailureMechanism.Sections.First();
using (new HydraRingCalculatorFactoryConfig())
{
var calculator = ((TestHydraRingCalculatorFactory) HydraRingCalculatorFactory.Instance).StructuresClosureCalculator;
// Call
new ClosingStructuresCalculationService().Calculate(calculation,
assessmentSectionStub,
failureMechanismSection,
closingStructuresFailureMechanism.GeneralInput,
closingStructuresFailureMechanism.Contribution,
testDataPath);
// Assert
StructuresClosureCalculationInput[] calculationInputs = calculator.ReceivedInputs.ToArray();
Assert.AreEqual(1, calculationInputs.Length);
Assert.AreEqual(testDataPath, calculator.HydraulicBoundaryDatabaseDirectory);
Assert.AreEqual(assessmentSectionStub.Id, calculator.RingId);
GeneralClosingStructuresInput generalInput = closingStructuresFailureMechanism.GeneralInput;
ClosingStructuresInput input = calculation.InputParameters;
var expectedInput = new StructuresClosureLowSillCalculationInput(
1300001,
new HydraRingSection(1, failureMechanismSection.GetSectionLength(), input.StructureNormalOrientation),
useForeshore ? input.ForeshoreGeometry.Select(c => new HydraRingForelandPoint(c.X, c.Y)) : new HydraRingForelandPoint[0],
useBreakWater ? new HydraRingBreakWater((int) input.BreakWater.Type, input.BreakWater.Height) : null,
generalInput.GravitationalAcceleration,
input.FactorStormDurationOpenStructure,
input.FailureProbabilityOpenStructure,
input.FailureProbabilityReparation,
input.IdenticalApertures,
input.AllowedLevelIncreaseStorage.Mean, input.AllowedLevelIncreaseStorage.StandardDeviation,
generalInput.ModelFactorStorageVolume.Mean, generalInput.ModelFactorStorageVolume.StandardDeviation,
input.StorageStructureArea.Mean, input.StorageStructureArea.CoefficientOfVariation,
generalInput.ModelFactorInflowVolume,
input.FlowWidthAtBottomProtection.Mean, input.FlowWidthAtBottomProtection.StandardDeviation,
input.CriticalOvertoppingDischarge.Mean, input.CriticalOvertoppingDischarge.CoefficientOfVariation,
input.FailureProbabilityStructureWithErosion,
input.StormDuration.Mean, input.StormDuration.CoefficientOfVariation,
input.ProbabilityOpenStructureBeforeFlooding,
input.ModelFactorSuperCriticalFlow.Mean, input.ModelFactorSuperCriticalFlow.StandardDeviation,
generalInput.ModelFactorSubCriticalFlow.Mean, generalInput.ModelFactorSubCriticalFlow.CoefficientOfVariation,
input.ThresholdHeightOpenWeir.Mean, input.ThresholdHeightOpenWeir.StandardDeviation,
input.InsideWaterLevel.Mean, input.InsideWaterLevel.StandardDeviation,
input.WidthFlowApertures.Mean, input.WidthFlowApertures.CoefficientOfVariation);
StructuresClosureLowSillCalculationInput actualInput = (StructuresClosureLowSillCalculationInput) calculationInputs[0];
HydraRingDataEqualityHelper.AreEqual(expectedInput, actualInput);
Assert.IsFalse(calculator.IsCanceled);
}
mockRepository.VerifyAll();
}
[Test]
[TestCase(true, false)]
[TestCase(true, true)]
[TestCase(false, false)]
public void Calculate_VariousFloodedCulvertCalculations_InputPropertiesCorrectlySentToCalculator(bool useForeshore, bool useBreakWater)
{
// Setup
var closingStructuresFailureMechanism = new ClosingStructuresFailureMechanism();
var mockRepository = new MockRepository();
var assessmentSectionStub = CreateAssessmentSectionStub(closingStructuresFailureMechanism, mockRepository);
mockRepository.ReplayAll();
closingStructuresFailureMechanism.AddSection(new FailureMechanismSection("test section", new[]
{
new Point2D(0, 0),
new Point2D(1, 1)
}));
ClosingStructuresCalculation calculation = new TestClosingStructuresCalculation()
{
InputParameters =
{
HydraulicBoundaryLocation = assessmentSectionStub.HydraulicBoundaryDatabase.Locations.First(hl => hl.Id == 1300001),
InflowModelType = ClosingStructureInflowModelType.LowSill
}
};
if (useForeshore)
{
calculation.InputParameters.ForeshoreProfile = new ForeshoreProfile(new Point2D(0, 0),
new[]
{
new Point2D(1, 1),
new Point2D(2, 2)
},
useBreakWater ? new BreakWater(BreakWaterType.Wall, 3.0) : null,
new ForeshoreProfile.ConstructionProperties());
}
FailureMechanismSection failureMechanismSection = closingStructuresFailureMechanism.Sections.First();
using (new HydraRingCalculatorFactoryConfig())
{
var calculator = ((TestHydraRingCalculatorFactory) HydraRingCalculatorFactory.Instance).StructuresClosureCalculator;
// Call
new ClosingStructuresCalculationService().Calculate(calculation,
assessmentSectionStub,
failureMechanismSection,
closingStructuresFailureMechanism.GeneralInput,
closingStructuresFailureMechanism.Contribution,
testDataPath);
// Assert
StructuresClosureCalculationInput[] calculationInputs = calculator.ReceivedInputs.ToArray();
Assert.AreEqual(1, calculationInputs.Length);
Assert.AreEqual(testDataPath, calculator.HydraulicBoundaryDatabaseDirectory);
Assert.AreEqual(assessmentSectionStub.Id, calculator.RingId);
GeneralClosingStructuresInput generalInput = closingStructuresFailureMechanism.GeneralInput;
ClosingStructuresInput input = calculation.InputParameters;
var expectedInput = new StructuresClosureFloodedCulvertCalculationInput(
1300001,
new HydraRingSection(1, failureMechanismSection.GetSectionLength(), input.StructureNormalOrientation),
useForeshore ? input.ForeshoreGeometry.Select(c => new HydraRingForelandPoint(c.X, c.Y)) : new HydraRingForelandPoint[0],
useBreakWater ? new HydraRingBreakWater((int) input.BreakWater.Type, input.BreakWater.Height) : null,
generalInput.GravitationalAcceleration,
input.FactorStormDurationOpenStructure,
input.FailureProbabilityOpenStructure,
input.FailureProbabilityReparation,
input.IdenticalApertures,
input.AllowedLevelIncreaseStorage.Mean, input.AllowedLevelIncreaseStorage.StandardDeviation,
generalInput.ModelFactorStorageVolume.Mean, generalInput.ModelFactorStorageVolume.StandardDeviation,
input.StorageStructureArea.Mean, input.StorageStructureArea.CoefficientOfVariation,
generalInput.ModelFactorInflowVolume,
input.FlowWidthAtBottomProtection.Mean, input.FlowWidthAtBottomProtection.StandardDeviation,
input.CriticalOvertoppingDischarge.Mean, input.CriticalOvertoppingDischarge.CoefficientOfVariation,
input.FailureProbabilityStructureWithErosion,
input.StormDuration.Mean, input.StormDuration.CoefficientOfVariation,
input.ProbabilityOpenStructureBeforeFlooding,
input.DrainCoefficient.Mean, input.DrainCoefficient.StandardDeviation,
input.AreaFlowApertures.Mean, input.AreaFlowApertures.StandardDeviation,
input.InsideWaterLevel.Mean, input.InsideWaterLevel.StandardDeviation);
StructuresClosureFloodedCulvertCalculationInput actualInput = (StructuresClosureFloodedCulvertCalculationInput) calculationInputs[0];
HydraRingDataEqualityHelper.AreEqual(expectedInput, actualInput);
Assert.IsFalse(calculator.IsCanceled);
}
mockRepository.VerifyAll();
}
private static IAssessmentSection CreateAssessmentSectionStub(IFailureMechanism failureMechanism, MockRepository mockRepository)
{
var assessmentSectionStub = mockRepository.Stub();
assessmentSectionStub.Stub(a => a.Id).Return("21");
assessmentSectionStub.Stub(a => a.FailureMechanismContribution).Return(new FailureMechanismContribution(new[]
{
failureMechanism
}, 1, 2));
assessmentSectionStub.HydraulicBoundaryDatabase = new HydraulicBoundaryDatabase
{
Locations =
{
new HydraulicBoundaryLocation(1300001, string.Empty, 0, 0)
}
};
return assessmentSectionStub;
}
}
}