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

pyqrack/qrack_simulator.py

@@ -60,6 +60,7 @@ class QrackSimulator:
         isHostPointer=(
             True if os.environ.get("PYQRACK_HOST_POINTER_DEFAULT_ON") else False
         ),
+        isSparse=False,
         noise=0,
         pyzxCircuit=None,
         qiskitCircuit=None,
@@ -99,6 +100,7 @@ class QrackSimulator:
                 isCpuGpuHybrid,
                 isOpenCL,
                 isHostPointer,
+                isSparse
             )
 
         self._throw_if_error()
@@ -116,36 +118,45 @@ class QrackSimulator:
             Qrack.qrack_lib.destroy(self.sid)
             self.sid = None
 
-    def _int_byref(self, a):
+    @staticmethod
+    def _int_byref(a):
         return (ctypes.c_int * len(a))(*a)
 
-    def _ulonglong_byref(self, a):
+    @staticmethod
+    def _ulonglong_byref(a):
         return (ctypes.c_ulonglong * len(a))(*a)
 
-    def _longlong_byref(self, a):
+    @staticmethod
+    def _longlong_byref(a):
         return (ctypes.c_longlong * 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 QrackSimulator._double_byref([float(item) for sublist in t for item in sublist])
 
-    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 _bool_byref(self, a):
+    @staticmethod
+    def _bool_byref(a):
         return (ctypes.c_bool * len(a))(*a)
 
-    def _qrack_complex_byref(self, a):
+    @staticmethod
+    def _qrack_complex_byref(a):
         t = [(c.real, c.imag) for c in a]
-        return self._real1_byref([float(item) for sublist in t for item in sublist])
+        return QrackSimulator._real1_byref([float(item) for sublist in t for item in sublist])
 
-    def _to_ubyte(self, nv, v):
+    @staticmethod
+    def _to_ubyte(nv, v):
         c = math.floor((nv - 1) / 8) + 1
         b = (ctypes.c_ubyte * (c * (1 << nv)))()
         n = 0
@@ -157,7 +168,8 @@ class QrackSimulator:
 
         return b
 
-    def _to_ulonglong(self, m, v):
+    @staticmethod
+    def _to_ulonglong(m, v):
         b = (ctypes.c_ulonglong * (m * len(v)))()
         n = 0
         for u in v:
@@ -169,7 +181,8 @@ class QrackSimulator:
         return b
 
     # See https://stackoverflow.com/questions/5389507/iterating-over-every-two-elements-in-a-list#answer-30426000
-    def _pairwise(self, it):
+    @staticmethod
+    def _pairwise(it):
         it = iter(it)
         while True:
             try:
@@ -195,7 +208,7 @@ class QrackSimulator:
 
     def set_device_list(self, d):
         """Set the GPU device ID"""
-        Qrack.qrack_lib.set_device_list(self.sid, len(d), self._longlong_byref(d))
+        Qrack.qrack_lib.set_device_list(self.sid, len(d), QrackSimulator._longlong_byref(d))
         self._throw_if_error()
 
     def clone(self):
@@ -359,7 +372,7 @@ class QrackSimulator:
             raise ValueError(
                 "2x2 matrix 'm' in QrackSimulator.mtrx() must contain at least 4 elements."
             )
-        Qrack.qrack_lib.Mtrx(self.sid, self._complex_byref(m), q)
+        Qrack.qrack_lib.Mtrx(self.sid, QrackSimulator._complex_byref(m), q)
         self._throw_if_error()
 
     def r(self, b, ph, q):
@@ -399,9 +412,9 @@ class QrackSimulator:
         Qrack.qrack_lib.Exp(
             self.sid,
             len(b),
-            self._ulonglong_byref(b),
+            QrackSimulator._ulonglong_byref(b),
             ctypes.c_double(ph),
-            self._ulonglong_byref(q),
+            QrackSimulator._ulonglong_byref(q),
         )
         self._throw_if_error()
 
@@ -418,7 +431,7 @@ class QrackSimulator:
         Raises:
             RuntimeError: QrackSimulator raised an exception.
         """
-        Qrack.qrack_lib.MCX(self.sid, len(c), self._ulonglong_byref(c), q)
+        Qrack.qrack_lib.MCX(self.sid, len(c), QrackSimulator._ulonglong_byref(c), q)
         self._throw_if_error()
 
     def mcy(self, c, q):
@@ -434,7 +447,7 @@ class QrackSimulator:
         Raises:
             RuntimeError: QrackSimulator raised an exception.
         """
-        Qrack.qrack_lib.MCY(self.sid, len(c), self._ulonglong_byref(c), q)
+        Qrack.qrack_lib.MCY(self.sid, len(c), QrackSimulator._ulonglong_byref(c), q)
         self._throw_if_error()
 
     def mcz(self, c, q):
@@ -450,7 +463,7 @@ class QrackSimulator:
         Raises:
             RuntimeError: QrackSimulator raised an exception.
         """
-        Qrack.qrack_lib.MCZ(self.sid, len(c), self._ulonglong_byref(c), q)
+        Qrack.qrack_lib.MCZ(self.sid, len(c), QrackSimulator._ulonglong_byref(c), q)
         self._throw_if_error()
 
     def mch(self, c, q):
@@ -466,7 +479,7 @@ class QrackSimulator:
         Raises:
             RuntimeError: QrackSimulator raised an exception.
         """
-        Qrack.qrack_lib.MCH(self.sid, len(c), self._ulonglong_byref(c), q)
+        Qrack.qrack_lib.MCH(self.sid, len(c), QrackSimulator._ulonglong_byref(c), q)
         self._throw_if_error()
 
     def mcs(self, c, q):
@@ -482,7 +495,7 @@ class QrackSimulator:
         Raises:
             RuntimeError: QrackSimulator raised an exception.
         """
-        Qrack.qrack_lib.MCS(self.sid, len(c), self._ulonglong_byref(c), q)
+        Qrack.qrack_lib.MCS(self.sid, len(c), QrackSimulator._ulonglong_byref(c), q)
         self._throw_if_error()
 
     def mct(self, c, q):
@@ -498,7 +511,7 @@ class QrackSimulator:
         Raises:
             RuntimeError: QrackSimulator raised an exception.
         """
-        Qrack.qrack_lib.MCT(self.sid, len(c), self._ulonglong_byref(c), q)
+        Qrack.qrack_lib.MCT(self.sid, len(c), QrackSimulator._ulonglong_byref(c), q)
         self._throw_if_error()
 
     def mcadjs(self, c, q):
@@ -514,7 +527,7 @@ class QrackSimulator:
         Raises:
             RuntimeError: QrackSimulator raised an exception.
         """
-        Qrack.qrack_lib.MCAdjS(self.sid, len(c), self._ulonglong_byref(c), q)
+        Qrack.qrack_lib.MCAdjS(self.sid, len(c), QrackSimulator._ulonglong_byref(c), q)
         self._throw_if_error()
 
     def mcadjt(self, c, q):
@@ -530,7 +543,7 @@ class QrackSimulator:
         Raises:
             RuntimeError: QrackSimulator raised an exception.
         """
-        Qrack.qrack_lib.MCAdjT(self.sid, len(c), self._ulonglong_byref(c), q)
+        Qrack.qrack_lib.MCAdjT(self.sid, len(c), QrackSimulator._ulonglong_byref(c), q)
         self._throw_if_error()
 
     def mcu(self, c, q, th, ph, la):
@@ -552,7 +565,7 @@ class QrackSimulator:
         Qrack.qrack_lib.MCU(
             self.sid,
             len(c),
-            self._ulonglong_byref(c),
+            QrackSimulator._ulonglong_byref(c),
             q,
             ctypes.c_double(th),
             ctypes.c_double(ph),
@@ -580,7 +593,7 @@ class QrackSimulator:
                 "2x2 matrix 'm' in QrackSimulator.mcmtrx() must contain at least 4 elements."
             )
         Qrack.qrack_lib.MCMtrx(
-            self.sid, len(c), self._ulonglong_byref(c), self._complex_byref(m), q
+            self.sid, len(c), QrackSimulator._ulonglong_byref(c), QrackSimulator._complex_byref(m), q
         )
         self._throw_if_error()
 
@@ -596,7 +609,7 @@ class QrackSimulator:
         Raises:
             RuntimeError: QrackSimulator raised an exception.
         """
-        Qrack.qrack_lib.MACX(self.sid, len(c), self._ulonglong_byref(c), q)
+        Qrack.qrack_lib.MACX(self.sid, len(c), QrackSimulator._ulonglong_byref(c), q)
         self._throw_if_error()
 
     def macy(self, c, q):
@@ -612,7 +625,7 @@ class QrackSimulator:
         Raises:
             RuntimeError: QrackSimulator raised an exception.
         """
-        Qrack.qrack_lib.MACY(self.sid, len(c), self._ulonglong_byref(c), q)
+        Qrack.qrack_lib.MACY(self.sid, len(c), QrackSimulator._ulonglong_byref(c), q)
         self._throw_if_error()
 
     def macz(self, c, q):
@@ -628,7 +641,7 @@ class QrackSimulator:
         Raises:
             RuntimeError: QrackSimulator raised an exception.
         """
-        Qrack.qrack_lib.MACZ(self.sid, len(c), self._ulonglong_byref(c), q)
+        Qrack.qrack_lib.MACZ(self.sid, len(c), QrackSimulator._ulonglong_byref(c), q)
         self._throw_if_error()
 
     def mach(self, c, q):
@@ -644,7 +657,7 @@ class QrackSimulator:
         Raises:
             RuntimeError: QrackSimulator raised an exception.
         """
-        Qrack.qrack_lib.MACH(self.sid, len(c), self._ulonglong_byref(c), q)
+        Qrack.qrack_lib.MACH(self.sid, len(c), QrackSimulator._ulonglong_byref(c), q)
         self._throw_if_error()
 
     def macs(self, c, q):
@@ -660,7 +673,7 @@ class QrackSimulator:
         Raises:
             RuntimeError: QrackSimulator raised an exception.
         """
-        Qrack.qrack_lib.MACS(self.sid, len(c), self._ulonglong_byref(c), q)
+        Qrack.qrack_lib.MACS(self.sid, len(c), QrackSimulator._ulonglong_byref(c), q)
         self._throw_if_error()
 
     def mact(self, c, q):
@@ -676,7 +689,7 @@ class QrackSimulator:
         Raises:
             RuntimeError: QrackSimulator raised an exception.
         """
-        Qrack.qrack_lib.MACT(self.sid, len(c), self._ulonglong_byref(c), q)
+        Qrack.qrack_lib.MACT(self.sid, len(c), QrackSimulator._ulonglong_byref(c), q)
         self._throw_if_error()
 
     def macadjs(self, c, q):
@@ -692,7 +705,7 @@ class QrackSimulator:
         Raises:
             RuntimeError: QrackSimulator raised an exception.
         """
-        Qrack.qrack_lib.MACAdjS(self.sid, len(c), self._ulonglong_byref(c), q)
+        Qrack.qrack_lib.MACAdjS(self.sid, len(c), QrackSimulator._ulonglong_byref(c), q)
         self._throw_if_error()
 
     def macadjt(self, c, q):
@@ -708,7 +721,7 @@ class QrackSimulator:
         Raises:
             RuntimeError: QrackSimulator raised an exception.
         """
-        Qrack.qrack_lib.MACAdjT(self.sid, len(c), self._ulonglong_byref(c), q)
+        Qrack.qrack_lib.MACAdjT(self.sid, len(c), QrackSimulator._ulonglong_byref(c), q)
         self._throw_if_error()
 
     def macu(self, c, q, th, ph, la):
@@ -730,7 +743,7 @@ class QrackSimulator:
         Qrack.qrack_lib.MACU(
             self.sid,
             len(c),
-            self._ulonglong_byref(c),
+            QrackSimulator._ulonglong_byref(c),
             q,
             ctypes.c_double(th),
             ctypes.c_double(ph),
@@ -758,7 +771,7 @@ class QrackSimulator:
                 "2x2 matrix 'm' in QrackSimulator.macmtrx() must contain at least 4 elements."
             )
         Qrack.qrack_lib.MACMtrx(
-            self.sid, len(c), self._ulonglong_byref(c), self._complex_byref(m), q
+            self.sid, len(c), QrackSimulator._ulonglong_byref(c), QrackSimulator._complex_byref(m), q
         )
         self._throw_if_error()
 
@@ -783,7 +796,7 @@ class QrackSimulator:
                 "2x2 matrix 'm' in QrackSimulator.ucmtrx() must contain at least 4 elements."
             )
         Qrack.qrack_lib.UCMtrx(
-            self.sid, len(c), self._ulonglong_byref(c), self._complex_byref(m), q, p
+            self.sid, len(c), QrackSimulator._ulonglong_byref(c), QrackSimulator._complex_byref(m), q, p
         )
         self._throw_if_error()
 
