ldiv for Mapped and unmapped (#88)

* ldiv for Mapped and unmapped

* Update bases.jl

* Update runtests.jl

* Update runtests.jl

* v0.7.1

Author: dlfivefifty

Project.toml

@@ -1,6 +1,6 @@
 name = "ContinuumArrays"
 uuid = "7ae1f121-cc2c-504b-ac30-9b923412ae5c"
-version = "0.7.0"
+version = "0.7.1"
 
 [deps]
 ArrayLayouts = "4c555306-a7a7-4459-81d9-ec55ddd5c99a"

src/bases/bases.jl

@@ -85,6 +85,12 @@ end
     demap(A)\demap(B)
 end
 
+function copy(P::Ldiv{<:MappedBasisLayouts,<:AbstractLazyLayout})
+    A,B = P.A, P.B
+    demap(A) \ B[invmap(basismap(A))]
+end
+copy(P::Ldiv{<:MappedBasisLayouts,ApplyLayout{typeof(*)}}) = copy(Ldiv{UnknownLayout,ApplyLayout{typeof(*)}}(P.A,P.B))
+
 @inline copy(L::Ldiv{<:AbstractBasisLayout,<:SubBasisLayouts}) = apply(\, L.A, ApplyQuasiArray(L.B))
 @inline function copy(L::Ldiv{<:SubBasisLayouts,<:AbstractBasisLayout}) 
     P = parent(L.A)
@@ -338,15 +344,14 @@ sublayout(::MappedWeightedBasisLayout, ::Type{<:Tuple{<:Inclusion,<:AbstractVect
 
 demap(x) = x
 demap(x::BroadcastQuasiArray) = BroadcastQuasiArray(x.f, map(demap, arguments(x))...)
-demap(V::SubQuasiArray{<:Any,2,<:Any,<:Tuple{<:AbstractAffineQuasiVector,<:Slice}}) = parent(V)
-demap(V::SubQuasiArray{<:Any,1,<:Any,<:Tuple{<:AbstractAffineQuasiVector}}) = parent(V)
+demap(V::SubQuasiArray{<:Any,2,<:Any,<:Tuple{Any,Slice}}) = parent(V)
+demap(V::SubQuasiArray{<:Any,1}) = parent(V)
 function demap(V::SubQuasiArray{<:Any,2}) 
     kr, jr = parentindices(V)
     demap(parent(V)[kr,:])[:,jr]
 end
 
-basismap(x::SubQuasiArray{<:Any,2,<:Any,<:Tuple{<:AbstractAffineQuasiVector,<:Any}}) = parentindices(x)[1]
-basismap(x::SubQuasiArray{<:Any,1,<:Any,<:Tuple{<:AbstractAffineQuasiVector}}) = parentindices(x)[1]
+basismap(x::SubQuasiArray) = parentindices(x)[1]
 basismap(x::BroadcastQuasiArray) = basismap(x.args[1])
 
 

test/runtests.jl

@@ -431,6 +431,11 @@ end
             @test f + g isa Expansion
             @test (f+g)[0.1] ≈ f[0.1] + g[0.1]
         end
+
+        @testset "vec demap" begin
+            @test L[y,:] \ exp.(axes(L,1))[y] ≈ L[y,:] \ exp.(y) ≈  factorize(L[y,:]) \ exp.(y)
+            @test ContinuumArrays.demap(view(axes(L,1),y)) == axes(L,1)
+        end
     end
 
     @testset "diff" begin