fix honeycomb examples

Author: J-C-Q

examples/Honeycomb/entanglement.jl

@@ -11,9 +11,9 @@ function simulate(path::String; depth=100, L=12, averaging=10, resolution=45, ty
     progressMeter = ProgressMeter.Progress(length(points) * averaging; dt=1.0)
     Threads.@threads for (px, py, pz) in points
         if type == :Kitaev
-            circuit = MeasurementOnlyKitaev(geometry, px, py, pz; depth)
+            circuit = compile(MeasurementOnlyKitaev(geometry, px, py, pz; depth))
         elseif type == :Kekule
-            circuit = MeasurementOnlyKekule(geometry, px, py, pz; depth)
+            circuit = compile(MeasurementOnlyKekule(geometry, px, py, pz; depth))
         else
             throw(ArgumentError("Unsupported type $type. Choose one of :Kitaev, :Kekule"))
         end
@@ -22,9 +22,9 @@ function simulate(path::String; depth=100, L=12, averaging=10, resolution=45, ty
         for _ in 1:averaging
             result = execute(circuit, sim)
             for l in 1:L+1
-                entropy[l] += QuantumClifford.entanglement_entropy(result, z_subsystems[l])
-                entropy[l] += QuantumClifford.entanglement_entropy(result, x_subsystems[l])
-                entropy[l] += QuantumClifford.entanglement_entropy(result, y_subsystems[l])
+                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

examples/Honeycomb/information.jl

@@ -11,19 +11,19 @@ function simulate(path::String; depth=100, L=12, averaging=10, resolution=45, ty
     progressMeter = ProgressMeter.Progress(length(points) * averaging; dt=1.0)
     Threads.@threads for (px, py, pz) in points
         if type == :Kitaev
-            circuit = MeasurementOnlyKitaev(geometry, px, py, pz; depth)
+            circuit = compile(MeasurementOnlyKitaev(geometry, px, py, pz; depth))
         elseif type == :Kekule
-            circuit = MeasurementOnlyKekule(geometry, px, py, pz; depth)
+            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))
         tmi = 0
         for _ in 1:averaging
             result = execute(circuit, sim)
-            tmi += QuantumClifford.tmi(result, z_subsystems)
-            tmi += QuantumClifford.tmi(result, x_subsystems)
-            tmi += QuantumClifford.tmi(result, y_subsystems)
+            tmi += QuantumClifford.tmi(result.stab, z_subsystems)
+            tmi += QuantumClifford.tmi(result.stab, x_subsystems)
+            tmi += QuantumClifford.tmi(result.stab, y_subsystems)
             ProgressMeter.next!(progressMeter)
         end
         tmi /= 3averaging