Merge branch 'dev'

Author: mihakralj

lib/averages/Mama.cs

@@ -111,7 +111,7 @@ private static double CalculateHilbertTransform(CircularBuffer buffer, double ad
     [MethodImpl(MethodImplOptions.AggressiveInlining)]
     private double CalculatePeriod(double im, double re)
     {
-        if (im == 0 || re == 0) return _pd[^2];
+        if (System.Math.Abs(im) <= double.Epsilon || System.Math.Abs(re) <= double.Epsilon) return _pd[^2];
         return _twoPi / System.Math.Atan(im / re);
     }
 
@@ -167,7 +167,7 @@ protected override double Calculation()
             _pd[^1] = AdjustPeriod(_pd[^1]);
 
             // Phase calculation
-            double phase = _i1[^1] != 0 ? System.Math.Atan(_q1[^1] / _i1[^1]) * _radToDeg : _ph[^2];
+            double phase = Math.Abs(_i1[^1]) >= double.Epsilon ? System.Math.Atan(_q1[^1] / _i1[^1]) * _radToDeg : _ph[^2];
             _ph.Add(phase, Input.IsNew);
 
             // Adaptive alpha

lib/errors/Mase.cs

@@ -120,7 +120,7 @@ private double CalculateMase()
         double sumAbsoluteError = CalculateSumAbsoluteError(actualValues, predictedValues);
         double naiveForecastError = CalculateNaiveForecastError(actualValues, naiveValues);
 
-        return naiveForecastError != 0 ? (sumAbsoluteError / _actualBuffer.Count) / naiveForecastError : double.PositiveInfinity;
+        return naiveForecastError >= double.Epsilon ? (sumAbsoluteError / _actualBuffer.Count) / naiveForecastError : double.PositiveInfinity;
     }
 
     /// <summary>