#region License Information /* HeuristicLab * Copyright (C) 2002-2019 Heuristic and Evolutionary Algorithms Laboratory (HEAL) * * This file is part of HeuristicLab. * * HeuristicLab 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. * * HeuristicLab 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 HeuristicLab. If not, see . */ #endregion using System; using System.Linq; using HeuristicLab.Analysis; using HeuristicLab.Common; using HeuristicLab.Core; using HeuristicLab.Data; using HeuristicLab.Operators; using HeuristicLab.Optimization; using HeuristicLab.Optimization.Operators; using HeuristicLab.Parameters; using HEAL.Attic; using HeuristicLab.PluginInfrastructure; using HeuristicLab.Random; namespace HeuristicLab.Algorithms.VariableNeighborhoodSearch { [Item("Variable Neighborhood Search (VNS)", "A variable neighborhood search algorithm based on the description in Mladenovic, N. and Hansen, P. (1997). Variable neighborhood search. Computers & Operations Research Volume 24, Issue 11, pp. 1097-1100.")] [Creatable(CreatableAttribute.Categories.SingleSolutionAlgorithms, Priority = 110)] [StorableType("E5FA1D1E-2611-4059-A836-4E89BB24342A")] public sealed class VariableNeighborhoodSearch : HeuristicOptimizationEngineAlgorithm, IStorableContent { public string Filename { get; set; } #region Problem Properties public override Type ProblemType { get { return typeof(ISingleObjectiveHeuristicOptimizationProblem); } } public new ISingleObjectiveHeuristicOptimizationProblem Problem { get { return (ISingleObjectiveHeuristicOptimizationProblem)base.Problem; } set { base.Problem = value; } } #endregion #region Parameter Properties private FixedValueParameter SeedParameter { get { return (FixedValueParameter)Parameters["Seed"]; } } private FixedValueParameter SetSeedRandomlyParameter { get { return (FixedValueParameter)Parameters["SetSeedRandomly"]; } } public IConstrainedValueParameter LocalImprovementParameter { get { return (IConstrainedValueParameter)Parameters["LocalImprovement"]; } } public IConstrainedValueParameter ShakingOperatorParameter { get { return (IConstrainedValueParameter)Parameters["ShakingOperator"]; } } private FixedValueParameter MaximumIterationsParameter { get { return (FixedValueParameter)Parameters["MaximumIterations"]; } } private ValueParameter AnalyzerParameter { get { return (ValueParameter)Parameters["Analyzer"]; } } public FixedValueParameter LocalImprovementMaximumIterationsParameter { get { return (FixedValueParameter)Parameters["LocalImprovementMaximumIterations"]; } } #endregion #region Properties private RandomCreator RandomCreator { get { return (RandomCreator)OperatorGraph.InitialOperator; } } public MultiAnalyzer Analyzer { get { return AnalyzerParameter.Value; } set { AnalyzerParameter.Value = value; } } private SolutionsCreator SolutionsCreator { get { return (SolutionsCreator)RandomCreator.Successor; } } private VariableNeighborhoodSearchMainLoop MainLoop { get { return FindMainLoop(SolutionsCreator.Successor); } } public int Seed { get { return SeedParameter.Value.Value; } set { SeedParameter.Value.Value = value; } } public bool SetSeedRandomly { get { return SetSeedRandomlyParameter.Value.Value; } set { SetSeedRandomlyParameter.Value.Value = value; } } public ILocalImprovementOperator LocalImprovement { get { return LocalImprovementParameter.Value; } set { LocalImprovementParameter.Value = value; } } public IMultiNeighborhoodShakingOperator ShakingOperator { get { return ShakingOperatorParameter.Value; } set { ShakingOperatorParameter.Value = value; } } public int MaximumIterations { get { return MaximumIterationsParameter.Value.Value; } set { MaximumIterationsParameter.Value.Value = value; } } public int LocalImprovementMaximumIterations { get { return LocalImprovementMaximumIterationsParameter.Value.Value; } set { LocalImprovementMaximumIterationsParameter.Value.Value = value; } } #endregion [Storable] private BestAverageWorstQualityAnalyzer qualityAnalyzer; [StorableConstructor] private VariableNeighborhoodSearch(StorableConstructorFlag _) : base(_) { } private VariableNeighborhoodSearch(VariableNeighborhoodSearch original, Cloner cloner) : base(original, cloner) { qualityAnalyzer = cloner.Clone(original.qualityAnalyzer); RegisterEventHandlers(); } public VariableNeighborhoodSearch() : base() { Parameters.Add(new FixedValueParameter("Seed", "The random seed used to initialize the new pseudo random number generator.", new IntValue(0))); Parameters.Add(new FixedValueParameter("SetSeedRandomly", "True if the random seed should be set to a random value, otherwise false.", new BoolValue(true))); Parameters.Add(new ConstrainedValueParameter("LocalImprovement", "The local improvement operation")); Parameters.Add(new ConstrainedValueParameter("ShakingOperator", "The operator that performs the shaking of solutions.")); Parameters.Add(new FixedValueParameter("MaximumIterations", "The maximum number of iterations which should be processed.", new IntValue(50))); Parameters.Add(new FixedValueParameter("LocalImprovementMaximumIterations", "The maximum number of iterations which should be performed in the local improvement phase.", new IntValue(50))); Parameters.Add(new ValueParameter("Analyzer", "The operator used to analyze the solution and moves.", new MultiAnalyzer())); RandomCreator randomCreator = new RandomCreator(); SolutionsCreator solutionsCreator = new SolutionsCreator(); VariableCreator variableCreator = new VariableCreator(); ResultsCollector resultsCollector = new ResultsCollector(); VariableNeighborhoodSearchMainLoop mainLoop = new VariableNeighborhoodSearchMainLoop(); OperatorGraph.InitialOperator = randomCreator; randomCreator.RandomParameter.ActualName = "Random"; randomCreator.SeedParameter.ActualName = SeedParameter.Name; randomCreator.SeedParameter.Value = null; randomCreator.SetSeedRandomlyParameter.ActualName = SetSeedRandomlyParameter.Name; randomCreator.SetSeedRandomlyParameter.Value = null; randomCreator.Successor = solutionsCreator; solutionsCreator.NumberOfSolutions = new IntValue(1); solutionsCreator.Successor = variableCreator; variableCreator.Name = "Initialize Evaluated Solutions"; variableCreator.CollectedValues.Add(new ValueParameter("Iterations", new IntValue(0))); variableCreator.CollectedValues.Add(new ValueParameter("EvaluatedSolutions", new IntValue(0))); variableCreator.CollectedValues.Add(new ValueParameter("CurrentNeighborhoodIndex", new IntValue(0))); variableCreator.CollectedValues.Add(new ValueParameter("NeighborhoodCount", new IntValue(0))); variableCreator.Successor = resultsCollector; resultsCollector.CollectedValues.Add(new LookupParameter("Evaluated Solutions", null, "EvaluatedSolutions")); resultsCollector.CollectedValues.Add(new LookupParameter("Iterations")); resultsCollector.ResultsParameter.ActualName = "Results"; resultsCollector.Successor = mainLoop; mainLoop.IterationsParameter.ActualName = "Iterations"; mainLoop.CurrentNeighborhoodIndexParameter.ActualName = "CurrentNeighborhoodIndex"; mainLoop.NeighborhoodCountParameter.ActualName = "NeighborhoodCount"; mainLoop.LocalImprovementParameter.ActualName = LocalImprovementParameter.Name; mainLoop.ShakingOperatorParameter.ActualName = ShakingOperatorParameter.Name; mainLoop.MaximumIterationsParameter.ActualName = MaximumIterationsParameter.Name; mainLoop.RandomParameter.ActualName = RandomCreator.RandomParameter.ActualName; mainLoop.ResultsParameter.ActualName = "Results"; mainLoop.AnalyzerParameter.ActualName = AnalyzerParameter.Name; mainLoop.EvaluatedSolutionsParameter.ActualName = "EvaluatedSolutions"; InitializeLocalImprovementOperators(); qualityAnalyzer = new BestAverageWorstQualityAnalyzer(); ParameterizeAnalyzers(); UpdateAnalyzers(); RegisterEventHandlers(); } public override IDeepCloneable Clone(Cloner cloner) { return new VariableNeighborhoodSearch(this, cloner); } [StorableHook(HookType.AfterDeserialization)] private void AfterDeserialization() { RegisterEventHandlers(); } public override void Prepare() { if (Problem != null) base.Prepare(); } private void RegisterEventHandlers() { LocalImprovementParameter.ValueChanged += new EventHandler(LocalImprovementParameter_ValueChanged); if (Problem != null) { Problem.Evaluator.QualityParameter.ActualNameChanged += new EventHandler(Evaluator_QualityParameter_ActualNameChanged); } } #region Events protected override void OnProblemChanged() { InitializeLocalImprovementOperators(); UpdateShakingOperators(); UpdateAnalyzers(); ParameterizeStochasticOperator(Problem.SolutionCreator); ParameterizeStochasticOperator(Problem.Evaluator); foreach (IOperator op in Problem.Operators.OfType()) ParameterizeStochasticOperator(op); ParameterizeSolutionsCreator(); ParameterizeMainLoop(); ParameterizeAnalyzers(); ParameterizeIterationBasedOperators(); ParameterizeLocalImprovementOperators(); Problem.Evaluator.QualityParameter.ActualNameChanged += new EventHandler(Evaluator_QualityParameter_ActualNameChanged); base.OnProblemChanged(); } protected override void Problem_SolutionCreatorChanged(object sender, EventArgs e) { ParameterizeStochasticOperator(Problem.SolutionCreator); ParameterizeSolutionsCreator(); base.Problem_SolutionCreatorChanged(sender, e); } protected override void Problem_EvaluatorChanged(object sender, EventArgs e) { ParameterizeStochasticOperator(Problem.Evaluator); ParameterizeSolutionsCreator(); ParameterizeMainLoop(); ParameterizeAnalyzers(); Problem.Evaluator.QualityParameter.ActualNameChanged += new EventHandler(Evaluator_QualityParameter_ActualNameChanged); base.Problem_EvaluatorChanged(sender, e); } protected override void Problem_OperatorsChanged(object sender, EventArgs e) { UpdateShakingOperators(); UpdateAnalyzers(); foreach (IOperator op in Problem.Operators.OfType()) ParameterizeStochasticOperator(op); ParameterizeIterationBasedOperators(); base.Problem_OperatorsChanged(sender, e); } private void Evaluator_QualityParameter_ActualNameChanged(object sender, EventArgs e) { ParameterizeMainLoop(); ParameterizeAnalyzers(); } private void LocalImprovementParameter_ValueChanged(object sender, EventArgs e) { ParameterizeLocalImprovementOperators(); } #endregion #region Helpers private void ParameterizeSolutionsCreator() { SolutionsCreator.EvaluatorParameter.ActualName = Problem.EvaluatorParameter.Name; SolutionsCreator.SolutionCreatorParameter.ActualName = Problem.SolutionCreatorParameter.Name; } private void ParameterizeStochasticOperator(IOperator op) { if (op is IStochasticOperator) ((IStochasticOperator)op).RandomParameter.ActualName = RandomCreator.RandomParameter.ActualName; } private void ParameterizeMainLoop() { MainLoop.EvaluatorParameter.ActualName = Problem.EvaluatorParameter.Name; MainLoop.MaximizationParameter.ActualName = Problem.MaximizationParameter.Name; MainLoop.QualityParameter.ActualName = Problem.Evaluator.QualityParameter.ActualName; } private void ParameterizeAnalyzers() { qualityAnalyzer.ResultsParameter.ActualName = "Results"; if (Problem != null) { qualityAnalyzer.MaximizationParameter.ActualName = Problem.MaximizationParameter.Name; qualityAnalyzer.QualityParameter.ActualName = Problem.Evaluator.QualityParameter.ActualName; qualityAnalyzer.QualityParameter.Depth = 1; qualityAnalyzer.BestKnownQualityParameter.ActualName = Problem.BestKnownQualityParameter.Name; } } private void ParameterizeIterationBasedOperators() { if (Problem != null) { foreach (IIterationBasedOperator op in Problem.Operators.OfType()) { op.IterationsParameter.ActualName = "Iterations"; op.MaximumIterationsParameter.ActualName = MaximumIterationsParameter.Name; } } } private void ParameterizeLocalImprovementOperators() { foreach (ILocalImprovementOperator op in LocalImprovementParameter.ValidValues) { op.MaximumIterationsParameter.Value = null; op.MaximumIterationsParameter.ActualName = LocalImprovementMaximumIterationsParameter.Name; var algOp = op as ILocalImprovementAlgorithmOperator; if (algOp != null && algOp != LocalImprovementParameter.Value) algOp.Problem = null; } if (LocalImprovementParameter.Value is ILocalImprovementAlgorithmOperator) ((ILocalImprovementAlgorithmOperator)LocalImprovementParameter.Value).Problem = Problem; } private void InitializeLocalImprovementOperators() { LocalImprovementParameter.ValidValues.Clear(); if (Problem != null) { // Regular ILocalImprovementOperators queried from Problem foreach (var op in Problem.Operators.OfType().Where(x => !(x is ILocalImprovementAlgorithmOperator))) { LocalImprovementParameter.ValidValues.Add(op); } // ILocalImprovementAlgorithmOperators queried from ApplicationManager var algOps = ApplicationManager.Manager.GetInstances() .Where(x => x.ProblemType.IsInstanceOfType(Problem)); foreach (var op in algOps) { if (LocalImprovementParameter.ValidValues.All(x => x.GetType() != op.GetType())) LocalImprovementParameter.ValidValues.Add(op); } } } private void UpdateShakingOperators() { IMultiNeighborhoodShakingOperator oldShakingOperator = ShakingOperator; IMultiNeighborhoodShakingOperator defaultShakingOperator = Problem.Operators.OfType().FirstOrDefault(); ShakingOperatorParameter.ValidValues.Clear(); foreach (IMultiNeighborhoodShakingOperator op in Problem.Operators.OfType().OrderBy(op => op.Name)) { ShakingOperatorParameter.ValidValues.Add(op); op.CurrentNeighborhoodIndexParameter.ActualName = "CurrentNeighborhoodIndex"; op.NeighborhoodCountParameter.ActualName = "NeighborhoodCount"; } if (oldShakingOperator != null) { IMultiNeighborhoodShakingOperator shakingOperator = ShakingOperatorParameter.ValidValues.FirstOrDefault(x => x.GetType() == oldShakingOperator.GetType()); if (shakingOperator != null) ShakingOperator = shakingOperator; else oldShakingOperator = null; } if (oldShakingOperator == null && defaultShakingOperator != null) { ShakingOperator = defaultShakingOperator; } } private void UpdateAnalyzers() { Analyzer.Operators.Clear(); if (Problem != null) { foreach (IAnalyzer analyzer in Problem.Operators.OfType()) { foreach (IScopeTreeLookupParameter param in analyzer.Parameters.OfType()) param.Depth = 1; Analyzer.Operators.Add(analyzer, analyzer.EnabledByDefault); } } Analyzer.Operators.Add(qualityAnalyzer, qualityAnalyzer.EnabledByDefault); } private VariableNeighborhoodSearchMainLoop FindMainLoop(IOperator start) { IOperator mainLoop = start; while (mainLoop != null && !(mainLoop is VariableNeighborhoodSearchMainLoop)) mainLoop = ((SingleSuccessorOperator)mainLoop).Successor; if (mainLoop == null) return null; else return (VariableNeighborhoodSearchMainLoop)mainLoop; } #endregion } }