qubasic 0.6.4__tar.gz → 0.6.5__tar.gz

{qubasic-0.6.4 → qubasic-0.6.5}/CHANGELOG.md

@@ -1,5 +1,10 @@
 # Changelog
 
+## 0.6.5 (2026-06-19)
+
+### Fixed
+- The Bloch sphere Y component had an inverted sign, so a qubit in |+i> was reported at -Y (and |-i> at +Y). `BLOCH`, `DRAW`, `PRINT QUBIT(n)`, and the memory-mapped Bloch-Y field (`$0100 + q*8 + 2`) were all affected, and `BLOCH` mislabeled |+i> as "|-i>". X and Z were already correct. Y now matches the standard convention (`<Y> = +2 Im(<1|rho|0>)`), verified against Qiskit expectation values.
+
 ## 0.6.4 (2026-06-19)
 
 ### Fixed

{qubasic-0.6.4/qubasic.egg-info → qubasic-0.6.5}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: qubasic
-Version: 0.6.4
+Version: 0.6.5
 Summary: Quantum BASIC Interactive Terminal
 Author-email: "Charles C. Norton" <machineelv@gmail.com>
 License-Expression: MIT

{qubasic-0.6.4 → qubasic-0.6.5}/pyproject.toml

@@ -4,7 +4,7 @@ build-backend = "setuptools.build_meta"
 
 [project]
 name = "qubasic"
-version = "0.6.4"
+version = "0.6.5"
 description = "Quantum BASIC Interactive Terminal"
 readme = "README.md"
 requires-python = ">=3.10"

{qubasic-0.6.4 → qubasic-0.6.5/qubasic.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: qubasic
-Version: 0.6.4
+Version: 0.6.5
 Summary: Quantum BASIC Interactive Terminal
 Author-email: "Charles C. Norton" <machineelv@gmail.com>
 License-Expression: MIT

{qubasic-0.6.4 → qubasic-0.6.5}/qubasic_core/__init__.py

@@ -28,7 +28,7 @@ __all__ = [
     'GATE_TABLE', 'GATE_ALIASES',
 ]
 
-__version__ = '0.6.4'
+__version__ = '0.6.5'
 
 def __getattr__(name):
     """Lazy import heavy modules on first access."""

{qubasic-0.6.4 → qubasic-0.6.5}/qubasic_core/display.py

@@ -292,10 +292,12 @@ class DisplayMixin:
 
         rho_00 = np.sum(np.abs(t0)**2)
         rho_11 = np.sum(np.abs(t1)**2)
-        rho_01 = np.sum(np.conj(t0) * t1)
+        # conj(t0)*t1 sums to rho_10 = <1|rho|0>. <X> = 2 Re(rho_10) (Re is
+        # symmetric in rho_01/rho_10), and <Y> = -2 Im(rho_01) = +2 Im(rho_10).
+        rho_10 = np.sum(np.conj(t0) * t1)
 
-        x = float(2 * rho_01.real)
-        y = float(-2 * rho_01.imag)
+        x = float(2 * rho_10.real)
+        y = float(2 * rho_10.imag)
         z = float(rho_00 - rho_11)
 
         return x, y, z

{qubasic-0.6.4 → qubasic-0.6.5}/qubasic_core/memory.py

@@ -112,17 +112,19 @@ class MemoryMixin:
         t0, t1 = t[0].flatten(), t[1].flatten()
         rho_00 = float(np.sum(np.abs(t0) ** 2))
         rho_11 = float(np.sum(np.abs(t1) ** 2))
-        rho_01 = complex(np.sum(np.conj(t0) * t1))
+        # conj(t0)*t1 sums to rho_10 = <1|rho|0>. Bloch X = 2 Re(rho_10) and
+        # Bloch Y = -2 Im(rho_01) = +2 Im(rho_10).
+        rho_10 = complex(np.sum(np.conj(t0) * t1))
         # Recover amplitudes with alpha as the real phase reference; beta then
         # carries the relative phase (real and imaginary parts both exposed).
         alpha = np.sqrt(max(0.0, rho_00))
         if alpha > 1e-12:
-            beta = rho_01 / alpha
+            beta = rho_10 / alpha
         else:
             beta = complex(np.sqrt(max(0.0, rho_11)))
         return [rho_11,
-                float(2 * rho_01.real),
-                float(-2 * rho_01.imag),
+                float(2 * rho_10.real),
+                float(2 * rho_10.imag),
                 float(rho_00 - rho_11),
                 float(alpha),
                 0.0,

{qubasic-0.6.4 → qubasic-0.6.5}/tests/test_features.py

@@ -2595,6 +2595,31 @@ class TestCommandSurfaceCoverage(unittest.TestCase):
         self.assertIsInstance(t._status_prompt(), str)
         capture(t._suggest_command, 'XYZ')
 
+    def test_bloch_vector_axes(self):
+        """_bloch_vector and the $0100 By field carry the right sign on every
+        axis. Regression: the Y component was inverted (|+i> read as -Y)."""
+        import numpy as np
+        from qubasic_core.terminal import QBasicTerminal
+        t = QBasicTerminal()
+        r2 = 1.0 / np.sqrt(2)
+        cases = [
+            (np.array([1, 0], dtype=complex),       (0.0, 0.0, +1.0)),   # |0>  +Z
+            (np.array([0, 1], dtype=complex),       (0.0, 0.0, -1.0)),   # |1>  -Z
+            (np.array([r2, r2], dtype=complex),     (+1.0, 0.0, 0.0)),   # |+>  +X
+            (np.array([r2, -r2], dtype=complex),    (-1.0, 0.0, 0.0)),   # |->  -X
+            (np.array([r2, 1j * r2], dtype=complex),  (0.0, +1.0, 0.0)), # |+i> +Y
+            (np.array([r2, -1j * r2], dtype=complex), (0.0, -1.0, 0.0)), # |-i> -Y
+        ]
+        for sv, (ex, ey, ez) in cases:
+            x, y, z = t._bloch_vector(sv, 0, 1)
+            self.assertAlmostEqual(x, ex, places=6)
+            self.assertAlmostEqual(y, ey, places=6)
+            self.assertAlmostEqual(z, ez, places=6)
+        # Memory-mapped Bloch Y for q0 ($0102) must agree for |+i>.
+        t.num_qubits = 1
+        t.last_sv = np.array([r2, 1j * r2], dtype=complex)
+        self.assertAlmostEqual(t._peek(0x0102), +1.0, places=6)
+
     def test_program_and_io_surface(self):
         import builtins
         t = self._ready()