pyqrack-cpu 1.76.0__py3-none-macosx_14_0_arm64.whl → 1.80.3__py3-none-macosx_14_0_arm64.whl

pyqrack/__init__.py

@@ -10,8 +10,8 @@ from .qrack_ace_backend import QrackAceBackend
 from .qrack_circuit import QrackCircuit
 from .qrack_neuron import QrackNeuron
 from .qrack_neuron_torch_layer import (
-    QrackTorchNeuron,
-    QrackNeuronFunction,
+    QrackNeuronTorch,
+    QrackNeuronTorchFunction,
     QrackNeuronTorchLayer,
 )
 from .qrack_simulator import QrackSimulator

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):
@@ -432,7 +426,8 @@ class QrackAceBackend:
     def _unpack(self, lq):
         return self._qubits[lq]
 
-    def _get_qb_lhv_indices(self, hq):
+    @staticmethod
+    def _get_qb_lhv_indices(hq):
         qb = []
         if len(hq) < 2:
             qb = [0]
@@ -446,7 +441,8 @@ class QrackAceBackend:
 
         return qb, lhv
 
-    def _get_lhv_bloch_angles(self, sim):
+    @staticmethod
+    def _get_lhv_bloch_angles(sim):
         # Z axis
         z = sim.bloch[2]
 
@@ -503,7 +499,8 @@ class QrackAceBackend:
 
         sim.mtrx([m00, m01, m10, m11], q)
 
-    def _rotate_lhv_to_bloch(self, sim, delta_azimuth, delta_inclination):
+    @staticmethod
+    def _rotate_lhv_to_bloch(sim, delta_azimuth, delta_inclination):
         # Apply rotation as "Azimuth, Inclination" (AI)
         cosA = math.cos(delta_azimuth)
         sinA = math.sin(delta_azimuth)
@@ -523,7 +520,7 @@ class QrackAceBackend:
         if len(hq) == 1:
             return
 
-        qb, lhv = self._get_qb_lhv_indices(hq)
+        qb, lhv = QrackAceBackend._get_qb_lhv_indices(hq)
 
         if phase:
             for q in qb:
@@ -542,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
@@ -572,7 +567,7 @@ class QrackAceBackend:
                 a, i = [0, 0, 0, 0, 0], [0, 0, 0, 0, 0]
                 a[0], i[0] = self._get_bloch_angles(hq[0])
                 a[1], i[1] = self._get_bloch_angles(hq[1])
-                a[2], i[2] = self._get_lhv_bloch_angles(hq[2])
+                a[2], i[2] = QrackAceBackend._get_lhv_bloch_angles(hq[2])
                 a[3], i[3] = self._get_bloch_angles(hq[3])
                 a[4], i[4] = self._get_bloch_angles(hq[4])
 
@@ -588,7 +583,7 @@ class QrackAceBackend:
                 i_target /= 5
                 for x in range(5):
                     if x == 2:
