Add Semiseparable

Author: dlfivefifty

Project.toml

@@ -13,17 +13,17 @@ LazyArrays = "5078a376-72f3-5289-bfd5-ec5146d43c02"
 LazyBandedMatrices = "d7e5e226-e90b-4449-9968-0f923699bf6f"
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
 MatrixFactorizations = "a3b82374-2e81-5b9e-98ce-41277c0e4c87"
+SemiseparableMatrices = "f8ebbe35-cbfb-4060-bf7f-b10e4670cf57"
 
 [compat]
-ArrayLayouts = "0.1"
-BandedMatrices = "0.14"
+ArrayLayouts = "0.1, 0.2"
+BandedMatrices = "0.14, 0.15"
 BlockArrays = "0.11"
 BlockBandedMatrices = "0.7.1"
 FillArrays = "0.8.4"
-InfiniteArrays = "0.6.1"
+InfiniteArrays = "0.6.1, 0.7"
 LazyArrays = "0.15"
 LazyBandedMatrices = "0.2"
-MatrixFactorizations = "0.2"
 julia = "1.3"
 
 [extras]

src/InfiniteLinearAlgebra.jl

@@ -1,5 +1,5 @@
 module InfiniteLinearAlgebra
-using BlockArrays, BlockBandedMatrices, BandedMatrices, LazyArrays, LazyBandedMatrices, 
+using BlockArrays, BlockBandedMatrices, BandedMatrices, LazyArrays, LazyBandedMatrices, SemiseparableMatrices,
         FillArrays, InfiniteArrays, MatrixFactorizations, ArrayLayouts, LinearAlgebra
 
 import Base: +, -, *, /, \, ^, OneTo, getindex, promote_op, _unsafe_getindex, print_matrix_row, size, axes,

test/test_infqr.jl

@@ -1,8 +1,9 @@
 using InfiniteLinearAlgebra, LinearAlgebra, BandedMatrices, InfiniteArrays, MatrixFactorizations, LazyArrays, FillArrays, SpecialFunctions, Test
-import LazyArrays: colsupport, rowsupport, MemoryLayout, DenseColumnMajor, TriangularLayout, resizedata!
+import LazyArrays: colsupport, rowsupport, MemoryLayout, DenseColumnMajor, TriangularLayout, resizedata!, arguments
 import LazyBandedMatrices: BroadcastBandedLayout
 import BandedMatrices: _BandedMatrix, _banded_qr!, BandedColumns
 import InfiniteLinearAlgebra: partialqr!, AdaptiveQRData, AdaptiveLayout
+import SemiseparableMatrices: AlmostBandedLayout
 
 
 @testset "Adaptive QR" begin
@@ -132,4 +133,11 @@ import InfiniteLinearAlgebra: partialqr!, AdaptiveQRData, AdaptiveLayout
         @test F.Q[1:10,1:10] == Eye(10)
         @test F.R[1:10,1:10] == A[1:10,1:10]
     end
+
+    @testset "almost-banded" begin
+        A = Vcat(Ones(1,∞), BandedMatrix(0 => -Ones(∞), 1 => 1:∞))
+        @test MemoryLayout(typeof(A)) == AlmostBandedLayout()
+        B,C = arguments(A)
+        
+    end
 end