@@ -807,7 +820,7 @@ class QrackSimulator:
                 "Multiplex matrix 'm' in QrackSimulator.multiplex1_mtrx() must contain at least (4 * 2 ** len(c)) elements."
             )
         Qrack.qrack_lib.Multiplex1Mtrx(
-            self.sid, len(c), self._ulonglong_byref(c), q, self._complex_byref(m)
+            self.sid, len(c), QrackSimulator._ulonglong_byref(c), q, QrackSimulator._complex_byref(m)
         )
         self._throw_if_error()
 
@@ -822,7 +835,7 @@ class QrackSimulator:
         Raises:
             RuntimeError: QrackSimulator raised an exception.
         """
-        Qrack.qrack_lib.MX(self.sid, len(q), self._ulonglong_byref(q))
+        Qrack.qrack_lib.MX(self.sid, len(q), QrackSimulator._ulonglong_byref(q))
         self._throw_if_error()
 
     def my(self, q):
@@ -836,7 +849,7 @@ class QrackSimulator:
         Raises:
             RuntimeError: QrackSimulator raised an exception.
         """
-        Qrack.qrack_lib.MY(self.sid, len(q), self._ulonglong_byref(q))
+        Qrack.qrack_lib.MY(self.sid, len(q), QrackSimulator._ulonglong_byref(q))
         self._throw_if_error()
 
     def mz(self, q):
