support padded coefficients (for ∞-dimensional case) (#93)

* support padded coefficients (for ∞-dimensional case)

* add tests

Author: dlfivefifty

src/ContinuumArrays.jl

@@ -8,7 +8,7 @@ import Base.Broadcast: materialize, BroadcastStyle, broadcasted
 import LazyArrays: MemoryLayout, Applied, ApplyStyle, flatten, _flatten, colsupport, most, combine_mul_styles, AbstractArrayApplyStyle,
                         adjointlayout, arguments, _mul_arguments, call, broadcastlayout, layout_getindex, UnknownLayout,
                         sublayout, sub_materialize, ApplyLayout, BroadcastLayout, combine_mul_styles, applylayout,
-                        simplifiable, _simplify, AbstractLazyLayout
+                        simplifiable, _simplify, AbstractLazyLayout, PaddedLayout
 import LinearAlgebra: pinv, dot, norm2
 import BandedMatrices: AbstractBandedLayout, _BandedMatrix
 import BlockArrays: block, blockindex, unblock, blockedrange, _BlockedUnitRange, _BlockArray

src/bases/bases.jl

@@ -122,11 +122,18 @@ copy(L::Ldiv{<:AbstractBasisLayout,BroadcastLayout{typeof(*)},<:Any,<:AbstractQu
 
 # expansion
 _grid(_, P) = error("Overload Grid")
-_grid(lay::ApplyLayout{typeof(*)}, P) = grid(first(arguments(lay,P)))
+
+_mul_grid(_, args) = grid(first(args))
+_mul_grid(::Tuple{Any,PaddedLayout}, (P,c)) = grid(P[:,colsupport(c)])
+
+function _grid(lay::ApplyLayout{typeof(*)}, P)
+    args = arguments(lay,P)
+    _mul_grid(map(MemoryLayout,args), args)
+end
 _grid(::MappedBasisLayout, P) = invmap(parentindices(P)[1])[grid(demap(P))]
 _grid(::SubBasisLayout, P) = grid(parent(P))
 _grid(::WeightedBasisLayouts, P) = grid(unweightedbasis(P))
-grid(P) = _grid(MemoryLayout(typeof(P)), P)
+grid(P) = _grid(MemoryLayout(P), P)
 
 
 struct TransformFactorization{T,Grid,Plan,IPlan} <: Factorization{T}

test/runtests.jl

@@ -605,4 +605,10 @@ include("test_basisconcat.jl")
     u = L*randn(6)
     rep = RecipesBase.apply_recipe(Dict{Symbol, Any}(), u)
     @test rep[1].args == (L.points,u[L.points])
+
+    @testset "padded" begin
+        u = L * Vcat(rand(3), Zeros(3))
+        rep = RecipesBase.apply_recipe(Dict{Symbol, Any}(), u)
+        @test rep[1].args == (L.points,u[L.points])
+    end
 end