SemiseparableMatrices v0.3

Author: dlfivefifty

Project.toml

@@ -27,7 +27,7 @@ InfiniteArrays = "0.12"
 LazyArrays = "0.21.8, 0.22"
 LazyBandedMatrices = "0.7.1"
 MatrixFactorizations = "0.8"
-SemiseparableMatrices = "0.2.7"
+SemiseparableMatrices = "0.2.7, 0.3"
 julia = "1.6"
 
 [extras]

test/test_infbanded.jl

@@ -1,4 +1,4 @@
-using InfiniteLinearAlgebra, InfiniteArrays, BandedMatrices, FillArrays, LazyBandedMatrices, Test
+using InfiniteLinearAlgebra, InfiniteArrays, BandedMatrices, FillArrays, LazyBandedMatrices, LazyArrays, Test
 import BandedMatrices: _BandedMatrix
 
 @testset "∞-banded" begin
@@ -140,4 +140,13 @@ import BandedMatrices: _BandedMatrix
         A, B = Diagonal(Fill(2,∞)) , LazyBandedMatrices.Bidiagonal(exp.(1:∞), exp.(1:∞), :L)
         @test (A*B)[1:10,1:10] ≈ (B*A)[1:10,1:10] ≈ 2B[1:10,1:10]
     end
+
+    @testset "concat" begin
+        H = ApplyArray(hvcat, 2, 1, [1 Zeros(1,∞)], [1; Zeros(∞)], Diagonal(1:∞))
+        @test bandwidths(H) == (1,1)
+        H = ApplyArray(hvcat, 2, 1, [0 Zeros(1,∞)], [0; Zeros(∞)], Diagonal(1:∞))
+        @test bandwidths(H) == (0,0)
+        H = ApplyArray(hvcat, (2,2), 1, [1 Zeros(1,∞)], [1; Zeros(∞)], Diagonal(1:∞))
+        @test_broken bandwidths(H) == (1,1)
+    end
 end