Changeset 15001 for branches/jschiess/HeuristicLab.Algorithms.SimulatedAnnealing/3.4/AcceptanceRatioReheatingOperator.cs

Timestamp:

05/30/17 07:50:37 (7 years ago)

Author:

jschiess

Message:

#1836 SA reheating strategies: Newest version of reheating

File:

• branches/jschiess/HeuristicLab.Algorithms.SimulatedAnnealing/3.4/AcceptanceRatioReheatingOperator.cs (modified) (6 diffs)

branches/jschiess/HeuristicLab.Algorithms.SimulatedAnnealing/3.4/AcceptanceRatioReheatingOperator.cs

@@ -37 +37 @@
         private const string IsAcceptedName = "IsAccepted";
         private const string AcceptanceMemoryName = "AcceptanceMemory";
+        private const string AverageAcceptanceRatioName = "AverageAcceptanceRatio";
 
 
@@ -73 +74 @@
             get { return (ILookupParameter<IntValue>)Parameters[MaximumIterationsName]; }
         }
-        private ValueParameter<DoubleValue> UpperTemperatureParameter
-        {
-            get { return (ValueParameter<DoubleValue>)Parameters[UpperTemperatureName]; }
+        private ValueLookupParameter<DoubleValue> UpperTemperatureParameter
+        {
+            get { return (ValueLookupParameter<DoubleValue>)Parameters[UpperTemperatureName]; }
         }
         private IConstrainedValueParameter<IDiscreteDoubleValueModifier> ReheatingOperatorParameter
@@ -100 +101 @@
         {
             get { return (ILookupParameter<ItemList<BoolValue>>)Parameters[AcceptanceMemoryName]; }
+        }
+        public ILookupParameter<DoubleValue> AverageAcceptanceRatioParameter
+        {
+            get { return (ILookupParameter<DoubleValue>)Parameters[AverageAcceptanceRatioName]; }
         }
 
@@ -118 +123 @@
             Parameters.Add(new LookupParameter<ItemList<BoolValue>>(AcceptanceMemoryName, "Memorizes the last N acceptance decisions."));
             Parameters.Add(new LookupParameter<BoolValue>(CoolingName, "True when the temperature should be cooled, false otherwise."));
-            Parameters.Add(new ValueParameter<DoubleValue>(UpperTemperatureName, "The upper bound of the temperature.", new DoubleValue(100d)));
+            Parameters.Add(new ValueLookupParameter<DoubleValue>(UpperTemperatureName, "The upper bound of the temperature.", "InitialTemperature"));
             Parameters.Add(new ConstrainedValueParameter<IDiscreteDoubleValueModifier>(ReheatingOperatorName, "The operator that reheats the temperature."));
-            Parameters.Add(new ValueParameter<DoubleRange>(ThresholdName, "The threshold controls the temperature. If the average ratio of accepted moves goes below the start of the range the temperature is heated. If the the average ratio of accepted moves goes beyond the end of the range the temperature is cooled again.", new DoubleRange(0.2, 0.1)));
+            Parameters.Add(new ValueParameter<DoubleRange>(ThresholdName, "The threshold controls the temperature. If the average ratio of accepted moves goes below the start of the range the temperature is heated. If the the average ratio of accepted moves goes beyond the end of the range the temperature is cooled again.", new DoubleRange(0.01, 0.1)));
             Parameters.Add(new ValueParameter<IntValue>(MemorySizeName, "The maximum size of the acceptance memory.", new IntValue(100)));
-
+            Parameters.Add(new LookupParameter<DoubleValue>(AverageAcceptanceRatioName, "Average acceptance over full acceptance memory."));
             #endregion
 
@@ -155 +160 @@
                 op.StartValueParameter.ActualName = StartTemperatureName;
                 op.StartValueParameter.Hidden = true;
-                op.EndValueParameter.ActualName = UpperTemperatureName;
-                op.EndValueParameter.Value = UpperTemperatureParameter.Value;
+                op.EndValueParameter.ActualName = EndTemperatureName;
                 op.EndValueParameter.Hidden = true;
             }
@@ -163 +167 @@
         public override IOperation Apply()
         {
-
+            var isAccepted = IsAcceptedParameter.ActualValue;
             var acceptances = AcceptanceMemoryParameter.ActualValue;
             var cooling = CoolingParameter.ActualValue.Value;
-
-            acceptances.Add(IsAcceptedParameter.ActualValue);
-
-            if (acceptances.Count > MemorySizeParameter.Value.Value) acceptances.RemoveAt(0);
-
-            var ratioStart = ThresholdParameter.Value.Start;
-            var ratioEnd = ThresholdParameter.Value.End;
-            var ratio = acceptances.Average(x => x.Value ? 1.0 : 0.0);
-
-            
-            if (!cooling && ratio >= ratioEnd)
+            var ratioAmount = AverageAcceptanceRatioParameter.ActualValue;
+
+            acceptances.Add(isAccepted);
+
+            ratioAmount.Value += isAccepted.Value ? 1 : 0;
+
+
+            if (acceptances.Count > MemorySizeParameter.Value.Value) {
+                ratioAmount.Value -= acceptances.ElementAt(0).Value ? 1 : 0;
+                acceptances.RemoveAt(0);
+            }
+
+            // only reheat when at least MemorySizeParameter.Value iterations have passed
+            if (acceptances.Count == MemorySizeParameter.Value.Value)
             {
-                // if we are heating, but should be cooling
-                cooling = true;
-                TemperatureStartIndexParameter.ActualValue.Value = Math.Max(0, IterationsParameter.ActualValue.Value - 1);
-                StartTemperatureParameter.ActualValue.Value = TemperatureParameter.ActualValue.Value;
-                EndTemperatureParameter.ActualValue.Value = LowerTemperatureParameter.ActualValue.Value;
-            }
-            else if (cooling && ratio <= ratioStart)
-            {
-                // if we are cooling but should be heating
-                cooling = false;
-                TemperatureStartIndexParameter.ActualValue.Value = Math.Max(0, IterationsParameter.ActualValue.Value - 1);
-                StartTemperatureParameter.ActualValue.Value = TemperatureParameter.ActualValue.Value;
-                EndTemperatureParameter.ActualValue.Value = UpperTemperatureParameter.Value.Value;
-            }
-
-            CoolingParameter.ActualValue.Value = cooling;
-
+                var ratio = ratioAmount.Value / MemorySizeParameter.Value.Value;
+                var ratioStart = ThresholdParameter.Value.Start;
+                var ratioEnd = ThresholdParameter.Value.End;
+
+                if (!cooling && ratio >= ratioEnd)
+                {
+                    // if we are heating, but should be cooling
+                    cooling = true;
+                    TemperatureStartIndexParameter.ActualValue.Value = Math.Max(0, IterationsParameter.ActualValue.Value - 1);
+                    StartTemperatureParameter.ActualValue.Value = TemperatureParameter.ActualValue.Value;
+                    EndTemperatureParameter.ActualValue.Value = LowerTemperatureParameter.ActualValue.Value;
+                }
+                else if (cooling && ratio <= ratioStart)
+                {
+                    // if we are cooling but should be heating
+                    cooling = false;
+                    TemperatureStartIndexParameter.ActualValue.Value = Math.Max(0, IterationsParameter.ActualValue.Value - 1);
+                    StartTemperatureParameter.ActualValue.Value = TemperatureParameter.ActualValue.Value;
+                    EndTemperatureParameter.ActualValue.Value = UpperTemperatureParameter.ActualValue.Value;
+                }
+
+                CoolingParameter.ActualValue.Value = cooling;
+            }
             return cooling ? Cool() : Heat();
         }