Changeset 17314 for branches/2994-AutoDiffForIntervals/HeuristicLab.Tests/HeuristicLab.Problems.DataAnalysis-3.4/AutoDiffIntervalTest.cs

Timestamp:

10/04/19 10:40:51 (5 years ago)

Author:

gkronber

Message:

#2994: fix bug in calculation of cos() and sin()

File:

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

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

@@ -142 +142 @@
       AssertAreEqualInterval(Cosine(new AlgebraicInterval(0, 2 * Math.PI)), new AlgebraicInterval(-1, 1));
       AssertAreEqualInterval(new AlgebraicInterval(-1, 1), Cosine(new AlgebraicInterval(Math.PI, 4 * Math.PI / 2)));
+
+      // some random expression that were problematic
+      AssertAreEqualInterval(new AlgebraicInterval(-0.0326026961803847, 1), new AlgebraicInterval(1, 2).Scale(1.73269159006613).Add(9.23027422279119).Cos());
+      AssertAreEqualInterval(new AlgebraicInterval(0.54030230586814, 1), new AlgebraicInterval(1, 2).Scale(3.39318455754522).Cos().IntPower(2).Cos());
+      AssertAreEqualInterval(new AlgebraicInterval(-0.825841274068331, 1), new AlgebraicInterval(1, 2).Scale(-16.7140105366707 * -0.223806041028346).Cos());
+      AssertAreEqualInterval(new AlgebraicInterval(0, 1), new AlgebraicInterval(1, 2).Scale(1.03513653909418).Exp().Cos().IntPower(2));
+      AssertAreEqualInterval(new AlgebraicInterval(-1, 1), new AlgebraicInterval(1, 2).Scale(0.535386838272755).Exp().IntPower(2).Cos());
     }
 
@@ -285 +292 @@
       };
       var r = eval.Evaluate(t, intervals, paramNodes, out double[] lg, out double[] ug);
-      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'))
+      // 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'))
     }
 
@@ -411 +419 @@
       Assert.AreEqual(3, r.UpperBound);
 
-      Assert.AreEqual(0.0, lg[0]); // x
-      Assert.AreEqual(0.0, ug[0]); XXXX
-
-      t = parser.Parse("sqrt(y)");
+      Assert.AreEqual(2.0 / 3.0, lg[0]); // x
+      Assert.AreEqual(1.0, ug[0]); 
+
+      t = parser.Parse("cuberoot(y)");
       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(0.0, r.UpperBound);
-
-      Assert.AreEqual(0.0, lg[0]); // y
-      Assert.AreEqual(0.0, ug[0], 1e-5);
+      Assert.AreEqual(Math.Pow(1, 1.0/3.0), r.LowerBound);
+      Assert.AreEqual(Math.Pow(2, 1.0 / 3.0), r.UpperBound);
+
+      Assert.AreEqual(1.0 / 3.0, lg[0]); // y
+      Assert.AreEqual(1.0 / 3.0 * Math.Pow(2, 1.0/3.0), ug[0], 1e-5);
 
       t = parser.Parse("sqrt(z)");
@@ -427 +435 @@
       r = eval.Evaluate(t, intervals, paramNodes, out lg, out ug);
       Assert.AreEqual(0.0, r.LowerBound);
-      Assert.AreEqual(0.0, r.UpperBound);
+      Assert.AreEqual(1.0, r.UpperBound);
 
       Assert.AreEqual(0.0, lg[0]); // z
-      Assert.AreEqual(0.0, ug[0], 1e-5);
+      Assert.AreEqual(1.0/3.0, ug[0], 1e-5);
     }