-                        self._rotate_lhv_to_bloch(
+                        QrackAceBackend._rotate_lhv_to_bloch(
                             hq[x], a_target - a[x], i_target - i[x]
                         )
                     else:
@@ -617,7 +612,7 @@ class QrackAceBackend:
                 a, i = [0, 0, 0], [0, 0, 0]
                 a[0], i[0] = self._get_bloch_angles(hq[0])
                 a[1], i[1] = self._get_bloch_angles(hq[1])
-                a[2], i[2] = self._get_lhv_bloch_angles(hq[2])
+                a[2], i[2] = QrackAceBackend._get_lhv_bloch_angles(hq[2])
 
                 a_target = 0
                 i_target = 0
@@ -631,7 +626,7 @@ class QrackAceBackend:
                 i_target /= 3
                 for x in range(3):
                     if x == 2:
-                        self._rotate_lhv_to_bloch(
+                        QrackAceBackend._rotate_lhv_to_bloch(
                             hq[x], a_target - a[x], i_target - i[x]
                         )
                     else:
@@ -653,8 +648,8 @@ class QrackAceBackend:
             for c in range(len(hq)):
                 h = hq[c]
                 if c == 2:
-                    a, i = self._get_lhv_bloch_angles(h)
-                    self._rotate_lhv_to_bloch(
+                    a, i = QrackAceBackend._get_lhv_bloch_angles(h)
+                    QrackAceBackend._rotate_lhv_to_bloch(
                         h,
                         math.atan(math.tan(a) * b) - a,
                         math.atan(math.tan(i) * b) - i,
@@ -674,7 +669,7 @@ class QrackAceBackend:
             self.sim[b[0]].u(b[1], th, ph, lm)
             return
 
-        qb, lhv = self._get_qb_lhv_indices(hq)
+        qb, lhv = QrackAceBackend._get_qb_lhv_indices(hq)
 
         for q in qb:
             b = hq[q]
@@ -693,7 +688,7 @@ class QrackAceBackend:
             self.sim[b[0]].r(p, th, b[1])
             return
 
-        qb, lhv = self._get_qb_lhv_indices(hq)
+        qb, lhv = QrackAceBackend._get_qb_lhv_indices(hq)
 
         for q in qb:
             b = hq[q]
@@ -721,7 +716,7 @@ class QrackAceBackend:
 
         self._correct(lq)
 
-        qb, lhv = self._get_qb_lhv_indices(hq)
+        qb, lhv = QrackAceBackend._get_qb_lhv_indices(hq)
 
         for q in qb:
             b = hq[q]
@@ -739,7 +734,7 @@ class QrackAceBackend:
             self.sim[b[0]].s(b[1])
             return
 
-        qb, lhv = self._get_qb_lhv_indices(hq)
+        qb, lhv = QrackAceBackend._get_qb_lhv_indices(hq)
 
         for q in qb:
             b = hq[q]
@@ -755,7 +750,7 @@ class QrackAceBackend:
             self.sim[b[0]].adjs(b[1])
             return
 
-        qb, lhv = self._get_qb_lhv_indices(hq)
+        qb, lhv = QrackAceBackend._get_qb_lhv_indices(hq)
 
         for q in qb:
             b = hq[q]
@@ -771,7 +766,7 @@ class QrackAceBackend:
             self.sim[b[0]].x(b[1])
             return
 
-        qb, lhv = self._get_qb_lhv_indices(hq)
+        qb, lhv = QrackAceBackend._get_qb_lhv_indices(hq)
 
         for q in qb:
             b = hq[q]
@@ -787,7 +782,7 @@ class QrackAceBackend:
             self.sim[b[0]].y(b[1])
             return
 
-        qb, lhv = self._get_qb_lhv_indices(hq)
+        qb, lhv = QrackAceBackend._get_qb_lhv_indices(hq)
 
         for q in qb:
             b = hq[q]
@@ -803,7 +798,7 @@ class QrackAceBackend:
             self.sim[b[0]].z(b[1])
             return
 
-        qb, lhv = self._get_qb_lhv_indices(hq)
+        qb, lhv = QrackAceBackend._get_qb_lhv_indices(hq)
 
         for q in qb:
             b = hq[q]
@@ -819,7 +814,7 @@ class QrackAceBackend:
             self.sim[b[0]].t(b[1])
             return
 
-        qb, lhv = self._get_qb_lhv_indices(hq)
+        qb, lhv = QrackAceBackend._get_qb_lhv_indices(hq)
 
         for q in qb:
             b = hq[q]
@@ -835,7 +830,7 @@ class QrackAceBackend:
             self.sim[b[0]].adjt(b[1])
             return
 
-        qb, lhv = self._get_qb_lhv_indices(hq)
+        qb, lhv = QrackAceBackend._get_qb_lhv_indices(hq)
 
         for q in qb:
             b = hq[q]
@@ -857,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
 
@@ -896,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):
@@ -917,8 +906,8 @@ class QrackAceBackend:
 
         self._correct(lq1)
 
-        qb1, lhv1 = self._get_qb_lhv_indices(hq1)
-        qb2, lhv2 = self._get_qb_lhv_indices(hq2)
+        qb1, lhv1 = QrackAceBackend._get_qb_lhv_indices(hq1)
+        qb2, lhv2 = QrackAceBackend._get_qb_lhv_indices(hq2)
         # Apply cross coupling on hardware qubits first
         self._apply_coupling(pauli, anti, qb1, lhv1, hq1, qb2, lhv2, hq2, lq1_lr)
         # Apply coupling to the local-hidden-variable target
@@ -1030,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
@@ -1065,7 +1052,7 @@ class QrackAceBackend:
         p = self.prob(lq)
         result = ((p + self._epsilon) >= 1) or (random.random() < p)
 
-        qb, lhv = self._get_qb_lhv_indices(hq)
+        qb, lhv = QrackAceBackend._get_qb_lhv_indices(hq)
 
         for q in qb:
             b = hq[q]
@@ -1091,7 +1078,7 @@ class QrackAceBackend:
 
         self._correct(lq)
 
-        qb, lhv = self._get_qb_lhv_indices(hq)
+        qb, lhv = QrackAceBackend._get_qb_lhv_indices(hq)
 
         for q in qb:
             b = hq[q]
@@ -1184,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":
@@ -1259,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)
@@ -1376,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
@@ -1436,14 +1415,13 @@ 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
 
         return _data
 
+    @staticmethod
     def get_qiskit_basis_gates():
         return [
             "id",
@@ -1528,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

@@ -59,7 +59,7 @@ class QrackCircuit:
             self.cid = Qrack.qrack_lib.qcircuit_inverse(clone_cid)
         elif len(past_light_cone) > 0:
             self.cid = Qrack.qrack_lib.qcircuit_past_light_cone(
-                clone_cid, len(past_light_cone), self._ulonglong_byref(past_light_cone)
+                clone_cid, len(past_light_cone), QrackCircuit._ulonglong_byref(past_light_cone)
             )
         else:
             self.cid = Qrack.qrack_lib.init_qcircuit_clone(clone_cid)
@@ -69,16 +69,20 @@ class QrackCircuit:
             Qrack.qrack_lib.destroy_qcircuit(self.cid)
             self.cid = None
 
-    def _ulonglong_byref(self, a):
+    @staticmethod
+    def _ulonglong_byref(a):
         return (ctypes.c_ulonglong * len(a))(*a)
 
-    def _double_byref(self, a):
+    @staticmethod
+    def _double_byref(a):
         return (ctypes.c_double * len(a))(*a)
 
-    def _complex_byref(self, a):
+    @staticmethod
+    def _complex_byref(a):
         t = [(c.real, c.imag) for c in a]
-        return self._double_byref([float(item) for sublist in t for item in sublist])
+        return QrackCircuit._double_byref([float(item) for sublist in t for item in sublist])
 
+    @staticmethod
     def _mtrx_to_u4(m):
         nrm = abs(m[0])
         if (nrm * nrm) < sys.float_info.epsilon:
@@ -103,6 +107,7 @@ class QrackCircuit:
 
         return th, ph, lm, np.angle(phase)
 
+    @staticmethod
     def _u3_to_mtrx(params):
         th = float(params[0])
         ph = float(params[1])
@@ -115,6 +120,7 @@ class QrackCircuit:
 
         return [c + 0j, -el * s, ep * s, ep * el * c]
 
+    @staticmethod
     def _u4_to_mtrx(params):
         m = QrackCircuit._u3_to_mtrx(params)
         g = np.exp(1j * float(params[3]))
@@ -123,6 +129,7 @@ class QrackCircuit:
 
         return m
 
+    @staticmethod
     def _make_mtrx_unitary(m):
         return QrackCircuit._u4_to_mtrx(QrackCircuit._mtrx_to_u4(m))
 
@@ -190,7 +197,7 @@ class QrackCircuit:
             raise ValueError(
                 "2x2 matrix 'm' in QrackCircuit.mtrx() must contain at least 4 elements."
             )
-        Qrack.qrack_lib.qcircuit_append_1qb(self.cid, self._complex_byref(m), q)
+        Qrack.qrack_lib.qcircuit_append_1qb(self.cid, QrackCircuit._complex_byref(m), q)
 
     def ucmtrx(self, c, m, q, p):
         """Multi-controlled single-target-qubit gate
@@ -214,7 +221,7 @@ class QrackCircuit:
                 "2x2 matrix 'm' in QrackCircuit.ucmtrx() must contain at least 4 elements."
             )
         Qrack.qrack_lib.qcircuit_append_mc(
-            self.cid, self._complex_byref(m), len(c), self._ulonglong_byref(c), q, p
+            self.cid, QrackCircuit._complex_byref(m), len(c), QrackCircuit._ulonglong_byref(c), q, p
         )
 
     def run(self, qsim):
@@ -248,20 +255,6 @@ class QrackCircuit:
         """
         Qrack.qrack_lib.qcircuit_out_to_file(self.cid, filename.encode("utf-8"))
 
-    def in_from_file(filename):
-        """Read in optimized circuit from file
-
-        Reads in an (optimized) circuit from a file named
-        according to the "filename" parameter.
-
-        Args:
-            filename: Name of file
-        """
-        out = QrackCircuit()
-        Qrack.qrack_lib.qcircuit_in_from_file(out.cid, filename.encode("utf-8"))
-
-        return out
-
     def out_to_string(self):
         """Output optimized circuit to string
 
@@ -273,19 +266,6 @@ class QrackCircuit:
 
         return out.value.decode("utf-8")
 
-    def file_gate_count(filename):
-        """File gate count
-
-        Return the count of gates in a QrackCircuit file
-
-        Args:
-            filename: Name of file
-        """
-        tokens = []
-        with open(filename, "r") as file:
-            tokens = file.read().split()
-        return int(tokens[1])
-
     def to_qiskit_circuit(self):
         """Convert to a Qiskit circuit
 
@@ -302,6 +282,36 @@ class QrackCircuit:
 
         return QrackCircuit.string_to_qiskit_circuit(self.out_to_string())
 
+    @staticmethod
+    def in_from_file(filename):
+        """Read in optimized circuit from file
+
+        Reads in an (optimized) circuit from a file named
+        according to the "filename" parameter.
+
+        Args:
+            filename: Name of file
+        """
+        out = QrackCircuit()
+        Qrack.qrack_lib.qcircuit_in_from_file(out.cid, filename.encode("utf-8"))
+
+        return out
+
+    @staticmethod
+    def file_gate_count(filename):
+        """File gate count
+
+        Return the count of gates in a QrackCircuit file
+
+        Args:
+            filename: Name of file
+        """
+        tokens = []
+        with open(filename, "r") as file:
+            tokens = file.read().split()
+        return int(tokens[1])
+
+    @staticmethod
     def file_to_qiskit_circuit(filename):
         """Convert an output file to a Qiskit circuit
 
@@ -325,6 +335,7 @@ class QrackCircuit:
         with open(filename, "r") as file:
             return QrackCircuit.string_to_qiskit_circuit(file.read())
 
+    @staticmethod
     def string_to_qiskit_circuit(circ_string):
         """Convert an output string to a Qiskit circuit
 
@@ -398,6 +409,7 @@ class QrackCircuit:
 
         return circ
 
+    @staticmethod
     def in_from_qiskit_circuit(circ):
         """Read a Qiskit circuit into a QrackCircuit
 
@@ -443,6 +455,7 @@ class QrackCircuit:
 
         return out
 
+    @staticmethod
     def file_to_quimb_circuit(
         filename,
         circuit_type=QuimbCircuitType.Circuit,
@@ -494,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,
@@ -526,9 +537,8 @@ class QrackCircuit:
 
         return tcirc
 
-    def file_to_tensorcircuit(
-        filename, inputs=None, circuit_params=None, binding_params=None
-    ):
+    @staticmethod
+    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
@@ -558,6 +568,7 @@ class QrackCircuit:
             qcirc, qcirc.num_qubits, inputs, circuit_params, binding_params
         )
 
+    @staticmethod
     def in_from_tensorcircuit(tcirc, enable_instruction=False, enable_inputs=False):
         """Convert a TensorCircuit circuit to a QrackCircuit
 

pyqrack/qrack_neuron.py

@@ -64,7 +64,7 @@ class QrackNeuron:
         self.nid = Qrack.qrack_lib.init_qneuron(
             simulator.sid,
             len(controls),
-            self._ulonglong_byref(controls),
+            QrackNeuron._ulonglong_byref(controls),
             target,
             activation_fn,
             alpha,
@@ -99,15 +99,32 @@ class QrackNeuron:
         self._throw_if_error()
         return result
 
-    def _ulonglong_byref(self, a):
+    @staticmethod
+    def _ulonglong_byref(a):
         return (ctypes.c_ulonglong * len(a))(*a)
 
-    def _real1_byref(self, a):
+    @staticmethod
+    def _real1_byref(a):
         # This needs to be c_double, if PyQrack is built with fp64.
         if Qrack.fppow < 6:
             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)
+        self._throw_if_error()
+        self.simulator = s
+
     def set_angles(self, a):
         """Directly sets the neuron parameters.
 
@@ -125,7 +142,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, self._real1_byref(a))
+        Qrack.qrack_lib.set_qneuron_angles(self.nid, QrackNeuron._real1_byref(a))
         self._throw_if_error()
 
     def get_angles(self):
@@ -137,7 +154,7 @@ class QrackNeuron:
         Raises:
             RuntimeError: QrackNeuron C++ library raised an exception.
         """
-        ket = self._real1_byref([0.0] * (1 << len(self.controls)))
+        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)
@@ -261,8 +278,9 @@ class QrackNeuron:
         Qrack.qrack_lib.qneuron_learn_permutation(self.nid, eta, e, r)
         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.
@@ -280,10 +298,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
@@ -311,8 +334,9 @@ class QrackNeuron:
 
         return vec_discrete
 
+    @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
@@ -326,8 +350,9 @@ class QrackNeuron:
         """
         return list(zip(*[item for sublist in arr for item in sublist]))
 
+    @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