Revert "Merge branch 'master' of https://github.com/JuliaMatrices/BlockBandedMatrices.jl"

This reverts commit f7b2f37, reversing
changes made to 4137e6b.

Author: dlfivefifty

Project.toml

@@ -1,6 +1,6 @@
 name = "BlockBandedMatrices"
 uuid = "ffab5731-97b5-5995-9138-79e8c1846df0"
-version = "0.8.2"
+version = "0.8.1"
 
 [deps]
 ArrayLayouts = "4c555306-a7a7-4459-81d9-ec55ddd5c99a"

src/AbstractBlockBandedMatrix.jl

@@ -1,24 +1,14 @@
+####
+# Matrix memory layout traits
+#
+# if MemoryLayout(A) returns BandedColumnMajor, you must override
+# pointer and leadingdimension
+# in addition to the banded matrix interface
+####
 
-"""
-   AbstractBlockBandedLayout
-
-isa a `MemoryLayout` that indicates that the array implements the block-banded
-interface.
-"""
 abstract type AbstractBlockBandedLayout <: AbstractBlockLayout end
-
-"""
-   AbstractBandedBlockBandedLayout
-
-isa a `MemoryLayout` that indicates that the array implements the banded-block-banded
-interface.
-"""
 abstract type AbstractBandedBlockBandedLayout <: AbstractBlockBandedLayout end
 
-
-struct BandedBlockBandedLayout <: AbstractBandedBlockBandedLayout end
-struct BlockBandedLayout <: AbstractBlockBandedLayout end
-
 struct BandedBlockBandedColumns{LAY} <: AbstractBandedBlockBandedLayout end
 struct BandedBlockBandedRows{LAY} <: AbstractBandedBlockBandedLayout end
 struct BlockBandedColumns{LAY} <: AbstractBlockBandedLayout end

src/BandedBlockBandedMatrix.jl

@@ -318,18 +318,13 @@ end
 end
 
 ## structured matrix methods ##
-function _bandedblockbanded_replace_in_print_matrix(A, i, j, s)
+function Base.replace_in_print_matrix(A::BandedBlockBandedMatrix, i::Integer, j::Integer, s::AbstractString)
     bi = findblockindex.(axes(A), (i,j))
     I,J = block.(bi)
     i,j = blockindex.(bi)
-    l,u = blockbandwidths(A)
-    λ,μ = subblockbandwidths(A)
-    -l ≤ Int(J-I) ≤ u && -λ ≤ j-i ≤ μ ? s : Base.replace_with_centered_mark(s)
+    -A.l ≤ Int(J-I) ≤ A.u && -A.λ ≤ j-i ≤ A.μ ? s : Base.replace_with_centered_mark(s)
 end
 
-Base.replace_in_print_matrix(A::BandedBlockBandedMatrix, i::Integer, j::Integer, s::AbstractString) = 
-    _bandedblockbanded_replace_in_print_matrix(A, i, j, s)
-
 
 ############
 # Indexing #
