From 0a54c1ff8bf8f9462029702985bdc9ae72ec1093 Mon Sep 17 00:00:00 2001 From: Thomas Christensen Date: Fri, 18 Oct 2024 10:41:03 +0200 Subject: [PATCH] fix-up to snf change --- .vendor/SmithNormalForm/src/snf.jl | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/.vendor/SmithNormalForm/src/snf.jl b/.vendor/SmithNormalForm/src/snf.jl index 0d9a5c4a..1865fbf2 100644 --- a/.vendor/SmithNormalForm/src/snf.jl +++ b/.vendor/SmithNormalForm/src/snf.jl @@ -221,7 +221,7 @@ function snf(M::AbstractMatrix{R}; inverse=true) where {R} # Λ′ = Λ*sign(Λ), T′ = sign(Λ)*T, and T⁻¹′ = T⁻¹*sign(Λ), # with the convention that sign(0) = 1. Then we still have that X = SΛT = SΛ′T′ # and also that Λ = S⁻¹XT⁻¹ ⇒ Λ′ = S⁻¹XT⁻¹′. - for j in axes(M, 2) # over column indices + for j in one(T):minimum(size(M)) Λⱼ = D[j,j] if Λⱼ < zero(R) @views V[j,:] .*= -one(R) # T′ = sign(Λ)*T [rows]