00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039 #ifdef _MSC_VER
00040
00041 #pragma warning(disable:4786)
00042
00043 #endif
00044
00045
00046
00047 #ifdef HAVE_CONFIG_H
00048
00049 #include <config.h>
00050
00051 #endif
00052
00053
00054
00055 #include <iostream>
00056
00057 #include "gp/eoParseTree.h"
00058
00059 #include "eo"
00060
00061
00062
00063 using namespace gp_parse_tree;
00064
00065 using namespace std;
00066
00067
00068
00069
00070
00071
00072
00073 #include "node.h"
00074
00075 #include "parameters.h"
00076
00077 #include "fitness.h"
00078
00079
00080
00081
00082
00083
00084
00085
00086
00087
00088
00089 typedef eoParseTree<FitnessType, Node > EoType;
00090
00091 typedef eoPop<EoType> Pop;
00092
00093
00094
00095
00096
00097
00098
00099 int main(int argc, char *argv[])
00100
00101 {
00102
00103
00104
00105
00106
00107 vector<Node> initSequence;
00108
00109
00110
00111
00112
00113 Parameters parameter(argc, argv);
00114
00115
00116
00117
00118
00119 rng.reseed(parameter.randomseed);
00120
00121
00122
00123
00124
00125 eoValueParam<unsigned> generationCounter(0, "Gen.");
00126
00127
00128
00129
00130
00131
00132
00133
00134
00135
00136
00137 eoIncrementor<unsigned> increment(generationCounter.value());
00138
00139
00140
00141
00142
00143
00144
00145
00146
00147
00148
00149 RegFitness eval(generationCounter, initSequence, parameter);
00150
00151
00152
00153
00154
00155 eoParseTreeDepthInit<FitnessType, Node> initializer(parameter.InitMaxDepth, initSequence, true, true);
00156
00157
00158
00159
00160
00161 Pop pop(parameter.population_size, initializer);
00162
00163
00164
00165
00166
00167 apply<EoType>(eval, pop);
00168
00169
00170
00171 generationCounter.value()++;
00172
00173
00174
00175
00176
00177
00178
00179
00180
00181 eoSubtreeXOver<FitnessType, Node> xover(parameter.MaxSize);
00182
00183
00184
00185
00186
00187 eoBranchMutation<FitnessType, Node> mutation(initializer, parameter.MaxSize);
00188
00189
00190
00191
00192
00193
00194
00195
00196
00197
00198
00199
00200
00201
00202
00203 eoSGATransform<EoType> transform(xover, parameter.xover_rate, mutation, parameter.mutation_rate);
00204
00205
00206
00207
00208
00209
00210
00211 eoDetTournamentSelect<EoType> selectOne(parameter.tournamentsize);
00212
00213
00214
00215 eoSelectMany<EoType> select(selectOne, parameter.offspring_size, eo_is_an_integer);
00216
00217
00218
00219
00220
00221
00222
00223
00224
00225
00226
00227
00228
00229 eoCommaReplacement<EoType> replace;
00230
00231
00232
00233
00234
00235 eoGenContinue<EoType> term(parameter.nGenerations);
00236
00237
00238
00239 eoCheckPoint<EoType> checkPoint(term);
00240
00241
00242
00243
00244
00245 eoAverageStat<EoType> avg;
00246
00247 eoBestFitnessStat<EoType> best;
00248
00249
00250
00251
00252
00253
00254
00255
00256
00257 checkPoint.add(increment);
00258
00259 checkPoint.add(avg);
00260
00261 checkPoint.add(best);
00262
00263
00264
00265 #ifdef HAVE_GNUPLOT
00266
00267 eoGnuplot1DMonitor gnuplotmonitor("gnuplotBestStats");
00268
00269 gnuplotmonitor.add(generationCounter);
00270
00271 gnuplotmonitor.add(best);
00272
00273
00274
00275 eoValueParam<string> dummy1("", "Smallest Tree Size");
00276
00277 gnuplotmonitor.add(dummy1);
00278
00279
00280
00281 eoGnuplot1DMonitor gnuplotAvgmonitor("gnuplotAvgStats");
00282
00283 gnuplotAvgmonitor.add(generationCounter);
00284
00285 gnuplotAvgmonitor.add(avg);
00286
00287
00288
00289 eoValueParam<string> dummy2("", "Average Tree Size");
00290
00291 gnuplotAvgmonitor.add(dummy2);
00292
00293
00294
00295 checkPoint.add(gnuplotmonitor);
00296
00297 checkPoint.add(gnuplotAvgmonitor);
00298
00299 #endif
00300
00301
00302
00303 eoEasyEA<EoType> gp(checkPoint, eval, select, transform, replace);
00304
00305
00306
00307 cout << "Initialization done" << endl;
00308
00309
00310
00311
00312
00313 try
00314
00315 {
00316
00317 gp(pop);
00318
00319 }
00320
00321 catch (exception& e)
00322
00323 {
00324
00325 cout << "exception: " << e.what() << endl;;
00326
00327 exit(EXIT_FAILURE);
00328
00329 }
00330
00331
00332
00333 return 1;
00334
00335
00336
00337 }
00338
00339
00340
00341
00342
00343
00344