@@ -418,43 +413,32 @@ sublayout(::BandedBlockBandedColumnMajor, ::Type{<:Tuple{BlockSlice{Block1},Bloc
 sublayout(::BandedBlockBandedColumnMajor, ::Type{<:Tuple{BlockSlice{BlockRange1},BlockSlice{Block1}}}) = BandedBlockBandedColumnMajor()
 sublayout(::BandedBlockBandedColumnMajor, ::Type{<:Tuple{BlockSlice{BlockRange1},BlockSlice{BlockIndexRange1}}}) = BandedBlockBandedColumnMajor()
 
-sublayout(::AbstractBandedBlockBandedLayout, ::Type{<:Tuple{BlockSlice{Block1},BlockSlice{Block1}}}) = BandedLayout()
-sublayout(::AbstractBandedBlockBandedLayout, ::Type{<:Tuple{BlockSlice{BlockIndexRange1},BlockSlice{BlockIndexRange1}}}) = BandedLayout()
-sublayout(::AbstractBandedBlockBandedLayout, ::Type{<:Tuple{BlockSlice{BlockIndexRange1},BlockSlice{Block1}}}) = BandedLayout()
-sublayout(::AbstractBandedBlockBandedLayout, ::Type{<:Tuple{BlockSlice{Block1},BlockSlice{BlockIndexRange1}}}) = BandedLayout()
-sublayout(::AbstractBandedBlockBandedLayout, ::Type{<:Tuple{BlockSlice{BlockRange1},BlockSlice{BlockRange1}}}) = BandedBlockBandedLayout()
-sublayout(::AbstractBandedBlockBandedLayout, ::Type{<:Tuple{BlockSlice{Block1},BlockSlice{BlockRange1}}}) = BandedBlockBandedLayout()
-sublayout(::AbstractBandedBlockBandedLayout, ::Type{<:Tuple{BlockSlice{BlockRange1},BlockSlice{Block1}}}) = BandedBlockBandedLayout()
-sublayout(::AbstractBandedBlockBandedLayout, ::Type{<:Tuple{BlockSlice{BlockRange1},BlockSlice{BlockIndexRange1}}}) = BandedBlockBandedLayout()
-
-
-sub_materialize(::AbstractBandedBlockBandedLayout, V, _) = BandedBlockBandedMatrix(V)
-
 isbanded(A::SubArray{<:Any,2,<:BandedBlockBandedMatrix}) = MemoryLayout(typeof(A)) == BandedColumnMajor()
 isbandedblockbanded(A::SubArray{<:Any,2,<:BandedBlockBandedMatrix}) = MemoryLayout(typeof(A)) == BandedBlockBandedColumnMajor()
 
 
-subblockbandwidths(V::SubArray{<:Any,2,<:Any,<:Tuple{<:BlockSlice{BlockRange1},<:BlockSlice{BlockRange1}}}) = 
-    subblockbandwidths(parent(V))
+subblockbandwidths(V::SubBandedBlockBandedMatrix) = subblockbandwidths(parent(V))
 
-function blockbandwidths(V::SubArray{<:Any,2,<:Any,<:Tuple{<:BlockSlice{BlockRange1},<:BlockSlice{Block1}}})
+function blockbandwidths(V::SubBandedBlockBandedMatrix{<:Any,BlockRange1,Block1})
     A = parent(V)
+
     KR = parentindices(V)[1].block.indices[1]
     J = parentindices(V)[2].block
     shift = Int(KR[1])-Int(J)
     blockbandwidth(A,1) - shift, blockbandwidth(A,2) + shift
 end
 
-function blockbandwidths(V::SubArray{<:Any,2,<:Any,<:Tuple{<:BlockSlice{<:Block1},<:BlockSlice{<:BlockRange1}}})
+function blockbandwidths(V::SubBandedBlockBandedMatrix{<:Any,Block1,BlockRange1})
     A = parent(V)
+
     K = parentindices(V)[1].block
     JR = parentindices(V)[2].block.indices[1]
     shift = Int(K)-Int(JR[1])
-    l,u = blockbandwidths(A)
-    l - shift, u + shift
+
+    blockbandwidth(A,1) - shift, blockbandwidth(A,2) + shift
 end
 
-function blockbandwidths(V::SubArray{<:Any,2,<:Any,<:Tuple{<:BlockSlice{<:BlockRange1},<:BlockSlice{<:BlockRange1}}})
+function blockbandwidths(V::SubBandedBlockBandedMatrix{<:Any,BlockRange1,BlockRange1})
     A = parent(V)
 
     KR = parentindices(V)[1].block.indices[1]

src/BlockBandedMatrices.jl

@@ -23,7 +23,7 @@ import SparseArrays: sparse
 
 import ArrayLayouts: BlasMatLmulVec,
                     triangularlayout, UpperTriangularLayout, TriangularLayout, MatLdivVec,
-                    triangulardata, sublayout, sub_materialize,
+                    triangulardata, sublayout, 
                     AbstractColumnMajor, DenseColumnMajor, ColumnMajor,
                     DiagonalLayout, MulAdd, mul, colsupport, rowsupport,
                     _qr, _factorize, _copyto!
@@ -38,7 +38,7 @@ import BandedMatrices: isbanded, bandwidths, bandwidth, banded_getindex, colrang
                         inbands_setindex!, inbands_getindex, banded_setindex!,
                         banded_generic_axpy!,
                         BlasFloat, banded_dense_axpy!, MemoryLayout,
-                        BandedLayout, BandedColumnMajor,
+                        BandedColumnMajor,
                         BandedSubBandedMatrix, bandeddata, tribandeddata,
                         _BandedMatrix, colstart, colstop, rowstart, rowstop,
                         BandedStyle, _fill_lmul!, bandshift,
@@ -50,7 +50,7 @@ export BandedBlockBandedMatrix, BlockBandedMatrix, BlockSkylineMatrix, blockband
 
 const Block1 = Block{1,Int}
 const BlockRange1 = BlockRange{1,Tuple{UnitRange{Int}}}
-const BlockIndexRange1 = BlockIndexRange{1,Tuple{UnitRange{Int}}}
+const BlockIndexRange1 = BlockIndexRange{1,Tuple{UnitRange{Int}}}    
 
 blockcolrange(A...) = blockcolsupport(A...)
 blockrowrange(A...) = blockrowsupport(A...)

src/linalg.jl

@@ -93,10 +93,6 @@ sublayout(::BlockBandedColumnMajor, ::Type{<:Tuple{<:BlockSlice{BlockRange1}, <:
 sublayout(::BlockBandedColumnMajor, ::Type{<:Tuple{<:BlockSlice{BlockRange1}, <:BlockSlice{BlockIndexRange1}}}) = ColumnMajor()
 sublayout(::BlockBandedColumnMajor, ::Type{<:Tuple{<:BlockSlice{BlockIndexRange1}, <:BlockSlice{BlockIndexRange1}}}) = ColumnMajor()
 
-isblockbanded(V::SubArray{<:Any,2,<:Any,<:Tuple{<:BlockSlice{BlockRange1}, <:BlockSlice{BlockRange1}}}) =
-    isblockbanded(parent(V))
-
-sub_materialize(::AbstractBlockBandedLayout, V, _) = BlockBandedMatrix(V)
 
 strides(V::SubBlockSkylineMatrix{<:Any,LL,UU,<:Union{BlockRange1,Block1},Block1}) where {LL,UU} =
     (1,parent(V).block_sizes.block_strides[Int(parentindices(V)[2].block)])
@@ -138,3 +134,5 @@ strides(V::SubBlockSkylineMatrix{T,LL,UU,BlockIndexRange1,BlockIndexRange1}) whe
     (1,parent(V).block_sizes.block_strides[Int(parentindices(V)[2].block.block)])
 
 MemoryLayout(V::SubBlockSkylineMatrix{T,LL,UU,BlockIndexRange1,BlockIndexRange1}) where {T,LL,UU} = ColumnMajor()
+
+

test/test_bandedblockbanded.jl

@@ -298,9 +298,7 @@ import BlockBandedMatrices: _BandedBlockBandedMatrix, blockcolsupport, blockrows
         @test B isa BandedBlockBandedMatrix
         @test blockbandwidths(V) == blockbandwidths(B) == (2,0)
         @test subblockbandwidths(V) == subblockbandwidths(B) == subblockbandwidths(A) == (λ,μ)
-        @test B == V == A[Block.(2:3), Block.(3:4)] 
-
-        @test A[Block.(2:3), Block.(3:4)] isa BandedBlockBandedMatrix
+        @test B == V
 
         x = randn(size(B,2))
         y = similar(x, size(B,1))

test/test_blockbanded.jl

@@ -11,7 +11,7 @@ import Base.Broadcast: materialize!
         cols = rows = 1:N
 
         @test Matrix(BlockBandedMatrix(Zeros(sum(rows),sum(cols)), rows, cols, (l,u))) ==
-            Array(BlockBandedMatrix(Zeros(sum(rows),sum(cols)), rows, cols, (l,u))) ==
+            Array(BlockBandedMatrix(Zeros(sum(rows),sum(cols)), rows, cols, (l,u))) == 
             zeros(Float64, 10, 10)
 
         @test Matrix(BlockBandedMatrix{Int}(Zeros(sum(rows),sum(cols)), rows,cols, (l,u))) ==
@@ -35,11 +35,13 @@ import Base.Broadcast: materialize!
         N = M = 4
         cols = rows = 1:N
         A = BlockBandedMatrix{Int}(undef, rows,cols, (l,u))
-        A.data .= 1:length(A.data)
+            A.data .= 1:length(A.data)
 
         @test A[1,1] == 1
         @test A[1,3] == 10
 
+
+
         @test blockbandwidth(A,1)  == 1
         @test blockbandwidths(A) == (l,u)
 
@@ -90,34 +92,12 @@ import Base.Broadcast: materialize!
         @test ret[1,2] == 0
     end
 
-    @testset "block-banded matrix interface for blockranges" begin
-        l , u = 1,1
-        N = M = 4
-        cols = rows = 1:N
-        A = BlockBandedMatrix{Int}(undef, rows,cols, (l,u))
-        A.data .= 1:length(A.data)
-
-        V = view(A, Block.(2:3), Block.(3:4))
-        @test isblockbanded(V)
-
-        B = BlockBandedMatrix(V)
-        @test B isa BlockBandedMatrix
-        @test blockbandwidths(V) == blockbandwidths(B) == (2,0)
-        @test B == V == A[Block.(2:3), Block.(3:4)]
-
-        @test A[Block.(2:3), Block.(3:4)] isa BlockBandedMatrix
-
-        x = randn(size(B,2))
-        y = similar(x, size(B,1))
-        @test all((similar(y) .= MulAdd(B, x)) .=== (similar(y) .= MulAdd(V,x)))
-    end
-
     @testset "BlockBandedMatrix indexing" begin
         l , u = 1,1
         N = M = 10
         cols = rows = 1:N
         A = BlockBandedMatrix{Float64}(undef, rows,cols, (l,u))
-        A.data .= 1:length(A.data)
+            A.data .= 1:length(A.data)
 
         A[1,1] = 5
         @test A[1,1] == 5