pyqrack-complex128 1.61.2__py3-none-manylinux_2_35_x86_64.whl → 1.82.1__py3-none-manylinux_2_35_x86_64.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

@@ -139,6 +139,12 @@ class LHVQubit:
         new_y = 2 * (rho_prime[0][1].imag - rho_prime[1][0].imag)
         new_z = 2 * rho_prime[0][0].real - 1  # since Tr(ρ') = 1
 
+        p = math.sqrt(new_x**2 + new_y**2 + new_z**2)
+
+        new_x /= p
+        new_y /= p
+        new_z /= p
+
         self.bloch = [new_x, new_y, new_z]
 
     def prob(self, basis=Pauli.PauliZ):
@@ -218,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,
     ):
@@ -273,7 +277,7 @@ class QrackAceBackend:
                 self._is_row_long_range[-1] = False
         sim_count = col_patch_count * row_patch_count
 
-        self._qubit_dict = {}
+        self._qubits = []
         sim_counts = [0] * sim_count
         sim_id = 0
         tot_qubits = 0
@@ -288,11 +292,7 @@ class QrackAceBackend:
                     sim_counts[t_sim_id] += 1
 
                     qubit.append(
-                        LHVQubit(
-                            toClone=(
-                                toClone._qubit_dict[tot_qubits][2] if toClone else None
-                            )
-                        )
+                        LHVQubit(toClone=(toClone._qubits[tot_qubits][2] if toClone else None))
                     )
 
                 if (not c) and (not r):
@@ -307,7 +307,7 @@ class QrackAceBackend:
                 if not c:
                     sim_id = (sim_id + 1) % sim_count
 
-                self._qubit_dict[tot_qubits] = qubit
+                self._qubits.append(qubit)
                 tot_qubits += 1
 
         self.sim = []
@@ -424,9 +424,10 @@ class QrackAceBackend:
         self._qec_x(c)
 
     def _unpack(self, lq):
-        return self._qubit_dict[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]
@@ -440,15 +441,16 @@ class QrackAceBackend:
 
         return qb, lhv
 
-    def _get_lhv_bloch_angles(self, sim):
+    @staticmethod
+    def _get_lhv_bloch_angles(sim):
         # Z axis
-        z = 1 - 2 * sim.prob(Pauli.PauliZ)
+        z = sim.bloch[2]
 
         # X axis
-        x = 1 - 2 * sim.prob(Pauli.PauliX)
+        x = sim.bloch[0]
 
         # Y axis
-        y = 1 - 2 * sim.prob(Pauli.PauliY)
+        y = sim.bloch[1]
 
         inclination = math.atan2(math.sqrt(x**2 + y**2), z)
         azimuth = math.atan2(y, x)
@@ -480,9 +482,7 @@ class QrackAceBackend:
 
         return azimuth, inclination
 
-    def _rotate_to_bloch(
-        self, hq, delta_azimuth, delta_inclination
-    ):
+    def _rotate_to_bloch(self, hq, delta_azimuth, delta_inclination):
         sim = self.sim[hq[0]]
         q = hq[1]
 
@@ -499,10 +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)
@@ -516,14 +514,13 @@ class QrackAceBackend:
 
         sim.mtrx([m00, m01, m10, m11])
 
-
     def _correct(self, lq, phase=False, skip_rotation=False):
         hq = self._unpack(lq)
 
         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,13 +567,15 @@ 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])
 
                 a_target = 0
                 i_target = 0
                 for x in range(5):
+                    if x == 2:
+                        continue
                     a_target += a[x]
                     i_target += i[x]
 
@@ -586,7 +583,9 @@ class QrackAceBackend:
                 i_target /= 5
                 for x in range(5):
                     if x == 2:
-                        self._rotate_lhv_to_bloch(hq[x], a_target - a[x], i_target - i[x])
+                        QrackAceBackend._rotate_lhv_to_bloch(
+                            hq[x], a_target - a[x], i_target - i[x]
+                        )
                     else:
                         self._rotate_to_bloch(hq[x], a_target - a[x], i_target - i[x])
 
@@ -613,11 +612,13 @@ 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
                 for x in range(3):
+                    if x == 2:
+                        continue
                     a_target += a[x]
                     i_target += i[x]
 
@@ -625,7 +626,9 @@ class QrackAceBackend:
                 i_target /= 3
                 for x in range(3):
                     if x == 2:
-                        self._rotate_lhv_to_bloch(hq[x], a_target - a[x], i_target - i[x])
+                        QrackAceBackend._rotate_lhv_to_bloch(
+                            hq[x], a_target - a[x], i_target - i[x]
+                        )
                     else:
                         self._rotate_to_bloch(hq[x], a_target - a[x], i_target - i[x])
 
@@ -636,6 +639,29 @@ class QrackAceBackend:
             b = hq[lhv]
             b.h()
 
+    def apply_magnetic_bias(self, q, b):
+        if b == 0:
+            return
+        b = math.exp(b)
+        for x in q:
+            hq = self._unpack(x)
+            for c in range(len(hq)):
+                h = hq[c]
+                if c == 2:
+                    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,
+                    )
+                else:
+                    a, i = self._get_bloch_angles(h)
+                    self._rotate_to_bloch(
+                        h,
+                        math.atan(math.tan(a) * b) - a,
+                        math.atan(math.tan(i) * b) - i,
+                    )
+
     def u(self, lq, th, ph, lm):
         hq = self._unpack(lq)
         if len(hq) < 2:
@@ -643,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]
@@ -662,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]
@@ -690,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]
@@ -708,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]
@@ -724,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]
@@ -740,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]
@@ -756,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]
@@ -772,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]
@@ -788,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]
@@ -804,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]
@@ -826,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
 
@@ -865,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):
@@ -886,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
@@ -999,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
@@ -1034,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]
@@ -1060,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]
@@ -1153,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":
@@ -1228,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)
@@ -1345,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
@@ -1405,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",
@@ -1497,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