pyqrack-complex128 1.78.3__py3-none-macosx_14_0_arm64.whl → 1.81.0__py3-none-macosx_14_0_arm64.whl

pyqrack/__init__.py

@@ -13,7 +13,6 @@ from .qrack_neuron_torch_layer import (
     QrackNeuronTorch,
     QrackNeuronTorchFunction,
     QrackNeuronTorchLayer,
-    QrackNeuronTorchLayerFunction,
 )
 from .qrack_simulator import QrackSimulator
 from .qrack_stabilizer import QrackStabilizer

pyqrack/qrack_ace_backend.py

@@ -224,9 +224,7 @@ class QrackAceBackend:
         isPaged=True,
         isCpuGpuHybrid=True,
         isOpenCL=True,
-        isHostPointer=(
-            True if os.environ.get("PYQRACK_HOST_POINTER_DEFAULT_ON") else False
-        ),
+        isHostPointer=(True if os.environ.get("PYQRACK_HOST_POINTER_DEFAULT_ON") else False),
         noise=0,
         toClone=None,
     ):
@@ -294,11 +292,7 @@ class QrackAceBackend:
                     sim_counts[t_sim_id] += 1
 
                     qubit.append(
-                        LHVQubit(
-                            toClone=(
-                                toClone._qubits[tot_qubits][2] if toClone else None
-                            )
-                        )
+                        LHVQubit(toClone=(toClone._qubits[tot_qubits][2] if toClone else None))
                     )
 
                 if (not c) and (not r):
@@ -545,9 +539,7 @@ class QrackAceBackend:
                 self.sim[hq[4][0]].prob(hq[4][1]),
             ]
             # Balancing suggestion from Elara (the custom OpenAI GPT)
