// Copyright (C) Stichting Deltares 2017. 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.Collections.Generic; using System.IO; using System.Linq; using Core.Common.Base.IO; using Core.Common.Utils; using log4net; using Ringtoets.Common.Data.Hydraulics; using Ringtoets.Common.Data.IllustrationPoints; using Ringtoets.Common.IO.HydraRing; using Ringtoets.Common.Service.IllustrationPoints; using Ringtoets.Common.Service.MessageProviders; using Ringtoets.Common.Service.Properties; using Ringtoets.HydraRing.Calculation.Calculator; using Ringtoets.HydraRing.Calculation.Calculator.Factory; using Ringtoets.HydraRing.Calculation.Data.Input.Hydraulics; using Ringtoets.HydraRing.Calculation.Exceptions; using HydraRingGeneralResult = Ringtoets.HydraRing.Calculation.Data.Output.IllustrationPoints.GeneralResult; namespace Ringtoets.Common.Service { /// /// Service that provides methods for performing Hydra-Ring calculations for design water level. /// public class DesignWaterLevelCalculationService { private static readonly ILog log = LogManager.GetLogger(typeof(DesignWaterLevelCalculationService)); private IDesignWaterLevelCalculator calculator; private bool canceled; /// /// Performs validation over the values on the given . /// Error and status information is logged during the execution of the operation. /// /// The file path of the hydraulic boundary database file which to validate. /// True if there were no validation errors; False otherwise. public static bool Validate(string hydraulicBoundaryDatabaseFilePath) { CalculationServiceHelper.LogValidationBegin(); string[] validationProblems = ValidateInput(hydraulicBoundaryDatabaseFilePath); CalculationServiceHelper.LogMessagesAsError(Resources.Hydraulic_boundary_database_connection_failed_0_, validationProblems); CalculationServiceHelper.LogValidationEnd(); return !validationProblems.Any(); } /// /// Performs a calculation for the design water level. /// /// The design water level calculation to use. /// The path which points to the hydraulic boundary database file. /// The norm of the assessment section. /// The object which is used to build log messages. /// Thrown when /// is null. /// Thrown when /// /// contains invalid characters. /// The target probability or the calculated probability falls outside the [0.0, 1.0] range and is not . /// /// Thrown when: /// /// No settings database file could be found at the location of /// with the same name. /// Unable to open settings database file. /// Unable to read required data from database file. /// /// Thrown when an error occurs while performing the calculation. public void Calculate(IHydraulicBoundaryWrapperCalculation designWaterLevelCalculation, string hydraulicBoundaryDatabaseFilePath, double norm, ICalculationMessageProvider messageProvider) { if (designWaterLevelCalculation == null) { throw new ArgumentNullException(nameof(designWaterLevelCalculation)); } string hlcdDirectory = Path.GetDirectoryName(hydraulicBoundaryDatabaseFilePath); CalculationServiceHelper.LogCalculationBegin(); calculator = HydraRingCalculatorFactory.Instance.CreateDesignWaterLevelCalculator(hlcdDirectory); var exceptionThrown = false; try { PerformCalculation(designWaterLevelCalculation, hydraulicBoundaryDatabaseFilePath, norm, messageProvider); } catch (HydraRingCalculationException e) { if (!canceled) { string lastErrorContent = calculator.LastErrorFileContent; log.Error(!string.IsNullOrEmpty(lastErrorContent) ? messageProvider.GetCalculationFailedMessage(designWaterLevelCalculation.Name, lastErrorContent) : messageProvider.GetCalculationFailedMessage(designWaterLevelCalculation.Name, e.Message), e); exceptionThrown = true; throw; } } finally { string lastErrorFileContent = calculator.LastErrorFileContent; bool errorOccurred = CalculationServiceHelper.HasErrorOccurred(canceled, exceptionThrown, lastErrorFileContent); if (errorOccurred) { log.Error(messageProvider.GetCalculationFailedMessage(designWaterLevelCalculation.Name, lastErrorFileContent)); } log.InfoFormat(Resources.DesignWaterLevelCalculationService_Calculate_Calculation_temporary_directory_can_be_found_on_location_0, calculator.OutputDirectory); CalculationServiceHelper.LogCalculationEnd(); if (errorOccurred) { throw new HydraRingCalculationException(lastErrorFileContent); } } } /// /// Cancels the currently running design water level calculation. /// public void Cancel() { calculator?.Cancel(); canceled = true; } /// /// Performs a calculation for the design water level. /// /// The design water level calculation to use. /// The path which points to the hydraulic boundary database file. /// The norm of the assessment section. /// The object which is used to build log messages. /// Thrown when: /// /// No settings database file could be found at the location of /// with the same name. /// Unable to open settings database file. /// Unable to read required data from database file. /// /// Thrown when an error occurs while performing the calculation. private void PerformCalculation(IHydraulicBoundaryWrapperCalculation designWaterLevelCalculation, string hydraulicBoundaryDatabaseFilePath, double norm, ICalculationMessageProvider messageProvider) { AssessmentLevelCalculationInput calculationInput = CreateInput(designWaterLevelCalculation, norm, hydraulicBoundaryDatabaseFilePath); calculator.Calculate(calculationInput); if (canceled || !string.IsNullOrEmpty(calculator.LastErrorFileContent)) { return; } HydraulicBoundaryLocationOutput hydraulicBoundaryLocationOutput = CreateHydraulicBoundaryLocationOutput( messageProvider, designWaterLevelCalculation.Name, calculationInput.Beta, norm, calculator.Converged); if (designWaterLevelCalculation.CalculateIllustrationPoints) { SetIllustrationPointsResult(hydraulicBoundaryLocationOutput, calculator.IllustrationPointsResult); } designWaterLevelCalculation.Output = hydraulicBoundaryLocationOutput; } /// /// Sets a based on the information /// of to the . /// /// The /// for which to set the . /// The to base the /// to create on. private static void SetIllustrationPointsResult(HydraulicBoundaryLocationOutput hydraulicBoundaryLocationOutput, HydraRingGeneralResult hydraRingGeneralResult) { if (hydraRingGeneralResult != null) { GeneralResultSubMechanismIllustrationPoint generalResult = GeneralResultSubmechanismIllustrationPointConverter.CreateGeneralResultSubmechanismIllustrationPoint(hydraRingGeneralResult); hydraulicBoundaryLocationOutput.SetIllustrationPoints(generalResult); } } /// /// Creates the output of the calculation. /// /// The object which is used to build log messages. /// The name of the hydraulic boundary location. /// The target reliability for the calculation. /// The target probability for the calculation. /// The value indicating whether the calculation converged. /// A . /// Thrown when /// or the calculated probability falls outside the [0.0, 1.0] range and is not . private HydraulicBoundaryLocationOutput CreateHydraulicBoundaryLocationOutput( ICalculationMessageProvider messageProvider, string hydraulicBoundaryLocationName, double targetReliability, double targetProbability, bool? calculatorConverged) { double designWaterLevel = calculator.DesignWaterLevel; double reliability = calculator.ReliabilityIndex; double probability = StatisticsConverter.ReliabilityToProbability(reliability); CalculationConvergence converged = RingtoetsCommonDataCalculationService.GetCalculationConvergence(calculatorConverged); if (converged != CalculationConvergence.CalculatedConverged) { log.Warn(messageProvider.GetCalculatedNotConvergedMessage(hydraulicBoundaryLocationName)); } return new HydraulicBoundaryLocationOutput(designWaterLevel, targetProbability, targetReliability, probability, reliability, converged); } /// /// Creates the input for a design water level calculation. /// /// The /// to create the input from. /// The norm to use during the calculation. /// The file path to the hydraulic /// boundary database. /// An . /// Thrown when the /// contains invalid characters. /// Thrown when: /// /// No settings database file could be found at the location of /// with the same name. /// Unable to open settings database file. /// Unable to read required data from database file. /// /// private static AssessmentLevelCalculationInput CreateInput(IHydraulicBoundaryWrapperCalculation designWaterLevelCalculation, double norm, string hydraulicBoundaryDatabaseFilePath) { var assessmentLevelCalculationInput = new AssessmentLevelCalculationInput(1, designWaterLevelCalculation.Id, norm); HydraRingSettingsDatabaseHelper.AssignSettingsFromDatabase(assessmentLevelCalculationInput, hydraulicBoundaryDatabaseFilePath); return assessmentLevelCalculationInput; } private static string[] ValidateInput(string hydraulicBoundaryDatabaseFilePath) { var validationResult = new List(); string validationProblem = HydraulicDatabaseHelper.ValidatePathForCalculation(hydraulicBoundaryDatabaseFilePath); if (!string.IsNullOrEmpty(validationProblem)) { validationResult.Add(validationProblem); } return validationResult.ToArray(); } } }