PaGMO  1.1.5
kur.cpp
1 /*****************************************************************************
2  * Copyright (C) 2004-2015 The PaGMO development team, *
3  * Advanced Concepts Team (ACT), European Space Agency (ESA) *
4  * *
5  * https://github.com/esa/pagmo *
6  * *
7  * act@esa.int *
8  * *
9  * This program is free software; you can redistribute it and/or modify *
10  * it under the terms of the GNU General Public License as published by *
11  * the Free Software Foundation; either version 2 of the License, or *
12  * (at your option) any later version. *
13  * *
14  * This program is distributed in the hope that it will be useful, *
15  * but WITHOUT ANY WARRANTY; without even the implied warranty of *
16  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
17  * GNU General Public License for more details. *
18  * *
19  * You should have received a copy of the GNU General Public License *
20  * along with this program; if not, write to the *
21  * Free Software Foundation, Inc., *
22  * 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. *
23  *****************************************************************************/
24 
25 #include <cmath>
26 
27 #include "../exceptions.h"
28 #include "../types.h"
29 #include "base.h"
30 #include "kur.h"
31 
32 namespace pagmo { namespace problem {
33 
41 kur::kur(size_type dim):base(dim,0,2)
42 {
43  // Set bounds.
44  set_lb(-5.0);
45  set_ub(5.0);
46 }
47 
50 {
51  return base_ptr(new kur(*this));
52 }
53 
56 {
57  unsigned int n = get_dimension();
58  pagmo_assert(f.size() == 2);
59  pagmo_assert(x.size() == n);
60 
61  f[0] = 0;
62  f[1] = 0;
63 
64  for(unsigned int i = 0; i < n-1; ++i) {
65  f[0] += -10 * exp( -0.2 * sqrt(x[i]*x[i] + x[i+1]*x[i+1]));
66  f[1] += pow(std::abs(x[i]),0.8) + 5 * sin(pow(x[i],3));
67  }
68  f[1] += pow(std::abs(x[n-1]),0.8) + 5 * sin(pow(x[n-1],3));
69 }
70 
71 std::string kur::get_name() const
72 {
73  return "Kursawe's study";
74 }
75 }}
76 
77 BOOST_CLASS_EXPORT_IMPLEMENT(pagmo::problem::kur)
Root PaGMO namespace.
boost::shared_ptr< base > base_ptr
Alias for shared pointer to base problem.
Definition: problem/base.h:62
std::vector< double > decision_vector
Decision vector type.
Definition: types.h:40
Kursawe's study.
Definition: kur.h:52
kur(size_type=10)
Definition: kur.cpp:41
Base problem class.
Definition: problem/base.h:148
size_type get_dimension() const
Return global dimension.
void set_lb(const decision_vector &)
Set lower bounds from pagmo::decision_vector.
std::string get_name() const
Get problem's name.
Definition: kur.cpp:71
void set_ub(const decision_vector &)
Set upper bounds from pagmo::decision_vector.
base_ptr clone() const
Clone method.
Definition: kur.cpp:49
std::vector< double > fitness_vector
Fitness vector type.
Definition: types.h:42
void objfun_impl(fitness_vector &, const decision_vector &) const
Implementation of the objective function.
Definition: kur.cpp:55
decision_vector::size_type size_type
Problem's size type: the same as pagmo::decision_vector's size type.
Definition: problem/base.h:160