pete
/
abp
1
1
Fork 0
Code Issues 0 Pull Requests 0 Releases 0 Wiki Activity
Anders and Briegel in Python
You can not select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
104 Commits
1 Branch
24MB
Tree: c1fe5f8b79
Branches Tags
${ item.name }
Create branch ${ searchTerm }
from 'c1fe5f8b79'
${ noResults }
abp/tests/test_against_anders_and_bri...

108 lines
2.2KB
Raw Blame History 

  1. from abp.graphstate import GraphState

  2. from anders_briegel import graphsim

  3. import random

  4. import difflib

  5. import re

  6. 

  7. def compare(a, b):

  8.     """ Sketchy as you like. Remove this abomination """

  9.     aa = a.get_adj_list()

  10.     bb = b.adj_list()

  11.     try:

  12.         assert re.sub("\\s", "", aa) == re.sub("\\s", "", bb)

  13.     except AssertionError:

  14.         print aa

  15.         print bb

  16. 

  17. def test_hadamard():

  18.     """ Test hadamards """

  19.     a = graphsim.GraphRegister(1)

  20.     b = GraphState()

  21.     b.add_vertex(0)

  22. 

  23.     compare(a, b)

  24.     a.hadamard(0)

  25.     b.act_hadamard(0)

  26.     compare(a, b)

  27.     a.hadamard(0)

  28.     b.act_hadamard(0)

  29.     compare(a, b)

  30. 

  31. 

  32. def test_local_rotations():

  33.     """ Test local rotations """

  34.     a = graphsim.GraphRegister(1)

  35.     b = GraphState()

  36.     b.add_vertex(0)

  37.     compare(a, b)

  38. 

  39.     for i in range(1000):

  40.         j = random.randint(0, 23)

  41.         a.local_op(0, graphsim.LocCliffOp(j))

  42.         b.act_local_rotation(0, j)

  43.         compare(a, b)

  44. 

  45. 

  46. def test_cz_table():

  47.     """ Test the CZ table """

  48.     for j in range(24):

  49.         a = graphsim.GraphRegister(2)

  50.         b = GraphState()

  51.         b.add_vertex(0)

  52.         b.add_vertex(1)

  53.         compare(a, b)

  54.         #a.local_op(0, graphsim.LocCliffOp(j))

  55.         #b.act_local_rotation(0, j)

  56. 

  57.         a.local_op(1, graphsim.LocCliffOp(j))

  58.         b.act_local_rotation(1, j)

  59.         a.cphase(0, 1)

  60.         b.act_cz(0, 1)

  61.         compare(a, b)

  62. 

  63. 

  64. 

  65. def _test_1():

  66.     """ TODO: this one always succeeds """

  67.     N=10

  68. 

  69.     a = graphsim.GraphRegister(N)

  70.     b = GraphState()

  71. 

  72.     for i in range(N):

  73.         a.hadamard(i)

  74.         b.add_vertex(i)

  75.         b.act_hadamard(i)

  76. 

  77.     for i in range(N-1):

  78.         a.cphase(i, i+1)

  79.         b.act_cz(i, i+1)

  80. 

  81.     compare(a, b)

  82. 

  83. 

  84. def _test_2():

  85.     """ TODO: This one fails at the moment """

  86.     N=10

  87. 

  88.     a = graphsim.GraphRegister(N)

  89.     b = GraphState()

  90. 

  91.     for i in range(N):

  92.         b.add_vertex(i)

  93. 

  94.     for i in range(100):

  95.         if random.random()>0.5:

  96.             j = random.randint(0, N-1)

  97.             a.hadamard(j)

  98.             b.act_hadamard(j)

  99.         else:

  100.             q = random.randint(0, N-2)

  101.             a.cphase(q, q+1)

  102.             b.act_cz(q, q+1)

  103. 

  104.     aa = a.get_adj_list()

  105.     bb = b.adj_list()

  106.     compare(a, b)