pennylane-qrack 0.10.0__tar.gz → 0.10.2__tar.gz

{pennylane_qrack-0.10.0 → pennylane_qrack-0.10.2}/Makefile

@@ -24,7 +24,7 @@ help:
 build-deps:
 ifneq ($(OS),Windows_NT)
 ifeq ($(QRACK_PRESENT),)
-	git clone https://github.com/unitaryfund/qrack.git; cd qrack; git checkout 61a9f1e4694eab612b5032e884e6f1696842e47d; cd ..
+	git clone https://github.com/unitaryfund/qrack.git; cd qrack; git checkout 0dd1e6ce371b69e21bfa33003bc678b3d42713d5; cd ..
 endif
 	mkdir -p qrack/build
 ifeq ($(UNAME_S),Linux)

{pennylane_qrack-0.10.0/pennylane_qrack.egg-info → pennylane_qrack-0.10.2}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: pennylane-qrack
-Version: 0.10.0
+Version: 0.10.2
 Summary: PennyLane plugin for Qrack.
 Home-page: http://github.com/vm6502q
 Maintainer: vm6502q

{pennylane_qrack-0.10.0 → pennylane_qrack-0.10.2}/catalyst/runtime/include/QuantumDevice.hpp

@@ -17,6 +17,7 @@
 #include <complex>
 #include <memory>
 #include <optional>
+#include <random>
 #include <vector>
 
 #include "DataView.hpp"
