Changeset 17348 for trunk

Timestamp:

11/04/19 14:40:48 (4 years ago)

Author:

chaider

Message:

#3032 Added analytical quotient for intervals

Location:

trunk

Files:

• HeuristicLab.Problems.DataAnalysis/3.4/Implementation/Interval.cs (modified) (1 diff)
• HeuristicLab.Tests/HeuristicLab.Problems.DataAnalysis-3.4/IntervalTest.cs (modified) (1 diff)

trunk/HeuristicLab.Problems.DataAnalysis/3.4/Implementation/Interval.cs

@@ -237 +237 @@
       return new Interval(Math.Min(absLower, absUpper), Math.Max(absLower, absUpper));
     }
+
+    public static Interval AnalyticalQuotient(Interval a, Interval b) {
+      if (b.LowerBound <= 0 && b.UpperBound >= 0) //Second interval goes over zero or is zero
+        return new Interval(a.LowerBound, a.UpperBound);
+
+      var values = new List<double> {
+        a.LowerBound / Math.Sqrt(Math.Pow(b.LowerBound, 2) + 1),
+        a.LowerBound / Math.Sqrt(Math.Pow(b.UpperBound, 2) + 1),
+        a.UpperBound / Math.Sqrt(Math.Pow(b.LowerBound, 2) + 1),
+        a.UpperBound / Math.Sqrt(Math.Pow(b.UpperBound, 2) + 1)
+      };
+
+      return new Interval(values.Min(), values.Max());
+    }
     #endregion
   }

trunk/HeuristicLab.Tests/HeuristicLab.Problems.DataAnalysis-3.4/IntervalTest.cs

@@ -182 +182 @@
       Assert.AreEqual(new Interval(0, 2), Interval.Absolute(new Interval(0, 2)));
     }
+
+    [TestMethod]
+    [TestCategory("Problems.DataAnalysis")]
+    [TestProperty("Time", "short")]
+    public void TestIntervalAnalyticalQuotient() {
+      //Second interval goes over zero
+      Assert.AreEqual(new Interval(-4, 2), Interval.AnalyticalQuotient(new Interval(-4, 2), new Interval(0, 5)));
+      Assert.AreEqual(new Interval(-4, 2), Interval.AnalyticalQuotient(new Interval(-4, 2), new Interval(-4, 0)));
+      Assert.AreEqual(new Interval(-4, 2), Interval.AnalyticalQuotient(new Interval(-4, 2), new Interval(-2, 4)));
+      Assert.AreEqual(new Interval(-4, -2), Interval.AnalyticalQuotient(new Interval(-4, -2), new Interval(0, 5)));
+      Assert.AreEqual(new Interval(-4, -2), Interval.AnalyticalQuotient(new Interval(-4, -2), new Interval(-4, 0)));
+      Assert.AreEqual(new Interval(-4, -2), Interval.AnalyticalQuotient(new Interval(-4, -2), new Interval(-2, 4)));
+      Assert.AreEqual(new Interval(3, 7), Interval.AnalyticalQuotient(new Interval(3, 7), new Interval(0, 5)));
+      Assert.AreEqual(new Interval(3, 7), Interval.AnalyticalQuotient(new Interval(3, 7), new Interval(-4, 0)));
+      Assert.AreEqual(new Interval(3, 7), Interval.AnalyticalQuotient(new Interval(3, 7), new Interval(-2, 4)));
+      Assert.AreEqual(new Interval(0, 0), Interval.AnalyticalQuotient(new Interval(0, 0), new Interval(0, 5)));
+      Assert.AreEqual(new Interval(-4, -4), Interval.AnalyticalQuotient(new Interval(-4, -4), new Interval(-4, 0)));
+      Assert.AreEqual(new Interval(4, 4), Interval.AnalyticalQuotient(new Interval(4, 4), new Interval(-2, 4)));
+
+      //First interval goes over zero
+      Assert.AreEqual(new Interval(-2.0/Math.Sqrt(5), 5.0 / Math.Sqrt(5)), Interval.AnalyticalQuotient(new Interval(-2, 5), new Interval(-5, -2)));
+      Assert.AreEqual(new Interval(-2.0 / Math.Sqrt(5), 5.0 / Math.Sqrt(5)), Interval.AnalyticalQuotient(new Interval(-2, 5), new Interval(2, 5)));
+      Assert.AreEqual(new Interval(2.0/Math.Sqrt(26), 5.0/Math.Sqrt(5)), Interval.AnalyticalQuotient(new Interval(2, 5), new Interval(-5, -2)));
+      Assert.AreEqual(new Interval(2.0 / Math.Sqrt(26), 5.0 / Math.Sqrt(5)), Interval.AnalyticalQuotient(new Interval(2, 5), new Interval(2, 5)));
+      Assert.AreEqual(new Interval(-5.0/Math.Sqrt(5), -2.0/Math.Sqrt(26)), Interval.AnalyticalQuotient(new Interval(-5, -2), new Interval(-5, -2)));
+      Assert.AreEqual(new Interval(-5.0 / Math.Sqrt(5), -2.0 / Math.Sqrt(26)), Interval.AnalyticalQuotient(new Interval(-5, -2), new Interval(2, 5)));
+    }
   }
 }