-            prms = math.sqrt(
-                (p[0] ** 2 + p[1] ** 2 + 3 * (p[2] ** 2) + p[3] ** 2 + p[4] ** 2) / 7
-            )
+            prms = math.sqrt((p[0] ** 2 + p[1] ** 2 + 3 * (p[2] ** 2) + p[3] ** 2 + p[4] ** 2) / 7)
             qrms = math.sqrt(
                 (
                     (1 - p[0]) ** 2
@@ -860,9 +852,7 @@ class QrackAceBackend:
             gate = self.sim[sim_id].macz if anti else self.sim[sim_id].mcz
             shadow = self._anti_cz_shadow if anti else self._cz_shadow
         else:
-            raise RuntimeError(
-                "QrackAceBackend._get_gate() should never return identity!"
-            )
+            raise RuntimeError("QrackAceBackend._get_gate() should never return identity!")
 
         return gate, shadow
 
@@ -899,11 +889,7 @@ class QrackAceBackend:
                 b2 = hq2[q2]
                 if b1[0] == b2[0]:
                     gate_fn([b1[1]], b2[1])
-                elif (
-                    lq1_lr
-                    or (b1[1] == b2[1])
-                    or ((len(qb1) == 2) and (b1[1] == (b2[1] & 1)))
-                ):
+                elif lq1_lr or (b1[1] == b2[1]) or ((len(qb1) == 2) and (b1[1] == (b2[1] & 1))):
                     shadow_fn(b1, b2)
 
     def _cpauli(self, lq1, lq2, anti, pauli):
@@ -1033,9 +1019,7 @@ class QrackAceBackend:
                 self.sim[hq[4][0]].prob(hq[4][1]),
             ]
             # Balancing suggestion from Elara (the custom OpenAI GPT)
-            prms = math.sqrt(
-                (p[0] ** 2 + p[1] ** 2 + 3 * (p[2] ** 2) + p[3] ** 2 + p[4] ** 2) / 7
-            )
+            prms = math.sqrt((p[0] ** 2 + p[1] ** 2 + 3 * (p[2] ** 2) + p[3] ** 2 + p[4] ** 2) / 7)
             qrms = math.sqrt(
                 (
                     (1 - p[0]) ** 2
@@ -1187,17 +1171,11 @@ class QrackAceBackend:
                 (-1 * float(operation.params[1])) + math.pi / 2,
             )
         elif name == "rx":
-            self._sim.r(
-                Pauli.PauliX, float(operation.params[0]), operation.qubits[0]._index
-            )
+            self._sim.r(Pauli.PauliX, float(operation.params[0]), operation.qubits[0]._index)
         elif name == "ry":
-            self._sim.r(
-                Pauli.PauliY, float(operation.params[0]), operation.qubits[0]._index
-            )
+            self._sim.r(Pauli.PauliY, float(operation.params[0]), operation.qubits[0]._index)
         elif name == "rz":
-            self._sim.r(
-                Pauli.PauliZ, float(operation.params[0]), operation.qubits[0]._index
-            )
+            self._sim.r(Pauli.PauliZ, float(operation.params[0]), operation.qubits[0]._index)
         elif name == "h":
             self._sim.h(operation.qubits[0]._index)
         elif name == "x":
@@ -1262,9 +1240,9 @@ class QrackAceBackend:
                         cregbit = clbit
 
                     regbit = 1 << cregbit
-                    self._classical_register = (
-                        self._classical_register & (~regbit)
-                    ) | (qubit_outcome << cregbit)
+                    self._classical_register = (self._classical_register & (~regbit)) | (
+                        qubit_outcome << cregbit
+                    )
 
         elif name == "bfunc":
             mask = int(operation.mask, 16)
@@ -1379,9 +1357,7 @@ class QrackAceBackend:
             if operation.name == "id" or operation.name == "barrier":
                 continue
 
-            if is_initializing and (
-                (operation.name == "measure") or (operation.name == "reset")
-            ):
+            if is_initializing and ((operation.name == "measure") or (operation.name == "reset")):
                 continue
 
             is_initializing = False
@@ -1439,9 +1415,7 @@ class QrackAceBackend:
                 self._sample_cregbits = []
 
         if self._sample_measure and (len(self._sample_qubits) > 0):
-            _data = self._add_sample_measure(
-                self._sample_qubits, self._sample_clbits, self._shots
-            )
+            _data = self._add_sample_measure(self._sample_qubits, self._sample_clbits, self._shots)
 
         del self._sim
 
@@ -1532,17 +1506,13 @@ class QrackAceBackend:
                 noise_model.add_quantum_error(depolarizing_error(y, 2), "cx", [a, b])
                 noise_model.add_quantum_error(depolarizing_error(y_cy, 2), "cy", [a, b])
                 noise_model.add_quantum_error(depolarizing_error(y_cy, 2), "cz", [a, b])
-                noise_model.add_quantum_error(
-                    depolarizing_error(y_swap, 2), "swap", [a, b]
-                )
+                noise_model.add_quantum_error(depolarizing_error(y_swap, 2), "swap", [a, b])
             else:
                 y_cy = 1 - (1 - y) ** 2
                 y_swap = 1 - (1 - y) ** 3
                 noise_model.add_quantum_error(depolarizing_error(y_cy, 2), "cx", [a, b])
                 noise_model.add_quantum_error(depolarizing_error(y_cy, 2), "cy", [a, b])
                 noise_model.add_quantum_error(depolarizing_error(y_cy, 2), "cz", [a, b])
-                noise_model.add_quantum_error(
-                    depolarizing_error(y_swap, 2), "swap", [a, b]
-                )
+                noise_model.add_quantum_error(depolarizing_error(y_swap, 2), "swap", [a, b])
 
         return noise_model

pyqrack/qrack_circuit.py

@@ -507,9 +507,7 @@ class QrackCircuit:
             )
             if circuit_type == QuimbCircuitType.Circuit
             else (
-                qtn.CircuitDense(
-                    N=qcirc.num_qubits, psi0=psi0, gate_opts=gate_opts, tags=tags
-                )
+                qtn.CircuitDense(N=qcirc.num_qubits, psi0=psi0, gate_opts=gate_opts, tags=tags)
                 if circuit_type == QuimbCircuitType.CircuitDense
                 else qtn.CircuitMPS(
                     N=qcirc.num_qubits,
@@ -540,9 +538,7 @@ class QrackCircuit:
         return tcirc
 
     @staticmethod
-    def file_to_tensorcircuit(
-        filename, inputs=None, circuit_params=None, binding_params=None
-    ):
+    def file_to_tensorcircuit(filename, inputs=None, circuit_params=None, binding_params=None):
         """Convert an output file to a TensorCircuit circuit
 
         Reads in an (optimized) circuit from a file named

pyqrack/qrack_neuron.py

@@ -31,7 +31,9 @@ class QrackNeuron:
         simulator(QrackSimulator): Simulator instance for all synaptic clefts of the neuron
         controls(list(int)): Indices of all "control" qubits, for neuron input
         target(int): Index of "target" qubit, for neuron output
-        tolerance(double): Rounding tolerance
+        activation_fn(NeuronActivationFn): Activation function choice
+        alpha(float): Activation function parameter, if required
+        angles(list[ctypes.c_float]): (or c_double) Memory for neuron prediction angles
     """
 
     def _get_error(self):
@@ -48,7 +50,6 @@ class QrackNeuron:
         target,
         activation_fn=NeuronActivationFn.Sigmoid,
         alpha=1.0,
-        tolerance=sys.float_info.epsilon,
         _init=True,
     ):
         self.simulator = simulator
@@ -56,7 +57,7 @@ class QrackNeuron:
         self.target = target
         self.activation_fn = activation_fn
         self.alpha = alpha
-        self.tolerance = tolerance
+        self.angles = QrackNeuron._real1_byref([0.0] * (1 << len(controls)))
 
         if not _init:
             return
@@ -66,9 +67,6 @@ class QrackNeuron:
             len(controls),
             QrackNeuron._ulonglong_byref(controls),
             target,
-            activation_fn,
-            alpha,
-            tolerance,
         )
 
         self._throw_if_error()
@@ -91,11 +89,9 @@ class QrackNeuron:
             self.simulator,
             self.controls,
             self.target,
-            self.activation_fn,
-            self.alpha,
-            self.tolerance,
         )
