Introduce GramMatrixPolynomial

dlfivefifty · dlfivefifty · commit 068784e4e863 · 2025-01-21T22:37:23.000Z
diff --git a/test/test_gram.jl b/test/test_gram.jl
@@ -0,0 +1,10 @@
+using ClassicalOrthogonalPolynomials, FastTransforms
+
+P = Legendre()
+x = axes(P,1)
+w = @.(1-x^2)
+
+μ = P'w
+X = jacobimatrix(P)
+n = 20
+@test GramMatrix(μ[1:2n], X[1:2n,1:2n]) ≈ (P' * (w .* P))[1:n,1:n]
