Changeset 968 for cpp/frams/genetics/genooperators.cpp

Timestamp:

06/30/20 00:30:39 (5 years ago)

Author:

Maciej Komosinski

Message:

Introduced a separate function to mutate neuron connection weight

File:

• cpp/frams/genetics/genooperators.cpp (modified) (8 diffs)

cpp/frams/genetics/genooperators.cpp ¶

@@ -9 +9 @@
 #include <frams/util/rndutil.h>
 
+//
+// custom distributions for mutations of various parameters
+//
 static double distrib_force[] =   // for '!'
 {
@@ -31 +34 @@
         -1, 1,         // ~linear
 };
-
+/*
+static double distrib_weight[] =
+{
+5,                 // distribution -999 _-^_^-_ +999
+-999, 999,         // each weight value may be useful, especially...
+-5, -0.3,          // ...little non-zero values
+-3, -0.6,
+0.6, 3,
+0.3, 5,
+};
+*/
 
 int GenoOperators::roulette(const double *probtab, const int count)
@@ -108 +121 @@
 double GenoOperators::getMutatedNeuroClassProperty(double current, Neuro *n, int i)
 {
-        if (i == -1) return mutateCreepNoLimit('f', current, 2, true); //i==-1: mutating weight of neural connection
+        if (i == -1)
+        {
+                logPrintf("GenoOperators", "getMutatedNeuroClassProperty", LOG_WARN, "Deprecated usage in C++ source: to mutate connection weight, use getMutatedNeuronConnectionWeight().");
+                return getMutatedNeuronConnectionWeight(current);
+        }
         Param p;
         if (i >= NEUROCLASS_PROP_OFFSET) { i -= NEUROCLASS_PROP_OFFSET; p = n->getClass()->getProperties(); }
@@ -117 +134 @@
 }
 
+double GenoOperators::getMutatedNeuronConnectionWeight(double current)
+{
+        return mutateCreepNoLimit('f', current, 2, true);
+}
+
 bool GenoOperators::mutatePropertyNaive(ParamInterface &p, int i)
 {
@@ -145 +167 @@
         if (p.type(i)[0] != 'f' && p.type(i)[0] != 'd') return false; //don't know how to mutate
         const char *n = p.id(i), *na = p.name(i);
-        if (strcmp(n, "si") == 0 && strcmp(na, "Sigmoid") == 0) newval = CustomRnd(distrib_sigmo); else
-                if (strcmp(n, "in") == 0 && strcmp(na, "Inertia") == 0) newval = CustomRnd(distrib_inertia); else
-                        if (strcmp(n, "fo") == 0 && strcmp(na, "Force") == 0) newval = CustomRnd(distrib_force); else
+        if (strcmp(n, "si") == 0 && strcmp(na, "Sigmoid") == 0) newval = round(CustomRnd(distrib_sigmo), 3); else
+                if (strcmp(n, "in") == 0 && strcmp(na, "Inertia") == 0) newval = round(CustomRnd(distrib_inertia), 3); else
+                        if (strcmp(n, "fo") == 0 && strcmp(na, "Force") == 0) newval = round(CustomRnd(distrib_force), 3); else
                         {
                                 double mn, mx, df;
@@ -167 +189 @@
         {
                 if (limit_precision_3digits)
-                        result = floor(result * 1000 + 0.5) / 1000.0; //round
+                        result = round(result, 3);
         }
         return result;
@@ -186 +208 @@
                 {
                         //reflect and wrap above may have changed the (limited) precision, so try to round again (maybe unnecessarily, because we don't know if reflect+wrap above were triggered)
-                        double result_try = floor(result * 1000 + 0.5) / 1000.0; //round
+                        double result_try = round(result, 3);
                         if (mn <= result_try && result_try <= mx) result = result_try; //after rounding still witin allowed range, so keep rounded value
                 }
@@ -375 +397 @@
                         int l = (int)strlen(n);
                         if (len >= l && l > Len && (strncmp(s, n, l) == 0)) { I = NEUROCLASS_PROP_OFFSET + i; Len = l; }
-                        if (also_v1_N_props) //recognize old properties symbols /=!
+                        if (also_v1_N_props) //recognize old symbols of properties:  /=!
                         {
                                 if (strcmp(n, "si") == 0) n = "/"; else