swap W indexing to make more sense

JuliaApproximation · Nov 4, 2024 · 7d91fc3 · 7d91fc3
1 parent 4947845
commit 7d91fc3
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/GramMatrix.jl b/src/GramMatrix.jl
@@ -59,7 +59,7 @@ function compute_skew_generators(W::GramMatrix{T}) where T
     n = size(W, 1)
     G = zeros(T, n, 2)
     G[n, 1] = one(T)
-    G[:, 2] .= W[n-1, :]*X[n-1, n] + W[n, :]*X[n, n] - X'W[:, n]
+    G[:, 2] .= W[:, n-1]*X[n-1, n] + W[:, n]*X[n, n] - X'W[:, n]
     return G
 end