Attempt to make dimensions consistent

Author: dlfivefifty

src/GramMatrix.jl

@@ -78,12 +78,12 @@ function GramMatrix(μ::AbstractVector{T}, X::XT, p0::T) where {T, XT <: Abstrac
             W[m, 2] = (X[m-1, m]*W[m-1, 1] + (X[m, m]-X[1, 1])*W[m, 1] + X[m+1, m]*W[m+1, 1])/X[2, 1]
         end
     end
-    @inbounds @simd for n in 3:n
-        for m in n:min(N-n+1, b+n)
-            W[m, n] = (X[m-1, m]*W[m-1, n-1] + (X[m, m]-X[n-1, n-1])*W[m, n-1] + X[m+1, m]*W[m+1, n-1] - X[n-2, n-1]*W[m, n-2])/X[n, n-1]
+    @inbounds @simd for j in 3:N
+        for m in j:min(N, b+j)
+            W[m, j] = (X[m-1, m]*W[m-1, j-1] + (X[m, m]-X[j-1, j-1])*W[m, j-1] + (m < N ? X[m+1, m]*W[m+1, j-1] : zero(T)) - X[j-2, j-1]*W[m, j-2])/X[j, j-1]
         end
     end
-    return GramMatrix(Symmetric(W[1:n, 1:n], :L), eval(XT.name.name)(view(X, 1:n, 1:n)))
+    return GramMatrix(Symmetric(W, :L), X)
 end
 
 #

test/grammatrixtests.jl

@@ -29,7 +29,7 @@ using FastTransforms, BandedMatrices, LazyArrays, LinearAlgebra, Test
         F = cholesky(G)
         @test F.L*F.L' ≈ W
 
-        @test G[1:66,1:66] ≈ GramMatrix(W[1:5, 1], X)
+        @test G ≈ GramMatrix(W[1:5, 1], X̃)
 
         X = BandedMatrix(SymTridiagonal(T[2n-1 for n in 1:n+b], T[-n for n in 1:n+b-1])) # Laguerre X, tests nonzero diagonal
         W = I+X^2+X^4