Merge branch 'master' into 1.0-dev

Author: Kolaru

.github/workflows/CI.yml

@@ -11,6 +11,7 @@ jobs:
       fail-fast: false
       matrix:
         version:
+          - '1.6'
           - '1'
           - 'nightly'
         os:

src/IntervalArithmetic.jl

@@ -1,7 +1,5 @@
 # This file is part of the IntervalArithmetic.jl package; MIT licensed
 
-__precompile__(true)
-
 module IntervalArithmetic
 
 import CRlibm
@@ -68,12 +66,7 @@ export
     pow, extended_div, nthroot,
     setformat
 
-if VERSION >= v"1.5.0-DEV.124"
-    import Base: isdisjoint
-else
-    export isdisjoint
-end
-
+import Base: isdisjoint
 
 export
     setindex   # re-export from StaticArrays for IntervalBox
@@ -124,4 +117,15 @@ include("plot_recipes/plot_recipes.jl")
 """
 const Region{T} = Union{Interval{T}, IntervalBox{T}}
 
+
+# These definitions has been put there because generated functions must be
+# defined after all methods they use.
+@generated function Interval{T}(x::AbstractIrrational) where T
+    res = atomic(Interval{T}, x())  # Precompute the interval
+    return :(return $res)  # Set body of the function to return the precomputed result
+end
+
+Interval{BigFloat}(x::AbstractIrrational) = atomic(Interval{BigFloat}, x)
+Interval(x::AbstractIrrational) = Interval{Float64}(x)
+
 end # module IntervalArithmetic

src/multidim/intervalbox.jl

@@ -18,14 +18,8 @@ IntervalBox(x) = IntervalBox(x...)
 IntervalBox(X::IntervalBox, n) = foldl(×, Iterators.repeated(X, n))
 
 # construct from two vectors giving bottom and top corners:
-function IntervalBox(los::AbstractVector, his::AbstractVector)
-    xs = map(los, his) do lo, hi
-        # Allow the creation of the intervals even when the bounds are in the
-        # wrong order
-        return lo <= hi ? interval(lo, hi) : interval(hi, lo)
-    end
-    return IntervalBox(xs...)
-end
+IntervalBox(lo::AbstractVector, hi::AbstractVector) = IntervalBox(force_interval.(lo, hi))
+IntervalBox(lo::SVector{N,T}, hi::SVector{N,T}) where {N,T} = IntervalBox(force_interval.(lo, hi))
 
 Base.@propagate_inbounds Base.getindex(X::IntervalBox, i) = X.v[i]
 

test/interval_tests/consistency.jl

@@ -264,8 +264,6 @@ import IntervalArithmetic: unsafe_interval
         @test isnan(mid(emptyinterval()))
         @test mid(entireinterval()) == 0.0
         @test isnan(mid(nai()))
-        # @test_throws InexactError nai(interval(1//2)) TODO move this test
-
         @test mid(2.125) == 2.125
         @test radius(2.125) == 0
     end

test/interval_tests/numeric.jl

@@ -329,10 +329,10 @@ end
     end
 
     @testset "Float32 intervals" begin
-
         a = interval(Float32, 1e38)
         b = interval(Float32, 1e2)
         @test a * b ≛ interval(Float32, floatmax(Float32), Inf)
+        @test Interval(1.0f0) ^ Interval(1.0f0) ≛ Interval(1.0f0) # test for PR #482
     end
 
 

test/runtests.jl

@@ -23,5 +23,7 @@ include_test("display_tests/display.jl")
 # Multidim tests
 include_test("multidim_tests/multidim.jl")
 
+include("multidim_tests/multidim.jl")
+
 # ITF1788 tests
-include_test("test_ITF1788/run_ITF1788.jl")  # TODO fix these tests
+include("test_ITF1788/run_ITF1788.jl")

test/test_ITF1788/abs_rev.jl

@@ -47,5 +47,4 @@
     @test abs_rev(interval(-Inf, -1.0), entireinterval())[2] === emptyinterval()
 
     @test abs_rev(interval(-Inf, 1.0), entireinterval())[2] === interval(-1.0, 1.0)
-
 end

test/test_ITF1788/atan2.jl

@@ -75,5 +75,4 @@
     @test atan(interval(0.0, 5.0), interval(0.0, 5.0)) === interval(0x0p0, 0x1.921FB54442D19p0)
 
     @test atan(interval(-5.0, 0.0), interval(0.0, 5.0)) === interval(-0x1.921FB54442D19p0, 0x0p0)
-
 end

test/test_ITF1788/c-xsc.jl

@@ -11,7 +11,6 @@
     @test -(interval(10.0, 20.0)) === interval(-20.0, -10.0)
 
     @test +(interval(10.0, 20.0)) === interval(10.0, 20.0)
-
 end
 
 @testset "cxsc.intervalmuldiv" begin
@@ -77,7 +76,6 @@ end
     @test /(interval(-2.0, 1.0), interval(-8.0, -4.0)) === interval(-0.25, 0.5)
 
     @test /(interval(-2.0, -1.0), interval(-8.0, -4.0)) === interval(0.125, 0.5)
-
 end
 
 @testset "cxsc.intervalsetops" begin