-
Notifications
You must be signed in to change notification settings - Fork 1
/
Species.js
197 lines (177 loc) · 7.67 KB
/
Species.js
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
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
class Species {
constructor(p) {
this.players = [];
this.bestFitness = 0;
this.champ;
this.averageFitness = 0;
this.staleness = 0; //how many generations the species has gone without an improvement
this.rep;
//--------------------------------------------
//coefficients for testing compatibility
this.excessCoeff = 1;
this.weightDiffCoeff = 0.5;
this.compatibilityThreshold = 3;
if (p) {
this.players.push(p);
//since it is the only one in the species it is by default the best
this.bestFitness = p.fitness;
this.rep = p.brain.clone();
this.champ = p.cloneForReplay();
}
}
//------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
//returns whether the parameter genome is in this species
sameSpecies(g) {
var compatibility;
var excessAndDisjoint = this.getExcessDisjoint(g, this.rep); //get the number of excess and disjoint genes between this player and the current species this.rep
var averageWeightDiff = this.averageWeightDiff(g, this.rep); //get the average weight difference between matching genes
var largeGenomeNormaliser = g.genes.length - 20;
if (largeGenomeNormaliser < 1) {
largeGenomeNormaliser = 1;
}
compatibility = (this.excessCoeff * excessAndDisjoint / largeGenomeNormaliser) + (this.weightDiffCoeff * averageWeightDiff); //compatibility formula
return (this.compatibilityThreshold > compatibility);
}
//------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
//add a player to the species
addToSpecies(p) {
this.players.push(p);
}
//------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
//returns the number of excess and disjoint genes between the 2 input genomes
//i.e. returns the number of genes which dont match
getExcessDisjoint(brain1, brain2) {
var matching = 0.0;
for (var i = 0; i < brain1.genes.length; i++) {
for (var j = 0; j < brain2.genes.length; j++) {
if (brain1.genes[i].innovationNo == brain2.genes[j].innovationNo) {
matching++;
break;
}
}
}
return (brain1.genes.length + brain2.genes.length - 2 * (matching)); //return no of excess and disjoint genes
}
//------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
//returns the avereage weight difference between matching genes in the input genomes
averageWeightDiff(brain1, brain2) {
if (brain1.genes.length == 0 || brain2.genes.length == 0) {
return 0;
}
var matching = 0;
var totalDiff = 0;
for (var i = 0; i < brain1.genes.length; i++) {
for (var j = 0; j < brain2.genes.length; j++) {
if (brain1.genes[i].innovationNo == brain2.genes[j].innovationNo) {
matching++;
totalDiff += abs(brain1.genes[i].weight - brain2.genes[j].weight);
break;
}
}
}
if (matching == 0) { //divide by 0 error
return 100;
}
return totalDiff / matching;
}
//-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
//sorts the species by fitness
sortSpecies() {
var temp = []; // new ArrayList < Player > ();
//selection short
for (var i = 0; i < this.players.length; i++) {
var max = 0;
var maxIndex = 0;
for (var j = 0; j < this.players.length; j++) {
if (this.players[j].fitness > max) {
max = this.players[j].fitness;
maxIndex = j;
}
}
temp.push(this.players[maxIndex]);
this.players.splice(maxIndex, 1);
// this.players.remove(maxIndex);
i--;
}
// this.players = (ArrayList) temp.clone();
arrayCopy(temp, this.players);
if (this.players.length == 0) {
this.staleness = 200;
return;
}
//if new best player
if (this.players[0].fitness > this.bestFitness) {
this.staleness = 0;
this.bestFitness = this.players[0].fitness;
this.rep = this.players[0].brain.clone();
this.champ = this.players[0].cloneForReplay();
} else { //if no new best player
this.staleness++;
}
}
//-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
//simple stuff
setAverage() {
var sum = 0;
for (var i = 0; i < this.players.length; i++) {
sum += this.players[i].fitness;
}
this.averageFitness = sum / this.players.length;
}
//-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
//gets baby from the this.players in this species
giveMeBaby(innovationHistory) {
var baby;
if (random(1) < 0.25) { //25% of the time there is no crossover and the child is simply a clone of a random(ish) player
baby = this.selectPlayer().clone();
} else { //75% of the time do crossover
//get 2 random(ish) parents
var parent1 = this.selectPlayer();
var parent2 = this.selectPlayer();
//the crossover function expects the highest fitness parent to be the object and the lowest as the argument
if (parent1.fitness < parent2.fitness) {
baby = parent2.crossover(parent1);
} else {
baby = parent1.crossover(parent2);
}
}
baby.brain.mutate(innovationHistory); //mutate that baby brain
return baby;
}
//-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
//selects a player based on it fitness
selectPlayer() {
var fitnessSum = 0;
for (var i = 0; i < this.players.length; i++) {
fitnessSum += this.players[i].fitness;
}
var rand = random(fitnessSum);
var runningSum = 0;
for (var i = 0; i < this.players.length; i++) {
runningSum += this.players[i].fitness;
if (runningSum > rand) {
return this.players[i];
}
}
//unreachable code to make the parser happy
return this.players[0];
}
//------------------------------------------------------------------------------------------------------------------------------------------
//kills off bottom half of the species
cull() {
if (this.players.length > 2) {
for (var i = this.players.length / 2; i < this.players.length; i++) {
// this.players.remove(i);
this.players.splice(i, 1);
i--;
}
}
}
//------------------------------------------------------------------------------------------------------------------------------------------
//in order to protect unique this.players, the fitnesses of each player is divided by the number of this.players in the species that that player belongs to
fitnessSharing() {
for (var i = 0; i < this.players.length; i++) {
this.players[i].fitness /= this.players.length;
}
}
}