Weighed - mapped OP grammatrix (#159)

Author: dlfivefifty

Project.toml

@@ -1,7 +1,7 @@
 name = "ClassicalOrthogonalPolynomials"
 uuid = "b30e2e7b-c4ee-47da-9d5f-2c5c27239acd"
 authors = ["Sheehan Olver <solver@mac.com>"]
-version = "0.11.8"
+version = "0.11.9"
 
 [deps]
 ArrayLayouts = "4c555306-a7a7-4459-81d9-ec55ddd5c99a"

src/ClassicalOrthogonalPolynomials.jl

@@ -210,6 +210,7 @@ function weightedgrammatrix_layout(::MappedOPLayout, P)
 end
 
 grammatrix_layout(::MappedOPLayout, P) = grammatrix_layout(MappedBasisLayout(), P)
+grammatrix_layout(::WeightedOPLayout{MappedOPLayout}, P) = grammatrix_layout(MappedBasisLayout(), P)
 
 OrthogonalPolynomial(w::Weight) =error("Override for $(typeof(w))")
 

test/test_jacobi.jl

@@ -172,6 +172,12 @@ import ClassicalOrthogonalPolynomials: recurrencecoefficients, basis, MulQuasiMa
             @test (D * (JacobiWeight(c,a) .* P))[0.1,1:5] ≈ ((JacobiWeight(c,a) .* P)[0.1+h,1:5]-(JacobiWeight(c,a) .* P)[0.1,1:5])/h atol=1E-5
             @test (D * (JacobiWeight(c,a) .* P))[0.1,1:5] ≈ (D * (JacobiWeight(c,a) .* Legendre()))[0.1,1:5]
         end
+
+        @testset "grammatrix" begin
+            W = Weighted(jacobi(1,1,0..1))
+            M = grammatrix(W)
+            @test M[1:10,1:10] == grammatrix(Weighted(Jacobi(1,1)))[1:10,1:10]/2
+        end
     end
 
     @testset "functions" begin