stim_ghz_state_circuit_simulation_and_detector_analysis.py

1import stim
2
3# Create a circuit that produces a GHZ state and measures it.
4circuit = stim.Circuit("""
5    H 0
6    CNOT 0 1
7    CNOT 1 2
8    M 0 1 2
9""")
10
11# Simulate the circuit 10 times.
12# Each sample is a bit-packed numpy array.
13sampler = circuit.compile_sampler()
14samples = sampler.sample(shots=10)
15print(samples)
16# [[0 0 0]
17#  [1 1 1]
18#  ...
19#  [0 0 0]]
20
21# We can also analyze the circuit's error properties.
22# For example, what are the possible detection events?
23# (In this case, there are no error sources, so no detection events.)
24dem = circuit.detector_error_model()
25print(dem)