@@ -850,7 +863,7 @@ class QrackSimulator:
         Raises:
             RuntimeError: QrackSimulator raised an exception.
         """
-        Qrack.qrack_lib.MZ(self.sid, len(q), self._ulonglong_byref(q))
+        Qrack.qrack_lib.MZ(self.sid, len(q), QrackSimulator._ulonglong_byref(q))
         self._throw_if_error()
 
     def mcr(self, b, ph, c, q):
@@ -873,7 +886,7 @@ class QrackSimulator:
             ctypes.c_ulonglong(b),
             ctypes.c_double(ph),
             len(c),
-            self._ulonglong_byref(c),
+            QrackSimulator._ulonglong_byref(c),
             q,
         )
         self._throw_if_error()
@@ -897,11 +910,11 @@ class QrackSimulator:
         Qrack.qrack_lib.MCExp(
             self.sid,
             len(b),
-            self._ulonglong_byref(b),
+            QrackSimulator._ulonglong_byref(b),
             ctypes.c_double(ph),
             len(cs),
-            self._ulonglong_byref(cs),
-            self._ulonglong_byref(q),
+            QrackSimulator._ulonglong_byref(cs),
+            QrackSimulator._ulonglong_byref(q),
         )
         self._throw_if_error()
 
@@ -983,7 +996,7 @@ class QrackSimulator:
         Raises:
             RuntimeError: QrackSimulator raised an exception.
         """
-        Qrack.qrack_lib.CSWAP(self.sid, len(c), self._ulonglong_byref(c), qi1, qi2)
+        Qrack.qrack_lib.CSWAP(self.sid, len(c), QrackSimulator._ulonglong_byref(c), qi1, qi2)
         self._throw_if_error()
 
     def acswap(self, c, qi1, qi2):
@@ -999,7 +1012,7 @@ class QrackSimulator:
         Raises:
             RuntimeError: QrackSimulator raised an exception.
         """
-        Qrack.qrack_lib.ACSWAP(self.sid, len(c), self._ulonglong_byref(c), qi1, qi2)
+        Qrack.qrack_lib.ACSWAP(self.sid, len(c), QrackSimulator._ulonglong_byref(c), qi1, qi2)
         self._throw_if_error()
 
     # standard operations
@@ -1053,9 +1066,16 @@ class QrackSimulator:
         Returns:
             Measurement result of all qubits.
         """
-        result = Qrack.qrack_lib.MAll(self.sid)
+        num_q = self.num_qubits()
+        num_words = (num_q + 63) // 64
+        _r = (ctypes.c_ulonglong * num_words)()
+        Qrack.qrack_lib.MAllLong(self.sid, _r)
         self._throw_if_error()
-        return result
+        r = 0
+        for w in range(num_words):
+            r <<= 64
+            r |= _r[w]
+        return r
 
     def measure_pauli(self, b, q):
         """Pauli Measurement gate
@@ -1076,7 +1096,7 @@ class QrackSimulator:
         if len(b) != len(q):
             raise RuntimeError("Lengths of list parameters are mismatched.")
         result = Qrack.qrack_lib.Measure(
-            self.sid, len(b), self._int_byref(b), self._ulonglong_byref(q)
+            self.sid, len(b), QrackSimulator._int_byref(b), QrackSimulator._ulonglong_byref(q)
         )
         self._throw_if_error()
         return result
@@ -1097,8 +1117,8 @@ class QrackSimulator:
         Returns:
             list of measurement result.
         """
-        m = self._ulonglong_byref([0] * s)
-        Qrack.qrack_lib.MeasureShots(self.sid, len(q), self._ulonglong_byref(q), s, m)
+        m = QrackSimulator._ulonglong_byref([0] * s)
+        Qrack.qrack_lib.MeasureShots(self.sid, len(q), QrackSimulator._ulonglong_byref(q), s, m)
         self._throw_if_error()
         return [m[i] for i in range(s)]
 
@@ -1114,7 +1134,8 @@ class QrackSimulator:
         self._throw_if_error()
 
     # arithmetic-logic-unit (ALU)
-    def _split_longs(self, a):
+    @staticmethod
+    def _split_longs(a):
         """Split operation
 
         Splits the given integer into 64 bit numbers.
@@ -1137,7 +1158,8 @@ class QrackSimulator:
             a = a >> 64
         return aParts
 
-    def _split_longs_2(self, a, m):
+    @staticmethod
+    def _split_longs_2(a, m):
         """Split simultanoues operation
 
         Splits 2 integers into same number of 64 bit numbers.
@@ -1176,13 +1198,13 @@ class QrackSimulator:
         Raises:
             RuntimeError: QrackSimulator raised an exception.
         """
-        aParts = self._split_longs(a)
+        aParts = QrackSimulator._split_longs(a)
         Qrack.qrack_lib.ADD(
             self.sid,
             len(aParts),
-            self._ulonglong_byref(aParts),
+            QrackSimulator._ulonglong_byref(aParts),
             len(q),
-            self._ulonglong_byref(q),
+            QrackSimulator._ulonglong_byref(q),
         )
         self._throw_if_error()
 
@@ -1198,13 +1220,13 @@ class QrackSimulator:
         Raises:
             RuntimeError: QrackSimulator raised an exception.
         """
-        aParts = self._split_longs(a)
+        aParts = QrackSimulator._split_longs(a)
         Qrack.qrack_lib.SUB(
             self.sid,
             len(aParts),
-            self._ulonglong_byref(aParts),
+            QrackSimulator._ulonglong_byref(aParts),
             len(q),
-            self._ulonglong_byref(q),
+            QrackSimulator._ulonglong_byref(q),
         )
         self._throw_if_error()
 
@@ -1222,14 +1244,14 @@ class QrackSimulator:
         Raises:
             RuntimeError: QrackSimulator raised an exception.
         """
-        aParts = self._split_longs(a)
+        aParts = QrackSimulator._split_longs(a)
         Qrack.qrack_lib.ADDS(
             self.sid,
             len(aParts),
-            self._ulonglong_byref(aParts),
+            QrackSimulator._ulonglong_byref(aParts),
             s,
             len(q),
-            self._ulonglong_byref(q),
+            QrackSimulator._ulonglong_byref(q),
         )
         self._throw_if_error()
 
@@ -1247,14 +1269,14 @@ class QrackSimulator:
         Raises:
             RuntimeError: QrackSimulator raised an exception.
         """
