From 3c796e77595cde44f652630367b117f7d432f3f8 Mon Sep 17 00:00:00 2001
From: "github-actions[bot]"
 <41898282+github-actions[bot]@users.noreply.github.com>
Date: Fri, 26 Jan 2024 15:27:12 -0500
Subject: [PATCH] CompatHelper: bump compat for Nemo to 0.40, (keep existing
 compat) (#228)

* CompatHelper: bump compat for Nemo to 0.40, (keep existing compat)


---------

Co-authored-by: CompatHelper Julia <compathelper_noreply@julialang.org>
Co-authored-by: Stefan Krastanov <github.acc@krastanov.org>
Co-authored-by: Stefan Krastanov <stefan@krastanov.org>
---
 Project.toml        | 6 +++---
 src/entanglement.jl | 2 +-
 src/randoms.jl      | 6 ++----
 test/test_jet.jl    | 2 +-
 4 files changed, 7 insertions(+), 9 deletions(-)

diff --git a/Project.toml b/Project.toml
index 0e2594548..9e92987d6 100644
--- a/Project.toml
+++ b/Project.toml
@@ -33,12 +33,12 @@ QuantumOpticsBase = "4f57444f-1401-5e15-980d-4471b28d5678"
 
 [extensions]
 QuantumCliffordGPUExt = "CUDA"
+QuantumCliffordLDPCDecodersExt = "LDPCDecoders"
 QuantumCliffordMakieExt = "Makie"
 QuantumCliffordPlotsExt = "Plots"
+QuantumCliffordPyQDecodersExt = "PyQDecoders"
 QuantumCliffordQOpticsExt = "QuantumOpticsBase"
 QuantumCliffordQuantikzExt = "Quantikz"
-QuantumCliffordLDPCDecodersExt = "LDPCDecoders"
-QuantumCliffordPyQDecodersExt = "PyQDecoders"
 
 [compat]
 CUDA = "4.4.0"
@@ -53,7 +53,7 @@ LDPCDecoders = "0.1.1"
 LinearAlgebra = "1.9"
 MacroTools = "0.5.9"
 Makie = "0.20"
-Nemo = "0.38"
+Nemo = "0.38, 0.39, 0.40"
 Plots = "1.38.0"
 PrecompileTools = "1.2"
 PyQDecoders = "0.1.1"
diff --git a/src/entanglement.jl b/src/entanglement.jl
index fafeec4dd..31e7c3e8c 100644
--- a/src/entanglement.jl
+++ b/src/entanglement.jl
@@ -201,7 +201,7 @@ function entanglement_entropy(state::AbstractStabilizer, subsystem::AbstractVect
     graph = Graphs.Graph(state)
     adjmat = Graphs.adjacency_matrix(graph)
     other_subsystem = filter(i->!(i in collect(subsystem)), 1:Graphs.nv(graph))
-    subadjmat = Nemo.matrix(Nemo.residue_ring(Nemo.ZZ, 2), collect(adjmat[subsystem,other_subsystem]))
+    subadjmat = Nemo.matrix(Nemo.GF(2), collect(adjmat[subsystem,other_subsystem]))
     LinearAlgebra.rank(subadjmat)
 end
 
diff --git a/src/randoms.jl b/src/randoms.jl
index 1af3dd392..0d1ec09a8 100644
--- a/src/randoms.jl
+++ b/src/randoms.jl
@@ -167,10 +167,8 @@ function precise_inv(a)::Matrix{UInt8}
 end
 
 function nemo_inv(a, n)::Matrix{UInt8}
-    binaryring = Nemo.residue_ring(Nemo.ZZ, 2) # TODO should I use GF(2) instead of ResidueRing(ZZ, 2)?
-    M = Nemo.matrix_space(binaryring, n, n)
-    inverted = inv(M(Matrix{Int}(a))) # Nemo is very picky about input data types
-    return (x->mod(UInt8(x.data),0x2)).(inverted)
+    inverted = inv(Nemo.matrix(Nemo.GF(2),a))
+    return collect(UInt8.(inverted.==1)) # maybe there is a better way to do the conversion
 end
 
 """Sample (h, S) from the distribution P_n(h, S) from Bravyi and Maslov Algorithm 1."""
diff --git a/test/test_jet.jl b/test/test_jet.jl
index bad6302f6..131ffa8eb 100644
--- a/test/test_jet.jl
+++ b/test/test_jet.jl
@@ -36,5 +36,5 @@ end
     )
     @show rep
     @test_broken length(JET.get_reports(rep)) == 0
-    @test length(JET.get_reports(rep)) <= 20
+    @test length(JET.get_reports(rep)) <= 25
 end