-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathIndividual.cpp
93 lines (69 loc) · 2.47 KB
/
Individual.cpp
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
#include "Individual.h"
using namespace std;
// Estandar para numeros flotantes: IEEE - 754
//Los Qubits de cada cromosoma son interpretados como numeros flotantes usando el estandar
float Individual::interpret_chromosome(size_t exponent_size, size_t mantissa_size) {
float exponent = 0.f;
float mantissa = 1.f;
int sign = chromosome[m - 1].measure();
for (int i = 0; i < exponent_size; ++i)
exponent += (chromosome[mantissa_size + i].measure() * std::pow(2, i));
exponent = std::pow(2, exponent - 127);
for (int i = 1; i <= mantissa_size; ++i)
mantissa += (chromosome[mantissa_size - i].measure() * std::pow(2, -i));
latest_interpretation = std::pow(-1, sign) * exponent * mantissa;
return latest_interpretation;
}
//Constructor copia de la clase individual
Individual::Individual(const Individual &o)
: m(o.m), fitness(o.fitness), latest_interpretation(o.latest_interpretation) {
chromosome = new Qubit[m];
for (int i = 0; i < m; ++i)
chromosome[i] = o.chromosome[i];
}
//Constructor de la clase individual
Individual::Individual(size_t size) {
m = size;
chromosome = new Qubit[m];
latest_interpretation = 0.f;
fitness = 0.f;
}
//destructor, elimina el array dinamico de cromosomas
Individual::~Individual() {
delete[] chromosome;
}
//Funcion que obtiene el numero flotante
const float Individual::getLatestInterpretation() const {
return latest_interpretation;
}
//funcion que como dice el nombre evalua el fitnes en cierto segmento de la funcion
void Individual::evaluate_fitness(const float a, const float b) {
float x = latest_interpretation;
fitness = (x < a || x > b) ? (float)INT_MIN
: 3 * std::cos(4 * x) + 4 * std::sin(3 * x) + x; //funcion a optimizar
if (isinf(fitness)) fitness = (float)INT_MIN;
}
//sobrecargas para tratar al individuo como un array para facilitar uso
Qubit& Individual::operator [] (size_t index) {
return chromosome[index];
}
const Qubit& Individual::operator [] (size_t index) const {
return chromosome[index];
}
//Un individuo es menor a otro si su fitnes es mayor
bool Individual::operator < (const Individual& o) const {
return (fitness > o.fitness);
}
//Sobrecarga del operador de asignacion, necesario para evitar conflictos posteriores
Individual& Individual::operator = (const Individual& o) {
if (this != &o) {
delete[] chromosome;
m = o.m;
fitness = o.fitness;
latest_interpretation = o.latest_interpretation;
chromosome = new Qubit[m];
for (int i = 0; i < m; ++i)
chromosome[i] = o.chromosome[i];
}
return *this;
}