support padded coefficients (for ∞-dimensional case)

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,7 +122,14 @@ 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))