Update runtests.jl

JuliaLinearAlgebra · Nov 29, 2024 · fa0f389 · fa0f389
1 parent a6e64c5
commit fa0f389
Showing 1 changed file with 0 additions and 111 deletions.
diff --git a/test/runtests.jl b/test/runtests.jl
@@ -284,56 +284,6 @@ end
 
 
 @testset "Algebra" begin
-    @testset "BandedMatrix" begin
-        A = BandedMatrix(-3 => Fill(7 / 10, ∞), -2 => 1:∞, 1 => Fill(2im, ∞))
-        @test A isa BandedMatrix{ComplexF64}
-        @test A[1:10, 1:10] == diagm(-3 => Fill(7 / 10, 7), -2 => 1:8, 1 => Fill(2im, 9))
-
-        A = BandedMatrix(0 => Vcat([1, 2, 3], Zeros(∞)), 1 => Vcat(1, Zeros(∞)))
-        @test A[1, 2] == 1
-
-        A = BandedMatrix(-3 => Fill(7 / 10, ∞), -2 => Fill(1, ∞), 1 => Fill(2im, ∞))
-        Ac = BandedMatrix(A')
-        At = BandedMatrix(transpose(A))
-        @test Ac[1:10, 1:10] ≈ (A')[1:10, 1:10] ≈ A[1:10, 1:10]'
-        @test At[1:10, 1:10] ≈ transpose(A)[1:10, 1:10] ≈ transpose(A[1:10, 1:10])
-
-        A = BandedMatrix(-1 => Vcat(Float64[], Fill(1 / 4, ∞)), 0 => Vcat([1.0 + im], Fill(0, ∞)), 1 => Vcat(Float64[], Fill(1, ∞)))
-        @test MemoryLayout(typeof(view(A.data, :, 1:10))) == ApplyLayout{typeof(hcat)}()
-        Ac = BandedMatrix(A')
-        At = BandedMatrix(transpose(A))
-        @test Ac[1:10, 1:10] ≈ (A')[1:10, 1:10] ≈ A[1:10, 1:10]'
-        @test At[1:10, 1:10] ≈ transpose(A)[1:10, 1:10] ≈ transpose(A[1:10, 1:10])
-
-        A = BandedMatrix(-2 => Vcat(Float64[], Fill(1 / 4, ∞)), 0 => Vcat([1.0 + im, 2, 3], Fill(0, ∞)), 1 => Vcat(Float64[], Fill(1, ∞)))
-        Ac = BandedMatrix(A')
-        At = BandedMatrix(transpose(A))
-        @test Ac[1:10, 1:10] ≈ (A')[1:10, 1:10] ≈ A[1:10, 1:10]'
-        @test At[1:10, 1:10] ≈ transpose(A)[1:10, 1:10] ≈ transpose(A[1:10, 1:10])
-
-        A = _BandedMatrix(Fill(1, 4, ∞), ℵ₀, 1, 2)
-        @test A^2 isa BandedMatrix
-        @test (A^2)[1:10, 1:10] == (A*A)[1:10, 1:10] == (A[1:100, 1:100]^2)[1:10, 1:10]
-        @test A^3 isa ApplyMatrix{<:Any,typeof(*)}
-        @test (A^3)[1:10, 1:10] == (A*A*A)[1:10, 1:10] == ((A*A)*A)[1:10, 1:10] == (A*(A*A))[1:10, 1:10] == (A[1:100, 1:100]^3)[1:10, 1:10]
-
-        @testset "∞ x finite" begin
-            A = BandedMatrix(1 => 1:∞) + BandedMatrix(-1 => Fill(2, ∞))
-            B = _BandedMatrix(randn(3, 5), ℵ₀, 1, 1)
-
-            @test lmul!(2.0, copy(B)')[:, 1:10] == (2B')[:, 1:10]
-
-            @test_throws ArgumentError BandedMatrix(A)
-            @test A * B isa MulMatrix
-            @test B'A isa MulMatrix
-
-            @test all(diag(A[1:6, 1:6]) .=== zeros(6))
-
-            @test (A*B)[1:7, 1:5] ≈ A[1:7, 1:6] * B[1:6, 1:5]
-            @test (B'A)[1:5, 1:7] ≈ (B')[1:5, 1:6] * A[1:6, 1:7]
-        end
-    end
-
     @testset "BlockTridiagonal" begin
         A = BlockTridiagonal(Vcat([fill(1.0, 2, 1), Matrix(1.0I, 2, 2), Matrix(1.0I, 2, 2), Matrix(1.0I, 2, 2)], Fill(Matrix(1.0I, 2, 2), ∞)),
             Vcat([zeros(1, 1)], Fill(zeros(2, 2), ∞)),
@@ -369,67 +319,6 @@ end
         @test_broken B[Block(1, 3)] == Zeros(2, 2)
     end
 
-    @testset "Fill" begin
-        A = _BandedMatrix(Ones(1, ∞), ℵ₀, -1, 1)
-        @test 1.0 .* A isa BandedMatrix{Float64,<:Fill}
-        @test Zeros(∞) .* A ≡ Zeros(∞, ∞) .* A ≡ A .* Zeros(1, ∞) ≡ A .* Zeros(∞, ∞) ≡ Zeros(∞, ∞)
-        @test Ones(∞) .* A isa BandedMatrix{Float64,<:Ones}
-        @test A .* Ones(1, ∞) isa BandedMatrix{Float64,<:Ones}
-        @test 2.0 .* A isa BandedMatrix{Float64,<:Fill}
-        @test A .* 2.0 isa BandedMatrix{Float64,<:Fill}
-        @test Eye(∞) * A isa BandedMatrix{Float64,<:Ones}
-        @test A * Eye(∞) isa BandedMatrix{Float64,<:Ones}
-
-        @test A * A isa BandedMatrix
-        @test (A*A)[1:10, 1:10] == BandedMatrix(2 => Ones(8))
-
-        Ã = _BandedMatrix(Fill(1, 1, ∞), ℵ₀, -1, 1)
-        @test A * Ã isa BandedMatrix
-        @test Ã * A isa BandedMatrix
-        @test Ã * Ã isa BandedMatrix
-
-        B = _BandedMatrix(Ones(1, 10), ℵ₀, -1, 1)
-        C = _BandedMatrix(Ones(1, 10), 10, -1, 1)
-        D = _BandedMatrix(Ones(1, ∞), 10, -1, 1)
-
-        @test (A*B)[1:10, 1:10] == (B*C)[1:10, 1:10] == (D*A)[1:10, 1:10] == D * B == (C*D)[1:10, 1:10] == BandedMatrix(2 => Ones(8))
-    end
-
-    @testset "Banded Broadcast" begin
-        A = _BandedMatrix((1:∞)', ℵ₀, -1, 1)
-        @test 2.0 .* A isa BandedMatrix{Float64,<:Adjoint}
-        @test A .* 2.0 isa BandedMatrix{Float64,<:Adjoint}
-        @test Eye(∞) * A isa BandedMatrix{Float64,<:Adjoint}
-        @test A * Eye(∞) isa BandedMatrix{Float64,<:Adjoint}
-        A = _BandedMatrix(Vcat((1:∞)', Ones(1, ∞)), ℵ₀, 0, 1)
-        @test 2.0 .* A isa BandedMatrix
-        @test A .* 2.0 isa BandedMatrix
-        @test Eye(∞) * A isa BandedMatrix
-        @test A * Eye(∞) isa BandedMatrix
-        b = 1:∞
-        @test bandwidths(b .* A) == (0, 1)
-
-        @test colsupport(b .* A, 1) == 1:1
-        @test Base.replace_in_print_matrix(b .* A, 2, 1, "0.0") == " ⋅ "
-        @test bandwidths(A .* b) == (0, 1)
-        @test A .* b' isa BroadcastArray
-        @test bandwidths(A .* b') == bandwidths(A .* b')
-        @test colsupport(A .* b', 3) == 2:3
-
-        A = _BandedMatrix(Ones{Int}(1, ∞), ℵ₀, 0, 0)'
-        B = _BandedMatrix((-2:-2:-∞)', ℵ₀, -1, 1)
-        C = Diagonal(2 ./ (1:2:∞))
-        @test bandwidths(A * (B * C)) == (-1, 1)
-        @test bandwidths((A * B) * C) == (-1, 1)
-
-        A = _BandedMatrix(Ones{Int}(1, ∞), ℵ₀, 0, 0)'
-        B = _BandedMatrix((-2:-2:-∞)', ℵ₀, -1, 1)
-        @test MemoryLayout(A + B) isa BroadcastBandedLayout{typeof(+)}
-        @test MemoryLayout(2 * (A + B)) isa BroadcastBandedLayout{typeof(*)}
-        @test bandwidths(A + B) == (0, 1)
-        @test bandwidths(2 * (A + B)) == (0, 1)
-    end
-
     @testset "Triangle OP recurrences" begin
         k = mortar(Base.OneTo.(1:∞))
         n = mortar(Fill.(1:∞, 1:∞))