pyqrack-cuda 1.27.0__tar.gz

pyqrack-cuda-1.27.0/LICENSE

@@ -0,0 +1,21 @@
+MIT License
+
+Copyright (c) 2021 vm6502q
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in all
+copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+SOFTWARE.

pyqrack-cuda-1.27.0/MANIFEST.in

@@ -0,0 +1,11 @@
+include MANIFEST.in
+include README.md
+recursive-include pyqrack *.so
+recursive-include pyqrack *.so.*
+recursive-include pyqrack *.a
+recursive-include pyqrack *.dylib
+recursive-include pyqrack *.lib
+recursive-include pyqrack *.dll
+recursive-include pyqrack *.exp
+recursive-include pyqrack qrack_cl_precompile
+recursive-include pyqrack qrack_cl_precompile.exe

pyqrack-cuda-1.27.0/PKG-INFO

@@ -0,0 +1,62 @@
+Metadata-Version: 2.1
+Name: pyqrack-cuda
+Version: 1.27.0
+Summary: pyqrack - Pure Python vm6502q/qrack Wrapper
+Home-page: https://github.com/vm6502q/pyqrack
+Author: Daniel Strano
+Author-email: dan@unitary.fund
+License: MIT
+Keywords: pyqrack qrack simulator quantum gpu
+Classifier: Environment :: Console
+Classifier: Intended Audience :: Developers
+Classifier: Intended Audience :: Science/Research
+Classifier: Operating System :: Microsoft :: Windows
+Classifier: Operating System :: MacOS
+Classifier: Operating System :: POSIX :: Linux
+Classifier: Programming Language :: C++
+Classifier: Programming Language :: Python :: 3.5
+Classifier: Programming Language :: Python :: 3.6
+Classifier: Programming Language :: Python :: 3.7
+Classifier: Programming Language :: Python :: 3.8
+Classifier: Programming Language :: Python :: 3.9
+Classifier: Programming Language :: Python :: 3.10
+Classifier: Topic :: Scientific/Engineering
+Description-Content-Type: text/markdown
+License-File: LICENSE
+Provides-Extra: dev
+Requires-Dist: pytest>=7.3.1; extra == "dev"
+
+# pyqrack
+[![Downloads](https://pepy.tech/badge/pyqrack)](https://pepy.tech/project/pyqrack) [![Downloads](https://pepy.tech/badge/pyqrack/month)](https://pepy.tech/project/pyqrack)
+
+Pure Python bindings for the pure C++11/OpenCL Qrack quantum computer simulator library
+
+(**PyQrack** is just pure Qrack.)
+
+To use this package, it's helpful to be familiar with [vm6502q/qrack](https://github.com/vm6502q/qrack). Users gain **much** more control over options by building **vm6502q/qrack** and [vm6502q/pyqrack](https://github.com/vm6502q/pyqrack) from source. For advanced users, building from source is the intended primary method of **PyQrack** package distribution and use.
+
+Import and instantiate [`QrackSimulator`](https://github.com/vm6502q/pyqrack/blob/main/pyqrack/qrack_simulator.py) instances. This simulator can perform arbitrary single qubit and controlled-single-qubit gates, as well as other specific gates like `SWAP`.
+
+Any 2x2 bit operator matrix is represented by a list of 4 `complex` floating point numbers, in [**row-major order**](https://en.wikipedia.org/wiki/Row-_and_column-major_order).
+
+Single and array "`b`" parameters represent [**Pauli operator bases**](https://en.wikipedia.org/wiki/Pauli_matrices). They are specified according to the enumeration of the [`Pauli`](https://github.com/vm6502q/pyqrack/blob/main/pyqrack/pauli.py) class.
+
+`MC[x]` and `MAC[x]` methods are controlled single bit gates, with as many control qubits as you specify via Python list `c` argument. `MCX` is multiply-controlled Pauli X, and `MACX` is "anti-"controlled Pauli X, i.e. "anti-control" activates the gate if all control bits are specifically **off**, as opposed to **on**.
+
+The package installation directory contains a `qrack_cl_precompile` folder with executables for supported platforms, to compile OpenCL kernels once, beforehand, avoiding the need to recompile "just-in-time" every time that you load this package in a Python interpreter. If you no longer want to use precompiled kernels, or if precompilation fails, just delete the `~/.qrack` directory, or the equivalent `.qrack` sub-directory in the user home folder of your operating system.
+
+To load the required **vm6502q/qrack** libraries from a different location, set the `PYQRACK_SHARED_LIB_PATH` environment variable. If an absolute path is set for this variable, it will be treated as an absolute path. If a relative path is used, the **pyqrack** package `qrack_system` directory is the root.
+
+PyQrack v0.4.6 adds experimental support for [PyZX](https://github.com/Quantomatic/pyzx) `Circuit` definitions as an intermediate representation for `QrackSimulator`. To try this, load a `Circuit` in PyZX, (use that module to optimize your circuit, as you like,) and create a `QrackSimulator()` instance using the `pyzxCircuit` named argument of the constructor, like so:
+
+```python
+sim = QrackSimulator(pyzxCircuit=c)
+```
+
+where `c` is a PyZX circuit object. The circuit will automatically be simulated in the constructed `QrackSimulator` instance. This also allows loading from QASM and other intermediate representations supported by PyZX.
+
+See [https://pyqrack.readthedocs.io/en/latest/](https://pyqrack.readthedocs.io/en/latest/) for an API reference.
+
+Please feel welcome to open an issue, if you'd like help. 😃
+
+**Special thanks go to Zeeshan Ahmed, for bug fixes and design suggestions, Ashish Panigrahi, for documentation and design suggestions, WingCode, for documentation, and to the broader community of Qrack contributors, for years of happy Qracking! You rock!**

pyqrack-cuda-1.27.0/README.md

@@ -0,0 +1,34 @@
+# pyqrack
+[![Downloads](https://pepy.tech/badge/pyqrack)](https://pepy.tech/project/pyqrack) [![Downloads](https://pepy.tech/badge/pyqrack/month)](https://pepy.tech/project/pyqrack)
+
+Pure Python bindings for the pure C++11/OpenCL Qrack quantum computer simulator library
+
+(**PyQrack** is just pure Qrack.)
+
+To use this package, it's helpful to be familiar with [vm6502q/qrack](https://github.com/vm6502q/qrack). Users gain **much** more control over options by building **vm6502q/qrack** and [vm6502q/pyqrack](https://github.com/vm6502q/pyqrack) from source. For advanced users, building from source is the intended primary method of **PyQrack** package distribution and use.
+
+Import and instantiate [`QrackSimulator`](https://github.com/vm6502q/pyqrack/blob/main/pyqrack/qrack_simulator.py) instances. This simulator can perform arbitrary single qubit and controlled-single-qubit gates, as well as other specific gates like `SWAP`.
+
+Any 2x2 bit operator matrix is represented by a list of 4 `complex` floating point numbers, in [**row-major order**](https://en.wikipedia.org/wiki/Row-_and_column-major_order).
+
+Single and array "`b`" parameters represent [**Pauli operator bases**](https://en.wikipedia.org/wiki/Pauli_matrices). They are specified according to the enumeration of the [`Pauli`](https://github.com/vm6502q/pyqrack/blob/main/pyqrack/pauli.py) class.
+
+`MC[x]` and `MAC[x]` methods are controlled single bit gates, with as many control qubits as you specify via Python list `c` argument. `MCX` is multiply-controlled Pauli X, and `MACX` is "anti-"controlled Pauli X, i.e. "anti-control" activates the gate if all control bits are specifically **off**, as opposed to **on**.
+
+The package installation directory contains a `qrack_cl_precompile` folder with executables for supported platforms, to compile OpenCL kernels once, beforehand, avoiding the need to recompile "just-in-time" every time that you load this package in a Python interpreter. If you no longer want to use precompiled kernels, or if precompilation fails, just delete the `~/.qrack` directory, or the equivalent `.qrack` sub-directory in the user home folder of your operating system.
+
+To load the required **vm6502q/qrack** libraries from a different location, set the `PYQRACK_SHARED_LIB_PATH` environment variable. If an absolute path is set for this variable, it will be treated as an absolute path. If a relative path is used, the **pyqrack** package `qrack_system` directory is the root.
+
+PyQrack v0.4.6 adds experimental support for [PyZX](https://github.com/Quantomatic/pyzx) `Circuit` definitions as an intermediate representation for `QrackSimulator`. To try this, load a `Circuit` in PyZX, (use that module to optimize your circuit, as you like,) and create a `QrackSimulator()` instance using the `pyzxCircuit` named argument of the constructor, like so:
+
+```python
+sim = QrackSimulator(pyzxCircuit=c)
+```
+
+where `c` is a PyZX circuit object. The circuit will automatically be simulated in the constructed `QrackSimulator` instance. This also allows loading from QASM and other intermediate representations supported by PyZX.
+
+See [https://pyqrack.readthedocs.io/en/latest/](https://pyqrack.readthedocs.io/en/latest/) for an API reference.
+
+Please feel welcome to open an issue, if you'd like help. 😃
+
+**Special thanks go to Zeeshan Ahmed, for bug fixes and design suggestions, Ashish Panigrahi, for documentation and design suggestions, WingCode, for documentation, and to the broader community of Qrack contributors, for years of happy Qracking! You rock!**

pyqrack-cuda-1.27.0/pyproject.toml

@@ -0,0 +1,6 @@
+[build-system]
+requires = [
+    "setuptools>=42",
+    "wheel"
+]
+build-backend = "setuptools.build_meta"

pyqrack-cuda-1.27.0/pyqrack/__init__.py

@@ -0,0 +1,13 @@
+# (C) Daniel Strano and the Qrack contributors 2017-2021. All rights reserved.
+#
+# Use of this source code is governed by an MIT-style license that can be
+# found in the LICENSE file or at https://opensource.org/licenses/MIT.
+
+from .qrack_system import QrackSystem, Qrack
+from .qrack_simulator import QrackSimulator
+from .qrack_neuron import QrackNeuron
+from .qrack_circuit import QrackCircuit
+from .pauli import Pauli
+from .neuron_activation_fn import NeuronActivationFn
+from .quimb_circuit_type import QuimbCircuitType
+from .util import convert_qiskit_circuit_to_qasm_experiment

pyqrack-cuda-1.27.0/pyqrack/neuron_activation_fn.py

@@ -0,0 +1,21 @@
+# (C) Daniel Strano and the Qrack contributors 2017-2021. All rights reserved.
+#
+# Pauli operators are specified for "b" (or "basis") parameters.
+#
+# Use of this source code is governed by an MIT-style license that can be
+# found in the LICENSE file or at https://opensource.org/licenses/MIT.
+
+from enum import IntEnum
+
+
+class NeuronActivationFn(IntEnum):
+    # Default
+    Sigmoid = 0
+    # Rectified linear 
+    ReLU = 1
+    # Gaussian linear
+    GeLU = 2
+    # Version of (default) "Sigmoid" with tunable sharpness
+    Generalized_Logistic = 3
+    # Leaky rectified linear
+    LeakyReLU = 4

pyqrack-cuda-1.27.0/pyqrack/pauli.py

@@ -0,0 +1,19 @@
+# (C) Daniel Strano and the Qrack contributors 2017-2021. All rights reserved.
+#
+# Pauli operators are specified for "b" (or "basis") parameters.
+#
+# Use of this source code is governed by an MIT-style license that can be
+# found in the LICENSE file or at https://opensource.org/licenses/MIT.
+
+from enum import IntEnum
+
+
+class Pauli(IntEnum):
+    # Pauli Identity operator. Corresponds to Q# constant "PauliI."
+    PauliI = 0
+    # Pauli X operator. Corresponds to Q# constant "PauliX."
+    PauliX = 1
+    # Pauli Y operator. Corresponds to Q# constant "PauliY."
+    PauliY = 3
+    # Pauli Z operator. Corresponds to Q# constant "PauliZ."
+    PauliZ = 2

pyqrack-cuda-1.27.0/pyqrack/qrack_circuit.py

@@ -0,0 +1,472 @@
+# (C) Daniel Strano and the Qrack contributors 2017-2023. All rights reserved.
+#
+# Use of this source code is governed by an MIT-style license that can be
+# found in the LICENSE file or at https://opensource.org/licenses/MIT.
+
+import ctypes
+
+from .qrack_system import Qrack
+from .quimb_circuit_type import QuimbCircuitType
+
+_IS_QISKIT_AVAILABLE = True
+try:
+    from qiskit.circuit.quantumcircuit import QuantumCircuit
+    from qiskit.compiler.transpiler import transpile
+    import numpy as np
+    import math
+except ImportError:
+    _IS_QISKIT_AVAILABLE = False
+
+_IS_QUIMB_AVAILABLE = True
+try:
+    import quimb as qu
+    import quimb.tensor as qtn
+except ImportError:
+    _IS_QUIMB_AVAILABLE = False
+
+_IS_TENSORCIRCUIT_AVAILABLE = True
+try:
+    import tensorcircuit as tc
+except ImportError:
+    _IS_TENSORCIRCUIT_AVAILABLE = False
+
+class QrackCircuit:
+    """Class that exposes the QCircuit class of Qrack
+
+    QrackCircuit allows the user to specify a unitary circuit, before running it.
+    Upon running the state, the result is a QrackSimulator state. Currently,
+    measurement is not supported, but measurement can be run on the resultant
+    QrackSimulator.
+
+    Attributes:
+        cid(int): Qrack ID of this circuit
+    """
+
+    def __init__(self, is_collapse = True, is_near_clifford = False, clone_cid = -1, is_inverse=False, past_light_cone = []):
+        if clone_cid < 0:
+            self.cid = Qrack.qrack_lib.init_qcircuit(is_collapse, is_near_clifford)
+        elif is_inverse:
+            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))
+        else:
+            self.cid = Qrack.qrack_lib.init_qcircuit_clone(clone_cid)
+
+    def __del__(self):
+        if self.cid is not None:
+            Qrack.qrack_lib.destroy_qcircuit(self.cid)
+            self.cid = None
+
+    def _ulonglong_byref(self, a):
+        return (ctypes.c_ulonglong * len(a))(*a)
+
+    def _double_byref(self, a):
+        return (ctypes.c_double * len(a))(*a)
+
+    def _complex_byref(self, a):
+        t = [(c.real, c.imag) for c in a]
+        return self._double_byref([float(item) for sublist in t for item in sublist])
+
+    def clone(self):
+        """Make a new circuit that is an exact clone of this circuit
+
+        Raises:
+            RuntimeError: QrackCircuit C++ library raised an exception.
+        """
+        return QrackCircuit(clone_cid = self.cid, is_inverse = False)
+
+    def inverse(self):
+        """Make a new circuit that is the exact inverse of this circuit
+
+        Raises:
+            RuntimeError: QrackCircuit C++ library raised an exception.
+        """
+        return QrackCircuit(clone_cid = self.cid, is_inverse = True)
+
+    def past_light_cone(self, q):
+        """Make a new circuit with just this circuits' past light cone for certain qubits.
+
+        Args:
+            q: list of qubit indices to include at beginning of past light cone
+
+        Raises:
+            RuntimeError: QrackCircuit C++ library raised an exception.
+        """
+        return QrackCircuit(clone_cid = self.cid, is_inverse = False, past_light_cone = q)
+
+    def get_qubit_count(self):
+        """Get count of qubits in circuit
+
+        Raises:
+            RuntimeError: QrackCircuit C++ library raised an exception.
+        """
+        return Qrack.qrack_lib.get_qcircuit_qubit_count(self.cid)
+
+    def swap(self, q1, q2):
+        """Add a 'Swap' gate to the circuit
+
+        Args:
+            q1: qubit index #1
+            q2: qubit index #2
+
+        Raises:
+            RuntimeError: QrackCircuit C++ library raised an exception.
+        """
+        Qrack.qrack_lib.qcircuit_swap(self.cid, q1, q2)
+
+    def mtrx(self, m, q):
+        """Operation from matrix.
+
+        Applies arbitrary operation defined by the given matrix.
+
+        Args:
+            m: row-major complex list representing the operator.
+            q: the qubit number on which the gate is applied to.
+
+        Raises:
+            ValueError: 2x2 matrix 'm' in QrackCircuit.mtrx() must contain at least 4 elements.
+            RuntimeError: QrackSimulator raised an exception.
+        """
+        if len(m) < 4:
+            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)
+
+    def ucmtrx(self, c, m, q, p):
+        """Multi-controlled single-target-qubit gate
+
+        Specify a controlled gate by its control qubits, its single-qubit
+        matrix "payload," the target qubit, and the permutation of qubits
+        that activates the gate.
+
+        Args:
+            c: list of controlled qubits
+            m: row-major complex list representing the operator.
+            q: target qubit
+            p: permutation of target qubits
+
+        Raises:
+            ValueError: 2x2 matrix 'm' in QrackCircuit.ucmtrx() must contain at least 4 elements.
+            RuntimeError: QrackSimulator raised an exception.
+        """
+        if len(m) < 4:
+            raise ValueError("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
+        )
+
+    def run(self, qsim):
+        """Run circuit on simulator
+
+        Run the encoded circuit on a specific simulator. The
+        result will remain in this simulator.
+
+        Args:
+            qsim: QrackSimulator on which to run circuit
+
+        Raises:
+            RuntimeError: QrackCircuit raised an exception.
+        """
+        qb_count = self.get_qubit_count()
+        sim_qb_count = qsim.num_qubits()
+        if sim_qb_count < qb_count:
+            for i in range(sim_qb_count, qb_count):
+                qsim.allocate_qubit(i)
+        Qrack.qrack_lib.qcircuit_run(self.cid, qsim.sid)
+        qsim._throw_if_error()
+
+    def out_to_file(self, filename):
+        """Output optimized circuit to file
+
+        Outputs the (optimized) circuit to a file named
+        according to the "filename" parameter.
+
+        Args:
+            filename: Name of file
+        """
+        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 file_to_qiskit_circuit(filename):
+        """Convert an output file to a Qiskit circuit
+
+        Reads in an (optimized) circuit from a file named
+        according to the "filename" parameter and outputs
+        a Qiskit circuit.
+
+        Args:
+            filename: Name of file
+
+        Raises:
+            RuntimeErorr: Before trying to file_to_qiskit_circuit() with
+                QrackCircuit, you must install Qiskit, numpy, and math!
+        """
+        if not _IS_QISKIT_AVAILABLE:
+            raise RuntimeError(
+                "Before trying to file_to_qiskit_circuit() with QrackCircuit, you must install Qiskit, numpy, and math!"
+            )
+
+        tokens = []
+        with open(filename, 'r') as file:
+            tokens = file.read().split()
+
+        i = 0
+        num_qubits = int(tokens[i])
+        i = i + 1
+        circ = QuantumCircuit(num_qubits)
+
+        num_gates = int(tokens[i])
+        i = i + 1
+
+        for g in range(num_gates):
+            target = int(tokens[i])
+            i = i + 1
+
+            control_count = int(tokens[i])
+            i = i + 1
+            controls = []
+            for j in range(control_count):
+                controls.append(int(tokens[i]))
+                i = i + 1
+
+            payload_count = int(tokens[i])
+            i = i + 1
+            payloads = {}
+            for j in range(payload_count):
+                key = int(tokens[i])
+                i = i + 1
+                op = np.zeros((2,2), dtype=complex)
+                row = []
+                for _ in range(2):
+                    amp = tokens[i].replace("(","").replace(")","").split(',')
+                    row.append(float(amp[0]) + float(amp[1])*1j)
+                    i = i + 1
+                l = math.sqrt(np.real(row[0] * np.conj(row[0]) + row[1] * np.conj(row[1])))
+                op[0][0] = row[0] / l
+                op[0][1] = row[1] / l
+
+                if np.abs(op[0][0] - row[0]) > 1e-5:
+                    print("Warning: gate ", str(g), ", payload ", str(j), " might not be unitary!")
+                if np.abs(op[0][1] - row[1]) > 1e-5:
+                    print("Warning: gate ", str(g), ", payload ", str(j), " might not be unitary!")
+
+                row = []
+                for _ in range(2):
+                    amp = tokens[i].replace("(","").replace(")","").split(',')
+                    row.append(float(amp[0]) + float(amp[1])*1j)
+                    i = i + 1
+                l = math.sqrt(np.real(row[0] * np.conj(row[0]) + row[1] * np.conj(row[1])))
+                op[1][0] = row[0] / l
+                op[1][1] = row[1] / l
+
+                ph = np.real(np.log(np.linalg.det(op)) / 1j)
+
+                op[1][0] = -np.exp(1j * ph) * np.conj(op[0][1])
+                op[1][1] = np.exp(1j * ph) * np.conj(op[0][0])
+
+                if np.abs(op[1][0] - row[0]) > 1e-5:
+                    print("Warning: gate ", str(g), ", payload ", str(j), " might not be unitary!")
+                if np.abs(op[1][1] - row[1]) > 1e-5:
+                    print("Warning: gate ", str(g), ", payload ", str(j), " might not be unitary!")
+
+                # Qiskit has a lower tolerance for deviation from numerically unitary.
+                payloads[key] = np.array(op)
+
+            gate_list = []
+            for j in range(1 << control_count):
+                if j in payloads:
+                    gate_list.append(payloads[j])
+                else:
+                    gate_list.append(np.array([[1, 0],[0, 1]]))
+
+            circ.uc(gate_list, controls, target)
+
+        return circ
+
+    def in_from_qiskit_circuit(circ):
+        """Read a Qiskit circuit into a QrackCircuit
+
+        Reads in a circuit from a Qiskit `QuantumCircuit`
+
+        Args:
+            circ: Qiskit circuit
+
+        Raises:
+            RuntimeErorr: Before trying to file_to_qiskit_circuit() with
+                QrackCircuit, you must install Qiskit, numpy, and math!
+        """
+        if not _IS_QISKIT_AVAILABLE:
+            raise RuntimeError(
+                "Before trying to file_to_qiskit_circuit() with QrackCircuit, you must install Qiskit, numpy, and math!"
+            )
+
+        out = QrackCircuit()
+
+        basis_gates = ["u", "cx"]
+        circ = transpile(circ, basis_gates=basis_gates, optimization_level=3)
+        for gate in circ.data:
+            o = gate.operation
+            if o.name == "u":
+                th = float(o.params[0])
+                ph = float(o.params[1])
+                lm = float(o.params[2])
+
+                c = math.cos(th / 2)
+                s = math.sin(th / 2)
+
+                op = [
+                    c + 0j,
+                    -np.exp(1j * lm) * s,
+                    np.exp(1j * ph) * s,
+                    np.exp(1j * (ph + lm)) * c
+                ]
+                out.mtrx(op, circ.find_bit(gate.qubits[0])[0])
+            else:
+                ctrls = []
+                for c in gate.qubits[0:1]:
+                    ctrls.append(circ.find_bit(c)[0])
+                out.ucmtrx(ctrls, [0, 1, 1, 0], circ.find_bit(gate.qubits[1])[0], 1)
+
+        return out
+
+    def file_to_quimb_circuit(
+        filename,
+        circuit_type=QuimbCircuitType.Circuit,
+        psi0=None,
+        gate_opts=None,
+        tags=None,
+        psi0_dtype='complex128',
+        psi0_tag='PSI0',
+        bra_site_ind_id='b{}'
+    ):
+        """Convert an output file to a Quimb circuit
+
+        Reads in an (optimized) circuit from a file named
+        according to the "filename" parameter and outputs
+        a Quimb circuit.
+
+        Args:
+            filename: Name of file
+            circuit_type: "QuimbCircuitType" enum value specifying type of Quimb circuit
+            psi0: The initial state, assumed to be |00000....0> if not given. The state is always copied and the tag PSI0 added
+            gate_opts: Default keyword arguments to supply to each gate_TN_1D() call during the circuit
+            tags: Tag(s) to add to the initial wavefunction tensors (whether these are propagated to the rest of the circuit’s tensors
+            psi0_dtype: Ensure the initial state has this dtype.
+            psi0_tag: Ensure the initial state has this tag.
+            bra_site_ind_id: Use this to label ‘bra’ site indices when creating certain (mostly internal) intermediate tensor networks.
+
+        Raises:
+            RuntimeErorr: Before trying to file_to_quimb_circuit() with
+                QrackCircuit, you must install quimb, Qiskit, numpy, and math!
+        """
+        if not _IS_QUIMB_AVAILABLE:
+            raise RuntimeError(
+                "Before trying to file_to_quimb_circuit() with QrackCircuit, you must install quimb, Qiskit, numpy, and math!"
+            )
+
+        qcirc = QrackCircuit.file_to_qiskit_circuit(filename)
+        basis_gates = ["u", "cx"]
+        qcirc = transpile(qcirc, basis_gates=basis_gates, optimization_level=3)
+
+        tcirc = qtn.Circuit(
+            N=qcirc.num_qubits,
+            psi0=psi0,
+            gate_opts=gate_opts,
+            tags=tags,
+            psi0_dtype=psi0_dtype,
+            psi0_tag=psi0_tag,
+            bra_site_ind_id=bra_site_ind_id
+        ) if circuit_type == QuimbCircuitType.Circuit else (
+            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,
+                    psi0=psi0,
+                    gate_opts=gate_opts,
+                    tags=tags,
+                    psi0_dtype=psi0_dtype,
+                    psi0_tag=psi0_tag,
+                    bra_site_ind_id=bra_site_ind_id
+                )
+        )
+        for gate in qcirc.data:
+            o = gate.operation
+            if o.name == "u":
+                th = float(o.params[0])
+                ph = float(o.params[1])
+                lm = float(o.params[2])
+
+                tcirc.apply_gate('U3', th, ph, lm, qcirc.find_bit(gate.qubits[0])[0])
+            else:
+                tcirc.apply_gate('CNOT', qcirc.find_bit(gate.qubits[0])[0], qcirc.find_bit(gate.qubits[1])[0])
+
+        return tcirc
+
+    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
+        according to the "filename" parameter and outputs
+        a TensorCircuit circuit.
+
+        Args:
+            filename: Name of file
+            inputs: pass-through to tensorcircuit.Circuit.from_qiskit
+            circuit_params: pass-through to tensorcircuit.Circuit.from_qiskit
+            binding_params: pass-through to tensorcircuit.Circuit.from_qiskit
+
+        Raises:
+            RuntimeErorr: Before trying to file_to_quimb_circuit() with
+                QrackCircuit, you must install TensorCircuit, Qiskit, numpy, and math!
+        """
+        if not _IS_TENSORCIRCUIT_AVAILABLE:
+            raise RuntimeError(
+                "Before trying to file_to_tensorcircuit() with QrackCircuit, you must install TensorCircuit, Qiskit, numpy, and math!"
+            )
+
+        qcirc = QrackCircuit.file_to_qiskit_circuit(filename)
+        basis_gates = ["u", "cx"]
+        qcirc = transpile(qcirc, basis_gates=basis_gates, optimization_level=3)
+
+        return tc.Circuit.from_qiskit(qcirc, qcirc.num_qubits, inputs, circuit_params, binding_params)
+
+    def in_from_tensorcircuit(tcirc, enable_instruction = False, enable_inputs = False):
+        """Convert a TensorCircuit circuit to a QrackCircuit
+
+        Accepts a TensorCircuit circuit and outputs an equivalent QrackCircuit
+
+        Args:
+            tcirc: TensorCircuit circuit
+            enable_instruction: whether to also export measurement and reset instructions
+            enable_inputs: whether to also export the inputs
+
+        Raises:
+            RuntimeErorr: Before trying to in_from_tensorcircuit() with
+                QrackCircuit, you must install TensorCircuit, Qiskit, numpy, and math!
+        """
+        if not _IS_TENSORCIRCUIT_AVAILABLE:
+            raise RuntimeError(
+                "Before trying to in_from_tensorcircuit() with QrackCircuit, you must install TensorCircuit, Qiskit, numpy, and math!"
+            )
+
+        # Convert from TensorCircuit to Qiskit
+        qcirc = tcirc.to_qiskit(enable_instruction, enable_inputs)
+
+        # Convert to QrackCircuit
+        return QrackCircuit.in_from_qiskit_circuit(qcirc)