Fix rounding in ^

JuliaIntervals · May 22, 2019 · 49c8ee4 · 49c8ee4
1 parent da493b5
commit 49c8ee4
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/intervals/functions.jl b/src/intervals/functions.jl
@@ -211,11 +211,11 @@ function ^(x::Interval, n::Integer)  # fast integer power
                         power_by_squaring(mag(x), n, RoundUp))
 
         elseif x.lo > 0
-            return Interval(power_by_squaring(x.lo, n, RoundUp),
+            return Interval(power_by_squaring(x.lo, n, RoundDown),
                             power_by_squaring(x.hi, n, RoundUp))
 
         else  # x.lo < x.hi < 0
-            return Interval(power_by_squaring(-x.hi, n, RoundUp),
+            return Interval(power_by_squaring(-x.hi, n, RoundDown),
                             power_by_squaring(-x.lo, n, RoundUp))
         end