singularities on onclusion

Author: dlfivefifty

src/ClassicalOrthogonalPolynomials.jl

@@ -167,7 +167,8 @@ singularities(::WeightedOPLayout, a) = singularities(weight(a))
 singularities(w) = singularities(MemoryLayout(w), w)
 singularities(::ExpansionLayout, f) = singularities(basis(f))
 
-singularities(S::SubQuasiArray) = singularities(parent(S))[parentindices(S)[1]]
+singularitiesview(w, ::Inclusion) = w # for now just assume it doesn't change
+singularities(S::SubQuasiArray) = singularitiesview(singularities(parent(S)), parentindices(S)[1])
 
 basis(lay, v, ax::Inclusion{<:Any,<:AbstractInterval}) = _basis(singularities(v))
 

test/test_legendre.jl

@@ -86,10 +86,17 @@ import QuasiArrays: MulQuasiArray
         @test W isa Clenshaw
         @test W * [1; 2; zeros(∞)] ≈ P \ (w .* (P[:,1:2] * [1,2]))
 
+        f = expand(P, exp)
+        @test (w .* f)[0.1] ≈ w[0.1]exp(0.1)
+
         M = P'P
         @test M isa Diagonal
         @test P'x ≈ [0; 2/3; zeros(∞)]
-        @test P'exp.(x) ≈ M * (P\exp.(x))
+        @test P'f ≈ M * (P\f)
+
+        v = f.args[2]
+        @test v'*M isa Adjoint
+        @test v'*M*v ≈ dot(f,f) ≈ f'f
     end
 
     @testset "test on functions" begin