Implement multithreaded stochastic swap in rust

.pylintrc

@@ -33,7 +33,7 @@ unsafe-load-any-extension=no
 # A comma-separated list of package or module names from where C extensions may
 # be loaded. Extensions are loading into the active Python interpreter and may
 # run arbitrary code
-extension-pkg-allow-list=retworkx, numpy, tweedledum
+extension-pkg-allow-list=retworkx, numpy, tweedledum, qiskit._accelerate
 
 
 [MESSAGES CONTROL]

Cargo.toml

@@ -0,0 +1,33 @@
+[package]
+name = "qiskit-terra"
+version = "0.20.0"
+edition = "2018"
+
+# See more keys and their definitions at https://doc.rust-lang.org/cargo/reference/manifest.html
+
+[lib]
+name = "qiskit_accelerate"
+crate-type = ["cdylib"]
+
+[dependencies]
+rayon = "1.5"
+numpy = "0.15.1"
+rand = "0.8"
+rand_pcg = "0.3"
+rand_distr = "0.4.3"
+
+[dependencies.pyo3]
+version = "0.15.1"
+features = ["extension-module", "hashbrown"]
+
+[dependencies.ndarray]
+version = "^0.15.0"
+features = ["rayon"]
+
+[dependencies.hashbrown]
+version = "0.11.2"
+features = ["rayon"]
+
+[profile.release]
+lto = 'fat'
+codegen-units = 1

MANIFEST.in

@@ -15,3 +15,6 @@ include test/python/pickles/*.pickle
 include test/python/qasm/*.qasm
 include test/python/visualization/references/*.png
 include test/python/notebooks/*.ipynb
+
+include Cargo.toml
+recursive-include src *

azure-pipelines.yml

@@ -286,6 +286,8 @@ stages:
             tools/verify_headers.py qiskit test
             python tools/find_optional_imports.py
             reno lint
+            cargo fmt --check
+            cargo clippy -- -D warnings
           displayName: 'Style and lint'
     - job: 'Docs'
       pool: {vmImage: 'ubuntu-latest'}

examples/python/stochastic_swap.py

@@ -73,23 +73,22 @@
 # Build the expected output to verify the pass worked
 expected = QuantumCircuit(qr, cr)
 expected.cx(qr[1], qr[2])
+expected.h(qr[2])
 expected.swap(qr[0], qr[1])
+expected.h(qr[0])
 expected.cx(qr[1], qr[3])
 expected.h(qr[3])
-expected.h(qr[2])
 expected.measure(qr[1], cr[0])
-expected.h(qr[0])
 expected.swap(qr[1], qr[3])
-expected.h(qr[3])
 expected.cx(qr[2], qr[1])
+expected.h(qr[3])
+expected.swap(qr[0], qr[1])
 expected.measure(qr[2], cr[2])
-expected.swap(qr[1], qr[3])
-expected.measure(qr[3], cr[3])
-expected.cx(qr[1], qr[0])
-expected.measure(qr[1], cr[0])
-expected.measure(qr[0], cr[1])
+expected.cx(qr[3], qr[1])
+expected.measure(qr[0], cr[3])
+expected.measure(qr[3], cr[0])
+expected.measure(qr[1], cr[1])
 expected_dag = circuit_to_dag(expected)
-
 # Run the pass on the dag from the input circuit
 pass_ = StochasticSwap(coupling, 20, 999)
 after = pass_.run(dag)

pyproject.toml

@@ -1,5 +1,6 @@
 [build-system]
-requires = ["Cython>=0.27.1", "setuptools", "wheel"]
+requires = ["Cython>=0.27.1", "setuptools", "wheel", "setuptools-rust"]
+build-backend = "setuptools.build_meta"
 
 [tool.black]
 line-length = 100