-        aParts = self._split_longs(a)
+        aParts = QrackSimulator._split_longs(a)
         Qrack.qrack_lib.SUBS(
             self.sid,
             len(aParts),
-            self._ulonglong_byref(aParts),
+            QrackSimulator._ulonglong_byref(aParts),
             s,
             len(q),
-            self._ulonglong_byref(q),
+            QrackSimulator._ulonglong_byref(q),
         )
         self._throw_if_error()
 
@@ -1285,14 +1307,14 @@ class QrackSimulator:
 
         if len(q) != len(o):
             raise RuntimeError("Lengths of list parameters are mismatched.")
-        aParts = self._split_longs(a)
+        aParts = QrackSimulator._split_longs(a)
         Qrack.qrack_lib.MUL(
             self.sid,
             len(aParts),
-            self._ulonglong_byref(aParts),
+            QrackSimulator._ulonglong_byref(aParts),
             len(q),
-            self._ulonglong_byref(q),
-            self._ulonglong_byref(o),
+            QrackSimulator._ulonglong_byref(q),
+            QrackSimulator._ulonglong_byref(o),
         )
         self._throw_if_error()
 
@@ -1324,14 +1346,14 @@ class QrackSimulator:
 
         if len(q) != len(o):
             raise RuntimeError("Lengths of list parameters are mismatched.")
-        aParts = self._split_longs(a)
+        aParts = QrackSimulator._split_longs(a)
         Qrack.qrack_lib.DIV(
             self.sid,
             len(aParts),
-            self._ulonglong_byref(aParts),
+            QrackSimulator._ulonglong_byref(aParts),
             len(q),
-            self._ulonglong_byref(q),
-            self._ulonglong_byref(o),
+            QrackSimulator._ulonglong_byref(q),
+            QrackSimulator._ulonglong_byref(o),
         )
         self._throw_if_error()
 
@@ -1352,15 +1374,15 @@ class QrackSimulator:
         """
         if len(q) != len(o):
             raise RuntimeError("Lengths of list parameters are mismatched.")
-        aParts, mParts = self._split_longs_2(a, m)
+        aParts, mParts = QrackSimulator._split_longs_2(a, m)
         Qrack.qrack_lib.MULN(
             self.sid,
             len(aParts),
-            self._ulonglong_byref(aParts),
-            self._ulonglong_byref(mParts),
+            QrackSimulator._ulonglong_byref(aParts),
+            QrackSimulator._ulonglong_byref(mParts),
             len(q),
-            self._ulonglong_byref(q),
-            self._ulonglong_byref(o),
+            QrackSimulator._ulonglong_byref(q),
+            QrackSimulator._ulonglong_byref(o),
         )
         self._throw_if_error()
 
@@ -1382,15 +1404,15 @@ class QrackSimulator:
         """
         if len(q) != len(o):
             raise RuntimeError("Lengths of list parameters are mismatched.")
-        aParts, mParts = self._split_longs_2(a, m)
+        aParts, mParts = QrackSimulator._split_longs_2(a, m)
         Qrack.qrack_lib.DIVN(
             self.sid,
             len(aParts),
-            self._ulonglong_byref(aParts),
-            self._ulonglong_byref(mParts),
+            QrackSimulator._ulonglong_byref(aParts),
+            QrackSimulator._ulonglong_byref(mParts),
             len(q),
-            self._ulonglong_byref(q),
-            self._ulonglong_byref(o),
+            QrackSimulator._ulonglong_byref(q),
+            QrackSimulator._ulonglong_byref(o),
         )
         self._throw_if_error()
 
@@ -1421,15 +1443,15 @@ class QrackSimulator:
 
         if len(q) != len(o):
             raise RuntimeError("Lengths of list parameters are mismatched.")
-        aParts, mParts = self._split_longs_2(a, m)
+        aParts, mParts = QrackSimulator._split_longs_2(a, m)
         Qrack.qrack_lib.POWN(
             self.sid,
             len(aParts),
-            self._ulonglong_byref(aParts),
-            self._ulonglong_byref(mParts),
+            QrackSimulator._ulonglong_byref(aParts),
+            QrackSimulator._ulonglong_byref(mParts),
             len(q),
-            self._ulonglong_byref(q),
-            self._ulonglong_byref(o),
+            QrackSimulator._ulonglong_byref(q),
+            QrackSimulator._ulonglong_byref(o),
         )
         self._throw_if_error()
 
@@ -1447,15 +1469,15 @@ class QrackSimulator:
         Raises:
             RuntimeError: QrackSimulator raised an exception.
         """
-        aParts = self._split_longs(a)
+        aParts = QrackSimulator._split_longs(a)
         Qrack.qrack_lib.MCADD(
             self.sid,
             len(aParts),
-            self._ulonglong_byref(aParts),
+            QrackSimulator._ulonglong_byref(aParts),
             len(c),
-            self._ulonglong_byref(c),
+            QrackSimulator._ulonglong_byref(c),
             len(q),
-            self._ulonglong_byref(q),
+            QrackSimulator._ulonglong_byref(q),
         )
         self._throw_if_error()
 
@@ -1473,15 +1495,15 @@ class QrackSimulator:
         Raises:
             RuntimeError: QrackSimulator raised an exception.
         """
-        aParts = self._split_longs(a)
+        aParts = QrackSimulator._split_longs(a)
         Qrack.qrack_lib.MCSUB(
             self.sid,
             len(aParts),
-            self._ulonglong_byref(aParts),
+            QrackSimulator._ulonglong_byref(aParts),
             len(c),
-            self._ulonglong_byref(c),
+            QrackSimulator._ulonglong_byref(c),
             len(q),
-            self._ulonglong_byref(q),
+            QrackSimulator._ulonglong_byref(q),
         )
         self._throw_if_error()
 
@@ -1514,15 +1536,15 @@ class QrackSimulator:
 
         if len(q) != len(o):
             raise RuntimeError("Lengths of list parameters are mismatched.")
-        aParts = self._split_longs(a)
+        aParts = QrackSimulator._split_longs(a)
         Qrack.qrack_lib.MCMUL(
             self.sid,
             len(aParts),
-            self._ulonglong_byref(aParts),
+            QrackSimulator._ulonglong_byref(aParts),
             len(c),
-            self._ulonglong_byref(c),
+            QrackSimulator._ulonglong_byref(c),
             len(q),
-            self._ulonglong_byref(q),
+            QrackSimulator._ulonglong_byref(q),
         )
         self._throw_if_error()
 
@@ -1556,15 +1578,15 @@ class QrackSimulator:
 
         if len(q) != len(o):
             raise RuntimeError("Lengths of list parameters are mismatched.")
-        aParts = self._split_longs(a)
+        aParts = QrackSimulator._split_longs(a)
         Qrack.qrack_lib.MCDIV(
             self.sid,
             len(aParts),
-            self._ulonglong_byref(aParts),
+            QrackSimulator._ulonglong_byref(aParts),
             len(c),
-            self._ulonglong_byref(c),
+            QrackSimulator._ulonglong_byref(c),
             len(q),
-            self._ulonglong_byref(q),
+            QrackSimulator._ulonglong_byref(q),
         )
         self._throw_if_error()
 
