Move out BandedMatrices code (#192)

* Drop Julia <v1.10

* Move out BandedMatrices code

* Update infblocktridiagonal.jl

* Update runtests.jl

* move out blockarrays

* move out block code

* compiles

* tests pass

* Update InfiniteLinearAlgebra.jl

Author: dlfivefifty

Project.toml

@@ -1,6 +1,6 @@
 name = "InfiniteLinearAlgebra"
 uuid = "cde9dba0-b1de-11e9-2c62-0bab9446c55c"
-version = "0.8.5"
+version = "0.9.0"
 
 [deps]
 ArrayLayouts = "4c555306-a7a7-4459-81d9-ec55ddd5c99a"
@@ -23,11 +23,11 @@ BandedMatrices = "1.7.2"
 BlockArrays = "1.0"
 BlockBandedMatrices = "0.13"
 FillArrays = "1.0"
-InfiniteArrays = "0.14"
+InfiniteArrays = "0.15"
 InfiniteRandomArrays = "0.2"
 Infinities = "0.1"
-LazyArrays = "2.0"
-LazyBandedMatrices = "0.10"
+LazyArrays = "2.3"
+LazyBandedMatrices = "0.11"
 LinearAlgebra = "1"
 MatrixFactorizations = "3.0"
 Random = "1"
@@ -39,8 +39,8 @@ julia = "1.10"
 
 [extras]
 Aqua = "4c88cf16-eb10-579e-8560-4a9242c79595"
-Random = "9a3f8284-a2c9-5f02-9a11-845980a1fd5c"
 InfiniteRandomArrays = "2bc77966-89c7-476d-a40f-269028fac4a9"
+Random = "9a3f8284-a2c9-5f02-9a11-845980a1fd5c"
 SpecialFunctions = "276daf66-3868-5448-9aa4-cd146d93841b"
 StaticArrays = "90137ffa-7385-5640-81b9-e52037218182"
 Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"

src/InfiniteLinearAlgebra.jl

@@ -14,17 +14,17 @@ import Base.Broadcast: BroadcastStyle, Broadcasted, broadcasted
 
 import ArrayLayouts: AbstractBandedLayout, AbstractQLayout, AdjQRPackedQLayout, CNoPivot, DenseColumnMajor, FillLayout,
                      MatLdivVec, MatLmulMat, MatLmulVec, MemoryLayout, QRPackedQLayout, RangeCumsum, TriangularLayout,
-                     TridiagonalLayout, __qr, _bidiag_forwardsub!, _factorize, _qr, check_mul_axes, colsupport,
+                     TridiagonalLayout, __qr, _factorize, _qr, check_mul_axes, colsupport,
                      diagonaldata, ldiv!, lmul!, mul, mulreduce, reflector!, reflectorApply!,
                      rowsupport, sub_materialize, subdiagonaldata, sublayout, supdiagonaldata, transposelayout,
                      triangulardata, triangularlayout, zero!, materialize!
 
 import BandedMatrices: BandedColumns, BandedMatrix, BandedMatrix, _BandedMatrix, AbstractBandedMatrix,
-                       _BandedMatrix, _BandedMatrix, _banded_qr, _banded_qr!, _default_banded_broadcast, banded_chol!,
+                       _BandedMatrix, _BandedMatrix, _banded_qr, _banded_qr!, banded_chol!,
                        banded_similar, bandedcolumns, bandeddata, bandwidths
 
 import BlockArrays: AbstractBlockLayout, BlockLayout, BlockSlice, BlockSlice1, BlockedOneTo,
-                    blockcolsupport, sizes_from_blocks, OneToCumsum, AbstractBlockedUnitRange
+                    blockcolsupport, sizes_from_blocks, AbstractBlockedUnitRange
 
 import BlockBandedMatrices: AbstractBlockBandedLayout, BlockBandedMatrix, BlockSkylineMatrix,
                             BlockSkylineSizes, BlockTridiagonal, _BlockBandedMatrix, _BlockSkylineMatrix,
@@ -53,35 +53,20 @@ import MatrixFactorizations: AdjQLPackedQLayout, LayoutQ, QL, QLPackedQ, QLPacke
                              ul!, ul_layout
 
 import SemiseparableMatrices: AbstractAlmostBandedLayout, _almostbanded_qr!
-
-
-LinearAlgebra._cut_B(x::AbstractVector, ::InfUnitRange) = x
-LinearAlgebra._cut_B(X::AbstractMatrix, ::InfUnitRange) = X
-
-
-if VERSION ≥ v"1.11.0-DEV.21"
-    using LinearAlgebra: UpperOrLowerTriangular
-else
-    const UpperOrLowerTriangular{T,S} = Union{LinearAlgebra.UpperTriangular{T,S},
-                                              LinearAlgebra.UnitUpperTriangular{T,S},
-                                              LinearAlgebra.LowerTriangular{T,S},
-                                              LinearAlgebra.UnitLowerTriangular{T,S}}
-end
-
+import InfiniteArrays: UpperOrLowerTriangular, TridiagonalToeplitzLayout, TriToeplitz, PertTridiagonalToeplitzLayout, PertConstRows, 
+                        subdiagonalconstant, diagonalconstant, supdiagonalconstant
 
 # BroadcastStyle(::Type{<:BandedMatrix{<:Any,<:Any,<:OneToInf}}) = LazyArrayStyle{2}()
 
-function ArrayLayouts._power_by_squaring(_, ::NTuple{2,InfiniteCardinal{0}}, A::AbstractMatrix{T}, p::Integer) where T
-    if p < 0
-        inv(A)^(-p)
-    elseif p == 0
-        Eye{T}(∞)
-    elseif p == 1
-        copy(A)
-    else
-        A*A^(p-1)
-    end
-end
+const InfiniteArraysBandedMatricesExt = Base.get_extension(InfiniteArrays, :InfiniteArraysBandedMatricesExt)
+const InfiniteArraysBlockArraysExt = Base.get_extension(InfiniteArrays, :InfiniteArraysBlockArraysExt)
+const InfBandedMatrix = InfiniteArraysBandedMatricesExt.InfBandedMatrix
+const InfToeplitz = InfiniteArraysBandedMatricesExt.InfToeplitz
+const PertToeplitzLayout = InfiniteArraysBandedMatricesExt.PertToeplitzLayout
+const PertToeplitz = InfiniteArraysBandedMatricesExt.PertToeplitz
+const BlockTriPertToeplitz = InfiniteArraysBlockArraysExt.BlockTriPertToeplitz
+const BlockTridiagonalToeplitzLayout = InfiniteArraysBlockArraysExt.BlockTridiagonalToeplitzLayout
+
 
 function choplength(c::AbstractVector, tol)
     @inbounds for k = length(c):-1:1
@@ -128,12 +113,10 @@ pad(c::Transpose, ax, bx) = transpose(pad(parent(c), bx, ax))
 pad(c::Adjoint, ax, bx) = adjoint(pad(parent(c), bx, ax))
 pad(c::BlockVec, ax::BlockedOneTo{Int,<:InfStepRange}) = BlockVec(pad(c.args[1], size(c.args[1],1), ∞))
 
-export Vcat, Fill, ql, ql!, ∞, ContinuousSpectrumError, BlockTridiagonal
+export ∞, ContinuousSpectrumError, BlockTridiagonal
 
 include("banded/hessenbergq.jl")
 
-include("banded/infbanded.jl")
-include("blockbanded/blockbanded.jl")
 include("banded/infqltoeplitz.jl")
 include("banded/infreversecholeskytoeplitz.jl")
 include("banded/infreversecholeskytridiagonal.jl")