Fix tests

Author: OlivierHnt

src/matmul.jl

@@ -50,7 +50,7 @@ Matrix multiplication mode type.
 
 Available mode types:
 - `:slow` (default): generic algorithm.
-- `:fast` : Rump's algorithm.
+- `:fast`: Rump's algorithm.
 """
 struct MatMulMode{T} end
 
@@ -122,18 +122,18 @@ end
 # Note: Rump's algorithm
 
 _mul!(::MatMulMode{:fast}, C, A::AbstractMatrix{<:Interval{<:Rational}}, B::AbstractVecOrMat{<:Interval{<:Rational}}, α, β) =
-    LinearAlgebra._mul!(C, A, B, α, β)
+    _mul!(MatMulMode{:slow}(), C, A, B, α, β)
 _mul!(::MatMulMode{:fast}, C, A::AbstractMatrix{<:Interval{<:Rational}}, B::AbstractVecOrMat, α, β) =
-    LinearAlgebra._mul!(C, A, B, α, β)
+    _mul!(MatMulMode{:slow}(), C, A, B, α, β)
 _mul!(::MatMulMode{:fast}, C, A::AbstractMatrix, B::AbstractVecOrMat{<:Interval{<:Rational}}, α, β) =
-    LinearAlgebra._mul!(C, A, B, α, β)
+    _mul!(MatMulMode{:slow}(), C, A, B, α, β)
 
 _mul!(::MatMulMode{:fast}, C, A::AbstractMatrix{<:Complex{<:Interval{<:Rational}}}, B::AbstractVecOrMat{<:Complex{<:Interval{<:Rational}}}, α, β) =
-    LinearAlgebra._mul!(C, A, B, α, β)
+    _mul!(MatMulMode{:slow}(), C, A, B, α, β)
 _mul!(::MatMulMode{:fast}, C, A::AbstractMatrix{<:Complex{<:Interval{<:Rational}}}, B::AbstractVecOrMat{<:Interval{<:Rational}}, α, β) =
-    LinearAlgebra._mul!(C, A, B, α, β)
+    _mul!(MatMulMode{:slow}(), C, A, B, α, β)
 _mul!(::MatMulMode{:fast}, C, A::AbstractMatrix{<:Interval{<:Rational}}, B::AbstractVecOrMat{<:Complex{<:Interval{<:Rational}}}, α, β) =
-    LinearAlgebra._mul!(C, A, B, α, β)
+    _mul!(MatMulMode{:slow}(), C, A, B, α, β)
 
 function _mul!(::MatMulMode{:fast}, C, A, B, α, β)
     Int == Int32 || return _fastmul!(C, A, B, α, β)

test/interval_tests/linearalgebra.jl

@@ -1,21 +1,23 @@
 import LinearAlgebra
 
 @testset "Matrix inversion" begin
+    IntervalArithmetic.configure(; matmul = :slow)
     A = [interval(2) interval(1, 2) ; interval(0) interval(1)]
     @test all(isequal_interval.(inv(A), [interval(0, 1) interval(-1.25, -0.25) ; interval(-0.5, 0.5) interval(0.5, 1.5)]))
     B = [interval(2) interval(1, 2) ; interval(0) interval(0, 1)]
     @test all(isnai, inv(B))
 end
 
 @testset "Matrix multiplication" begin
+    IntervalArithmetic.configure(; matmul = :fast)
     A = [interval(2, 4) interval(-2, 1) ; interval(-1, 2) interval(2, 4)]
     imA = interval(im) * A
 
-    @test all(isequal_interval.(A * A, [interval(-2, 19.5) interval(-16, 10) ; interval(-10, 16) interval(-2, 19.5)]))
-    @test all(isequal_interval.(A * mid.(A), [interval(5, 12.5) interval(-8, 2) ; interval(-2, 8) interval(5, 12.5)]))
-    @test all(isequal_interval.(mid.(A) * A, [interval(5, 12.5) interval(-8, 2) ; interval(-2, 8) interval(5, 12.5)]))
+    @test all(issubset_interval.([interval(0, 18) interval(-16, 8) ; interval(-8, 16) interval(0, 18)], A * A))
+    @test all(issubset_interval.([interval(5, 12.5) interval(-8, 2) ; interval(-2, 8) interval(5, 12.5)], A * mid.(A)))
+    @test all(issubset_interval.([interval(5, 12.5) interval(-8, 2) ; interval(-2, 8) interval(5, 12.5)], mid.(A) * A))
 
-    @test all(isequal_interval.(imA * imA, -[interval(-2, 19.5) interval(-16, 10) ; interval(-10, 16) interval(-2, 19.5)]))
-    @test all(isequal_interval.(mid.(A) * imA, interval(im)*[interval(5, 12.5) interval(-8, 2) ; interval(-2, 8) interval(5, 12.5)]))
-    @test all(isequal_interval.(imA * mid.(A), interval(im)*[interval(5, 12.5) interval(-8, 2) ; interval(-2, 8) interval(5, 12.5)]))
+    @test all(issubset_interval.([interval(-18, 0) interval(-8, 16) ; interval(-16, 8) interval(-18, 0)], imA * imA))
+    @test all(issubset_interval.(interval(im)*[interval(5, 12.5) interval(-8, 2) ; interval(-2, 8) interval(5, 12.5)], mid.(A) * imA))
+    @test all(issubset_interval.(interval(im)*[interval(5, 12.5) interval(-8, 2) ; interval(-2, 8) interval(5, 12.5)], imA * mid.(A)))
 end