decoding

Author: J-C-Q

examples/Nishimori/NishimorisCat.jl

@@ -1,33 +1,7 @@
 using MonitoredQuantumCircuits
 import JLD2
 import ProgressMeter
-function simulate(path::String; L=12, averaging=10)
-    geometry = HoneycombGeometry(Periodic, L, L)
-
-    progressMeter = ProgressMeter.Progress(length(points) * averaging; dt=1.0)
-    Threads.@threads for (px, py, pz) in points
-        if type == :Kitaev
-            circuit = compile(MeasurementOnlyKitaev(geometry, px, py, pz; depth))
-        elseif type == :Kekule
-            circuit = compile(MeasurementOnlyKekule(geometry, px, py, pz; depth))
-        else
-            throw(ArgumentError("Unsupported type $type. Choose one of :Kitaev, :Kekule"))
-        end
-        sim = QuantumClifford.TableauSimulator(nQubits(geometry))
-        entropy = zeros(Float64, L+1)
-        for _ in 1:averaging
-            result = execute(circuit, sim)
-            for l in 1:L+1
-                entropy[l] += QuantumClifford.entanglement_entropy(result.stab, z_subsystems[l])
-                entropy[l] += QuantumClifford.entanglement_entropy(result.stab, x_subsystems[l])
-                entropy[l] += QuantumClifford.entanglement_entropy(result.stab, y_subsystems[l])
-            end
-            ProgressMeter.next!(progressMeter)
-        end
-        entropy ./= 3averaging
-        JLD2.save(
-            "$path/ENT_L=$(L)_px=$(px)_py=$(py)_pz=$(pz)_averaging=$(averaging)_depth=$(depth).jld2",
-            "entropy", entropy,
-            "probs", (px, py, pz))
-    end
+function simulate(; path="", L=12, averaging=10)
+    geometry = HoneycombGeometry(Open, L, L)
+    circuit = compile(NishimorisCat(geometry))
 end

src/circuit.jl

@@ -194,11 +194,18 @@ struct CompiledCircuit{Ops<:Tuple}
                 pos.ancilla[i] = ancilla[h]
             end
         end
+        nMeasurements = 0
+        for p in pointer
+            instruction = instructions[p]
+            if typeof(instruction) <: MeasurementOperation
+                nMeasurements += 1
+            end
+        end
 
 
 
 
-        new{Ops}(operations, positions, instructions, n_qubits, maximum([a.second for a in ancilla]) - n_qubits, pointer, qubit_map_compiled_to_geometry, qubit_map_geometry_to_compiled, 0)
+        new{Ops}(operations, positions, instructions, n_qubits, maximum([a.second for a in ancilla]) - n_qubits, pointer, qubit_map_compiled_to_geometry, qubit_map_geometry_to_compiled, nMeasurements)
     end
 end
 

src/circuits/nishimoris_cat.jl

@@ -1,4 +1,4 @@
-function NishimorisCat(geometry::HoneycombGeometry{Open}; depth::Integer=100)
+function NishimorisCat(geometry::HoneycombGeometry{Open})
     circuit = Circuit(geometry)
 
     for position in 1:nQubits(geometry)