add panqec codes and decoders

CondaPkg.toml

@@ -7,4 +7,5 @@ scipy = ""
 pymatching = ""
 ldpc = ">=2.2.8"
 sinter = ">=1.14"
-stim = ">=1.14"
+stim = ">=1.14"
+panqec = ">=0.1.7"

README.md

@@ -99,4 +99,46 @@ julia> syndrome = PyArray(H)*error .% 2
 
 julia> decoding = bpd.decode(np.array(syndrome))
 Python: array([0, 1, 0])
-```
+```
+
+## panqec
+
+The python panqec module is immediately available:
+
+```
+julia> using PyQDecoders
+
+julia> PyQDecoders.panqec
+Python: <module 'panqec' from ...>
+```
+
+Running the example from `panqec` [original tutorial](https://panqec.readthedocs.io/en/latest/tutorials/Panqec%20basics.html):
+
+```
+julia> using PyQDecoders: panqec, panqecdecoders
+
+julia> Toric2DCode = panqec.codes.Toric2DCode;
+
+julia> PauliErrorModel = panqec.error_models.PauliErrorModel;
+
+julia> MatchingDecoder = panqecdecoders.MatchingDecoder;
+
+julia> code = Toric2DCode(4)
+
+julia> error_model = PauliErrorModel(0.2, 0.3, 0.5)
+
+julia> p = 0.1
+
+julia> decoder = MatchingDecoder(code, error_model, p)
+
+julia> errors = error_model.generate(code, p)
+
+julia> syndrome = code.measure_syndrome(errors)
+
+julia> println("Errors: ", errors)
+Errors: [0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0
+ 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]
+
+julia> println("Syndrome: ", syndrome)
+Syndrome: [0 1 0 0 0 1 0 0 0 0 1 0 0 0 1 0 1 1 0 0 0 0 0 0 0 0 0 0 1 1 0 0]
+```

src/PyQDecoders.jl

@@ -7,6 +7,9 @@ const np = PythonCall.pynew()
 const pm = PythonCall.pynew()
 const ldpc = PythonCall.pynew()
 const ldpccodes = PythonCall.pynew()
+const panqec = PythonCall.pynew()
+const panqeccodes = PythonCall.pynew()
+const panqecdecoders = PythonCall.pynew()
 
 function __init__()
     PythonCall.pycopy!(sp, PythonCall.pyimport("scipy"))
@@ -15,6 +18,9 @@ function __init__()
     PythonCall.pycopy!(pm, PythonCall.pyimport("pymatching"))
     PythonCall.pycopy!(ldpc, PythonCall.pyimport("ldpc"))
     PythonCall.pycopy!(ldpccodes, PythonCall.pyimport("ldpc.codes"))
+    PythonCall.pycopy!(panqec, PythonCall.pyimport("panqec"))
+    PythonCall.pycopy!(panqeccodes, PythonCall.pyimport("panqec.codes"))
+    PythonCall.pycopy!(panqecdecoders, PythonCall.pyimport("panqec.decoders"))
 end
 
 end # module