#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.Collections.Generic; using System.Linq; using HeuristicLab.Common; using HeuristicLab.Core; using HeuristicLab.Data; using HeuristicLab.Optimization; using HeuristicLab.Parameters; using HEAL.Attic; namespace HeuristicLab.Problems.TestFunctions.MultiObjective { [StorableType("13D363E4-76FF-4A5A-9B2C-767D9E880E4B")] [Item("HypervolumeAnalyzer", "Computes the enclosed Hypervolume between the current front and a given reference Point")] public class HypervolumeAnalyzer : MOTFAnalyzer { public ILookupParameter ReferencePointParameter { get { return (ILookupParameter)Parameters["ReferencePoint"]; } } public IResultParameter HypervolumeResultParameter { get { return (IResultParameter)Parameters["Hypervolume"]; } } public IResultParameter BestKnownHypervolumeResultParameter { get { return (IResultParameter)Parameters["Best known hypervolume"]; } } public IResultParameter HypervolumeDistanceResultParameter { get { return (IResultParameter)Parameters["Absolute Distance to BestKnownHypervolume"]; } } [StorableConstructor] protected HypervolumeAnalyzer(StorableConstructorFlag _) : base(_) { } protected HypervolumeAnalyzer(HypervolumeAnalyzer original, Cloner cloner) : base(original, cloner) { } public override IDeepCloneable Clone(Cloner cloner) { return new HypervolumeAnalyzer(this, cloner); } public HypervolumeAnalyzer() { Parameters.Add(new LookupParameter("ReferencePoint", "The reference point for hypervolume calculation")); Parameters.Add(new ResultParameter("Hypervolume", "The hypervolume of the current generation")); Parameters.Add(new ResultParameter("Best known hypervolume", "The optimal hypervolume")); Parameters.Add(new ResultParameter("Absolute Distance to BestKnownHypervolume", "The difference between the best known and the current hypervolume")); HypervolumeResultParameter.DefaultValue = new DoubleValue(0); BestKnownHypervolumeResultParameter.DefaultValue = new DoubleValue(0); HypervolumeDistanceResultParameter.DefaultValue = new DoubleValue(0); } public override IOperation Apply() { var qualities = QualitiesParameter.ActualValue; var testFunction = TestFunctionParameter.ActualValue; int objectives = qualities[0].Length; var referencePoint = ReferencePointParameter.ActualValue; double best = BestKnownHypervolumeResultParameter.ActualValue.Value; if (referencePoint.SequenceEqual(testFunction.ReferencePoint(objectives))) { best = Math.Max(best, testFunction.OptimalHypervolume(objectives)); } IEnumerable front = NonDominatedSelect.SelectNonDominatedVectors(qualities.Select(q => q.ToArray()), testFunction.Maximization(objectives), true); double hv = Hypervolume.Calculate(front, referencePoint.ToArray(), testFunction.Maximization(objectives)); if (hv > best) { best = hv; } HypervolumeResultParameter.ActualValue.Value = hv; BestKnownHypervolumeResultParameter.ActualValue.Value = best; HypervolumeDistanceResultParameter.ActualValue.Value = best - hv; return base.Apply(); } } }