Changeset 17792 for branches/3076_IA_evaluators_analyzers/HeuristicLab.Problems.DataAnalysis.Symbolic.Regression/3.4/SymbolicRegressionConstraintAnalyzer.cs

Timestamp:

12/15/20 14:52:22 (3 years ago)

Author:

dpiringe

Message:

#3076

• updated SymbolicRegressionSingleObjectiveMetaModelAnalyzer to find a meta model with and without constant optimization
• added a new result parameter of type DataTable in SymbolicRegressionConstraintAnalyzer to show the absolute number of unsatisfied individuals
• added soft constraints handling for SymbolicRegressionSingleObjectiveConstraintEvaluator

File:

• branches/3076_IA_evaluators_analyzers/HeuristicLab.Problems.DataAnalysis.Symbolic.Regression/3.4/SymbolicRegressionConstraintAnalyzer.cs (modified) (7 diffs)

branches/3076_IA_evaluators_analyzers/HeuristicLab.Problems.DataAnalysis.Symbolic.Regression/3.4/SymbolicRegressionConstraintAnalyzer.cs

@@ -25 +25 @@
     private const string MaximumStepsParameterName = "Maximum Steps";
     private const string StartPenalityParameterName = "Start Penality";
+    private const string ConstraintUnsatisfiedSolutionsParameterName = "Constraint Unsatisfied Solutions";
 
     #region parameter properties
@@ -49 +50 @@
       (IValueParameter<DoubleValue>) Parameters[StartPenalityParameterName];
 
+    public IResultParameter<DataTable> ConstraintUnsatisfiedSolutionsParameter =>
+      (IResultParameter<DataTable>)Parameters[ConstraintUnsatisfiedSolutionsParameterName];
+
     #endregion
 
     public override bool EnabledByDefault => false;
     public static int Iterations { get; set; } = 0;
+
+    public IBoundsEstimator BoundsEstimator { get; set; } = new IABoundsEstimator();
 
     [StorableConstructor]
@@ -79 +85 @@
       Parameters.Add(new ValueParameter<DoubleValue>(StartPenalityParameterName,
         "The start value for the penality multiplier.", new DoubleValue(0.5)));
+      Parameters.Add(new ResultParameter<DataTable>(ConstraintUnsatisfiedSolutionsParameterName,
+        "Shows the number of solutions with unsatisfied constraints."));
 
 
@@ -92 +100 @@
           XAxisTitle = "Generations",
           YAxisTitle = "Penality Multiplier"
+        }
+      };
+
+      ConstraintUnsatisfiedSolutionsParameter.DefaultValue = new DataTable(ConstraintUnsatisfiedSolutionsParameterName) {
+        VisualProperties = {
+          XAxisTitle = "Generations",
+          YAxisTitle = "Constraint Unsatisfied Solutions"
         }
       };
@@ -137 +152 @@
         Parameters.Add(new ValueParameter<DoubleValue>(StartPenalityParameterName,
           "The start value for the penality multiplier.", new DoubleValue(0.5)));
+
+
+      if(!Parameters.ContainsKey(ConstraintUnsatisfiedSolutionsParameterName)) {
+        Parameters.Add(new ResultParameter<DataTable>(ConstraintUnsatisfiedSolutionsParameterName,
+          "Shows the number of solutions with unsatisfied constraints."));
+
+        ConstraintUnsatisfiedSolutionsParameter.DefaultValue = new DataTable(ConstraintUnsatisfiedSolutionsParameterName) {
+            VisualProperties = {
+            XAxisTitle = "Generations",
+            YAxisTitle = "Constraint Unsatisfied Solutions"
+          }
+        };
+      }
     }
 
@@ -158 +186 @@
       var constraints = problemData.IntervalConstraints.EnabledConstraints;
       var variableRanges = problemData.VariableRanges.GetReadonlyDictionary();
+      var intervalCollection = problemData.VariableRanges;
       var newDataTable = ConstraintViolationParameter.ActualValue;
+      var solutions = this.SymbolicExpressionTree.ToArray();
 
       if (newDataTable.Rows.Count == 0)
@@ -176 +206 @@
         penalityDataTable.Rows.Add(new DataRow(rowName));
       penalityDataTable.Rows[rowName].Values.Add(penalityMultiplier.Value);
+
+      var constraintUnsatisfiedSolutionsDataTable = ConstraintUnsatisfiedSolutionsParameter.ActualValue;
+      if (constraintUnsatisfiedSolutionsDataTable.Rows.Count == 0)
+        constraintUnsatisfiedSolutionsDataTable.Rows.Add(new DataRow(ConstraintUnsatisfiedSolutionsParameterName));
+
+      constraintUnsatisfiedSolutionsDataTable.Rows[ConstraintUnsatisfiedSolutionsParameterName]
+        .Values
+        .Add(
+          solutions
+          .Where(s => IntervalUtil.IntervalConstraintsViolation(constraints, BoundsEstimator, intervalCollection, s).Any(x => x != 0.0))
+          .Count());
 
       return base.Apply();