Merge branch 'dl/clenshawargorder' of https://github.com/JuliaApproximation/FastTransforms.jl into dl/clenshawargorder

Author: dlfivefifty

Project.toml

@@ -2,6 +2,7 @@ name = "FastTransforms"
 uuid = "057dd010-8810-581a-b7be-e3fc3b93f78c"
 version = "0.17"
 
+
 [deps]
 AbstractFFTs = "621f4979-c628-5d54-868e-fcf4e3e8185c"
 ArrayLayouts = "4c555306-a7a7-4459-81d9-ec55ddd5c99a"

src/GramMatrix.jl

@@ -64,7 +64,7 @@ function GramMatrix(μ::AbstractVector{T}, X::XT, p0::T) where {T, XT <: Abstrac
     n = (N+1)÷2
     @assert N == size(X, 1) == size(X, 2)
     @assert bandwidths(X) == (1, 1)
-    W = Matrix{T}(undef, N, N)
+    W = LowerTriangular(Matrix{T}(undef, N, N))
     if n > 0
         @inbounds for m in 1:N
             W[m, 1] = p0*μ[m]