pyqrack-complex128 1.80.6__tar.gz → 1.83.0__tar.gz

{pyqrack_complex128-1.80.6 → pyqrack_complex128-1.83.0}/Makefile

@@ -30,7 +30,7 @@ build-deps:
 	rm -rf pyqrack/qrack_system/qrack_cl_precompile
 ifneq ($(OS),Windows_NT)
 ifeq ($(QRACK_PRESENT),)
-	git clone https://github.com/unitaryfund/qrack.git; cd qrack; git checkout df0d945f5c6c37be5bb2f619587f76ae2e2187ec; cd ..
+	git clone https://github.com/unitaryfund/qrack.git; cd qrack; git checkout 0033dfb996eb198a05cf34bcc2cf513b52298272; cd ..
 endif
 	mkdir -p qrack/build
 ifeq ($(UNAME_S),Linux)

{pyqrack_complex128-1.80.6/pyqrack_complex128.egg-info → pyqrack_complex128-1.83.0}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: pyqrack-complex128
-Version: 1.80.6
+Version: 1.83.0
 Summary: pyqrack - Pure Python vm6502q/qrack Wrapper
 Home-page: https://github.com/vm6502q/pyqrack
 Author: Daniel Strano

{pyqrack_complex128-1.80.6 → pyqrack_complex128-1.83.0}/pyqrack/qrack_neuron.py

@@ -31,7 +31,9 @@ class QrackNeuron:
         simulator(QrackSimulator): Simulator instance for all synaptic clefts of the neuron
         controls(list(int)): Indices of all "control" qubits, for neuron input
         target(int): Index of "target" qubit, for neuron output
-        tolerance(double): Rounding tolerance
+        activation_fn(NeuronActivationFn): Activation function choice
+        alpha(float): Activation function parameter, if required
+        angles(list[ctypes.c_float]): (or c_double) Memory for neuron prediction angles
     """
 
     def _get_error(self):
@@ -48,15 +50,15 @@ class QrackNeuron:
         target,
         activation_fn=NeuronActivationFn.Sigmoid,
         alpha=1.0,
-        tolerance=sys.float_info.epsilon,
         _init=True,
+        _isTorch=False
     ):
         self.simulator = simulator
         self.controls = controls
         self.target = target
         self.activation_fn = activation_fn
         self.alpha = alpha
-        self.tolerance = tolerance
+        self.angles = None if _isTorch else QrackNeuron._real1_byref([0.0] * (1 << len(controls)))
 
         if not _init:
             return
@@ -66,9 +68,6 @@ class QrackNeuron:
             len(controls),
             QrackNeuron._ulonglong_byref(controls),
             target,
-            activation_fn,
-            alpha,
-            tolerance,
         )
 
         self._throw_if_error()
@@ -91,11 +90,9 @@ class QrackNeuron:
             self.simulator,
             self.controls,
             self.target,
-            self.activation_fn,
-            self.alpha,
-            self.tolerance,
         )
-        self.nid = Qrack.qrack_lib.clone_qneuron(self.simulator.sid)
+        result.nid = Qrack.qrack_lib.clone_qneuron(self.simulator.sid)
+        result.angles = self.angles[:]
         self._throw_if_error()
         return result
 
@@ -110,20 +107,50 @@ class QrackNeuron:
             return (ctypes.c_float * len(a))(*a)
         return (ctypes.c_double * len(a))(*a)
 
-    def set_simulator(self, s):
+    def set_simulator(self, s, controls=None, target=None):
         """Set the neuron simulator
 
         Set the simulator used by this neuron
 
         Args:
             s(QrackSimulator): The simulator to use
+            controls(list[int]): The control qubit IDs to use
+            target(int): The output qubit ID to use
 
         Raises:
             RuntimeError: QrackSimulator raised an exception.
         """
-        Qrack.qrack_lib.set_qneuron_sim(self.nid, s.sid)
+        if controls is None:
+            controls = self.controls
+        if target is None:
+            target = self.target
+        Qrack.qrack_lib.set_qneuron_sim(
+            self.nid,
+            s.sid,
+            len(controls),
+            QrackNeuron._ulonglong_byref(controls),
+            target,
+        )
         self._throw_if_error()
         self.simulator = s
+        self.controls = controls
+        self.target = target
+
+    def set_qubit_ids(self, controls, target=None):
+        """Set the neuron qubit identifiers
+
+        Set the control and target qubits within the simulator
+
+        Args:
+            controls(list[int]): The control qubit IDs to use
+            target(int): The output qubit ID to use
+
+        Raises:
+            RuntimeError: QrackSimulator raised an exception.
+        """
+        if target is None:
+            target = self.target
+        self.set_simulator(self.simulator, controls, target)
 
     def set_angles(self, a):
         """Directly sets the neuron parameters.
@@ -142,8 +169,7 @@ class QrackNeuron:
             raise ValueError(
                 "Angles 'a' in QrackNeuron.set_angles() must contain at least (2 ** len(self.controls)) elements."
             )