@@ -1587,17 +1609,17 @@ class QrackSimulator:
         """
         if len(q) != len(o):
             raise RuntimeError("Lengths of list parameters are mismatched.")
-        aParts, mParts = self._split_longs_2(a, m)
+        aParts, mParts = QrackSimulator._split_longs_2(a, m)
         Qrack.qrack_lib.MCMULN(
             self.sid,
             len(aParts),
-            self._ulonglong_byref(aParts),
+            QrackSimulator._ulonglong_byref(aParts),
             len(c),
-            self._ulonglong_byref(c),
-            self._ulonglong_byref(mParts),
+            QrackSimulator._ulonglong_byref(c),
+            QrackSimulator._ulonglong_byref(mParts),
             len(q),
-            self._ulonglong_byref(q),
-            self._ulonglong_byref(o),
+            QrackSimulator._ulonglong_byref(q),
+            QrackSimulator._ulonglong_byref(o),
         )
         self._throw_if_error()
 
@@ -1621,17 +1643,17 @@ class QrackSimulator:
         """
         if len(q) != len(o):
             raise RuntimeError("Lengths of list parameters are mismatched.")
-        aParts, mParts = self._split_longs_2(a, m)
+        aParts, mParts = QrackSimulator._split_longs_2(a, m)
         Qrack.qrack_lib.MCDIVN(
             self.sid,
             len(aParts),
-            self._ulonglong_byref(aParts),
+            QrackSimulator._ulonglong_byref(aParts),
             len(c),
-            self._ulonglong_byref(c),
-            self._ulonglong_byref(mParts),
+            QrackSimulator._ulonglong_byref(c),
+            QrackSimulator._ulonglong_byref(mParts),
             len(q),
-            self._ulonglong_byref(q),
-            self._ulonglong_byref(o),
+            QrackSimulator._ulonglong_byref(q),
+            QrackSimulator._ulonglong_byref(o),
         )
         self._throw_if_error()
 
@@ -1664,17 +1686,17 @@ class QrackSimulator:
 
         if len(q) != len(o):
             raise RuntimeError("Lengths of list parameters are mismatched.")
-        aParts, mParts = self._split_longs_2(a, m)
+        aParts, mParts = QrackSimulator._split_longs_2(a, m)
         Qrack.qrack_lib.MCPOWN(
             self.sid,
             len(aParts),
-            self._ulonglong_byref(aParts),
+            QrackSimulator._ulonglong_byref(aParts),
             len(c),
-            self._ulonglong_byref(c),
-            self._ulonglong_byref(mParts),
+            QrackSimulator._ulonglong_byref(c),
+            QrackSimulator._ulonglong_byref(mParts),
             len(q),
-            self._ulonglong_byref(q),
-            self._ulonglong_byref(o),
+            QrackSimulator._ulonglong_byref(q),
+            QrackSimulator._ulonglong_byref(o),
         )
         self._throw_if_error()
 
@@ -1706,10 +1728,10 @@ class QrackSimulator:
         Qrack.qrack_lib.LDA(
             self.sid,
             len(qi),
-            self._ulonglong_byref(qi),
+            QrackSimulator._ulonglong_byref(qi),
             len(qv),
-            self._ulonglong_byref(qv),
-            self._to_ubyte(len(qv), t),
+            QrackSimulator._ulonglong_byref(qv),
+            QrackSimulator._to_ubyte(len(qv), t),
         )
         self._throw_if_error()
 
@@ -1741,10 +1763,10 @@ class QrackSimulator:
             self.sid,
             s,
             len(qi),
-            self._ulonglong_byref(qi),
+            QrackSimulator._ulonglong_byref(qi),
             len(qv),
-            self._ulonglong_byref(qv),
-            self._to_ubyte(len(qv), t),
+            QrackSimulator._ulonglong_byref(qv),
+            QrackSimulator._to_ubyte(len(qv), t),
         )
         self._throw_if_error()
 
@@ -1776,10 +1798,10 @@ class QrackSimulator:
             self.sid,
             s,
             len(qi),
-            self._ulonglong_byref(qi),
+            QrackSimulator._ulonglong_byref(qi),
             len(qv),
-            self._ulonglong_byref(qv),
-            self._to_ubyte(len(qv), t),
+            QrackSimulator._ulonglong_byref(qv),
+            QrackSimulator._to_ubyte(len(qv), t),
         )
         self._throw_if_error()
 
@@ -1809,7 +1831,7 @@ class QrackSimulator:
             )
 
         Qrack.qrack_lib.Hash(
-            self.sid, len(q), self._ulonglong_byref(q), self._to_ubyte(len(q), t)
+            self.sid, len(q), QrackSimulator._ulonglong_byref(q), QrackSimulator._to_ubyte(len(q), t)
         )
         self._throw_if_error()
 
@@ -2030,7 +2052,7 @@ class QrackSimulator:
         Raises:
             RuntimeError: QrackSimulator raised an exception.
         """
-        Qrack.qrack_lib.QFT(self.sid, len(qs), self._ulonglong_byref(qs))
+        Qrack.qrack_lib.QFT(self.sid, len(qs), QrackSimulator._ulonglong_byref(qs))
         self._throw_if_error()
 
     def iqft(self, qs):
@@ -2045,7 +2067,7 @@ class QrackSimulator:
         Raises:
             RuntimeError: QrackSimulator raised an exception.
         """
-        Qrack.qrack_lib.IQFT(self.sid, len(qs), self._ulonglong_byref(qs))
+        Qrack.qrack_lib.IQFT(self.sid, len(qs), QrackSimulator._ulonglong_byref(qs))
         self._throw_if_error()
 
     # pseudo-quantum
@@ -2119,15 +2141,13 @@ class QrackSimulator:
                 "QrackSimulator with isTensorNetwork=True option cannot compose()! (Turn off just this option, in the constructor.)"
             )
 
