1  #region License Information


2  /* HeuristicLab


3  * Copyright (C) 20022019 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  using System;


22  using System.Collections.Generic;


23  using HeuristicLab.Common;


24  using HeuristicLab.Core;


25  using HeuristicLab.Encodings.RealVectorEncoding;


26  using HEAL.Fossil;


27 


28  namespace HeuristicLab.Problems.TestFunctions.MultiObjective {


29  [Item("IHR3", "Testfunction as defined as IHR3 in \"Igel, C., Hansen, N., & Roth, S. (2007). Covariance matrix adaptation for multiobjective optimization. Evolutionary computation, 15(1), 128.\" [24.06.16]")]


30  [StorableType("316D5351762D4883ACEF06F4EAFA73AE")]


31  public class IHR3 : IHR {


32  protected override IEnumerable<double[]> GetOptimalParetoFront(int objectives) {


33  List<double[]> res = new List<double[]>();


34  for (int i = 0; i <= 500; i++) {


35  RealVector r = new RealVector(objectives);


36  r[0] = 1 / 500.0 * i;


37 


38  res.Add(this.Evaluate(r, objectives));


39  }


40  return res;


41  }


42 


43  protected override double GetBestKnownHypervolume(int objectives) {


44  return Hypervolume.Calculate(GetOptimalParetoFront(objectives), GetReferencePoint(objectives), GetMaximization(objectives));


45  }


46 


47  [StorableConstructor]


48  protected IHR3(StorableConstructorFlag _) : base(_) { }


49  protected IHR3(IHR3 original, Cloner cloner) : base(original, cloner) { }


50  public IHR3() : base() {


51  }


52 


53  public override IDeepCloneable Clone(Cloner cloner) {


54  return new IHR3(this, cloner);


55  }


56 


57  protected override double F1(RealVector y) {


58  return Math.Abs(y[0]);


59  }


60 


61  protected override double F2(RealVector y) {


62  var g = G(y);


63  var d = H(y[0], y) / g;


64  return g * HF(1  Math.Sqrt(d)  d * Math.Sin(10 * Math.PI * y[0]), y);


65  }


66 


67  protected override double G(RealVector y) {


68  double sum = 0.0;


69  for (int i = 1; i < y.Length; i++) {


70  sum += HG(y[i]);


71  }


72  return 1 + 9 * sum / (y.Length  1);


73  }


74  }


75  }