-        self.nid = Qrack.qrack_lib.clone_qneuron(self.simulator.sid)
+        result.nid = Qrack.qrack_lib.clone_qneuron(self.simulator.sid)
+        result.angles = self.angles[:]
         self._throw_if_error()
         return result
 
@@ -110,6 +106,27 @@ class QrackNeuron:
             return (ctypes.c_float * len(a))(*a)
         return (ctypes.c_double * len(a))(*a)
 
+    def set_simulator(self, s):
+        """Set the neuron simulator
+
+        Set the simulator used by this neuron
+
+        Args:
+            s(QrackSimulator): The simulator to use
+
+        Raises:
+            RuntimeError: QrackSimulator raised an exception.
+        """
+        Qrack.qrack_lib.set_qneuron_sim(
+            self.nid,
+            s.sid,
+            len(self.controls),
+            QrackNeuron._ulonglong_byref(self.controls),
+            self.target,
+        )
+        self._throw_if_error()
+        self.simulator = s
+
     def set_angles(self, a):
         """Directly sets the neuron parameters.
 
@@ -127,8 +144,7 @@ class QrackNeuron:
             raise ValueError(
                 "Angles 'a' in QrackNeuron.set_angles() must contain at least (2 ** len(self.controls)) elements."
             )
-        Qrack.qrack_lib.set_qneuron_angles(self.nid, QrackNeuron._real1_byref(a))
-        self._throw_if_error()
+        self.angles = QrackNeuron._real1_byref(a)
 
     def get_angles(self):
         """Directly gets the neuron parameters.
@@ -139,10 +155,7 @@ class QrackNeuron:
         Raises:
             RuntimeError: QrackNeuron C++ library raised an exception.
         """
-        ket = QrackNeuron._real1_byref([0.0] * (1 << len(self.controls)))
-        Qrack.qrack_lib.get_qneuron_angles(self.nid, ket)
-        self._throw_if_error()
-        return list(ket)
+        return list(self.angles)
 
     def set_alpha(self, a):
         """Set the neuron 'alpha' parameter.
@@ -151,13 +164,8 @@ class QrackNeuron:
         parameter that is applied as a power to its angles, before
         learning and prediction. This makes the activation function
         sharper (or less sharp).
-
-        Raises:
-            RuntimeError: QrackNeuron C++ library raised an exception.
         """
         self.alpha = a
-        Qrack.qrack_lib.set_qneuron_alpha(self.nid, a)
-        self._throw_if_error()
 
     def set_activation_fn(self, f):
         """Sets the activation function of this QrackNeuron
