Stieltjes points legendre (#15)

* Update stieltjes.jl * Update test_hilbert.jl
JuliaApproximation · Jun 13, 2023 · 253297c · 253297c
1 parent 1cb02f5
commit 253297c
Show file tree

Hide file tree

Showing 2 changed files with 12 additions and 0 deletions.
diff --git a/src/stieltjes.jl b/src/stieltjes.jl
@@ -196,6 +196,11 @@ end
     transpose(RecurrenceArray(z, (A,B,C), data))
 end
 
+@simplify function *(S::StieltjesPoints, P::Legendre)
+    S * Weighted(P)
+end
+
+
 
 ##
 # mapped

diff --git a/test/test_hilbert.jl b/test/test_hilbert.jl
@@ -138,6 +138,13 @@ end
         f = wT * (T \ exp.(x))
         @test inv.(z .- x') * f ≈ [2.8826861116458593, 1.6307809018753612]
     end
+
+    @testset "StieltjesPoints Legendre" begin
+        z = [2.,3.]
+        P = Legendre()
+        x = axes(P,1)
+        @test (inv.(z .- x') * P)[:,1:5] ≈ [(inv.(2 .- x')*P)[1:5]'; (inv.(3 .- x')*P)[1:5]']
+    end
 end