avoid simplifying Inv * Lazy

JuliaLinearAlgebra · Jul 20, 2023 · 5d10bd9 · 5d10bd9
1 parent 9960f3d
commit 5d10bd9
Show file tree

Hide file tree

Showing 2 changed files with 5 additions and 1 deletion.
diff --git a/src/LazyBandedMatrices.jl b/src/LazyBandedMatrices.jl
@@ -935,8 +935,11 @@ copy(M::Mul{<:BroadcastLayouts,ApplyLayout{typeof(*)}}) = simplify(M)
 copy(M::Mul{<:AbstractInvLayout{<:BandedLazyLayouts},<:StructuredLazyLayouts}) = ArrayLayouts.ldiv(pinv(M.A), M.B)
 copy(M::Mul{<:StructuredApplyLayouts{typeof(*)},<:AbstractInvLayout{<:BandedLazyLayouts}}) = simplify(M)
 copy(M::Mul{ApplyLayout{typeof(*)},<:AbstractInvLayout{<:BandedLazyLayouts}}) = simplify(M)
+
 copy(M::Mul{<:AbstractInvLayout{<:BandedLazyLayouts},<:StructuredApplyLayouts{typeof(*)}}) = simplify(M)
 copy(M::Mul{<:AbstractInvLayout, <:StructuredApplyLayouts{typeof(*)}}) = simplify(M)
+simplifiable(::Mul{<:AbstractInvLayout, <:StructuredLazyLayouts}) = Val(false)
+copy(M::Mul{<:AbstractInvLayout, <:StructuredLazyLayouts}) = simplify(M)
 
 
 copy(L::Ldiv{<:StructuredLazyLayouts, <:StructuredLazyLayouts}) = lazymaterialize(\, L.A, L.B)

diff --git a/test/test_banded.jl b/test/test_banded.jl
@@ -1,4 +1,5 @@
 using LazyBandedMatrices, LazyArrays, Test
+using BandedMatrices: _BandedMatrix
 
 struct PseudoBandedMatrix{T} <: AbstractMatrix{T}
     data::Array{T}
@@ -710,7 +711,7 @@ LinearAlgebra.lmul!(β::Number, A::PseudoBandedMatrix) = (lmul!(β, A.data); A)
         M = ApplyArray(*,L,L)
         @test Bi * M ≈ inv(B) * L *L
         @test M* Bi ≈ L *L * inv(B)
-        @test Li * M ≈ M * Li ≈ L 
+        @test Li * M ≈ M * Li ≈ L
 
         x = Vcat([1,2], Zeros(8))
         @test Li * x ≈ L \ x