extl/intelligence/ga/fopti/test/fopti_ga_test.h

   1 /* ///////////////////////////////////////////////////////////////////////
   2  * File:                fopti_ga_test.h
   3  *
   4  * Created:             08.11.13
   5  * Updated:             08.11.13
   6  *
   7  * Brief: Unit-testing
   8  *
   9  * [<Home>]
  10  * Copyright (c) 2008-2020, Waruqi All rights reserved.
  11  * //////////////////////////////////////////////////////////////////// */
  12 #ifndef EXTL_INTELLIGENCE_GA_FOPTI_GA_TEST_H
  13 #define EXTL_INTELLIGENCE_GA_FOPTI_GA_TEST_H
  14 #include "../fopti.h"
  15 /* ///////////////////////////////////////////////////////////////////////
  16  * unit_test namespace
  17  */
  18 EXTL_INTELLIGENCE_BEGIN_WHOLE_NAMESPACE
  19 EXTL_TEST_NAME_BEGIN_NAMESPACE(fopti_ga_test)
  20
  21 /* ///////////////////////////////////////////////////////////////////////
  22  * Unit-testing
  23  */
  24 struct fopti_ga_test
  25 {
  26         // for maximum optimization
  27         // [-1, 20]
  28         // x=19.850052
  29         // max f(x)=21.850026
  30         static double ga1_func(double x)
  31         {
  32                 return x * sin(10 * 3.1415926 * x) + 2.0;
  33         }
  34
  35         // x >= -1 => -1 * (x + 1) <= 0
  36         static double ga1_cons_1(double x)
  37         {
  38                 return -1 *(x + 1);
  39         }
  40
  41         // x <= 20 => (x - 20) <= 0
  42         static double ga1_cons_2(double x)
  43         {
  44                 return (x - 20);
  45         }
  46
  47         // for minimum optimization
  48         // [-100, 100]
  49         // (x, y): (-0.0898, 0.7126) or (0.0898, -0.7126)
  50         // min f(x, y)=-1.031628
  51         static double ga2_func(double x, double y)
  52         {
  53                 return (4 - 2.1 * x * x + x * x * x * x / 3) * x * x + x * y + (-4 + 4 * y * y) * y * y;
  54         }
  55
  56         // x >= -100 => -1 * (x + 100) <= 0
  57         static double ga2_cons_1(double x, double /*y*/)
  58         {
  59                 return -1 *(x + 100);
  60         }
  61
  62         // x <= 100 => (x - 100) <= 0
  63         static double ga2_cons_2(double x, double /*y*/)
  64         {
  65                 return (x - 100);
  66         }
  67
  68         // y >= -100 => -1 * (y + 100) <= 0
  69         static double ga2_cons_3(double /*x*/, double y)
  70         {
  71                 return -1 *(y + 100);
  72         }
  73
  74         // y <= 100 => (y - 100) <= 0
  75         static double ga2_cons_4(double /*x*/, double y)
  76         {
  77                 return (y - 100);
  78         }
  79
  80         // for maximum optimization
  81         // [-100, 100]
  82         // (x, y): (0, 0)
  83         // min f(x, y)=1
  84         static double ga3_func(double x, double y)
  85         {
  86                 double a = sin(sqrt(x * x + y * y));
  87                 double b = (1 + 0.001 * (x * x + y * y));
  88                 return 0.5 - (a * a - 0.5) / (b * b);
  89         }
  90
  91         // x >= -100 => -1 * (x + 100) <= 0
  92         static double ga3_cons_1(double x, double /*y*/)
  93         {
  94                 return -1 *(x + 100);
  95         }
  96
  97         // x <= 100 => (x - 100) <= 0
  98         static double ga3_cons_2(double x, double /*y*/)
  99         {
 100                 return (x - 100);
 101         }
 102
 103         // y >= -100 => -1 * (y + 100) <= 0
 104         static double ga3_cons_3(double /*x*/, double y)
 105         {
 106                 return -1 *(y + 100);
 107         }
 108
 109         // y <= 100 => (y - 100) <= 0
 110         static double ga3_cons_4(double /*x*/, double y)
 111         {
 112                 return (y - 100);
 113         }
 114
 115         // for minimum optimization
 116         // [0, ...]
 117         // (x, y): (1, 1)
 118         // min f(x, y)=10
 119         static double ga4_func(double x, double y)
 120         {
 121                 return x * x + y * y + 8;
 122         }
 123         // x >= 0 => -1 * (x) <= 0
 124         static double ga4_cons_1(double x, double /*y*/)
 125         {
 126                 return -1 * x;
 127         }
 128
 129         // y >= 0 => -1 * (y) <= 0
 130         static double ga4_cons_2(double /*x*/, double y)
 131         {
 132                 return -1 * y ;
 133         }
 134
 135         // x * x - y >= 0 => -1 * (x * x - y) <= 0
 136         static double ga4_cons_3(double x, double y)
 137         {
 138                 return -1 * (x * x - y);
 139         }
 140
 141         // -x - y * y + 2 = 0 => (-x - y * y + 2)^2 = 0
 142         static double ga4_cons_4(double x, double y)
 143         {
 144                 return (2 - x - y * y) * (2 - x - y * y);
 145         }
 146
 147         fopti_ga_test()
 148         {
 149                 EXTL_TEST_TRACE(_T("\nfopti_ga_test:"));
 150
 151                 typedef func_ptr<double (*)(double)>                    ga1_func_type;
 152                 typedef fopti_ga<ga1_func_type>                                 ga1_ga_type;
 153                 typedef ga1_ga_type::scopes_type                                ga1_scopes_type;
 154                 typedef ga1_scopes_type::scope_type                             ga1_scope_type;
 155                 typedef ga1_ga_type::constraints_type                   ga1_constraints_type;
 156
 157                 ga1_ga_type ga1_ga      (       ga1_func_type(&ga1_func)
 158                                                         ,       ga1_scopes_type(ga1_scope_type(-1, 20))
 159                                                         ,       ga1_constraints_type(ga1_func_type(&ga1_cons_1), ga1_func_type(&ga1_cons_2))
 160                                                         );
 161                 ga1_ga.generate();
 162
 163                 typedef func_ptr<double (*)(double, double)>    ga2_func_type;
 164                 typedef fopti_ga<ga2_func_type, e_false_v>              ga2_ga_type;
 165                 typedef ga2_ga_type::scopes_type                                ga2_scopes_type;
 166                 typedef ga2_scopes_type::scope_type                             ga2_scope_type;
 167                 typedef ga2_ga_type::constraints_type                   ga2_constraints_type;
 168
 169                 ga2_ga_type ga2_ga      (       ga2_func_type(&ga2_func)
 170                                                         ,       ga2_scopes_type(ga2_scope_type(-100, 100), ga2_scope_type(-100, 100))
 171                                                         ,       ga2_constraints_type(ga2_func_type(&ga2_cons_1), ga2_func_type(&ga2_cons_2), ga2_func_type(&ga2_cons_3), ga2_func_type(&ga2_cons_4))
 172                                                         );
 173                 ga2_ga.generate();
 174
 175                 typedef func_ptr<double (*)(double, double)>    ga3_func_type;
 176                 typedef fopti_ga<ga3_func_type>                                 ga3_ga_type;
 177                 typedef ga3_ga_type::scopes_type                                ga3_scopes_type;
 178                 typedef ga3_scopes_type::scope_type                             ga3_scope_type;
 179                 typedef ga3_ga_type::constraints_type                   ga3_constraints_type;
 180
 181                 ga3_ga_type ga3_ga      (       ga3_func_type(&ga3_func)
 182                                                         ,       ga3_scopes_type(ga3_scope_type(-100, 100), ga3_scope_type(-100, 100))
 183                                                         ,       ga3_constraints_type(ga3_func_type(&ga3_cons_1), ga3_func_type(&ga3_cons_2), ga3_func_type(&ga3_cons_3), ga3_func_type(&ga3_cons_4))
 184                                                         );
 185                 ga3_ga.generate();
 186
 187                 typedef func_ptr<double (*)(double, double)>    ga4_func_type;
 188                 typedef fopti_ga<ga4_func_type, e_false_v>              ga4_ga_type;
 189                 typedef ga4_ga_type::scopes_type                                ga4_scopes_type;
 190                 typedef ga4_scopes_type::scope_type                             ga4_scope_type;
 191                 typedef ga4_ga_type::constraints_type                   ga4_constraints_type;
 192
 193                 ga4_ga_type ga4_ga      (       ga4_func_type(&ga4_func)
 194                                                         ,       ga4_scopes_type(ga4_scope_type(0, 100), ga4_scope_type(0, 100))
 195                                                         ,       ga4_constraints_type(ga4_func_type(&ga4_cons_1), ga4_func_type(&ga4_cons_2), ga4_func_type(&ga4_cons_3), ga4_func_type(&ga4_cons_4))
 196                                                         );
 197                 ga4_ga.generate();
 198         }
 199 };
 200
 201 fopti_ga_test g_fopti_ga_test;
 202
 203 /* ///////////////////////////////////////////////////////////////////////
 204  * unit_test namespace
 205  */
 206 EXTL_TEST_NAME_END_NAMESPACE(fopti_ga_test)
 207 EXTL_INTELLIGENCE_END_WHOLE_NAMESPACE
 208
 209 /* //////////////////////////////////////////////////////////////////// */
 210 #endif /* EXTL_INTELLIGENCE_GA_FOPTI_GA_TEST_H */
 211 /* //////////////////////////////////////////////////////////////////// */
 212