Changeset 751 for cpp/frams/genetics/oper_fx.cpp

Timestamp:

02/28/18 23:41:58 (5 years ago)

Author:

Maciej Komosinski

Message:

More versatile mutation function for numbers

File:

• cpp/frams/genetics/oper_fx.cpp (modified) (4 diffs)

cpp/frams/genetics/oper_fx.cpp

@@ -80 +80 @@
 double GenoOperators::mutateNeuProperty(double current, Neuro *n, int i)
 {
-        if (i == -1) return mutateCreepNoLimit('f', current, -10, 10); //i==-1: mutating weight of neural connection
+        if (i == -1) return mutateCreepNoLimit('f', current, 2, true); //i==-1: mutating weight of neural connection
         Param p;
         if (i >= 100) { i -= 100; p = n->getClass()->getProperties(); }
@@ -97 +97 @@
         ExtValue ev;
         p.get(i, ev);
-        ev.setDouble(mutateCreep(p.type(i)[0], ev.getDouble(), mn, mx));
+        ev.setDouble(mutateCreep(p.type(i)[0], ev.getDouble(), mn, mx, true));
         p.set(i, ev);
         return true;
@@ -123 +123 @@
                 double mn, mx, df;
                 getMinMaxDef(&p, i, mn, mx, df);
-                newval = mutateCreep(p.type(i)[0], oldval, mn, mx);
+                newval = mutateCreep(p.type(i)[0], oldval, mn, mx, true);
                         }
         return true;
 }
 
-double GenoOperators::mutateCreepNoLimit(char type, double current, double mn, double mx)
-{
-        double result = RndGen.Gauss(current, (mx - mn) / 2 / 5); // /halfinterval, 5 times narrower
-        if (type == 'd') { result = int(result + 0.5); if (result == current) result += randomN(2) * 2 - 1; }
-        else result = floor(result * 1000 + 0.5) / 1000.0; //round
+double GenoOperators::mutateCreepNoLimit(char type, double current, double stddev, bool limit_precision_3digits)
+{
+        double result = RndGen.Gauss(current, stddev);
+        if (type == 'd')
+        {
+                result = int(result + 0.5);
+                if (result == current) result += randomN(2) * 2 - 1; //force some change
+        }
+        else
+        {
+                if (limit_precision_3digits)
+                        result = floor(result * 1000 + 0.5) / 1000.0; //round
+        }
         return result;
 }
 
-double GenoOperators::mutateCreep(char type, double current, double mn, double mx)
-{
-        double result = mutateCreepNoLimit(type, current, mn, mx); //TODO consider that when boundary is touched (reflect/absorb below), the default precision (3 digits) may change. Is it good or bad?
+double GenoOperators::mutateCreep(char type, double current, double mn, double mx, double stddev, bool limit_precision_3digits)
+{
+        double result = mutateCreepNoLimit(type, current, stddev, limit_precision_3digits);
+        //TODO consider that when boundary is touched (reflect+absorb below), the requested precision (3 digits) may change. Is it good or bad?
         //reflect:
         if (result > mx) result = mx - (result - mx); else
@@ -146 +155 @@
                 if (result < mn) result = mn;
         return result;
+}
+
+double GenoOperators::mutateCreep(char type, double current, double mn, double mx, bool limit_precision_3digits)
+{
+        double stddev = (mx - mn) / 2 / 5; // magic arbitrary formula for stddev, which becomes /halfinterval, 5 times narrower
+        return mutateCreep(type, current, mn, mx, stddev, limit_precision_3digits);
 }