-        Qrack.qrack_lib.Compose(self.sid, other.sid, self._ulonglong_byref(q))
+        Qrack.qrack_lib.Compose(self.sid, other.sid, QrackSimulator._ulonglong_byref(q))
         self._throw_if_error()
 
     def decompose(self, q):
         """Decompose system
 
-        Decompose the given qubit out of the system.
-        Warning: The qubit subsystem state must be separable, or the behavior
-        of this method is undefined.
+        Factorize a set of contiguous bits with minimal fidelity loss.
 
         Args:
             q: qubit id
@@ -2137,7 +2157,7 @@ class QrackSimulator:
             RuntimeError: QrackSimulator with isTensorNetwork=True option cannot decompose()! (Turn off just this option, in the constructor.)
 
         Returns:
-            State of the systems.
+            Decomposed subsystem simulator.
         """
         if self.is_tensor_network:
             raise RuntimeError(
@@ -2147,17 +2167,15 @@ class QrackSimulator:
         other = QrackSimulator()
         Qrack.qrack_lib.destroy(other.sid)
         l = len(q)
-        other.sid = Qrack.qrack_lib.Decompose(self.sid, l, self._ulonglong_byref(q))
+        other.sid = Qrack.qrack_lib.Decompose(self.sid, l, QrackSimulator._ulonglong_byref(q))
         self._throw_if_error()
         return other
 
     def dispose(self, q):
         """Dispose qubits
 
-        Minimally decompose a set of contiguous bits from the separably
-        composed unit, and discard the separable bits.
-        Warning: The qubit subsystem state must be separable, or the behavior
-        of this method is undefined.
+        Factorize a set of contiguous bits with minimal fidelity loss,
+        and discard the separable bits.
 
         Args:
             q: qubit
@@ -2165,9 +2183,6 @@ class QrackSimulator:
         Raises:
             RuntimeError: QrackSimulator raised an exception.
             RuntimeError: QrackSimulator with isTensorNetwork=True option cannot dispose()! (Turn off just this option, in the constructor.)
-
-        Returns:
-            State of the systems.
         """
         if self.is_tensor_network:
             raise RuntimeError(
@@ -2175,7 +2190,7 @@ class QrackSimulator:
             )
 
         l = len(q)
-        Qrack.qrack_lib.Dispose(self.sid, l, self._ulonglong_byref(q))
+        Qrack.qrack_lib.Dispose(self.sid, l, QrackSimulator._ulonglong_byref(q))
         self._throw_if_error()
 
     ## miscellaneous
@@ -2245,7 +2260,7 @@ class QrackSimulator:
         Raises:
             RuntimeError: QrackSimulator raised an exception.
         """
-        Qrack.qrack_lib.InKet(self.sid, self._qrack_complex_byref(ket))
+        Qrack.qrack_lib.InKet(self.sid, QrackSimulator._qrack_complex_byref(ket))
         self._throw_if_error()
 
     def out_ket(self):
@@ -2262,10 +2277,10 @@ class QrackSimulator:
             list representing the state vector.
         """
         amp_count = 1 << self.num_qubits()
-        ket = self._qrack_complex_byref([complex(0, 0)] * amp_count)
+        ket = QrackSimulator._qrack_complex_byref([complex(0, 0)] * amp_count)
         Qrack.qrack_lib.OutKet(self.sid, ket)
         self._throw_if_error()
-        return [complex(r, i) for r, i in self._pairwise(ket)]
+        return [complex(r, i) for r, i in QrackSimulator._pairwise(ket)]
 
     def out_probs(self):
         """Get basis dimension probabilities
@@ -2280,11 +2295,77 @@ class QrackSimulator:
             list representing the basis dimension probabilities.
         """
         prob_count = 1 << self.num_qubits()
-        probs = self._real1_byref([0.0] * prob_count)
+        probs = QrackSimulator._real1_byref([0.0] * prob_count)
         Qrack.qrack_lib.OutProbs(self.sid, probs)
         self._throw_if_error()
         return list(probs)
 
+    def out_rdm(self, q):
+        """Get reduced density matrix
+
+        Returns the raw reduced density matrix of the simulator, for the qubit list.
+        Warning: State vector or is not always the internal representation leading
+        to sub-optimal performance of the method.
+
+        Raises:
+            RuntimeError: QrackSimulator raised an exception.
+
+        Returns:
+            flat list structure representing the reduced density matrix.
+        """
+        amp_count = 1 << len(q)
+        sqr_amp_count = amp_count * amp_count
+        flat_rdm = QrackSimulator._qrack_complex_byref([complex(0, 0)] * sqr_amp_count)
+        Qrack.qrack_lib.OutReducedDensityMatrix(self.sid, len(q), QrackSimulator._ulonglong_byref(q), flat_rdm)
+        self._throw_if_error()
+        return [complex(r, i) for r, i in QrackSimulator._pairwise(flat_rdm)]
+
+    def highest_prob_perm(self):
+        """Get the permutation (bit string) with the highest probability
+
+        Returns the single highest-probability bit string in the Hilbert space
+
+        Raises:
+            RuntimeError: QrackSimulator raised an exception.
+
+        Returns:
+            Highest probability dimension index
+        """
+        num_q = self.num_qubits()
+        num_words = (num_q + 63) // 64
+        _r = (ctypes.c_ulonglong * num_words)()
+        Qrack.qrack_lib.HighestProbAll(self.sid, _r)
+        self._throw_if_error()
+        r = 0
+        for w in range(num_words):
+            r <<= 64
+            r |= _r[w]
+        return r
+
+    def highest_n_prob_perm(self, n):
+        """Get the top n permutations (bit strings) with the highest probability
+
+        Returns the top n highest-probability bit strings in the Hilbert space
+
+        Raises:
+            RuntimeError: QrackSimulator raised an exception.
+
+        Returns:
+            Top n highest probability dimension indices
+        """
+        num_q = self.num_qubits()
+        num_words = (num_q + 63) // 64
+        _r = (ctypes.c_ulonglong * (num_words * n))()
+        Qrack.qrack_lib.HighestProbAllN(self.sid, n, _r)
+        self._throw_if_error()
+        r = [0] * n
+        for i in range(n):
+            r[i] = 0
+            for w in range(num_words):
+                r[i] <<= 64
+                r[i] |= _r[(i * num_words) + w]
+        return r
+
     def prob_all(self, q):
         """Probabilities of all subset permutations
 
@@ -2299,8 +2380,8 @@ class QrackSimulator:
         Returns:
             list representing the state vector.
         """
-        probs = self._real1_byref([0.0] * (1 << len(q)))
-        Qrack.qrack_lib.ProbAll(self.sid, len(q), self._ulonglong_byref(q), probs)
+        probs = QrackSimulator._real1_byref([0.0] * (1 << len(q)))
+        Qrack.qrack_lib.ProbAll(self.sid, len(q), QrackSimulator._ulonglong_byref(q), probs)
         self._throw_if_error()
         return list(probs)
 
@@ -2362,7 +2443,7 @@ class QrackSimulator:
         if len(q) != len(c):
             raise RuntimeError("prob_perm argument lengths do not match.")
         result = Qrack.qrack_lib.PermutationProb(
-            self.sid, len(q), self._ulonglong_byref(q), self._bool_byref(c)
+            self.sid, len(q), QrackSimulator._ulonglong_byref(q), QrackSimulator._bool_byref(c)
         )
         self._throw_if_error()
         return result
@@ -2390,7 +2471,7 @@ class QrackSimulator:
         if len(q) != len(c):
             raise RuntimeError("prob_perm argument lengths do not match.")
         result = Qrack.qrack_lib.PermutationProbRdm(
-            self.sid, len(q), self._ulonglong_byref(q), self._bool_byref(c), r
+            self.sid, len(q), QrackSimulator._ulonglong_byref(q), QrackSimulator._bool_byref(c), r
         )
         self._throw_if_error()
         return result
@@ -2411,7 +2492,7 @@ class QrackSimulator:
             Expectation value
         """
         result = Qrack.qrack_lib.PermutationExpectation(
-            self.sid, len(q), self._ulonglong_byref(q)
+            self.sid, len(q), QrackSimulator._ulonglong_byref(q)
         )
         self._throw_if_error()
         return result
@@ -2434,7 +2515,7 @@ class QrackSimulator:
             Expectation value
         """
         result = Qrack.qrack_lib.PermutationExpectationRdm(
-            self.sid, len(q), self._ulonglong_byref(q), r
+            self.sid, len(q), QrackSimulator._ulonglong_byref(q), r
         )
         self._throw_if_error()
         return result
@@ -2460,7 +2541,7 @@ class QrackSimulator:
             raise RuntimeError("factorized_expectation argument lengths do not match.")
         m = max([(x.bit_length() + 63) // 64 for x in c])
         result = Qrack.qrack_lib.FactorizedExpectation(
-            self.sid, len(q), self._ulonglong_byref(q), m, self._to_ulonglong(m, c)
+            self.sid, len(q), QrackSimulator._ulonglong_byref(q), m, QrackSimulator._to_ulonglong(m, c)
         )
         self._throw_if_error()
         return result
@@ -2490,7 +2571,7 @@ class QrackSimulator:
             )
         m = max([(x.bit_length() + 63) // 64 for x in c])
         result = Qrack.qrack_lib.FactorizedExpectationRdm(
-            self.sid, len(q), self._ulonglong_byref(q), m, self._to_ulonglong(m, c), r
+            self.sid, len(q), QrackSimulator._ulonglong_byref(q), m, QrackSimulator._to_ulonglong(m, c), r
         )
         self._throw_if_error()
         return result
@@ -2517,7 +2598,7 @@ class QrackSimulator:
                 "factorized_expectation_rdm argument lengths do not match."
             )
         result = Qrack.qrack_lib.FactorizedExpectationFp(
-            self.sid, len(q), self._ulonglong_byref(q), self._real1_byref(c)
+            self.sid, len(q), QrackSimulator._ulonglong_byref(q), QrackSimulator._real1_byref(c)
         )
         self._throw_if_error()
         return result
@@ -2546,7 +2627,7 @@ class QrackSimulator:
                 "factorized_expectation_fp_rdm argument lengths do not match."
             )
         result = Qrack.qrack_lib.FactorizedExpectationFpRdm(
-            self.sid, len(q), self._ulonglong_byref(q), self._real1_byref(c), r
+            self.sid, len(q), QrackSimulator._ulonglong_byref(q), QrackSimulator._real1_byref(c), r
         )
         self._throw_if_error()
         return result
@@ -2570,7 +2651,7 @@ class QrackSimulator:
         if (3 * len(q)) != len(b):
             raise RuntimeError("unitary_expectation argument lengths do not match.")
         result = Qrack.qrack_lib.UnitaryExpectation(
-            self.sid, len(q), self._ulonglong_byref(q), self._real1_byref(b)
+            self.sid, len(q), QrackSimulator._ulonglong_byref(q), QrackSimulator._real1_byref(b)
         )
         self._throw_if_error()
         return result
@@ -2594,7 +2675,7 @@ class QrackSimulator:
         if (len(q) << 2) != len(b):
             raise RuntimeError("matrix_expectation argument lengths do not match.")
         result = Qrack.qrack_lib.MatrixExpectation(
-            self.sid, len(q), self._ulonglong_byref(q), self._complex_byref(b)
+            self.sid, len(q), QrackSimulator._ulonglong_byref(q), QrackSimulator._complex_byref(b)
         )
         self._throw_if_error()
         return result
@@ -2626,9 +2707,9 @@ class QrackSimulator:
         result = Qrack.qrack_lib.UnitaryExpectationEigenVal(
             self.sid,
             len(q),
-            self._ulonglong_byref(q),
-            self._real1_byref(b),
-            self._real1_byref(e),
+            QrackSimulator._ulonglong_byref(q),
+            QrackSimulator._real1_byref(b),
+            QrackSimulator._real1_byref(e),
         )
         self._throw_if_error()
         return result
@@ -2660,9 +2741,9 @@ class QrackSimulator:
         result = Qrack.qrack_lib.MatrixExpectationEigenVal(
             self.sid,
             len(q),
-            self._ulonglong_byref(q),
-            self._complex_byref(b),
-            self._real1_byref(e),
+            QrackSimulator._ulonglong_byref(q),
+            QrackSimulator._complex_byref(b),
+            QrackSimulator._real1_byref(e),
         )
         self._throw_if_error()
         return result
@@ -2687,7 +2768,7 @@ class QrackSimulator:
         if len(q) != len(b):
             raise RuntimeError("pauli_expectation argument lengths do not match.")
         result = Qrack.qrack_lib.PauliExpectation(
-            self.sid, len(q), self._ulonglong_byref(q), self._ulonglong_byref(b)
+            self.sid, len(q), QrackSimulator._ulonglong_byref(q), QrackSimulator._ulonglong_byref(b)
         )
         self._throw_if_error()
         return result
@@ -2707,7 +2788,7 @@ class QrackSimulator:
         Returns:
             float variance
         """
-        result = Qrack.qrack_lib.Variance(self.sid, len(q), self._ulonglong_byref(q))
+        result = Qrack.qrack_lib.Variance(self.sid, len(q), QrackSimulator._ulonglong_byref(q))
         self._throw_if_error()
         return result
 
@@ -2729,7 +2810,7 @@ class QrackSimulator:
             variance
         """
         result = Qrack.qrack_lib.VarianceRdm(
-            self.sid, len(q), self._ulonglong_byref(q), r
+            self.sid, len(q), QrackSimulator._ulonglong_byref(q), r
         )
         self._throw_if_error()
         return result
@@ -2755,7 +2836,7 @@ class QrackSimulator:
             raise RuntimeError("factorized_variance argument lengths do not match.")
         m = max([(x.bit_length() + 63) // 64 for x in c])
         result = Qrack.qrack_lib.FactorizedVariance(
-            self.sid, len(q), self._ulonglong_byref(q), m, self._to_ulonglong(m, c)
+            self.sid, len(q), QrackSimulator._ulonglong_byref(q), m, QrackSimulator._to_ulonglong(m, c)
         )
         self._throw_if_error()
         return result
@@ -2783,7 +2864,7 @@ class QrackSimulator:
             raise RuntimeError("factorized_variance_rdm argument lengths do not match.")
         m = max([(x.bit_length() + 63) // 64 for x in c])
         result = Qrack.qrack_lib.FactorizedVarianceRdm(
-            self.sid, len(q), self._ulonglong_byref(q), m, self._to_ulonglong(m, c), r
+            self.sid, len(q), QrackSimulator._ulonglong_byref(q), m, QrackSimulator._to_ulonglong(m, c), r
         )
         self._throw_if_error()
         return result
@@ -2808,7 +2889,7 @@ class QrackSimulator:
         if (len(q) << 1) != len(c):
             raise RuntimeError("factorized_variance_rdm argument lengths do not match.")
         result = Qrack.qrack_lib.FactorizedVarianceFp(
-            self.sid, len(q), self._ulonglong_byref(q), self._real1_byref(c)
+            self.sid, len(q), QrackSimulator._ulonglong_byref(q), QrackSimulator._real1_byref(c)
         )
         self._throw_if_error()
         return result
@@ -2837,7 +2918,7 @@ class QrackSimulator:
                 "factorized_variance_fp_rdm argument lengths do not match."
             )
         result = Qrack.qrack_lib.FactorizedVarianceFpRdm(
-            self.sid, len(q), self._ulonglong_byref(q), self._real1_byref(c), r
+            self.sid, len(q), QrackSimulator._ulonglong_byref(q), QrackSimulator._real1_byref(c), r
         )
         self._throw_if_error()
         return result
@@ -2861,7 +2942,7 @@ class QrackSimulator:
         if (3 * len(q)) != len(b):
             raise RuntimeError("unitary_variance argument lengths do not match.")
         result = Qrack.qrack_lib.UnitaryVariance(
-            self.sid, len(q), self._ulonglong_byref(q), self._real1_byref(b)
+            self.sid, len(q), QrackSimulator._ulonglong_byref(q), QrackSimulator._real1_byref(b)
         )
         self._throw_if_error()
         return result
@@ -2885,7 +2966,7 @@ class QrackSimulator:
         if (len(q) << 2) != len(b):
             raise RuntimeError("matrix_variance argument lengths do not match.")
         result = Qrack.qrack_lib.MatrixVariance(
-            self.sid, len(q), self._ulonglong_byref(q), self._complex_byref(b)
+            self.sid, len(q), QrackSimulator._ulonglong_byref(q), QrackSimulator._complex_byref(b)
         )
         self._throw_if_error()
         return result
@@ -2917,9 +2998,9 @@ class QrackSimulator:
         result = Qrack.qrack_lib.UnitaryVarianceEigenVal(
             self.sid,
             len(q),
-            self._ulonglong_byref(q),
-            self._real1_byref(b),
-            self._real1_byref(e),
+            QrackSimulator._ulonglong_byref(q),
+            QrackSimulator._real1_byref(b),
+            QrackSimulator._real1_byref(e),
         )
         self._throw_if_error()
         return result
@@ -2951,9 +3032,9 @@ class QrackSimulator:
         result = Qrack.qrack_lib.MatrixVarianceEigenVal(
             self.sid,
             len(q),
-            self._ulonglong_byref(q),
-            self._complex_byref(b),
-            self._real1_byref(e),
+            QrackSimulator._ulonglong_byref(q),
+            QrackSimulator._complex_byref(b),
+            QrackSimulator._real1_byref(e),
         )
         self._throw_if_error()
         return result
@@ -2978,7 +3059,7 @@ class QrackSimulator:
         if len(q) != len(b):
             raise RuntimeError("pauli_variance argument lengths do not match.")
         result = Qrack.qrack_lib.PauliVariance(
-            self.sid, len(q), self._ulonglong_byref(q), self._ulonglong_byref(b)
+            self.sid, len(q), QrackSimulator._ulonglong_byref(q), QrackSimulator._ulonglong_byref(b)
         )
         self._throw_if_error()
         return result
@@ -3002,7 +3083,7 @@ class QrackSimulator:
         if len(b) != len(q):
             raise RuntimeError("Lengths of list parameters are mismatched.")
         result = Qrack.qrack_lib.JointEnsembleProbability(
-            self.sid, len(b), self._ulonglong_byref(b), q
+            self.sid, len(b), QrackSimulator._ulonglong_byref(b), q
         )
         self._throw_if_error()
         return result
@@ -3031,7 +3112,7 @@ class QrackSimulator:
             )
 
         Qrack.qrack_lib.PhaseParity(
-            self.sid, ctypes.c_double(la), len(q), self._ulonglong_byref(q)
+            self.sid, ctypes.c_double(la), len(q), QrackSimulator._ulonglong_byref(q)
         )
         self._throw_if_error()
 
@@ -3057,7 +3138,7 @@ class QrackSimulator:
                 "QrackStabilizer cannot phase_root_n()! (Create a QrackSimulator instead, also with isTensorNetwork=False.)"
             )
 
-        Qrack.qrack_lib.PhaseRootN(self.sid, n, len(q), self._ulonglong_byref(q))
+        Qrack.qrack_lib.PhaseRootN(self.sid, n, len(q), QrackSimulator._ulonglong_byref(q))
         self._throw_if_error()
 
     def try_separate_1qb(self, qi1):
@@ -3114,7 +3195,7 @@ class QrackSimulator:
             State of all the qubits.
         """
         result = Qrack.qrack_lib.TrySeparateTol(
-            self.sid, len(qs), self._ulonglong_byref(qs), t
+            self.sid, len(qs), QrackSimulator._ulonglong_byref(qs), t
         )
         self._throw_if_error()
         return result
@@ -3130,7 +3211,7 @@ class QrackSimulator:
         Raises:
             Runtimeerror: QrackSimulator raised an exception.
         """
-        result = Qrack.qrack_lib.Separate(self.sid, len(qs), self._ulonglong_byref(qs))
+        result = Qrack.qrack_lib.Separate(self.sid, len(qs), QrackSimulator._ulonglong_byref(qs))
         self._throw_if_error()
 
     def get_unitary_fidelity(self):
@@ -3254,6 +3335,41 @@ class QrackSimulator:
         Qrack.qrack_lib.SetNoiseParameter(self.sid, np)
         self._throw_if_error()
 
+    def set_ace_max_qb(self, qb):
+        """Set "automatic circuit elision" (ACE) max qubits
+
+        If isSchmidtDecompose=True, maximum entangled subsytem size
+        of this simulator will be capped to 'qb', and entangling
+        gates that would exceed that size are replaced with gate
+        shadows.
+
+        Args:
+            qb: maximum subsystem qubits
+
+        Raises:
+            RuntimeError: QrackSimulator raised an exception.
+        """
+        Qrack.qrack_lib.SetAceMaxQb(self.sid, qb)
+        self._throw_if_error()
+
+    def set_sparse_ace_max_mb(self, mb):
+        """Set sparse "automatic circuit elision" (ACE) max memory
+
+        If isSchmidtDecompose=True, isSparse=True, and
+        isOpenCL=False, maximum subsytem size memory MB of this
+        simulator will be capped to 'mb', and entangling gates
+        that would exceed that size are replaced with gate
+        shadows.
+
+        Args:
+            mb: maximum subsystem memory in MB
+
+        Raises:
+            RuntimeError: QrackSimulator raised an exception.
+        """
+        Qrack.qrack_lib.SetSparseAceMaxMb(self.sid, mb)
+        self._throw_if_error()
+
     def normalize(self):
         """Normalize the state
 
@@ -4389,6 +4505,7 @@ class QrackSimulator:
 
         return _data
 
+    @staticmethod
     def get_qiskit_basis_gates():
         return [
             "id",