@@ -103,6 +104,26 @@ struct QuantumDevice {
      */
     [[nodiscard]] virtual auto GetDeviceShots() const -> size_t = 0;
 
+    /**
+     * @brief Set the PRNG of the device.
+     *
+     * The Catalyst runtime enables seeded program execution on non-hardware devices.
+     * A random number generator instance is managed by the runtime to predictably
+     * generate results for non-deterministic programs, such as those involving `Measure`
+     * calls.
+     * Devices implementing support for this feature do not need to use the provided
+     * PRNG instance as their sole source of random numbers, but it is expected that the
+     * the same instance state will predictable and reproducibly generate the same
+     * program results. It is also expected that the provided PRNG state is evolved
+     * sufficiently so that two device executions sharing the same instance do not produce
+     * identical results.
+     * The provided PRNG instance is not thread-locked, and devices wishing to share it
+     * across threads will need to provide their own thread-safety.
+     *
+     * @param gen The std::mt19937 PRNG object.
+     */
+    virtual void SetDevicePRNG([[maybe_unused]] std::mt19937 *gen){};
+
     /**
      * @brief Start recording a quantum tape if provided.
      *
@@ -138,6 +159,28 @@ struct QuantumDevice {
      */
     virtual void PrintState() = 0;
 
+    /**
+     * @brief Prepare subsystems using the given ket vector in the computational basis.
+     *
+     * @param state A state vector of size 2**len(wires)
+     * @param wires The wire(s) the operation acts on
+     */
+    virtual void SetState(DataView<std::complex<double>, 1> &, std::vector<QubitIdType> &)
+    {
+        RT_FAIL("Unsupported functionality");
+    }
+
+    /**
+     * @brief Prepares a single computational basis state.
+     *
+     * @param n Prepares the basis state |n>, where n is an array of integers from the set {0, 1}
+     * @param wires The wire(s) the operation acts on
+     */
+    virtual void SetBasisState(DataView<int8_t, 1> &, std::vector<QubitIdType> &)
+    {
+        RT_FAIL("Unsupported functionality");
+    }
+
     /**
      * @brief Apply a single gate to the state vector of a device with its name if this is
      * supported.

{pennylane_qrack-0.10.0 → pennylane_qrack-0.10.2}/catalyst/runtime/include/RuntimeCAPI.h

@@ -24,7 +24,7 @@ extern "C" {
 
 // Quantum Runtime Instructions
 void __catalyst__rt__fail_cstr(const char *);
-void __catalyst__rt__initialize();
+void __catalyst__rt__initialize(uint32_t *seed);
 void __catalyst__rt__device_init(int8_t *, int8_t *, int8_t *);
 void __catalyst__rt__device_release();
 void __catalyst__rt__finalize();
@@ -32,6 +32,7 @@ void __catalyst__rt__toggle_recorder(bool);
 void __catalyst__rt__print_state();
 void __catalyst__rt__print_tensor(OpaqueMemRefT *, bool);
 void __catalyst__rt__print_string(char *);
+void __catalyst__rt__assert_bool(bool, char *);
 int64_t __catalyst__rt__array_get_size_1d(QirArray *);
 int8_t *__catalyst__rt__array_get_element_ptr_1d(QirArray *, int64_t);
 
@@ -47,6 +48,8 @@ RESULT *__catalyst__rt__result_get_one();
 RESULT *__catalyst__rt__result_get_zero();
 
 // Quantum Gate Set Instructions
+void __catalyst__qis__SetState(MemRefT_CplxT_double_1d *, uint64_t, ...);
+void __catalyst__qis__SetBasisState(MemRefT_int8_1d *, uint64_t, ...);
 void __catalyst__qis__Identity(QUBIT *, const Modifiers *);
 void __catalyst__qis__PauliX(QUBIT *, const Modifiers *);
 void __catalyst__qis__PauliY(QUBIT *, const Modifiers *);

{pennylane_qrack-0.10.0 → pennylane_qrack-0.10.2}/catalyst/runtime/include/Types.h

@@ -126,6 +126,15 @@ struct MemRefT_int64_1d {
     size_t strides[1];
 };
 
+// MemRefT<int64_t, dimension=1> type
+struct MemRefT_int8_1d {
+    int8_t *data_allocated;
+    int8_t *data_aligned;
+    size_t offset;
+    size_t sizes[1];
+    size_t strides[1];
+};
+
 // PairT<MemRefT<double, dimension=1>, MemRefT<int64, dimension=2>> type
 struct PairT_MemRefT_double_int64_1d {
     struct MemRefT_double_1d first;

{pennylane_qrack-0.10.0 → pennylane_qrack-0.10.2}/pennylane_qrack/_version.py

@@ -16,4 +16,4 @@
    Version number (major.minor.patch[-label])
 """
 
-__version__ = "0.10.0"
+__version__ = "0.10.2"

{pennylane_qrack-0.10.0 → pennylane_qrack-0.10.2}/pennylane_qrack/qrack_device.py

@@ -24,7 +24,7 @@ import itertools as it
 
 import numpy as np
 
-from pennylane import QubitDevice, DeviceError
+from pennylane import QubitDevice, DeviceError, QuantumFunctionError
 from pennylane.ops import (
     QubitStateVector,
     BasisState,
@@ -175,31 +175,34 @@ class QrackDevice(QubitDevice):
 
     def __init__(self, wires=0, shots=None, **kwargs):
         options = dict(kwargs)
-        if 'isStabilizerHybrid' in options:
-            self.isStabilizerHybrid = options['isStabilizerHybrid']
-        if 'isTensorNetwork' in options:
-            self.isTensorNetwork = options['isTensorNetwork']
-        if 'isSchmidtDecompose' in options:
-            self.isSchmidtDecompose = options['isSchmidtDecompose']
-        if 'isBinaryDecisionTree' in options:
-            self.isBinaryDecisionTree = options['isBinaryDecisionTree']
-        if 'isOpenCL' in options:
-            self.isOpenCL = options['isOpenCL']
-        if 'isHostPointer' in options:
-            self.isHostPointer = options['isHostPointer']
-        if 'noise' in options:
-            self.noise = options['noise']
+        if "isStabilizerHybrid" in options:
+            self.isStabilizerHybrid = options["isStabilizerHybrid"]
+        if "isTensorNetwork" in options:
+            self.isTensorNetwork = options["isTensorNetwork"]
+        if "isSchmidtDecompose" in options:
+            self.isSchmidtDecompose = options["isSchmidtDecompose"]
+        if "isBinaryDecisionTree" in options:
+            self.isBinaryDecisionTree = options["isBinaryDecisionTree"]
+        if "isOpenCL" in options:
+            self.isOpenCL = options["isOpenCL"]
+        if "isHostPointer" in options:
+            self.isHostPointer = options["isHostPointer"]
+        if "noise" in options:
+            self.noise = options["noise"]
+            if (self.noise != 0) and (shots is None):
+                raise ValueError("Shots must be finite for noisy simulation (not analytical mode).")
         super().__init__(wires=wires, shots=shots)
         self._state = QrackSimulator(
             self.num_wires,
-            isStabilizerHybrid = self.isStabilizerHybrid,
+            isStabilizerHybrid=self.isStabilizerHybrid,
             isTensorNetwork=self.isTensorNetwork,
             isSchmidtDecompose=self.isSchmidtDecompose,
             isBinaryDecisionTree=self.isBinaryDecisionTree,
             isOpenCL=self.isOpenCL,
             isHostPointer=self.isHostPointer,
-            noise = self.noise
+            noise=self.noise,
         )
+        self._circuit = []
 
     def _reverse_state(self):
         end = self.num_wires - 1
@@ -216,7 +219,14 @@ class QrackDevice(QubitDevice):
             operations (List[pennylane.Operation]): operations to be applied
         """
 
-        for i, op in enumerate(operations):
+        self._circuit = self._circuit + operations
+        if self.noise == 0:
+            self._apply()
+            self._circuit = []
+        # else: Defer application until shots or expectation values are requested
+
+    def _apply(self):
+        for op in self._circuit:
             if isinstance(op, QubitStateVector):
                 self._apply_qubit_state_vector(op)
             elif isinstance(op, BasisState):
@@ -647,20 +657,36 @@ class QrackDevice(QubitDevice):
         # estimate the ev
         return np.mean(self.sample(observable))
 
+    def _generate_sample(self):
+        rev_sample = self._state.m_all()
+        sample = 0
+        for i in range(self.num_wires):
+            if (rev_sample & (1 << i)) > 0:
+                sample |= 1 << (self.num_wires - (i + 1))
+        return sample
+
     def generate_samples(self):
         if self.shots is None:
-            raise qml.QuantumFunctionError(
+            raise QuantumFunctionError(
                 "The number of shots has to be explicitly set on the device "
                 "when using sample-based measurements."
             )
 
+        if self.noise != 0:
+            samples = []
+            for _ in range(self.shots):
+                self._state.reset_all()
+                self._apply()
+                samples.append(self._generate_sample())
+            self._samples = QubitDevice.states_to_binary(np.array(samples), self.num_wires)
+            self._circuit = []
+
+            return self._samples
+
         if self.shots == 1:
-            rev_sample = self._state.m_all()
-            sample = 0
-            for i in range(self.num_wires):
-                if (rev_sample & (1 << i)) > 0:
-                    sample |= 1 << (self.num_wires - (i + 1))
-            self._samples = QubitDevice.states_to_binary(np.array([sample]), self.num_wires)
+            self._samples = QubitDevice.states_to_binary(
+                np.array([self._generate_sample()]), self.num_wires
+            )
 
             return self._samples
 

{pennylane_qrack-0.10.0 → pennylane_qrack-0.10.2/pennylane_qrack.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: pennylane-qrack
-Version: 0.10.0
+Version: 0.10.2
 Summary: PennyLane plugin for Qrack.
 Home-page: http://github.com/vm6502q
 Maintainer: vm6502q

{pennylane_qrack-0.10.0 → pennylane_qrack-0.10.2}/setup.py

@@ -19,10 +19,11 @@ import subprocess
 from setuptools import setup
 from setuptools.command.build_py import build_py
 
+
 class Build(build_py):
     def run(self):
         protoc_command = ["make", "build-deps"]
-        if os.name != 'nt':
+        if os.name != "nt":
             if subprocess.call(protoc_command) != 0:
                 sys.exit(-1)
         super().run()
@@ -31,11 +32,7 @@ class Build(build_py):
 with open("./pennylane_qrack/_version.py") as f:
     (version,) = re.findall('__version__ = "(.*)"', f.read())
 
-requirements = [
-    "pennylane>=0.32",
-    "pyqrack>=1.30.0",
-    "numpy>=1.16"
-]
+requirements = ["pennylane>=0.32", "pyqrack>=1.30.0", "numpy>=1.16"]
 
 info = {
     "name": "pennylane-qrack",
@@ -45,11 +42,9 @@ info = {
     "url": "http://github.com/vm6502q",
     "license": "Apache License 2.0",
     "packages": ["pennylane_qrack"],
-    "cmdclass": { 'build_py': Build },
+    "cmdclass": {"build_py": Build},
     "entry_points": {
-        "pennylane.plugins": [
-            "qrack.simulator = pennylane_qrack.qrack_device:QrackDevice"
-        ]
+        "pennylane.plugins": ["qrack.simulator = pennylane_qrack.qrack_device:QrackDevice"]
     },
     "description": "PennyLane plugin for Qrack.",
     "long_description": open("README.rst").read(),

{pennylane_qrack-0.10.0 → pennylane_qrack-0.10.2}/tests/test_apply.py

@@ -260,9 +260,7 @@ class TestStateApply:
 
     @pytest.mark.parametrize("device_wires", [3, 4, 5])
     @pytest.mark.parametrize("op_wires", [[0], [2], [0, 1], [1, 0], [2, 0]])
-    def test_qubit_state_vector_on_wires_subset(
-        self, init_state, device_wires, op_wires, tol
-    ):
+    def test_qubit_state_vector_on_wires_subset(self, init_state, device_wires, op_wires, tol):
         """Test QubitStateVector application on a subset of device wires"""
         dev = QrackDevice(device_wires, isOpenCL=False)
         state = init_state(len(op_wires))
@@ -314,7 +312,7 @@ class TestStateApply:
             qml.StatePrep(state, wires=[0])
             qml.apply(op)
             return qml.probs()
-        
+
         dev._obs_queue = []
 
         res = circuit()
@@ -352,7 +350,7 @@ class TestStateApply:
             qml.StatePrep(state, wires=[0])
             qml.apply(op)
             return qml.probs()
-        
+
         dev._obs_queue = []
 
         res = circuit()
@@ -364,9 +362,7 @@ class TestStateApply:
     @pytest.mark.parametrize("theta", [0.5432, -0.232])
     @pytest.mark.parametrize("omega", [1.213, -0.221])
     @pytest.mark.parametrize("op,func", single_qubit_three_param)
-    def test_single_qubit_three_parameters(
-        self, init_state, op, func, phi, theta, omega, tol
-    ):
+    def test_single_qubit_three_parameters(self, init_state, op, func, phi, theta, omega, tol):
         """Test PauliX application"""
         dev = QrackDevice(1, isOpenCL=False)
         state = init_state(1)
@@ -383,14 +379,13 @@ class TestStateApply:
     @pytest.mark.parametrize("theta", [0.5432, -0.232])
     @pytest.mark.parametrize("omega", [1.213, -0.221])
     @pytest.mark.parametrize("op,func", single_qubit_three_param)
-    def test_single_qubit_three_parameters_qjit(
-        self, init_state, op, func, phi, theta, omega, tol
-    ):
+    def test_single_qubit_three_parameters_qjit(self, init_state, op, func, phi, theta, omega, tol):
         """Test PauliX application"""
         dev = QrackDevice(1, isOpenCL=False)
         state = init_state(1)
 
         op.data = [phi, theta, omega]
+
         @qjit
         @qml.qnode(dev)
         def circuit():
@@ -430,6 +425,7 @@ class TestStateApply:
             qml.StatePrep(state, wires=[0, 1])
             qml.apply(op)
             return qml.probs()
+
         dev._obs_queue = []
 
         res = circuit()
@@ -462,12 +458,14 @@ class TestStateApply:
         state = init_state(N)
 
         op = qml.QubitUnitary(mat, wires=list(range(N)))
+
         @qjit
         @qml.qnode(dev)
         def circuit():
             qml.StatePrep(state, wires=list(range(N)))
             qml.apply(op)
             return qml.probs()
+
         dev._obs_queue = []
 
         res = circuit()
@@ -506,6 +504,7 @@ class TestStateApply:
             qml.StatePrep(state, wires=[0, 1, 2])
             qml.apply(op)
             return qml.probs()
+
         dev._obs_queue = []
 
         res = circuit()
@@ -549,6 +548,7 @@ class TestStateApply:
         state = init_state(2)
 
         op.data = [theta]
+
         @qjit
         @qml.qnode(dev)
         def circuit():
@@ -567,9 +567,7 @@ class TestStateApply:
     @pytest.mark.parametrize("theta", [0.5432, -0.232])
     @pytest.mark.parametrize("omega", [1.213, -0.221])
     @pytest.mark.parametrize("op,func", two_qubit_three_param)
-    def test_two_qubit_three_parameters(
-        self, init_state, op, func, phi, theta, omega, tol
-    ):
+    def test_two_qubit_three_parameters(self, init_state, op, func, phi, theta, omega, tol):
         """Test parametrized two qubit gates application"""
         dev = QrackDevice(2, isOpenCL=False)
         state = init_state(2)
@@ -588,9 +586,7 @@ class TestStateApply:
         """Test that apply fails for incorrect state preparation."""
         dev = QrackDevice(2, isOpenCL=False)
 
-        with pytest.raises(
-            ValueError, match="Sum of amplitudes-squared does not equal one."
-        ):
+        with pytest.raises(ValueError, match="Sum of amplitudes-squared does not equal one."):
             dev.apply([qml.QubitStateVector(np.array([1, -1]), wires=[0])])
 
         with pytest.raises(

{pennylane_qrack-0.10.0 → pennylane_qrack-0.10.2}/tests/test_integration.py

@@ -17,6 +17,7 @@ import pennylane as qml
 from pennylane_qrack.qrack_device import QrackDevice
 from catalyst import qjit
 
+
 class TestIntegration:
     """Some basic integration tests."""