PaGMO  1.1.5
pol.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 <boost/math/constants/constants.hpp>
26 #include <cmath>
27 
28 #include "../exceptions.h"
29 #include "../types.h"
30 #include "base.h"
31 #include "pol.h"
32 
33 namespace pagmo { namespace problem {
34 
40 pol::pol():base(2,0,2)
41 {
42  // Set bounds.
43  set_lb(-std::atan(1.0f) * 4.0f); //-PI
44  set_ub(std::atan(1.0f) * 4.0f); //PI
45 }
46 
49 {
50  return base_ptr(new pol(*this));
51 }
52 
55 {
56  pagmo_assert(f.size() == 2);
57  pagmo_assert(x.size() == 2);
58 
59  f[0] = 0;
60  f[1] = 0;
61 
62  double a1 = 0.5 * sin(1.0) - 2 * cos(1.0) + sin(2.0) - 1.5 * cos(2.0);
63 
64  double a2 = 1.5 * sin(1.0) - cos(1.0) + 2 * sin(2.0) - 0.5 * cos(2.0);
65 
66  double b1 = 0.5 * sin(x[0]) - 2 * cos(x[0]) + sin(x[1]) - 1.5 * cos(x[1]);
67 
68  double b2 = 1.5 * sin(x[0]) - cos(x[0]) + 2 * sin(x[1]) - 0.5 * cos(x[1]);
69 
70  f[0] = 1 + (a1-b1)*(a1-b1) + (a2-b2)*(a2-b2);
71  f[1] = (x[0]+3)*(x[0]+3) + (x[1]+1)*(x[1]+1);
72 }
73 
74 std::string pol::get_name() const
75 {
76  return "Poloni's study";
77 }
78 
79 
80 }} //namespaces
81 
82 BOOST_CLASS_EXPORT_IMPLEMENT(pagmo::problem::pol)
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
std::string get_name() const
Get problem's name.
Definition: pol.cpp:74
void objfun_impl(fitness_vector &, const decision_vector &) const
Implementation of the objective function.
Definition: pol.cpp:54
Base problem class.
Definition: problem/base.h:148
void set_lb(const decision_vector &)
Set lower bounds from pagmo::decision_vector.
base_ptr clone() const
Clone method.
Definition: pol.cpp:48
void set_ub(const decision_vector &)
Set upper bounds from pagmo::decision_vector.
std::vector< double > fitness_vector
Fitness vector type.
Definition: types.h:42
Poloni's study.
Definition: pol.h:63