Support Symmetric(::BandedBlockBandedMatrix)

Author: dlfivefifty

Project.toml

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

src/AbstractBlockBandedMatrix.jl

@@ -39,6 +39,21 @@ conjlayout(::Type{T}, ::BandedBlockBandedRows{Lay}) where {T<:Complex,Lay} = Ban
 conjlayout(::Type{T}, ::BlockBandedColumns{Lay}) where {T<:Complex,Lay} = BlockBandedColumns{typeof(conjlayout(T,Lay))}()
 conjlayout(::Type{T}, ::BlockBandedRows{Lay}) where {T<:Complex,Lay} = BlockBandedRows{typeof(conjlayout(T,Lay))}()
 
+symmetriclayout(::AbstractBandedBlockBandedLayout) = BandedBlockBandedLayout()
+hermitianlayout(::AbstractBandedBlockBandedLayout) = BandedBlockBandedLayout()
+
+function blockbandwidths(S::HermOrSym)
+    l, u = blockbandwidths(parent(S))
+    if S.uplo == 'U'
+        (u, u)
+    else
+        (l, l)
+    end
+end
+function subblockbandwidths(S::HermOrSym)
+    m = max(subblockbandwidths(parent(S))...)
+    (m,m)
+end
 
 # AbstractBandedMatrix must implement
 

src/BlockBandedMatrices.jl

@@ -27,7 +27,7 @@ import ArrayLayouts: BlasMatLmulVec, MatLmulVec, MatLmulMat,
                     AbstractColumnMajor, DenseColumnMajor, ColumnMajor,
                     DiagonalLayout, MulAdd, mul, colsupport, rowsupport,
                     _qr, _factorize, _copyto!, zero!, layout_replace_in_print_matrix,
-                    transposelayout, conjlayout
+                    transposelayout, conjlayout, symmetriclayout, hermitianlayout
 
 import BlockArrays: blocksize, blockcheckbounds, BlockedUnitRange, blockisequal, DefaultBlockAxis,
                         Block, BlockSlice, unblock, block, blockindex,

test/test_bandedblockbanded.jl

@@ -516,6 +516,20 @@ import ArrayLayouts: RangeCumsum
         @test convert(AbstractArray{Float64}, A) == convert(AbstractMatrix{Float64}, A) == A
         @test convert(AbstractArray{Int}, A) ≡ convert(AbstractMatrix{Int}, A) ≡ A
     end
+
+    @testset "Symmetric" begin
+        l, u = 2, 1
+        λ, μ = 2, 1
+        N = M = 4
+        cols = rows = 1:N
+        data = reshape(collect(1:(λ+μ+1)*(l+u+1)*sum(cols)), ((λ + μ + 1) * (l + u + 1), sum(cols)))
+        A = _BandedBlockBandedMatrix(data, rows, cols, (l, u), (λ, μ))
+
+        @test BandedBlockBandedMatrix(Symmetric(A)) == Symmetric(Matrix(A))
+        @test BandedBlockBandedMatrix(Symmetric(A,:L)) == Symmetric(Matrix(A),:L)
+        @test BandedBlockBandedMatrix(Hermitian(A)) == Hermitian(Matrix(A))
+        @test BandedBlockBandedMatrix(Hermitian(A,:L)) == Hermitian(Matrix(A),:L)
+    end
 end
 
 if false # turned off since tests have check-bounds=yes