Changeset 17350 for trunk/HeuristicLab.Problems.DataAnalysis

Timestamp:

11/13/19 14:50:26 (4 years ago)

Author:

chaider

Message:

#3032

• Added IsAlmost method with an epsilon parameter
• Added check for floating point issues at interval constructor
• Updated absolute operator (jump over zero between lower and upper bound)
• Updated analytical quotient to use already implemented operators to build it
• Updated test cases for operators
• Refactored test method names

File:

• trunk/HeuristicLab.Problems.DataAnalysis/3.4/Implementation/Interval.cs (modified) (6 diffs)

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

@@ -37 +37 @@
     protected Interval(StorableConstructorFlag _) { }
 
+    /// <summary>
+    /// Creates an interval with given bounds, where lower bound must be smaller than
+    /// the upper bound. Floating point precision errors trough calculations are fixed by,
+    /// checking if the intervals are almost the same (E-12). If this is the case, the bounds 
+    /// will be set to the bound closer to zero.
+    /// </summary>
+    /// <param name="lowerBound">Lower bound of the interval</param>
+    /// <param name="upperBound">Upper bound of the interval</param>
     public Interval(double lowerBound, double upperBound) {
+      if (lowerBound.IsAlmost(upperBound)) {
+        //If the bounds go over zero
+        if (lowerBound <= 0 && upperBound >= 0) {
+          lowerBound = 0.0;
+          upperBound = 0.0;
+          //Interval is negative
+        } else if (upperBound < 0) {
+          lowerBound = upperBound;
+          //Interval is positive
+        } else {
+          upperBound = lowerBound;
+        }
+      }
+
       if (lowerBound > upperBound)
         throw new ArgumentException("LowerBound must be smaller than UpperBound.");
@@ -58 +80 @@
                 double.IsNaN(LowerBound) || double.IsNaN(UpperBound);
       }
+    }
+
+    /// <summary>
+    /// True if the interval is positive without zero
+    /// </summary>
+    public bool IsPositive {
+      get => LowerBound > 0.0; 
+    }
+
+    /// <summary>
+    /// True if the interval is negative without zero
+    /// </summary>
+    public bool IsNegative {
+      get => UpperBound < 0.0;
     }
 
@@ -128 +164 @@
     }
 
-    //mkommend: Division by intervals containing 0 is implemented as defined in
+    //Division by intervals containing 0 is implemented as defined in
     //http://en.wikipedia.org/wiki/Interval_arithmetic
     public static Interval Divide(Interval a, Interval b) {
@@ -171 +207 @@
     public static Interval Tangens(Interval a) {
       return Interval.Divide(Interval.Sine(a), Interval.Cosine(a));
-    }  
+    }
     public static Interval HyperbolicTangent(Interval a) {
       return new Interval(Math.Tanh(a.LowerBound), Math.Tanh(a.UpperBound));
@@ -206 +242 @@
       if (a.UpperBound <= 0) return new Interval(a.UpperBound * a.UpperBound, a.LowerBound * a.LowerBound);     // interval is negative
       else if (a.LowerBound >= 0) return new Interval(a.LowerBound * a.LowerBound, a.UpperBound * a.UpperBound); // interval is positive
-      else return new Interval(0, Math.Max(a.LowerBound*a.LowerBound, a.UpperBound*a.UpperBound)); // interval goes over zero
+      else return new Interval(0, Math.Max(a.LowerBound * a.LowerBound, a.UpperBound * a.UpperBound)); // interval goes over zero
     }
 
@@ -235 +271 @@
       var absLower = Math.Abs(a.LowerBound);
       var absUpper = Math.Abs(a.UpperBound);
-      return new Interval(Math.Min(absLower, absUpper), Math.Max(absLower, absUpper));
+      var min = Math.Min(absLower, absUpper);
+      var max = Math.Max(absLower, absUpper);
+
+      if (a.Contains(0.0)) {
+        min = 0.0;
+      }
+
+      return new Interval(min, max);
     }
 
     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());
+      var dividend = a;
+      var divisor = Add(Square(b), new Interval(1.0, 1.0));
+      divisor = SquareRoot(divisor);
+
+      var quotient = Divide(dividend, divisor);
+      return quotient;
     }
     #endregion