Changeset 17325 for branches/2994-AutoDiffForIntervals/HeuristicLab.Tests/HeuristicLab.Problems.DataAnalysis.Symbolic-3.4/IntervalEvaluatorAutoDiffTest.cs

Timestamp:

10/09/19 11:13:11 (5 years ago)

Author:

gkronber

Message:

#2994: worked on ConstrainedNLS

File:

• branches/2994-AutoDiffForIntervals/HeuristicLab.Tests/HeuristicLab.Problems.DataAnalysis.Symbolic-3.4/IntervalEvaluatorAutoDiffTest.cs (modified) (4 diffs)

branches/2994-AutoDiffForIntervals/HeuristicLab.Tests/HeuristicLab.Problems.DataAnalysis.Symbolic-3.4/IntervalEvaluatorAutoDiffTest.cs

@@ -61 +61 @@
       var eval = new IntervalEvaluator();
       var parser = new InfixExpressionParser();
+      var intervals = new Dictionary<string, Interval>() {
+        { "x", new Interval(1, 2) },
+        { "unit", new Interval(0, 1) },
+        { "neg", new Interval(-1, 0) },
+      };
       var t = parser.Parse("sqr(x)");
       var paramNodes = t.IterateNodesPostfix().Where(n => n.SubtreeCount == 0).ToArray();
-      var intervals = new Dictionary<string, Interval>() {
-        { "x", new Interval(1, 2) },
-        { "y", new Interval(0, 1) }
-      };
-      var r = eval.Evaluate(t, intervals, paramNodes, out double[] lg, out double[] ug);
-      // TODO
-      // Assert.AreEqual(XXX, r.LowerBound);
-      // Assert.AreEqual(XXX, r.UpperBound);
-      // 
-      // Assert.AreEqual(XXX, lg[0]); // x
-      // Assert.AreEqual(XXX, ug[0]);
-      // 
-      // for  { "x", new Interval(1, 2) },
-      //   { "y", new Interval(0, 1) },
-      // 
-      // 0 <> -2,50012500572888E-05 for y in SQR(LOG('y'))
-      // 0 <> 2, 49987500573946E-05 for x in SQR(LOG('x'))
+      var r = eval.Evaluate(t, intervals, paramNodes, out double[] lg, out double[] ug);
+      Assert.AreEqual(1, r.LowerBound);
+      Assert.AreEqual(4, r.UpperBound);
+
+      Assert.AreEqual(2.0, lg[0]); // x
+      Assert.AreEqual(8.0, ug[0]);
+
+      t = parser.Parse("sqr(log(unit))");
+      paramNodes = t.IterateNodesPostfix().Where(n => n.SubtreeCount == 0).ToArray();
+      r = eval.Evaluate(t, intervals, paramNodes, out lg, out ug);
+      Assert.AreEqual(0.0, r.LowerBound);
+      Assert.AreEqual(double.PositiveInfinity, r.UpperBound);
+
+      Assert.AreEqual(0.0, lg[0]); // x
+      Assert.AreEqual(double.NaN, ug[0]);
+
     }
 
@@ -126 +130 @@
       var eval = new IntervalEvaluator();
       var parser = new InfixExpressionParser();
+      var intervals = new Dictionary<string, Interval>() {
+        { "x", new Interval(3, 4) },
+        { "z", new Interval(1, 2) }
+      };
       var t = parser.Parse("cos(x)");
       var paramNodes = t.IterateNodesPostfix().Where(n => n.SubtreeCount == 0).ToArray();
-      var intervals = new Dictionary<string, Interval>() {
-        { "x", new Interval(3, 4) },
-      };
       var r = eval.Evaluate(t, intervals, paramNodes, out double[] lg, out double[] ug);
       Assert.AreEqual(-1, r.LowerBound); //  3..4 crosses pi and cos(pi) == -1
@@ -137 +142 @@
       Assert.AreEqual(0, lg[0]); // x
       Assert.AreEqual(-4 * Math.Sin(4), ug[0]);
+
+      t = parser.Parse("LOG(COS('z'))");
+      paramNodes = t.IterateNodesPostfix().Where(n => n.SubtreeCount == 0).ToArray();
+      r = eval.Evaluate(t, intervals, paramNodes, out  lg, out  ug);
+      Assert.AreEqual(double.NaN, r.LowerBound); 
+      Assert.AreEqual(Math.Log(Math.Cos(1)), r.UpperBound); 
+
+      Assert.AreEqual(-2 * Math.Sin(2) / Math.Cos(2), lg[0], 1e-5); // x
+      Assert.AreEqual(-1 * Math.Sin(1) / Math.Cos(1), ug[0], 1e-5);
+      
     }
 
@@ -182 +197 @@
       r = eval.Evaluate(t, intervals, paramNodes, out lg, out ug);
 
-      Assert.AreEqual(0.5 * Math.Sqrt(1e-10), lg[0], 1e-6); // z          --> lim x -> 0 (sqrt(x)) = 0
+      Assert.AreEqual(0.5 * Math.Sqrt(1e-10), lg[0], 1e-5); // --> lim x -> 0 (sqrt(x)) = 0
       Assert.AreEqual(0.5, ug[0], 1e-5);
+
+      t = parser.Parse("sqrt(y - z)"); // 1..2 - 0..1
+      paramNodes = t.IterateNodesPostfix().Where(n => n.SubtreeCount == 0).ToArray();
+      r = eval.Evaluate(t, intervals, paramNodes, out lg, out ug);
+      Assert.AreEqual(0, r.LowerBound);
+      Assert.AreEqual(Math.Sqrt(2), r.UpperBound);
+
+      Assert.AreEqual(double.PositiveInfinity, lg[0], 1e-5); // y
+      Assert.AreEqual(1/ Math.Sqrt(2)  , ug[0], 1e-5);
+      Assert.AreEqual(double.NegativeInfinity, lg[1], 1e-5); // z
+      Assert.AreEqual(0.0   , ug[1], 1e-5);
     }