Overload blocbanded_col/rowsupport for infinite case (#155)

Author: dlfivefifty

Project.toml

@@ -1,6 +1,6 @@
 name = "InfiniteLinearAlgebra"
 uuid = "cde9dba0-b1de-11e9-2c62-0bab9446c55c"
-version = "0.7.2"
+version = "0.7.3"
 
 [deps]
 ArrayLayouts = "4c555306-a7a7-4459-81d9-ec55ddd5c99a"

src/blockbanded/blockbanded.jl

@@ -12,6 +12,9 @@ function BlockArrays.sortedunion(b, ::AbstractVector{<:PosInfinity})
 end
 BlockArrays.sortedunion(a::OneToInfCumsum, ::OneToInfCumsum) = a
 
+BlockBandedMatrices.blockbanded_colstop(A, x::InfiniteCardinal{0}) = x
+BlockBandedMatrices.blockbanded_rowstop(A, x::InfiniteCardinal{0}) = x
+
 
 
 function BlockArrays.sortedunion(a::Vcat{Int,1,<:Tuple{Union{Int,AbstractVector{Int}},<:AbstractRange}},

test/runtests.jl

@@ -231,6 +231,8 @@ include("test_infbanded.jl")
             Dy = BlockBandedMatrices._BandedBlockBandedMatrix((k .+ (b + c))', axes(k, 1), (-1, 1), (-1, 1))
             N = 100
             @test Dy[Block.(1:N), Block.(1:N)] == BlockBandedMatrices._BandedBlockBandedMatrix((k.+(b+c))[Block.(1:N)]', axes(k, 1)[Block.(1:N)], (-1, 1), (-1, 1))
+            @test colsupport(Dy, axes(Dy,2)) == 1:∞
+            @test rowsupport(Dy, axes(Dy,1)) == 2:∞
         end
 
         @testset "Symmetric" begin