InfBandedMatrices -> InfLinearAlgebra

Author: dlfivefifty

Project.toml

@@ -1,5 +1,6 @@
-name = "InfiniteBandedMatrices"
-uuid = "75c36928-c531-5f82-ab28-1f08c0e381e4"
+name = "InfiniteLinearAlgebra"
+uuid = "cde9dba0-b1de-11e9-2c62-0bab9446c55c"
+version = "0.0.1"
 
 [deps]
 BandedMatrices = "aae01518-5342-5314-be14-df237901396f"

README.md

@@ -1,4 +1,4 @@
-# InfiniteBandedMatrices.jl
+# InfiniteLinearAlgebra.jl
 
 A Julia repository for linear algebra with infinite banded and block-banded matrices
 
@@ -10,6 +10,6 @@ ql(A - 5*I)
 ```
 The infinite-dimensional QL decomposition is a subtly thing: its defined when the operator has non-positive Fredholm index, and if the Fredholm index is not zero, it may not be unique. For the Bull head matrix `A`, here are plots of `ql(A-λ*I).L[1,1]` alongside the image of the symbol `A`, which depicts the essential spectrum of `A` and where the Fredholm index changes. Note we have two plots as the regions with negative Fredholm index  have multiple QL decompositions. Where the Fredholm index is positive, the QL decomposition doesn't exist and is depected in black.
 
-<img src=https://github.com/JuliaMatrices/InfiniteBandedMatrices.jl/raw/master/images/ql1.png width=500 height=400>
-<img src=https://github.com/JuliaMatrices/InfiniteBandedMatrices.jl/raw/master/images/ql2.png width=500 height=400>
+<img src=https://github.com/JuliaMatrices/InfiniteLinearAlgebra.jl/raw/master/images/ql1.png width=500 height=400>
+<img src=https://github.com/JuliaMatrices/InfiniteLinearAlgebra.jl/raw/master/images/ql2.png width=500 height=400>
 

examples/blocktridiagonal.jl

@@ -1,4 +1,4 @@
-using Revise, InfiniteBandedMatrices, BlockBandedMatrices, BandedMatrices, BlockArrays, InfiniteArrays, FillArrays, LazyArrays, Test
+using Revise, InfiniteLinearAlgebra, BlockBandedMatrices, BandedMatrices, BlockArrays, InfiniteArrays, FillArrays, LazyArrays, Test
 
 A = BlockTridiagonal(Vcat([fill(1.0,2,1),Matrix(1.0I,2,2),Matrix(1.0I,2,2),Matrix(1.0I,2,2)],Fill(Matrix(1.0I,2,2), ∞)), 
                        Vcat([zeros(1,1)], Fill(zeros(2,2), ∞)), 

examples/intervalarithmetic.jl

@@ -1,2 +1,2 @@
-using InfiniteBandedMatrices
+using InfiniteLinearAlgebra
 import IntervalArithmetic

examples/jacobi.jl

@@ -1,4 +1,4 @@
-using Revise, InfiniteBandedMatrices, MatrixFactorizations
+using Revise, InfiniteLinearAlgebra, MatrixFactorizations
 import IntervalArithmetic
 import IntervalArithmetic: Interval, emptyinterval
 import MatrixFactorizations: add_mul_signed, reflector!
@@ -26,12 +26,12 @@ ql(Δ - Interval(4.2)I).L[1,1]
 rig_qltail(0.5,-1.98,2.0)
 
 A = (Δ-λ*I)
-import InfiniteBandedMatrices: _qlfactUnblocked!
-InfiniteBandedMatrices.qltail(0.5,Interval(-4.1,-4.0), 0.5)
+import InfiniteLinearAlgebra: _qlfactUnblocked!
+InfiniteLinearAlgebra.qltail(0.5,Interval(-4.1,-4.0), 0.5)
 
-InfiniteBandedMatrices.qltail(0.5,-4.1, 0.5)
+InfiniteLinearAlgebra.qltail(0.5,-4.1, 0.5)
 
-InfiniteBandedMatrices.qltail(0.5,-4.0, 0.5)
+InfiniteLinearAlgebra.qltail(0.5,-4.0, 0.5)
 
 Z,A,B = 0.5,-4.0,0.5
 
@@ -47,7 +47,7 @@ d,e = d2,e2
 rig_qltail(Z,A,B,d,e)
 
 
-InfiniteBandedMatrices.qltail(Z,A,B)
+InfiniteLinearAlgebra.qltail(Z,A,B)
 
 @which reflector!(X[2:-1:1,end])
 
@@ -80,7 +80,7 @@ sign(Interval(-10,10))
 
 length(d)
 
-InfiniteBandedMatrices.qltail(Z,A,B)
+InfiniteLinearAlgebra.qltail(Z,A,B)
 
 
 

examples/periodicschrodinger.jl

@@ -1,4 +1,4 @@
-using InfiniteBandedMatrices, Plots
+using InfiniteLinearAlgebra, Plots
 
 # Consider a rank-1 perturbation of periodic semi-infinite Schrödinger (Jacobi) operator.
 # We can construct it as a block matrix as follows:

examples/toeplitz.jl

@@ -1,6 +1,6 @@
-using Revise, InfiniteBandedMatrices, BlockBandedMatrices, BlockArrays, BandedMatrices, LazyArrays, FillArrays, MatrixFactorizations, Plots
+using Revise, InfiniteLinearAlgebra, BlockBandedMatrices, BlockArrays, BandedMatrices, LazyArrays, FillArrays, MatrixFactorizations, Plots
 import MatrixFactorizations: reflectorApply!, QLPackedQ
-import InfiniteBandedMatrices: blocktailiterate, _ql, qltail, rightasymptotics
+import InfiniteLinearAlgebra: blocktailiterate, _ql, qltail, rightasymptotics
 import BandedMatrices: bandeddata,_BandedMatrix
 
 function ℓ11(A,λ; kwds...) 
@@ -82,7 +82,7 @@ a = z -> 4z^10 + 1/z
 θ = range(0,2π; length=1000)
 plot!(a.(exp.(im.*θ)); linewidth=2.0, linecolor=:blue, legend=false)
 
-import InfiniteBandedMatrices: tail_de
+import InfiniteLinearAlgebra: tail_de
 a = reverse(A.data.applied.args[1]) .+ 0im
 
 a = randn(10) .+ im*randn(10); a[1] += 10; a[end] = 1;

src/InfiniteBandedMatrices.jl renamed to src/InfiniteLinearAlgebra.jl

@@ -1,4 +1,4 @@
-module InfiniteBandedMatrices
+module InfiniteLinearAlgebra
 using BlockArrays, BlockBandedMatrices, BandedMatrices, LazyArrays, FillArrays, InfiniteArrays, MatrixFactorizations, LinearAlgebra
 
 import Base: +, -, *, /, \, OneTo, getindex, promote_op, _unsafe_getindex, print_matrix_row, size,
@@ -32,12 +32,12 @@ end
 
 export Vcat, Fill, ql, ql!, ∞, ContinuousSpectrumError, BlockTridiagonal
 
-include("hessenbergq.jl")
+include("banded/hessenbergq.jl")
 
-include("infbanded.jl")
-include("infblocktridiagonal.jl")
-include("infqltoeplitz.jl")
+include("banded/infbanded.jl")
+include("blockbanded/infblocktridiagonal.jl")
+include("banded/infqltoeplitz.jl")
 include("infql.jl")
-include("rigorous.jl")
+include("banded/rigorous.jl")
 
 end # module