-        Qrack.qrack_lib.set_qneuron_angles(self.nid, QrackNeuron._real1_byref(a))
-        self._throw_if_error()
+        self.angles = QrackNeuron._real1_byref(a)
 
     def get_angles(self):
         """Directly gets the neuron parameters.
@@ -154,10 +180,7 @@ class QrackNeuron:
         Raises:
             RuntimeError: QrackNeuron C++ library raised an exception.
         """
-        ket = QrackNeuron._real1_byref([0.0] * (1 << len(self.controls)))
-        Qrack.qrack_lib.get_qneuron_angles(self.nid, ket)
-        self._throw_if_error()
-        return list(ket)
+        return list(self.angles)
 
     def set_alpha(self, a):
         """Set the neuron 'alpha' parameter.
@@ -166,13 +189,8 @@ class QrackNeuron:
         parameter that is applied as a power to its angles, before
         learning and prediction. This makes the activation function
         sharper (or less sharp).
-
-        Raises:
-            RuntimeError: QrackNeuron C++ library raised an exception.
         """
         self.alpha = a
-        Qrack.qrack_lib.set_qneuron_alpha(self.nid, a)
-        self._throw_if_error()
 
     def set_activation_fn(self, f):
         """Sets the activation function of this QrackNeuron
@@ -180,13 +198,8 @@ class QrackNeuron:
         Nonlinear activation functions can be important to neural net
         applications, like DNN. The available activation functions are
         enumerated in `NeuronActivationFn`.
-
-        Raises:
-            RuntimeError: QrackNeuron C++ library raised an exception.
         """
         self.activation_fn = f
-        Qrack.qrack_lib.set_qneuron_activation_fn(self.nid, f)
-        self._throw_if_error()
 
     def predict(self, e=True, r=True):
         """Predict based on training
@@ -205,7 +218,7 @@ class QrackNeuron:
         Raises:
             RuntimeError: QrackNeuron C++ library raised an exception.
         """
-        result = Qrack.qrack_lib.qneuron_predict(self.nid, e, r)
+        result = Qrack.qrack_lib.qneuron_predict(self.nid, self.angles, e, r, self.activation_fn, self.alpha)
         self._throw_if_error()
         return result
 
@@ -221,7 +234,7 @@ class QrackNeuron:
         Raises:
             RuntimeError: QrackNeuron C++ library raised an exception.
         """
-        result = Qrack.qrack_lib.qneuron_unpredict(self.nid, e)
+        result = Qrack.qrack_lib.qneuron_unpredict(self.nid, self.angles, e, self.activation_fn, self.alpha)
         self._throw_if_error()
         return result
 
@@ -237,7 +250,7 @@ class QrackNeuron:
         Raises:
             RuntimeError: QrackNeuron C++ library raised an exception.
         """
-        Qrack.qrack_lib.qneuron_learn_cycle(self.nid, e)
+        Qrack.qrack_lib.qneuron_learn_cycle(self.nid, self.angles, e, self.activation_fn, self.alpha)
         self._throw_if_error()
 
     def learn(self, eta, e=True, r=True):
@@ -256,7 +269,7 @@ class QrackNeuron:
         Raises:
             RuntimeError: QrackNeuron C++ library raised an exception.
         """
-        Qrack.qrack_lib.qneuron_learn(self.nid, eta, e, r)
+        Qrack.qrack_lib.qneuron_learn(self.nid, self.angles, eta, e, r, self.activation_fn, self.alpha)
         self._throw_if_error()
 
     def learn_permutation(self, eta, e=True, r=True):
@@ -275,7 +288,7 @@ class QrackNeuron:
         Raises:
             RuntimeError: QrackNeuron C++ library raised an exception.
         """
-        Qrack.qrack_lib.qneuron_learn_permutation(self.nid, eta, e, r)
+        Qrack.qrack_lib.qneuron_learn_permutation(self.nid, self.angles, eta, e, r, self.activation_fn, self.alpha)
         self._throw_if_error()
 
     @staticmethod

{pyqrack_complex128-1.80.6 → pyqrack_complex128-1.83.0}/pyqrack/qrack_neuron_torch_layer.py

@@ -6,6 +6,7 @@
 # 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
 import itertools
 import math
 import random
@@ -16,12 +17,14 @@ try:
     import torch
     import torch.nn as nn
     from torch.autograd import Function
+    import numpy as np
 except ImportError:
     _IS_TORCH_AVAILABLE = False
 
 from .pauli import Pauli
 from .qrack_neuron import QrackNeuron
 from .qrack_simulator import QrackSimulator
+from .qrack_system import Qrack
 from .neuron_activation_fn import NeuronActivationFn
 
 
@@ -29,66 +32,66 @@ from .neuron_activation_fn import NeuronActivationFn
 param_shift_eps = math.pi / 2
 # Neuron angle initialization
 init_phi = math.asin(0.5)
+# Systemic floating-point type
+fp_type = ctypes.c_float if Qrack.fppow <= 5 else ctypes.c_double
 
 
 class QrackNeuronTorchFunction(Function if _IS_TORCH_AVAILABLE else object):
     """Static forward/backward/apply functions for QrackNeuronTorch"""
 
     @staticmethod
