LZ2.h
Go to the documentation of this file.
1 //===========================================================================
2 /*!
3  *
4  *
5  * \brief Multi-objective optimization benchmark function LZ2.
6  *
7  * The function is described in
8  *
9  * H. Li and Q. Zhang.
10  * Multiobjective Optimization Problems with Complicated Pareto Sets, MOEA/D and NSGA-II,
11  * IEEE Trans on Evolutionary Computation, 2(12):284-302, April 2009.
12  *
13  *
14  *
15  * \author -
16  * \date -
17  *
18  *
19  * \par Copyright 1995-2017 Shark Development Team
20  *
21  * <BR><HR>
22  * This file is part of Shark.
23  * <http://shark-ml.org/>
24  *
25  * Shark is free software: you can redistribute it and/or modify
26  * it under the terms of the GNU Lesser General Public License as published
27  * by the Free Software Foundation, either version 3 of the License, or
28  * (at your option) any later version.
29  *
30  * Shark is distributed in the hope that it will be useful,
31  * but WITHOUT ANY WARRANTY; without even the implied warranty of
32  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
33  * GNU Lesser General Public License for more details.
34  *
35  * You should have received a copy of the GNU Lesser General Public License
36  * along with Shark. If not, see <http://www.gnu.org/licenses/>.
37  *
38  */
39 //===========================================================================
40 #ifndef SHARK_OBJECTIVEFUNCTIONS_BENCHMARK_LZ2_H
41 #define SHARK_OBJECTIVEFUNCTIONS_BENCHMARK_LZ2_H
42 
45 
46 namespace shark {
47 /*! \brief Multi-objective optimization benchmark function LZ2.
48 *
49 * The function is described in
50 *
51 * H. Li and Q. Zhang.
52 * Multiobjective Optimization Problems with Complicated Pareto Sets, MOEA/D and NSGA-II,
53 * IEEE Trans on Evolutionary Computation, 2(12):284-302, April 2009.
54 */
55 struct LZ2 : public MultiObjectiveFunction
56 {
57 
58  LZ2(std::size_t numVariables = 0) : m_handler(SearchPointType(numVariables,-1),SearchPointType(numVariables,1) ){
59  announceConstraintHandler(&m_handler);
60  }
61 
62  /// \brief From INameable: return the class name.
63  std::string name() const
64  { return "LZ2"; }
65 
66  std::size_t numberOfObjectives()const{
67  return 2;
68  }
69 
70  std::size_t numberOfVariables()const{
71  return m_handler.dimensions();
72  }
73 
75  return true;
76  }
77 
78  /// \brief Adjusts the number of variables if the function is scalable.
79  /// \param [in] numberOfVariables The new dimension.
81  SearchPointType lb(numberOfVariables,-1);
82  SearchPointType ub(numberOfVariables, 1);
83  lb(0) = 0;
84  ub(0) = 1;
85  m_handler.setBounds(lb, ub);
86  }
87 
88  ResultType eval( const SearchPointType & x ) const {
90 
91  ResultType value( 2, 0 );
92 
93  std::size_t counter1 = 0, counter2 = 0;
94  for( std::size_t i = 1; i < x.size(); i++ ) {
95  if( i % 2 == 0 ) {
96  counter2++;
97  value[1] += sqr( x(i) - ::sin( 6 * M_PI * x( 0 ) + i*M_PI/x.size() ) );
98  } else {
99  counter1++;
100  value[0] += sqr( x(i) - ::sin( 6 * M_PI * x( 0 ) + i*M_PI/x.size() ) );
101  }
102  }
103 
104  value[0] *= 2./counter1;
105  value[0] += x( 0 );
106 
107  value[1] *= 2./counter2;
108  value[1] += 1 - ::sqrt( x( 0 ) );
109 
110  return value;
111  }
112 private:
114 };
115 
116 }
117 #endif