restore sub/diagonaldata

Author: dlfivefifty

Project.toml

@@ -17,7 +17,7 @@ MatrixFactorizations = "a3b82374-2e81-5b9e-98ce-41277c0e4c87"
 SemiseparableMatrices = "f8ebbe35-cbfb-4060-bf7f-b10e4670cf57"
 
 [compat]
-ArrayLayouts = "0.6"
+ArrayLayouts = "0.5.1, 0.6"
 BandedMatrices = "0.16"
 BlockArrays = "0.14.5"
 BlockBandedMatrices = "0.10"

src/InfiniteLinearAlgebra.jl

@@ -8,7 +8,7 @@ import Base: +, -, *, /, \, ^, OneTo, getindex, promote_op, _unsafe_getindex, si
 import Base.Broadcast: BroadcastStyle, Broadcasted, broadcasted
 
 import ArrayLayouts: colsupport, rowsupport, triangularlayout, MatLdivVec, triangulardata, TriangularLayout, TridiagonalLayout, 
-                        sublayout, _qr, __qr, MatLmulVec, MatLmulMat, AbstractQLayout, materialize!, subdiag, supdiag,
+                        sublayout, _qr, __qr, MatLmulVec, MatLmulMat, AbstractQLayout, materialize!, diagonaldata, subdiagonaldata, supdiagonaldata,
                         _bidiag_forwardsub!, mulreduce, RangeCumsum
 import BandedMatrices: BandedMatrix, _BandedMatrix, AbstractBandedMatrix, bandeddata, bandwidths, BandedColumns, bandedcolumns,
                         _default_banded_broadcast, banded_similar

src/banded/infbanded.jl

@@ -1,3 +1,5 @@
+_BandedMatrix(data::AbstractMatrix{T}, ::Infinity, l, u) where T = _BandedMatrix(data, ℵ₀, l, u)
+
 ###
 # BandIndexing
 ###
@@ -439,8 +441,8 @@ mulreduce(M::Mul{<:AbstractQLayout, <:PertToeplitzLayout}) = ApplyArray(M)
 
 function _bidiag_forwardsub!(M::Ldiv{<:Any,<:PaddedLayout})
     A, b_in = M.A, M.B
-    dv = diag(A)
-    ev = subdiag(A)
+    dv = diagonaldata(A)
+    ev = subdiagonaldata(A)
     b = paddeddata(b_in)
     N = length(b)
     b[1] = bj1 = dv[1]\b[1]

src/inful.jl

@@ -35,26 +35,26 @@ function _ultailL1(C::AbstractMatrix, A::AbstractMatrix, B::AbstractMatrix)
 end
 
 function _ul(::TridiagonalToeplitzLayout, J::AbstractMatrix, ::Val{false}; check::Bool = true)
-    C = getindex_value(subdiag(J))
-    A = getindex_value(diag(J))
-    B = getindex_value(supdiag(J))
+    C = getindex_value(subdiagonaldata(J))
+    A = getindex_value(diagonaldata(J))
+    B = getindex_value(supdiagonaldata(J))
     L = _ultailL1(C, A, B)
     U = B/L
     UL(Tridiagonal(Fill(convert(typeof(L),C),∞), Fill(L,∞), Fill(U,∞)), OneToInf(), 0)
 end
 
 function _ul(::TridiagonalToeplitzLayout, J::AbstractMatrix, ::Val{true}; check::Bool = true)
-    C = getindex_value(subdiag(J))
-    A = getindex_value(diag(J))
-    B = getindex_value(supdiag(J))
+    C = getindex_value(subdiagonaldata(J))
+    A = getindex_value(diagonaldata(J))
+    B = getindex_value(supdiagonaldata(J))
     A^2 ≥ 4B*C || error("Pivotting not implemented")
     ul(J, Val(false))
 end
 
 function _ul(::BlockTridiagonalToeplitzLayout, J::AbstractMatrix, ::Val{false}; check::Bool = true)
-    C = getindex_value(subdiag(blocks(J)))
-    A = getindex_value(diag(blocks(J)))
-    B = getindex_value(supdiag(blocks(J)))
+    C = getindex_value(subdiagonaldata(blocks(J)))
+    A = getindex_value(diagonaldata(blocks(J)))
+    B = getindex_value(supdiagonaldata(blocks(J)))
     L = _ultailL1(C, A, B)
     U = B/L
     UL(mortar(Tridiagonal(Fill(convert(typeof(L),C),∞), Fill(L,∞), Fill(U,∞))), OneToInf(), 0)