add tests and use expand

Author: dlfivefifty

Project.toml

@@ -8,6 +8,7 @@ ArrayLayouts = "4c555306-a7a7-4459-81d9-ec55ddd5c99a"
 BandedMatrices = "aae01518-5342-5314-be14-df237901396f"
 BlockArrays = "8e7c35d0-a365-5155-bbbb-fb81a777f24e"
 DomainSets = "5b8099bc-c8ec-5219-889f-1d9e522a28bf"
+FastTransforms = "057dd010-8810-581a-b7be-e3fc3b93f78c"
 FillArrays = "1a297f60-69ca-5386-bcde-b61e274b549b"
 InfiniteArrays = "4858937d-0d70-526a-a4dd-2d5cb5dd786c"
 Infinities = "e1ba4f0e-776d-440f-acd9-e1d2e9742647"

src/bases/bases.jl

@@ -122,7 +122,7 @@ copy(L::Ldiv{<:AbstractBasisLayout,<:AbstractWeightLayout}) = transform_ldiv(L.A
 @inline function _broadcast_mul_ldiv(::Tuple{Any,AbstractBasisLayout}, A, B)
     a,b = arguments(B)
     @assert a isa AbstractQuasiVector # Only works for vec .* mat
-    ab = (A * (A \ a)) .* b # broadcasted should be overloaded
+    ab = expand(a) .* b # broadcasted should be overloaded
     MemoryLayout(ab) isa BroadcastLayout && return transform_ldiv(A, ab)
     A \ ab
 end
@@ -158,7 +158,7 @@ copy(L::Ldiv{<:MappedBasisLayouts,BroadcastLayout{typeof(*)}}) = _broadcast_mul_
     a,b = arguments(B)
     @assert a isa AbstractQuasiVector # Only works for vec .* mat
     A = Ac'
-    ab = (A * (A \ a)) .* b # broadcasted should be overloaded
+    ab = expand(a) .* b # broadcasted should be overloaded
     MemoryLayout(ab) isa BroadcastLayout && return Ac*transform_ldiv(A, ab)
     Ac*ab
 end

test/Project.toml

@@ -0,0 +1 @@
+[deps]

test/test_chebyshev.jl

@@ -1,4 +1,5 @@
-using ContinuumArrays, LinearAlgebra, FastTransforms, QuasiArrays, ArrayLayouts, Base64, LazyArrays, Test
+using ContinuumArrays, LinearAlgebra, QuasiArrays, ArrayLayouts, Base64, LazyArrays, Test
+using FastTransforms
 import ContinuumArrays: Basis, Weight, Map, LazyQuasiArrayStyle, TransformFactorization,
                         ExpansionLayout, checkpoints, MappedBasisLayout, MappedWeightedBasisLayout,
                         SubWeightedBasisLayout, WeightedBasisLayout, WeightLayout, basis, grammatrix
@@ -159,6 +160,8 @@ Base.:(==)(::FooBasis, ::FooBasis) = true
 
         ã = T * (T \ a)
         @test T \ (ã .* ã) ≈ [1.5,1,0.5,0,0]
+
+        @test T'*(a .* T) isa Matrix
     end
 
     @testset "sum/dot/diff" begin