working

Author: J-C-Q

Project.toml

@@ -4,6 +4,7 @@ authors = ["Quinten Preiss"]
 version = "0.1.0-DEV"
 
 [deps]
+AllocCheck = "9b6a8646-10ed-4001-bbdc-1d2f46dfbb1a"
 CSV = "336ed68f-0bac-5ca0-87d4-7b16caf5d00b"
 Colors = "5ae59095-9a9b-59fe-a467-6f913c188581"
 Combinatorics = "861a8166-3701-5b0c-9a16-15d98fcdc6aa"
@@ -34,6 +35,7 @@ MonitoredQuantumCircuitsMPIExt = "MPI"
 MonitoredQuantumCircuitsMakieExt = "Makie"
 
 [compat]
+AllocCheck = "0.2.2"
 CSV = "0.10.14"
 CUDA = "5.7.2"
 DataFrames = "1.6.1"

examples/2D_GHZ/simulation.jl

@@ -0,0 +1,29 @@
+function simulate_QuantumClifford(;shots=10000,postselect=false)
+    g = IBMQ_Falcon()
+    backend = QuantumClifford.TableauSimulator(nQubits(g);
+        mixed=false, basis=:Z)
+
+    magnetization = Int64[]
+    p = (s) -> begin
+        if postselect && all(i->i==0, (@view backend.measurements[1:11])) || !postselect
+            push!(magnetization, sum(i->2*i-1, (@view backend.measurements[12:end])))
+        end
+    end
+    execute!(()->monitoredGHZ!(backend, g; tApi=1/4), backend, p; shots=shots)
+    return magnetization
+end
+
+function simulate_Qiskit(;shots=10000,postselect=false)
+    g = IBMQ_Falcon()
+    backend = Qiskit.CliffordSimulator(nQubits(g);ancillas=nControlQubits(g))
+
+    magnetization = Int64[]
+    p = (s) -> begin
+        measurements = Qiskit.get_measurements(backend, s)
+        if postselect && all(i->i==0, (@view measurements[1:11])) || !postselect
+            push!(magnetization, sum(i->2*i-1, measurements[12:end]))
+        end
+    end
+    execute!(()->monitoredGHZ!(backend, g; tApi=1/4), backend, p; shots=shots)
+    return magnetization
+end

examples/Honeycomb/simulation.jl

@@ -0,0 +1,36 @@
+function simulate_QuantumClifford(L,n;shots=10000,depth=100)
+    g = HoneycombGeometry(Periodic, L, L)
+    ps = generateProbs(;n)
+    backend = QuantumClifford.TableauSimulator(nQubits(g); mixed=true)
+
+    partitionsX = subsystems(g, 4; cutType=:X)
+    partitionsY = subsystems(g, 4; cutType=:Y)
+    partitionsZ = subsystems(g, 4; cutType=:Z)
+
+    tmi = zeros(Float64, length(ps))
+    for (i, (px,py,pz)) in enumerate(ps)
+        post = (s) -> begin
+            tmi[i] += QuantumClifford.tmi(backend.state, partitionsX)
+            tmi[i] += QuantumClifford.tmi(backend.state, partitionsY)
+            tmi[i] += QuantumClifford.tmi(backend.state, partitionsZ)
+        end
+        execute!(()->measurementOnlyKitaev!(backend, g, px,py,pz; depth), backend, post; shots=shots)
+    end
+    return tmi ./= 3shots
+end
+
+function generateProbs(; n=45)
+    points = Vector{NTuple{3,Float64}}(undef, n*(n + 1) ÷ 2)
+    m = 1
+    for (k, i) in enumerate(range(0, 1, n))
+        for j in range(i, 1, n - k + 1)
+            px = i
+            py = j - i
+            pz = 1 - j
+            points[m] = (px, py, pz)
+            m += 1
+        end
+    end
+
+    return points
+end

examples/TransverseFieldIsing/simulation.jl

@@ -0,0 +1,14 @@
+function simulate_QuantumClifford(L,ps;shots=10000,depth=100)
+    g = ChainGeometry(Periodic, L)
+    backend = QuantumClifford.TableauSimulator(nQubits(g); mixed=false, basis=:X)
+
+    entanglement = zeros(Float64, length(ps))
+    for (i, p) in enumerate(ps)
+        post = (s) -> begin
+            entanglement[i] += QuantumClifford.entanglement_entropy(backend.state, 1:div(L,2))
+        end
+        execute!(()->monitoredTransverseFieldIsing!(backend, g, p; depth), backend, post; shots=shots)
+    end
+
+    return entanglement ./= shots
+end