inf QR works!

Author: dlfivefifty

src/infqr.jl

@@ -152,4 +152,25 @@ function lmul!(adjA::Adjoint{<:Any,<:QRPackedQ{<:Any,<:AdaptiveQRFactors}}, B::C
         end
     end
     B
-end
+end
+
+function (*)(A::QRPackedQ{T,<:AdaptiveQRFactors}, x::AbstractVector{S}) where {T,S}
+    TS = promote_op(matprod, T, S)
+    lmul!(A, Base.copymutable(convert(AbstractVector{TS},x)))
+end
+
+function (*)(A::Adjoint{T,<:QRPackedQ{T,<:AdaptiveQRFactors}}, x::AbstractVector{S}) where {T,S}
+    TS = promote_op(matprod, T, S)
+    lmul!(A, Base.copymutable(convert(AbstractVector{TS},x)))
+end
+
+ldiv!(R::UpperTriangular{<:Any,<:AdaptiveQRFactors}, B::AbstractVector) = materialize!(Ldiv(R, B))
+
+
+
+ldiv!(dest::AbstractVector, F::QR{<:Any,<:AdaptiveQRFactors}, b::AbstractVector) = 
+    ldiv!(F, copyto!(dest, b))
+ldiv!(F::QR{<:Any,<:AdaptiveQRFactors}, b::AbstractVector) = ldiv!(F.R, lmul!(F.Q',b))
+\(F::QR{<:Any,<:AdaptiveQRFactors}, B::AbstractVector) = ldiv!(F.R, F.Q'B)
+
+factorize(A::BandedMatrix{<:Any,<:Any,<:OneToInf}) = qr(A)

test/test_infqr.jl

@@ -63,15 +63,21 @@ import InfiniteLinearAlgebra: partialqr!, AdaptiveQRData, AdaptiveLayout
         Q,R = qr(A);
         b = Vcat([1.,2,3],Zeros(∞))
         @test lmul!(Q, Base.copymutable(b)).data ≈ qr(A[1:4,1:3]).Q*[1,2,3]
+
         @test Q[1,1] ≈ -1/sqrt(2)
         @test Q[200_000,200_000] ≈ -1.0
         @test Q[1:101,1:100] ≈ qr(A[1:101,1:100]).Q[:,1:100]
 
+        @test lmul!(Q, Base.copymutable(b)) == Q*b
+
         r = lmul!(Q', Base.copymutable(b))
         nr = length(r.data)
         @test qr(A[1:nr+1,1:nr]).Q'b[1:nr+1] ≈ r[1:nr+1]
 
-        materialize!(Ldiv(R, r))
+        @test Q'*b == r
+
+        @test qr(A)\b == A\b
+        @test (A*(A\b))[1:100] ≈ [1:3; Zeros(97)]
     end
 end