-    def forward(ctx, x, neuron):
-        ctx.neuron = neuron
-        ctx.simulator = neuron.simulator
-        ctx.save_for_backward(x)
-
+    def _apply(angles, neuron):
         # Baseline probability BEFORE applying this neuron's unitary
         pre_prob = neuron.simulator.prob(neuron.target)
 
-        angles = x.detach().cpu().numpy() if x.requires_grad else x.numpy()
-        neuron.set_angles(angles)
-        neuron.predict(True, False)
+        neuron.angles = angles.ctypes.data_as(ctypes.POINTER(fp_type))
 
         # Probability AFTER applying this neuron's unitary
-        post_prob = neuron.simulator.prob(neuron.target)
-        ctx.post_prob = post_prob
+        post_prob = neuron.predict(True, False)
 
+        neuron.angles = None
+
+        # Angle difference
         delta = math.asin(post_prob) - math.asin(pre_prob)
-        ctx.delta = delta
 
         # Return shape: (1,)
-        return x.new_tensor([delta])
+        return delta, max(math.sqrt(max(1.0 - post_prob * post_prob, 0.0)), 1e-6)
 
     @staticmethod
-    def backward(ctx, grad_output):
-        (x,) = ctx.saved_tensors
-        neuron = ctx.neuron
-        neuron.set_simulator(ctx.simulator)
-        post_prob = ctx.post_prob
+    def forward(ctx, x, neuron):
+        ctx.neuron = neuron
+        ctx.simulator = neuron.simulator
+        ctx.save_for_backward(x)
 
         angles = x.detach().cpu().numpy() if x.requires_grad else x.numpy()
+        delta, denom = QrackNeuronTorchFunction._apply(angles, neuron)
 
+        ctx.denom = denom
+
+        return x.new_tensor([delta])
+
+    @staticmethod
+    def _backward(angles, neuron, simulator):
         # Restore simulator to state BEFORE this neuron's unitary
-        neuron.set_angles(angles)
+        neuron.set_simulator(simulator)
+        neuron.angles = angles.ctypes.data_as(ctypes.POINTER(fp_type))
         neuron.unpredict()
         pre_sim = neuron.simulator
 
-        grad_x = torch.zeros_like(x)
+        grad_x = np.zeros(angles.shape[0], dtype=fp_type)
 
-        for i in range(x.shape[0]):
+        for i in range(angles.shape[0]):
             angle = angles[i]
 
             # θ + π/2
             angles[i] = angle + param_shift_eps
-            neuron.set_angles(angles)
             neuron.simulator = pre_sim.clone()
-            neuron.predict(True, False)
-            p_plus = neuron.simulator.prob(neuron.target)
+            p_plus = neuron.predict(True, False)
 
             # θ − π/2
             angles[i] = angle - param_shift_eps
-            neuron.set_angles(angles)
             neuron.simulator = pre_sim.clone()
-            neuron.predict(True, False)
-            p_minus = neuron.simulator.prob(neuron.target)
+            p_minus = neuron.predict(True, False)
 
             # Parameter-shift gradient
             grad_x[i] = 0.5 * (p_plus - p_minus)
@@ -97,9 +100,24 @@ class QrackNeuronTorchFunction(Function if _IS_TORCH_AVAILABLE else object):
 
         # Restore simulator
         neuron.set_simulator(pre_sim)
+        neuron.angles = None
+
+        return grad_x
+
+    @staticmethod
+    def backward(ctx, grad_output):
+        (x,) = ctx.saved_tensors
+        neuron = ctx.neuron
+        simulator = ctx.simulator
+        denom = ctx.denom
+
+        angles = x.detach().cpu().numpy() if x.requires_grad else x.numpy()
+
+        grad_x = QrackNeuronTorchFunction._backward(angles, neuron, simulator)
+        grad_x = torch.tensor(grad_x, device=x.device, dtype=x.dtype)
 
         # Apply chain rule and upstream gradient
-        grad_x *= grad_output[0] / math.sqrt(max(1.0 - post_prob * post_prob, 1e-6))
+        grad_x *= grad_output[0] / denom
 
         return grad_x, None
 
@@ -120,6 +138,10 @@ class QrackNeuronTorch(nn.Module if _IS_TORCH_AVAILABLE else object):
         return QrackNeuronTorchFunction.apply(self.weights, self.neuron)
 
 
+def dummy_post_init_fn(simulator):
+    pass
+
+
 class QrackNeuronTorchLayer(nn.Module if _IS_TORCH_AVAILABLE else object):
     """Torch layer wrapper for QrackNeuron (with maximally expressive set of neurons between inputs and outputs)
 
@@ -143,9 +165,8 @@ class QrackNeuronTorchLayer(nn.Module if _IS_TORCH_AVAILABLE else object):
         lowest_combo_count=0,
         highest_combo_count=2,
         activation=int(NeuronActivationFn.Generalized_Logistic),
-        dtype=torch.float if _IS_TORCH_AVAILABLE else float,
         parameters=None,
-        post_init_fn=lambda simulator: None,
+        post_init_fn=dummy_post_init_fn,
         **kwargs
     ):
         """