@@ -165,13 +173,8 @@ class QrackNeuron:
         Nonlinear activation functions can be important to neural net
         applications, like DNN. The available activation functions are
         enumerated in `NeuronActivationFn`.
-
-        Raises:
-            RuntimeError: QrackNeuron C++ library raised an exception.
         """
         self.activation_fn = f
-        Qrack.qrack_lib.set_qneuron_activation_fn(self.nid, f)
-        self._throw_if_error()
 
     def predict(self, e=True, r=True):
         """Predict based on training
@@ -190,7 +193,7 @@ class QrackNeuron:
         Raises:
             RuntimeError: QrackNeuron C++ library raised an exception.
         """
-        result = Qrack.qrack_lib.qneuron_predict(self.nid, e, r)
+        result = Qrack.qrack_lib.qneuron_predict(self.nid, self.angles, e, r, self.activation_fn, self.alpha)
         self._throw_if_error()
         return result
 
@@ -206,7 +209,7 @@ class QrackNeuron:
         Raises:
             RuntimeError: QrackNeuron C++ library raised an exception.
         """
-        result = Qrack.qrack_lib.qneuron_unpredict(self.nid, e)
+        result = Qrack.qrack_lib.qneuron_unpredict(self.nid, self.angles, e, self.activation_fn, self.alpha)
         self._throw_if_error()
         return result
 
@@ -222,7 +225,7 @@ class QrackNeuron:
         Raises:
             RuntimeError: QrackNeuron C++ library raised an exception.
         """
-        Qrack.qrack_lib.qneuron_learn_cycle(self.nid, e)
+        Qrack.qrack_lib.qneuron_learn_cycle(self.nid, self.angles, e, self.activation_fn, self.alpha)
         self._throw_if_error()
 
     def learn(self, eta, e=True, r=True):
@@ -241,7 +244,7 @@ class QrackNeuron:
         Raises:
             RuntimeError: QrackNeuron C++ library raised an exception.
         """
-        Qrack.qrack_lib.qneuron_learn(self.nid, eta, e, r)
+        Qrack.qrack_lib.qneuron_learn(self.nid, self.angles, eta, e, r, self.activation_fn, self.alpha)
         self._throw_if_error()
 
     def learn_permutation(self, eta, e=True, r=True):
@@ -260,12 +263,12 @@ class QrackNeuron:
         Raises:
             RuntimeError: QrackNeuron C++ library raised an exception.
         """
-        Qrack.qrack_lib.qneuron_learn_permutation(self.nid, eta, e, r)
+        Qrack.qrack_lib.qneuron_learn_permutation(self.nid, self.angles, eta, e, r, self.activation_fn, self.alpha)
         self._throw_if_error()
 
     @staticmethod
     def quantile_bounds(vec, bits):
-        """ Calculate vector quantile bounds
+        """Calculate vector quantile bounds
 
         This is a static helper method to calculate the quantile
         bounds of 2 ** bits worth of quantiles.
@@ -283,11 +286,15 @@ class QrackNeuron:
         n = len(vec)
         vec_sorted = sorted(vec)
 
-        return [vec_sorted[0]] + [vec_sorted[(k * n) // bins] for k in range(1, bins)] + [vec_sorted[-1]]
+        return (
+            [vec_sorted[0]]
+            + [vec_sorted[(k * n) // bins] for k in range(1, bins)]
+            + [vec_sorted[-1]]
+        )
 
     @staticmethod
     def discretize(vec, bounds):
-        """ Discretize vector by quantile bounds
+        """Discretize vector by quantile bounds
 
         This is a static helper method to discretize a numerical
         vector according to quantile bounds calculated by the
@@ -317,7 +324,7 @@ class QrackNeuron:
 
     @staticmethod
     def flatten_and_transpose(arr):
-        """ Flatten and transpose feature matrix
+        """Flatten and transpose feature matrix
 
         This is a static helper method to convert a multi-feature
         bit-row matrix to an observation-row matrix with flat
@@ -333,7 +340,7 @@ class QrackNeuron:
 
     @staticmethod
     def bin_endpoints_average(bounds):
-        """ Bin endpoints average
+        """Bin endpoints average
 
         This is a static helper method that accepts the output
         bins from quantile_bounds() and returns the average points