Example #1
0
"""
Test routines for the functions in 'surface_code.py'

created on: 20/07/17

@author: eduardo
"""
import unittest
import numpy as np
import numpy.testing as nptest
import surface_code

# initialize the surface code object
distance = 6
surface = "toric"
sc = surface_code.SurfaceCode(distance, surface)

sc.select_measurement_protocol(0, 0, "single_rounds")

print("-------------------Testing random noise------------------------")
sc.apply_qubit_error(0.5, 0.5)
print(sc.qubits)

print("-------------------Testing stabilizer measurements-------------")
sc.reset()
sc.apply_qubit_error(.3, .3)
sc.measure_stabilizer_type("star")
sc.measure_stabilizer_type("plaq")
print(sc.tags)
print(sc.qubits)
Example #2
0
pg = 0.008
eta = 0.0
a0 = 0.0
a1 = 0.0
protocol = "LOCAL"
theta = .0
PERFECT_LAST_ROUND = False

p = 0.029
q = 0.029
iterations = 1
cycles = 1

# Initialize objects
fail_rate = 0
sc = surface_code.SurfaceCode(distance, topology)
lc = layers.Layers(sc)
# sc.init_error_obj(topology, ps, pm, pg, eta, a0, a1, theta, protocol)

# Choose a measurement protocol
# sc.select_measurement_protocol(0, 0, "single", 0.0)

# Perform measurements
for i in range(iterations):

    # Errors and measurements
    # Random errors
    if q != 0:
        for t in range(cycles):
            sc.apply_qubit_error(p, 0)
            sc.measure_all_stabilizers()