@@ -160,7 +181,6 @@ class QrackNeuronTorchLayer(nn.Module if _IS_TORCH_AVAILABLE else object):
             lowest_combo_count (int): (Optional) Lowest combination count of input qubits iterated (0 is bias)
             highest_combo_count (int): (Optional) Highest combination count of input qubits iterated
             activation (int): (Optional) Integer corresponding to choice of activation function from NeuronActivationFn
-            dtype (type): (Optional) dtype of tensor objects used
             parameters (list[float]): (Optional) Flat list of initial neuron parameters, corresponding to little-endian basis states of input + hidden qubits, repeated for ascending combo count, repeated for each output index
             post_init_fn (Callable[QrackSimulator]): (Optional) Function that is applied after forward(x) state initialization, before inference. (As the function depends on nothing but the simulator, it's differentiable.)
         """
@@ -175,7 +195,7 @@ class QrackNeuronTorchLayer(nn.Module if _IS_TORCH_AVAILABLE else object):
             range(input_qubits + hidden_qubits, input_qubits + hidden_qubits + output_qubits)
         )
         self.activation = NeuronActivationFn(activation)
-        self.dtype = dtype
+        self.dtype = torch.float if Qrack.fppow <= 5 else torch.double
         self.apply_fn = QrackNeuronTorchFunction.apply
         self.post_init_fn = post_init_fn
 
@@ -189,15 +209,15 @@ class QrackNeuronTorchLayer(nn.Module if _IS_TORCH_AVAILABLE else object):
                     angles = (
                         (
                             torch.tensor(
-                                parameters[param_count : (param_count + p_count)], dtype=dtype
+                                parameters[param_count : (param_count + p_count)], dtype=self.dtype
                             )
                             if parameters
-                            else torch.zeros(p_count, dtype=dtype)
+                            else torch.zeros(p_count, dtype=self.dtype)
                         )
                     )
                     neurons.append(
                         QrackNeuronTorch(
-                            QrackNeuron(self.simulator, input_subset, output_id, activation), angles
+                            QrackNeuron(self.simulator, input_subset, output_id, activation, _isTorch=True), angles
                         )
                     )
                     param_count += p_count

{pyqrack_complex128-1.80.6 → pyqrack_complex128-1.83.0}/pyqrack/qrack_simulator.py

@@ -77,7 +77,6 @@ class QrackSimulator:
                 "Cannot clone a QrackSimulator and specify its qubit length at the same time, in QrackSimulator constructor!"
             )
 
-        self.is_tensor_network = isTensorNetwork
         self.is_pure_stabilizer = False
 
         if cloneSid > -1:
@@ -1307,15 +1306,11 @@ class QrackSimulator:
 
         Raises:
             RuntimeError: QrackSimulator raised an exception.
-            RuntimeError: QrackSimulator with isTensorNetwork=True option cannot mul()! (Turn off just this option, in the constructor.)
+            RuntimeError: QrackStabilizer cannot mul()! (Create a QrackSimulator instead.)
         """
-        if self.is_tensor_network:
-            raise RuntimeError(
-                "QrackSimulator with isTensorNetwork=True option cannot mul()! (Turn off just this option, in the constructor.)"
-            )
         if self.is_pure_stabilizer:
             raise RuntimeError(
-                "QrackStabilizer cannot mul()! (Create a QrackSimulator instead, also with isTensorNetwork=False.)"
+                "QrackStabilizer cannot mul()! (Create a QrackSimulator instead.)"
             )
 
         if len(q) != len(o):
@@ -1346,15 +1341,11 @@ class QrackSimulator:
 
         Raises:
             RuntimeError: QrackSimulator raised an exception.
-            RuntimeError: QrackSimulator with isTensorNetwork=True option cannot div()! (Turn off just this option, in the constructor.)
+            RuntimeError: QrackStabilizer cannot div()! (Create a QrackSimulator instead.)
         """
-        if self.is_tensor_network:
-            raise RuntimeError(
-                "QrackSimulator with isTensorNetwork=True option cannot div()! (Turn off just this option, in the constructor.)"
-            )
         if self.is_pure_stabilizer:
             raise RuntimeError(
-                "QrackStabilizer cannot div()! (Create a QrackSimulator instead, also with isTensorNetwork=False.)"
+                "QrackStabilizer cannot div()! (Create a QrackSimulator instead.)"
             )
 
         if len(q) != len(o):
@@ -1443,15 +1434,11 @@ class QrackSimulator:
 
         Raises:
             RuntimeError: QrackSimulator raised an exception.
-            RuntimeError: QrackSimulator with isTensorNetwork=True option cannot pown()! (Turn off just this option, in the constructor.)
+            RuntimeError: QrackStabilizer cannot pown()! (Create a QrackSimulator instead.)
         """
-        if self.is_tensor_network:
-            raise RuntimeError(
-                "QrackSimulator with isTensorNetwork=True option cannot pown()! (Turn off just this option, in the constructor.)"
-            )
         if self.is_pure_stabilizer:
             raise RuntimeError(
-                "QrackStabilizer cannot pown()! (Create a QrackSimulator instead, also with isTensorNetwork=False.)"
+                "QrackStabilizer cannot pown()! (Create a QrackSimulator instead.)"
             )
 
         if len(q) != len(o):
@@ -1536,15 +1523,11 @@ class QrackSimulator:
 
         Raises:
             RuntimeError: QrackSimulator raised an exception.
-            RuntimeError: QrackSimulator with isTensorNetwork=True option cannot mcmul()! (Turn off just this option, in the constructor.)
+            RuntimeError: QrackStabilizer cannot mcmul()! (Create a QrackSimulator instead.)
         """
-        if self.is_tensor_network:
-            raise RuntimeError(
-                "QrackSimulator with isTensorNetwork=True option cannot mcmul()! (Turn off just this option, in the constructor.)"
-            )
         if self.is_pure_stabilizer:
             raise RuntimeError(
-                "QrackStabilizer cannot mcmul()! (Create a QrackSimulator instead, also with isTensorNetwork=False.)"
+                "QrackStabilizer cannot mcmul()! (Create a QrackSimulator instead.)"
             )
 
         if len(q) != len(o):
@@ -1578,15 +1561,11 @@ class QrackSimulator:
 
         Raises:
             RuntimeError: QrackSimulator raised an exception.
-            RuntimeError: QrackSimulator with isTensorNetwork=True option cannot mcdiv()! (Turn off just this option, in the constructor.)
+            RuntimeError: QrackStabilizer cannot mcdiv()! (Create a QrackSimulator instead.)
         """
-        if self.is_tensor_network:
-            raise RuntimeError(
-                "QrackSimulator with isTensorNetwork=True option cannot mcdiv()! (Turn off just this option, in the constructor.)"
-            )
         if self.is_pure_stabilizer:
             raise RuntimeError(
-                "QrackStabilizer cannot mcdiv()! (Create a QrackSimulator instead, also with isTensorNetwork=False.)"
+                "QrackStabilizer cannot mcdiv()! (Create a QrackSimulator instead.)"
             )
 
         if len(q) != len(o):
@@ -1686,15 +1665,11 @@ class QrackSimulator:
 
         Raises:
             RuntimeError: QrackSimulator raised an exception.
-            RuntimeError: QrackSimulator with isTensorNetwork=True option cannot mcpown()! (Turn off just this option, in the constructor.)
+            RuntimeError: QrackStabilizer cannot mcpown()! (Create a QrackSimulator instead.)
         """
-        if self.is_tensor_network:
-            raise RuntimeError(
-                "QrackSimulator with isTensorNetwork=True option cannot mcpown()! (Turn off just this option, in the constructor.)"
-            )
         if self.is_pure_stabilizer:
             raise RuntimeError(
-                "QrackStabilizer cannot mcpown()! (Create a QrackSimulator instead, also with isTensorNetwork=False.)"
+                "QrackStabilizer cannot mcpown()! (Create a QrackSimulator instead.)"
             )
 
         if len(q) != len(o):
@@ -1727,15 +1702,11 @@ class QrackSimulator:
 
         Raises:
             RuntimeError: QrackSimulator raised an exception.
-            RuntimeError: QrackSimulator with isTensorNetwork=True option cannot lda()! (Turn off just this option, in the constructor.)
+            RuntimeError: QrackStabilizer cannot lda()! (Create a QrackSimulator instead.)
         """
-        if self.is_tensor_network:
-            raise RuntimeError(
-                "QrackSimulator with isTensorNetwork=True option cannot lda()! (Turn off just this option, in the constructor.)"
-            )
         if self.is_pure_stabilizer:
             raise RuntimeError(
-                "QrackStabilizer cannot lda()! (Create a QrackSimulator instead, also with isTensorNetwork=False.)"
+                "QrackStabilizer cannot lda()! (Create a QrackSimulator instead.)"
             )
 
         Qrack.qrack_lib.LDA(
@@ -1761,15 +1732,11 @@ class QrackSimulator:
 
         Raises:
             RuntimeError: QrackSimulator raised an exception.
-            RuntimeError: QrackSimulator with isTensorNetwork=True option cannot adc()! (Turn off just this option, in the constructor.)
+            RuntimeError: QrackStabilizer cannot adc()! (Create a QrackSimulator instead.)
         """
-        if self.is_tensor_network:
-            raise RuntimeError(
-                "QrackSimulator with isTensorNetwork=True option cannot adc()! (Turn off just this option, in the constructor.)"
-            )
         if self.is_pure_stabilizer:
             raise RuntimeError(
-                "QrackStabilizer cannot adc()! (Create a QrackSimulator instead, also with isTensorNetwork=False.)"
+                "QrackStabilizer cannot adc()! (Create a QrackSimulator instead.)"
             )
 
         Qrack.qrack_lib.ADC(
@@ -1796,15 +1763,11 @@ class QrackSimulator:
 
         Raises:
             RuntimeError: QrackSimulator raised an exception.
-            RuntimeError: QrackSimulator with isTensorNetwork=True option cannot sbc()! (Turn off just this option, in the constructor.)
+            RuntimeError: QrackStabilizer cannot sbc()! (Create a QrackSimulator instead.)
         """
-        if self.is_tensor_network:
-            raise RuntimeError(
-                "QrackSimulator with isTensorNetwork=True option cannot sbc()! (Turn off just this option, in the constructor.)"
-            )
         if self.is_pure_stabilizer:
             raise RuntimeError(
-                "QrackStabilizer cannot sbc()! (Create a QrackSimulator instead, also with isTensorNetwork=False.)"
+                "QrackStabilizer cannot sbc()! (Create a QrackSimulator instead.)"
             )
 
         Qrack.qrack_lib.SBC(
@@ -1832,15 +1795,11 @@ class QrackSimulator:
 
         Raises:
             RuntimeError: QrackSimulator raised an exception.
-            RuntimeError: QrackSimulator with isTensorNetwork=True option cannot hash()! (Turn off just this option, in the constructor.)
+            RuntimeError: QrackStabilizer cannot hash()! (Create a QrackSimulator instead.)
         """
-        if self.is_tensor_network:
-            raise RuntimeError(
-                "QrackSimulator with isTensorNetwork=True option cannot hash()! (Turn off just this option, in the constructor.)"
-            )
         if self.is_pure_stabilizer:
             raise RuntimeError(
-                "QrackStabilizer cannot hash()! (Create a QrackSimulator instead, also with isTensorNetwork=False.)"
+                "QrackStabilizer cannot hash()! (Create a QrackSimulator instead.)"
             )
 
         Qrack.qrack_lib.Hash(
@@ -2150,13 +2109,7 @@ class QrackSimulator:
 
         Raises:
             RuntimeError: QrackSimulator raised an exception.
-            RuntimeError: QrackSimulator with isTensorNetwork=True option cannot compose()! (Turn off just this option, in the constructor.)
         """
-        if self.is_tensor_network:
-            raise RuntimeError(
-                "QrackSimulator with isTensorNetwork=True option cannot compose()! (Turn off just this option, in the constructor.)"
-            )
-
         Qrack.qrack_lib.Compose(self.sid, other.sid, QrackSimulator._ulonglong_byref(q))
         self._throw_if_error()
 
@@ -2170,16 +2123,10 @@ class QrackSimulator:
 
         Raises:
             RuntimeError: QrackSimulator raised an exception.
-            RuntimeError: QrackSimulator with isTensorNetwork=True option cannot decompose()! (Turn off just this option, in the constructor.)
 
         Returns:
             Decomposed subsystem simulator.
         """
-        if self.is_tensor_network:
-            raise RuntimeError(
-                "QrackSimulator with isTensorNetwork=True option cannot decompose()! (Turn off just this option, in the constructor.)"
-            )
-
         other = QrackSimulator()
         Qrack.qrack_lib.destroy(other.sid)
         l = len(q)
@@ -2198,13 +2145,7 @@ class QrackSimulator:
 
         Raises:
             RuntimeError: QrackSimulator raised an exception.
-            RuntimeError: QrackSimulator with isTensorNetwork=True option cannot dispose()! (Turn off just this option, in the constructor.)
         """
-        if self.is_tensor_network:
-            raise RuntimeError(
-                "QrackSimulator with isTensorNetwork=True option cannot dispose()! (Turn off just this option, in the constructor.)"
-            )
-
         l = len(q)
         Qrack.qrack_lib.Dispose(self.sid, l, QrackSimulator._ulonglong_byref(q))
         self._throw_if_error()
@@ -3128,15 +3069,11 @@ class QrackSimulator:
 
         Raises:
             RuntimeError: QrackSimulator raised an exception.
-            RuntimeError: QrackSimulator with isTensorNetwork=True option cannot phase_parity()! (Turn off just this option, in the constructor.)
+            RuntimeError: QrackStabilizer cannot phase_parity()! (Create a QrackSimulator instead.)
         """
-        if self.is_tensor_network:
-            raise RuntimeError(
-                "QrackSimulator with isTensorNetwork=True option cannot phase_parity()! (Turn off just this option, in the constructor.)"
-            )
         if self.is_pure_stabilizer:
             raise RuntimeError(
-                "QrackStabilizer cannot phase_parity()! (Create a QrackSimulator instead, also with isTensorNetwork=False.)"
+                "QrackStabilizer cannot phase_parity()! (Create a QrackSimulator instead.)"
             )
 
         Qrack.qrack_lib.PhaseParity(
@@ -3155,15 +3092,11 @@ class QrackSimulator:
 
         Raises:
             RuntimeError: QrackSimulator raised an exception.
-            RuntimeError: QrackSimulator with isTensorNetwork=True option cannot phase_root_n()! (Turn off just this option, in the constructor.)
+            RuntimeError: QrackStabilizer cannot phase_root_n()! (Create a QrackSimulator instead.)
         """
-        if self.is_tensor_network:
-            raise RuntimeError(
-                "QrackSimulator with isTensorNetwork=True option cannot phase_root_n()! (Turn off just this option, in the constructor.)"
-            )
         if self.is_pure_stabilizer:
             raise RuntimeError(
-                "QrackStabilizer cannot phase_root_n()! (Create a QrackSimulator instead, also with isTensorNetwork=False.)"
+                "QrackStabilizer cannot phase_root_n()! (Create a QrackSimulator instead.)"
             )
 
         Qrack.qrack_lib.PhaseRootN(self.sid, n, len(q), QrackSimulator._ulonglong_byref(q))
@@ -3348,6 +3281,21 @@ class QrackSimulator:
         Qrack.qrack_lib.SetTInjection(self.sid, iti)
         self._throw_if_error()
 
+    def set_use_exact_near_clifford(self, enc):
+        """Set option to use exact near-Clifford simulation
+
+        If t-injection is available, this method turns exact near-Clifford
+        simulation off/on. (Default is on.) Approximate methods are faster.
+
+        Args:
+            enc: "on/off" for exact near-Clifford simulation
+
+        Raises:
+            RuntimeError: QrackSimulator raised an exception.
+        """
+        Qrack.qrack_lib.SetUseExactNearClifford(self.sid, enc)
+        self._throw_if_error()
+
     def set_noise_parameter(self, np):
         """Set noise parameter option
 

{pyqrack_complex128-1.80.6 → pyqrack_complex128-1.83.0}/pyqrack/qrack_system/qrack_system.py

@@ -170,14 +170,14 @@ class QrackSystem:
         self.qrack_lib.HighestProbAllN.argtypes = [c_ulonglong, c_ulonglong, POINTER(c_ulonglong)]
 
         self.qrack_lib.ProbAll.restype = None
-        if self.fppow == 5:
+        if self.fppow <= 5:
             self.qrack_lib.ProbAll.argtypes = [
                 c_ulonglong,
                 c_ulonglong,
                 POINTER(c_ulonglong),
                 POINTER(c_float),
             ]
-        elif self.fppow == 6:
+        else:
             self.qrack_lib.ProbAll.argtypes = [
                 c_ulonglong,
                 c_ulonglong,
@@ -242,7 +242,7 @@ class QrackSystem:
             c_bool,
         ]
 
-        if self.fppow == 5:
+        if self.fppow <= 5:
             self.qrack_lib.FactorizedExpectationFp.restype = c_double
             self.qrack_lib.FactorizedExpectationFp.argtypes = [
                 c_ulonglong,
@@ -288,7 +288,7 @@ class QrackSystem:
                 POINTER(c_float),
                 POINTER(c_float),
             ]
-        elif self.fppow == 6:
+        else:
             self.qrack_lib.FactorizedExpectationFp.restype = c_double
             self.qrack_lib.FactorizedExpectationFp.argtypes = [
                 c_ulonglong,
@@ -377,7 +377,7 @@ class QrackSystem:
             c_bool,
         ]
 
-        if self.fppow == 5:
+        if self.fppow <= 5:
             self.qrack_lib.FactorizedVarianceFp.restype = c_double
             self.qrack_lib.FactorizedVarianceFp.argtypes = [
                 c_ulonglong,
@@ -423,7 +423,7 @@ class QrackSystem:
                 POINTER(c_float),
                 POINTER(c_float),
             ]
-        elif self.fppow == 6:
+        else:
             self.qrack_lib.FactorizedVarianceFp.restype = c_double
             self.qrack_lib.FactorizedVarianceFp.argtypes = [
                 c_ulonglong,
@@ -1220,6 +1220,9 @@ class QrackSystem:
         self.qrack_lib.SetTInjection.restype = None
         self.qrack_lib.SetTInjection.argtypes = [c_ulonglong, c_bool]
 
+        self.qrack_lib.SetUseExactNearClifford.restype = None
+        self.qrack_lib.SetUseExactNearClifford.argtypes = [c_ulonglong, c_bool]
+
         self.qrack_lib.SetNoiseParameter.restype = None
         self.qrack_lib.SetNoiseParameter.argtypes = [c_ulonglong, c_double]
 
@@ -1244,9 +1247,6 @@ class QrackSystem:
             c_ulonglong,
             POINTER(c_ulonglong),
             c_ulonglong,
-            c_ulonglong,
-            c_double,
-            c_double,
         ]
 
         self.qrack_lib.clone_qneuron.restype = c_ulonglong
@@ -1256,55 +1256,60 @@ class QrackSystem:
         self.qrack_lib.destroy_qneuron.argtypes = [c_ulonglong]
 
         self.qrack_lib.set_qneuron_sim.restype = None
-        self.qrack_lib.set_qneuron_sim.argtypes = [c_ulonglong, c_ulonglong]
-
-        self.qrack_lib.set_qneuron_angles.restype = None
-        self.qrack_lib.get_qneuron_angles.restype = None
-
-        if self.fppow == 5:
-            self.qrack_lib.set_qneuron_angles.argtypes = [c_ulonglong, POINTER(c_float)]
-            self.qrack_lib.get_qneuron_angles.argtypes = [c_ulonglong, POINTER(c_float)]
-        elif self.fppow == 6:
-            self.qrack_lib.set_qneuron_angles.argtypes = [
-                c_ulonglong,
-                POINTER(c_double),
-            ]
-            self.qrack_lib.get_qneuron_angles.argtypes = [
-                c_ulonglong,
-                POINTER(c_double),
-            ]
+        self.qrack_lib.set_qneuron_sim.argtypes = [
+            c_ulonglong,
+            c_ulonglong,
+            c_ulonglong,
+            POINTER(c_ulonglong),
+            c_ulonglong
+        ]
 
-        self.qrack_lib.set_qneuron_alpha.restype = None
-        self.qrack_lib.set_qneuron_alpha.argtypes = [c_ulonglong, c_double]
+        if self.fppow <= 5:
+            self.qrack_lib.qneuron_predict.restype = c_double
+            self.qrack_lib.qneuron_predict.argtypes = [c_ulonglong, POINTER(c_float), c_bool, c_bool, c_ulonglong, c_double]
 
-        self.qrack_lib.get_qneuron_alpha.restype = c_double
-        self.qrack_lib.get_qneuron_alpha.argtypes = [c_ulonglong]
+            self.qrack_lib.qneuron_unpredict.restype = c_double
+            self.qrack_lib.qneuron_unpredict.argtypes = [c_ulonglong, POINTER(c_float), c_bool, c_ulonglong, c_double]
 
-        self.qrack_lib.set_qneuron_activation_fn.restype = None
-        self.qrack_lib.set_qneuron_activation_fn.argtypes = [c_ulonglong, c_ulonglong]
+            self.qrack_lib.qneuron_learn_cycle.restype = c_double
+            self.qrack_lib.qneuron_learn_cycle.argtypes = [c_ulonglong, POINTER(c_float), c_bool, c_ulonglong, c_double]
 
-        self.qrack_lib.get_qneuron_activation_fn.restype = c_ulonglong
-        self.qrack_lib.get_qneuron_activation_fn.argtypes = [c_ulonglong]
+            self.qrack_lib.qneuron_learn.restype = None
+            self.qrack_lib.qneuron_learn.argtypes = [c_ulonglong, POINTER(c_float), c_double, c_bool, c_bool, c_ulonglong, c_double]
 
-        self.qrack_lib.qneuron_predict.restype = c_double
-        self.qrack_lib.qneuron_predict.argtypes = [c_ulonglong, c_bool, c_bool]
+            self.qrack_lib.qneuron_learn_permutation.restype = None
+            self.qrack_lib.qneuron_learn_permutation.argtypes = [
+                c_ulonglong,
+                POINTER(c_float),
+                c_double,
+                c_bool,
+                c_bool,
+                c_ulonglong,
+                c_double,
+            ]
+        else:
+            self.qrack_lib.qneuron_predict.restype = c_double
+            self.qrack_lib.qneuron_predict.argtypes = [c_ulonglong, POINTER(c_double), c_bool, c_bool, c_ulonglong, c_double]
 
-        self.qrack_lib.qneuron_unpredict.restype = c_double
-        self.qrack_lib.qneuron_unpredict.argtypes = [c_ulonglong, c_bool]
+            self.qrack_lib.qneuron_unpredict.restype = c_double
+            self.qrack_lib.qneuron_unpredict.argtypes = [c_ulonglong, POINTER(c_double), c_bool, c_ulonglong, c_double]
 
-        self.qrack_lib.qneuron_learn_cycle.restype = c_double
-        self.qrack_lib.qneuron_learn_cycle.argtypes = [c_ulonglong, c_bool]
+            self.qrack_lib.qneuron_learn_cycle.restype = c_double
+            self.qrack_lib.qneuron_learn_cycle.argtypes = [c_ulonglong, POINTER(c_double), c_bool, c_ulonglong, c_double]
 
-        self.qrack_lib.qneuron_learn.restype = None
-        self.qrack_lib.qneuron_learn.argtypes = [c_ulonglong, c_double, c_bool, c_bool]
+            self.qrack_lib.qneuron_learn.restype = None
+            self.qrack_lib.qneuron_learn.argtypes = [c_ulonglong, POINTER(c_double), c_double, c_bool, c_bool, c_ulonglong, c_double]
 
-        self.qrack_lib.qneuron_learn_permutation.restype = None
-        self.qrack_lib.qneuron_learn_permutation.argtypes = [
-            c_ulonglong,
-            c_double,
-            c_bool,
-            c_bool,
-        ]
+            self.qrack_lib.qneuron_learn_permutation.restype = None
+            self.qrack_lib.qneuron_learn_permutation.argtypes = [
+                c_ulonglong,
+                POINTER(c_double),
+                c_double,
+                c_bool,
+                c_bool,
+                c_ulonglong,
+                c_double,
+            ]
 
         self.qrack_lib.init_qcircuit.restype = c_ulonglong
         self.qrack_lib.init_qcircuit.argtypes = [c_bool, c_bool]

{pyqrack_complex128-1.80.6 → pyqrack_complex128-1.83.0/pyqrack_complex128.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: pyqrack-complex128
-Version: 1.80.6
+Version: 1.83.0
 Summary: pyqrack - Pure Python vm6502q/qrack Wrapper
 Home-page: https://github.com/vm6502q/pyqrack
 Author: Daniel Strano

{pyqrack_complex128-1.80.6 → pyqrack_complex128-1.83.0}/setup.py

@@ -7,7 +7,7 @@ from setuptools import setup
 from setuptools.command.build_py import build_py
 
 
-VERSION = "1.80.6"
+VERSION = "1.83.0"
 
 # Read long description from README.
 README_PATH = os.path.join(os.path.abspath(os.path.dirname(__file__)), 'README.md')