Free cookie consent management tool by TermsFeed Policy Generator

source: trunk/sources/HeuristicLab.Problems.TestFunctions/3.3/SingleObjectiveTestFunctionProblem.cs @ 3751

Last change on this file since 3751 was 3751, checked in by abeham, 14 years ago

#893

  • fixed parameterization of strategy parameter operators in the test functions
File size: 18.5 KB
Line 
1#region License Information
2/* HeuristicLab
3 * Copyright (C) 2002-2010 Heuristic and Evolutionary Algorithms Laboratory (HEAL)
4 *
5 * This file is part of HeuristicLab.
6 *
7 * HeuristicLab is free software: you can redistribute it and/or modify
8 * it under the terms of the GNU General Public License as published by
9 * the Free Software Foundation, either version 3 of the License, or
10 * (at your option) any later version.
11 *
12 * HeuristicLab is distributed in the hope that it will be useful,
13 * but WITHOUT ANY WARRANTY; without even the implied warranty of
14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
15 * GNU General Public License for more details.
16 *
17 * You should have received a copy of the GNU General Public License
18 * along with HeuristicLab. If not, see <http://www.gnu.org/licenses/>.
19 */
20#endregion
21
22using System;
23using System.Collections.Generic;
24using System.Collections.ObjectModel;
25using System.Drawing;
26using System.Linq;
27using HeuristicLab.Common;
28using HeuristicLab.Core;
29using HeuristicLab.Data;
30using HeuristicLab.Encodings.RealVectorEncoding;
31using HeuristicLab.Optimization;
32using HeuristicLab.Parameters;
33using HeuristicLab.Persistence.Default.CompositeSerializers.Storable;
34using HeuristicLab.PluginInfrastructure;
35using HeuristicLab.Operators;
36using HeuristicLab.Problems.TestFunctions.Analyzers;
37
38namespace HeuristicLab.Problems.TestFunctions {
39  [Item("Single Objective Test Function", "Test function with real valued inputs and a single objective.")]
40  [StorableClass]
41  [Creatable("Problems")]
42  public sealed class SingleObjectiveTestFunctionProblem : ParameterizedNamedItem, ISingleObjectiveProblem {
43    [Storable]
44    private StdDevStrategyVectorCreator strategyVectorCreator;
45    [Storable]
46    private StdDevStrategyVectorCrossover strategyVectorCrossover;
47    [Storable]
48    private StdDevStrategyVectorManipulator strategyVectorManipulator;
49
50    public override Image ItemImage {
51      get { return HeuristicLab.Common.Resources.VS2008ImageLibrary.Type; }
52    }
53
54    #region Parameter Properties
55    public ValueParameter<BoolValue> MaximizationParameter {
56      get { return (ValueParameter<BoolValue>)Parameters["Maximization"]; }
57    }
58    IParameter ISingleObjectiveProblem.MaximizationParameter {
59      get { return MaximizationParameter; }
60    }
61    public ValueParameter<DoubleMatrix> BoundsParameter {
62      get { return (ValueParameter<DoubleMatrix>)Parameters["Bounds"]; }
63    }
64    public ValueParameter<IntValue> ProblemSizeParameter {
65      get { return (ValueParameter<IntValue>)Parameters["ProblemSize"]; }
66    }
67    public ValueParameter<IRealVectorCreator> SolutionCreatorParameter {
68      get { return (ValueParameter<IRealVectorCreator>)Parameters["SolutionCreator"]; }
69    }
70    IParameter IProblem.SolutionCreatorParameter {
71      get { return SolutionCreatorParameter; }
72    }
73    public ValueParameter<ISingleObjectiveTestFunctionProblemEvaluator> EvaluatorParameter {
74      get { return (ValueParameter<ISingleObjectiveTestFunctionProblemEvaluator>)Parameters["Evaluator"]; }
75    }
76    IParameter IProblem.EvaluatorParameter {
77      get { return EvaluatorParameter; }
78    }
79    public OptionalValueParameter<DoubleValue> BestKnownQualityParameter {
80      get { return (OptionalValueParameter<DoubleValue>)Parameters["BestKnownQuality"]; }
81    }
82    IParameter ISingleObjectiveProblem.BestKnownQualityParameter {
83      get { return BestKnownQualityParameter; }
84    }
85    #endregion
86
87    #region Properties
88    public BoolValue Maximization {
89      get { return MaximizationParameter.Value; }
90      set { MaximizationParameter.Value = value; }
91    }
92    public DoubleMatrix Bounds {
93      get { return BoundsParameter.Value; }
94      set { BoundsParameter.Value = value; }
95    }
96    public IntValue ProblemSize {
97      get { return ProblemSizeParameter.Value; }
98      set { ProblemSizeParameter.Value = value; }
99    }
100    public IRealVectorCreator SolutionCreator {
101      get { return SolutionCreatorParameter.Value; }
102      set { SolutionCreatorParameter.Value = value; }
103    }
104    ISolutionCreator IProblem.SolutionCreator {
105      get { return SolutionCreatorParameter.Value; }
106    }
107    public ISingleObjectiveTestFunctionProblemEvaluator Evaluator {
108      get { return EvaluatorParameter.Value; }
109      set { EvaluatorParameter.Value = value; }
110    }
111    ISingleObjectiveEvaluator ISingleObjectiveProblem.Evaluator {
112      get { return EvaluatorParameter.Value; }
113    }
114    IEvaluator IProblem.Evaluator {
115      get { return EvaluatorParameter.Value; }
116    }
117    public DoubleValue BestKnownQuality {
118      get { return BestKnownQualityParameter.Value; }
119      set { BestKnownQualityParameter.Value = value; }
120    }
121    private List<IOperator> operators;
122    public IEnumerable<IOperator> Operators {
123      get { return operators; }
124    }
125    private IEnumerable<IBestSingleObjectiveTestFunctionSolutionAnalyzer> BestSingleObjectiveTestFunctionSolutionAnalyzers {
126      get { return operators.OfType<IBestSingleObjectiveTestFunctionSolutionAnalyzer>(); }
127    }
128    #endregion
129
130    [StorableConstructor]
131    private SingleObjectiveTestFunctionProblem(bool deserializing) : base() { }
132    public SingleObjectiveTestFunctionProblem()
133      : base() {
134      UniformRandomRealVectorCreator creator = new UniformRandomRealVectorCreator();
135      AckleyEvaluator evaluator = new AckleyEvaluator();
136
137      Parameters.Add(new ValueParameter<BoolValue>("Maximization", "Set to false as most test functions are minimization problems.", new BoolValue(evaluator.Maximization)));
138      Parameters.Add(new ValueParameter<DoubleMatrix>("Bounds", "The lower and upper bounds in each dimension.", evaluator.Bounds));
139      Parameters.Add(new ValueParameter<IntValue>("ProblemSize", "The dimension of the problem.", new IntValue(2)));
140      Parameters.Add(new ValueParameter<IRealVectorCreator>("SolutionCreator", "The operator which should be used to create new test function solutions.", creator));
141      Parameters.Add(new ValueParameter<ISingleObjectiveTestFunctionProblemEvaluator>("Evaluator", "The operator which should be used to evaluate test function solutions.", evaluator));
142      Parameters.Add(new OptionalValueParameter<DoubleValue>("BestKnownQuality", "The quality of the best known solution of this test function.", new DoubleValue(evaluator.BestKnownQuality)));
143
144      strategyVectorCreator = new StdDevStrategyVectorCreator();
145      strategyVectorCreator.LengthParameter.ActualName = ProblemSizeParameter.Name;
146      strategyVectorCrossover = new StdDevStrategyVectorCrossover();
147      strategyVectorManipulator = new StdDevStrategyVectorManipulator();
148      strategyVectorManipulator.LearningRateParameter.Value = new DoubleValue(0.5);
149      strategyVectorManipulator.GeneralLearningRateParameter.Value = new DoubleValue(0.5);
150
151      creator.RealVectorParameter.ActualName = "Point";
152      ParameterizeSolutionCreator();
153      ParameterizeEvaluator();
154
155      Initialize();
156    }
157
158    public override IDeepCloneable Clone(Cloner cloner) {
159      SingleObjectiveTestFunctionProblem clone = (SingleObjectiveTestFunctionProblem)base.Clone(cloner);
160      clone.strategyVectorCreator = (StdDevStrategyVectorCreator)cloner.Clone(strategyVectorCreator);
161      clone.strategyVectorCrossover = (StdDevStrategyVectorCrossover)cloner.Clone(strategyVectorCrossover);
162      clone.strategyVectorManipulator = (StdDevStrategyVectorManipulator)cloner.Clone(strategyVectorManipulator);
163      clone.Initialize();
164      return clone;
165    }
166
167    #region Events
168    public event EventHandler SolutionCreatorChanged;
169    private void OnSolutionCreatorChanged() {
170      EventHandler handler = SolutionCreatorChanged;
171      if (handler != null) handler(this, EventArgs.Empty);
172    }
173    public event EventHandler EvaluatorChanged;
174    private void OnEvaluatorChanged() {
175      EventHandler handler = EvaluatorChanged;
176      if (handler != null) handler(this, EventArgs.Empty);
177    }
178    public event EventHandler OperatorsChanged;
179    private void OnOperatorsChanged() {
180      EventHandler handler = OperatorsChanged;
181      if (handler != null) handler(this, EventArgs.Empty);
182    }
183    public event EventHandler Reset;
184    private void OnReset() {
185      EventHandler handler = Reset;
186      if (handler != null) handler(this, EventArgs.Empty);
187    }
188
189    private void ProblemSizeParameter_ValueChanged(object sender, EventArgs e) {
190      ProblemSize.ValueChanged += new EventHandler(ProblemSize_ValueChanged);
191      ProblemSize_ValueChanged(null, EventArgs.Empty);
192    }
193    private void ProblemSize_ValueChanged(object sender, EventArgs e) {
194      if (ProblemSize.Value < 1) ProblemSize.Value = 1;
195      ParameterizeSolutionCreator();
196      strategyVectorManipulator.GeneralLearningRateParameter.Value = new DoubleValue(1.0 / Math.Sqrt(2 * ProblemSize.Value));
197      strategyVectorManipulator.LearningRateParameter.Value = new DoubleValue(1.0 / Math.Sqrt(2 * Math.Sqrt(ProblemSize.Value)));
198    }
199    private void SolutionCreatorParameter_ValueChanged(object sender, EventArgs e) {
200      ParameterizeSolutionCreator();
201      ParameterizeAnalyzers();
202      SolutionCreator_RealVectorParameter_ActualNameChanged(null, EventArgs.Empty);
203    }
204    private void SolutionCreator_RealVectorParameter_ActualNameChanged(object sender, EventArgs e) {
205      ParameterizeEvaluator();
206      ParameterizeOperators();
207      ParameterizeAnalyzers();
208    }
209    private void EvaluatorParameter_ValueChanged(object sender, EventArgs e) {
210      ParameterizeEvaluator();
211      UpdateMoveEvaluators();
212      ParameterizeAnalyzers();
213      Maximization.Value = Evaluator.Maximization;
214      BoundsParameter.Value = Evaluator.Bounds;
215      if (ProblemSize.Value < Evaluator.MinimumProblemSize)
216        ProblemSize.Value = Evaluator.MinimumProblemSize;
217      else if (ProblemSize.Value > Evaluator.MaximumProblemSize)
218        ProblemSize.Value = Evaluator.MaximumProblemSize;
219      BestKnownQuality = new DoubleValue(Evaluator.BestKnownQuality);
220      Evaluator_QualityParameter_ActualNameChanged(null, EventArgs.Empty);
221    }
222    private void Evaluator_QualityParameter_ActualNameChanged(object sender, EventArgs e) {
223      ParameterizeOperators();
224    }
225    private void BoundsParameter_ValueChanged(object sender, EventArgs e) {
226      Bounds.ToStringChanged += new EventHandler(Bounds_ToStringChanged);
227      Bounds_ToStringChanged(null, EventArgs.Empty);
228    }
229    private void Bounds_ToStringChanged(object sender, EventArgs e) {
230      if (Bounds.Columns != 2 || Bounds.Rows < 1)
231        Bounds = new DoubleMatrix(1, 2);
232      ParameterizeOperators();
233    }
234    private void Bounds_ItemChanged(object sender, EventArgs<int, int> e) {
235      if (e.Value2 == 0 && Bounds[e.Value, 1] <= Bounds[e.Value, 0])
236        Bounds[e.Value, 1] = Bounds[e.Value, 0] + 0.1;
237      if (e.Value2 == 1 && Bounds[e.Value, 0] >= Bounds[e.Value, 1])
238        Bounds[e.Value, 0] = Bounds[e.Value, 1] - 0.1;
239      ParameterizeOperators();
240    }
241    private void MoveGenerator_AdditiveMoveParameter_ActualNameChanged(object sender, EventArgs e) {
242      string name = ((ILookupParameter<AdditiveMove>)sender).ActualName;
243      foreach (IAdditiveRealVectorMoveOperator op in Operators.OfType<IAdditiveRealVectorMoveOperator>()) {
244        op.AdditiveMoveParameter.ActualName = name;
245      }
246    }
247    private void SphereEvaluator_Parameter_ValueChanged(object sender, EventArgs e) {
248      SphereEvaluator eval = (Evaluator as SphereEvaluator);
249      if (eval != null) {
250        foreach (ISphereMoveEvaluator op in Operators.OfType<ISphereMoveEvaluator>()) {
251          op.C = eval.C;
252          op.Alpha = eval.Alpha;
253        }
254      }
255    }
256    private void RastriginEvaluator_Parameter_ValueChanged(object sender, EventArgs e) {
257      RastriginEvaluator eval = (Evaluator as RastriginEvaluator);
258      if (eval != null) {
259        foreach (IRastriginMoveEvaluator op in Operators.OfType<IRastriginMoveEvaluator>()) {
260          op.A = eval.A;
261        }
262      }
263    }
264    private void strategyVectorCreator_BoundsParameter_ValueChanged(object sender, EventArgs e) {
265      strategyVectorManipulator.BoundsParameter.Value = strategyVectorCreator.BoundsParameter.Value;
266    }
267    private void strategyVectorCreator_StrategyParameterParameter_ActualNameChanged(object sender, EventArgs e) {
268      string name = strategyVectorCreator.StrategyParameterParameter.ActualName;
269      strategyVectorCrossover.ParentsParameter.ActualName = name;
270      strategyVectorCrossover.StrategyParameterParameter.ActualName = name;
271      strategyVectorManipulator.StrategyParameterParameter.ActualName = name;
272    }
273    #endregion
274
275    #region Helpers
276    [StorableHook(HookType.AfterDeserialization)]
277    private void Initialize() {
278      InitializeOperators();
279      ProblemSizeParameter.ValueChanged += new EventHandler(ProblemSizeParameter_ValueChanged);
280      ProblemSize.ValueChanged += new EventHandler(ProblemSize_ValueChanged);
281      BoundsParameter.ValueChanged += new EventHandler(BoundsParameter_ValueChanged);
282      Bounds.ToStringChanged += new EventHandler(Bounds_ToStringChanged);
283      Bounds.ItemChanged += new EventHandler<EventArgs<int, int>>(Bounds_ItemChanged);
284      SolutionCreatorParameter.ValueChanged += new EventHandler(SolutionCreatorParameter_ValueChanged);
285      SolutionCreator.RealVectorParameter.ActualNameChanged += new EventHandler(SolutionCreator_RealVectorParameter_ActualNameChanged);
286      EvaluatorParameter.ValueChanged += new EventHandler(EvaluatorParameter_ValueChanged);
287      Evaluator.QualityParameter.ActualNameChanged += new EventHandler(Evaluator_QualityParameter_ActualNameChanged);
288      strategyVectorCreator.BoundsParameter.ValueChanged += new EventHandler(strategyVectorCreator_BoundsParameter_ValueChanged);
289      strategyVectorCreator.StrategyParameterParameter.ActualNameChanged += new EventHandler(strategyVectorCreator_StrategyParameterParameter_ActualNameChanged);
290    }
291    private void ParameterizeAnalyzers() {
292      foreach (IBestSingleObjectiveTestFunctionSolutionAnalyzer analyzer in BestSingleObjectiveTestFunctionSolutionAnalyzers) {
293        analyzer.RealVectorParameter.ActualName = SolutionCreator.RealVectorParameter.ActualName;
294        analyzer.ResultsParameter.ActualName = "Results";
295      }
296    }
297    private void InitializeOperators() {
298      operators = new List<IOperator>();
299      operators.Add(new BestSingleObjectiveTestFunctionSolutionAnalyzer());
300      ParameterizeAnalyzers();
301      operators.AddRange(ApplicationManager.Manager.GetInstances<IRealVectorOperator>().Cast<IOperator>());
302      operators.Add(strategyVectorCreator);
303      operators.Add(strategyVectorCrossover);
304      operators.Add(strategyVectorManipulator);
305      UpdateMoveEvaluators();
306      ParameterizeOperators();
307      InitializeMoveGenerators();
308    }
309    private void InitializeMoveGenerators() {
310      foreach (IAdditiveRealVectorMoveOperator op in Operators.OfType<IAdditiveRealVectorMoveOperator>()) {
311        if (op is IMoveGenerator) {
312          op.AdditiveMoveParameter.ActualNameChanged += new EventHandler(MoveGenerator_AdditiveMoveParameter_ActualNameChanged);
313        }
314      }
315    }
316    private void UpdateMoveEvaluators() {
317      foreach (ISingleObjectiveTestFunctionMoveEvaluator op in Operators.OfType<ISingleObjectiveTestFunctionMoveEvaluator>().ToList())
318        operators.Remove(op);
319      foreach (ISingleObjectiveTestFunctionMoveEvaluator op in ApplicationManager.Manager.GetInstances<ISingleObjectiveTestFunctionMoveEvaluator>())
320        if (op.EvaluatorType == Evaluator.GetType()) {
321          operators.Add(op);
322          #region Synchronize evaluator specific parameters with the parameters of the corresponding move evaluators
323          if (op is ISphereMoveEvaluator) {
324            SphereEvaluator e = (Evaluator as SphereEvaluator);
325            e.AlphaParameter.ValueChanged += new EventHandler(SphereEvaluator_Parameter_ValueChanged);
326            e.CParameter.ValueChanged += new EventHandler(SphereEvaluator_Parameter_ValueChanged);
327            ISphereMoveEvaluator em = (op as ISphereMoveEvaluator);
328            em.C = e.C;
329            em.Alpha = e.Alpha;
330          } else if (op is IRastriginMoveEvaluator) {
331            RastriginEvaluator e = (Evaluator as RastriginEvaluator);
332            e.AParameter.ValueChanged += new EventHandler(RastriginEvaluator_Parameter_ValueChanged);
333            IRastriginMoveEvaluator em = (op as IRastriginMoveEvaluator);
334            em.A = e.A;
335          }
336          #endregion
337        }
338      ParameterizeOperators();
339      OnOperatorsChanged();
340    }
341    private void ParameterizeSolutionCreator() {
342      SolutionCreator.LengthParameter.Value = new IntValue(ProblemSize.Value);
343    }
344    private void ParameterizeEvaluator() {
345      Evaluator.PointParameter.ActualName = SolutionCreator.RealVectorParameter.ActualName;
346    }
347    private void ParameterizeOperators() {
348      foreach (IRealVectorCrossover op in Operators.OfType<IRealVectorCrossover>()) {
349        op.ParentsParameter.ActualName = SolutionCreator.RealVectorParameter.ActualName;
350        op.ChildParameter.ActualName = SolutionCreator.RealVectorParameter.ActualName;
351        op.BoundsParameter.ActualName = BoundsParameter.Name;
352      }
353      foreach (IRealVectorManipulator op in Operators.OfType<IRealVectorManipulator>()) {
354        op.RealVectorParameter.ActualName = SolutionCreator.RealVectorParameter.ActualName;
355        op.BoundsParameter.ActualName = BoundsParameter.Name;
356      }
357      foreach (IRealVectorMoveOperator op in Operators.OfType<IRealVectorMoveOperator>()) {
358        op.RealVectorParameter.ActualName = SolutionCreator.RealVectorParameter.ActualName;
359      }
360      foreach (IRealVectorMoveGenerator op in Operators.OfType<IRealVectorMoveGenerator>()) {
361        op.BoundsParameter.ActualName = BoundsParameter.Name;
362      }
363      foreach (ISingleObjectiveTestFunctionAdditiveMoveEvaluator op in Operators.OfType<ISingleObjectiveTestFunctionAdditiveMoveEvaluator>()) {
364        op.QualityParameter.ActualName = Evaluator.QualityParameter.ActualName;
365        op.RealVectorParameter.ActualName = SolutionCreator.RealVectorParameter.ActualName;
366      }
367      foreach (IRealVectorPSOEncoder op in Operators.OfType<IRealVectorPSOEncoder>()) {
368        ((ILookupParameter)op.OriginalRealVectorParameter).ActualName = SolutionCreator.RealVectorParameter.ActualName;
369        op.BoundsParameter.Value = (DoubleMatrix)BoundsParameter.Value.Clone();
370        op.BoundsParameter.ActualName = "ParticleBounds";
371      }
372    }
373    #endregion
374  }
375}
Note: See TracBrowser for help on using the repository browser.